Groups keyboard shortcuts have been updated
Dismiss
See shortcuts

Error while using v1.0.15 on GSNAP aligned bam files

39 views
Skip to first unread message

napr

unread,
Feb 17, 2017, 10:25:27 AM2/17/17
to strelka-discuss
I am using strelka v1.0.15 on GSNAP aligned RNA-Seq data and i'm getting this error

ERROR: Failed system call: '/data2/bsi/staff_analysis/m081429/Somatic_workflow_paper/install_tools1/STRELKA/strelka_workflow-1.0.15/libexec/strelka2 -clobber -filter-unanchored -min-paired-align-score 20 -min-single-align-score 10 -min-qscore 0 -report-range-begin 25000001 -report-range-end 50000000 -samtools-reference '/data2/bsi/staff_analysis/m081429/Somatic_workflow_paper/ref_files/GRCh37_gencode_v19_CTAT_lib_July272016/ref_genome.fa' -max-window-mismatch 3 20 -print-used-allele-counts -bam-seq-name 'chr6' -genome-size 2861343699 -max-indel-size 50 -indel-nonsite-match-prob 0.5 --min-contig-open-end-support 35 --somatic-snv-rate 0.000001 --shared-site-error-rate 0.0000005 --shared-site-error-strand-bias-fraction 0.5 --somatic-indel-rate 0.000001 --shared-indel-error-rate 0.000001 --tier2-min-single-align-score 5 --tier2-min-paired-align-score 5 --tier2-single-align-score-rescue-mode --tier2-mismatch-density-filter-count 10 --tier2-no-filter-unanchored --tier2-indel-nonsite-match-prob 0.25 --tier2-include-singleton --tier2-include-anomalous -bam-file /data2/external_data/Sun_Zhifu_zxs01/GSE40419/ERP001058_analysis_run_tum_pair_with_gsnap/ALIGNMENT/GSNAP/ERR164479_1.GSNAP.gatkin.splitNC.realign.recaliber.bam --tumor-bam-file /data2/external_data/Sun_Zhifu_zxs01/GSE40419/ERP001058_analysis_run_tum_pair_with_gsnap/ALIGNMENT/GSNAP/ERR164559_1.tumor.GSNAP.gatkin.splitNC.realign.recaliber.bam --somatic-snv-file /data2/external_data/Sun_Zhifu_zxs01/GSE40419/ERP001058_analysis_run_tum_pair_with_gsnap/SOMATIC_SNVINDEL/GSNAP/ERR164559_1.tumor.GSNAP.gatkin.splitNC.realign.recaliber.strelka/chromosomes/chr6/bins/0001/somatic.snvs.unfiltered.vcf --somatic-indel-file /data2/external_data/Sun_Zhifu_zxs01/GSE40419/ERP001058_analysis_run_tum_pair_with_gsnap/SOMATIC_SNVINDEL/GSNAP/ERR164559_1.tumor.GSNAP.gatkin.splitNC.realign.recaliber.strelka/chromosomes/chr6/bins/0001/somatic.indels.unfiltered.vcf --variant-window-flank-file 50 /data2/external_data/Sun_Zhifu_zxs01/GSE40419/ERP001058_analysis_run_tum_pair_with_gsnap/SOMATIC_SNVINDEL/GSNAP/ERR164559_1.tumor.GSNAP.gatkin.splitNC.realign.recaliber.strelka/chromosomes/chr6/bins/0001/somatic.indels.unfiltered.vcf.window --max-input-depth 10000 --ignore-conflicting-read-names --report-file /data2/external_data/Sun_Zhifu_zxs01/GSE40419/ERP001058_analysis_run_tum_pair_with_gsnap/SOMATIC_SNVINDEL/GSNAP/ERR164559_1.tumor.GSNAP.gatkin.splitNC.realign.recaliber.strelka/chromosomes/chr6/bins/0001/strelka.stats >| /data2/external_data/Sun_Zhifu_zxs01/GSE40419/ERP001058_analysis_run_tum_pair_with_gsnap/SOMATIC_SNVINDEL/GSNAP/ERR164559_1.tumor.GSNAP.gatkin.splitNC.realign.recaliber.strelka/chromosomes/chr6/bins/0001/strelka.stdout 2>| /data2/external_data/Sun_Zhifu_zxs01/GSE40419/ERP001058_analysis_run_tum_pair_with_gsnap/SOMATIC_SNVINDEL/GSNAP/ERR164559_1.tumor.GSNAP.gatkin.splitNC.realign.recaliber.strelka/chromosomes/chr6/bins/0001/strelka.stderr'

 at /data2/bsi/staff_analysis/m081429/Somatic_workflow_paper/install_tools1/STRELKA/strelka_workflow-1.0.15/libexec/../lib/Utils.pm line 44
        Utils::errorX('Failed system call: \'/data2/bsi/staff_analysis/m081429/Somat...') called at /data2/bsi/staff_analysis/m081429/Somatic_workflow_paper/install_tools1/STRELKA/strelka_workflow-1.0.15/libexec/../lib/Utils.pm line 57
        Utils::executeCmd('/data2/bsi/staff_analysis/m081429/Somatic_workflow_paper/inst...', 0) called at /data2/bsi/staff_analysis/m081429/Somatic_workflow_paper/install_tools1/STRELKA/strelka_workflow-1.0.15/libexec/callSomaticVariants.pl line 225
make: *** [/data2/external_data/Sun_Zhifu_zxs01/GSE40419/ERP001058_analysis_run_tum_pair_with_gsnap/SOMATIC_SNVINDEL/GSNAP/ERR164559_1.tumor.GSNAP.gatkin.splitNC.realign.recaliber.strelka/chromosomes/chr6/bins/0001/task.complete] Error 1

ERROR: Exception caught in align_pos() while realigning segment: 2 of read: STARLING_READ id: 718369 genomic_mapping?: tier1-mapped
full_segment_info:
key: ERR164479.10781396/1
id: 718369
seq:  GGGCTGTCTCTCCACCTCCTCACATTATGCTAACAGGAACGCAGACAGGCGTGACAGGAAAGATGATCGGGGAGGGAACACAGGTCAGTGTGGGGACAGGG
qual: CCCFFFFFGHHHHJJIIIHIIIJIIJGJJJJIJJJJJGIHGIGIGJJGEGG';FHGEGHFFEEDBECDBA>@BDD<B??????8@>ACCCCC@>@?BB?@B
GENOME ALIGNMENT pos: 31236721 strand: F path: 15S3M6N35I48M
is_realigned? 0
buffer_pos: 0
partial_segment 1 :
GENOME ALIGNMENT pos: 31236721 strand: F path: 15S3M
is_realigned? 0
buffer_pos: 31236706
partial_segment 2 :
GENOME ALIGNMENT pos: 31236730 strand: F path: 35I48M
is_realigned? 0
buffer_pos: 31236695

FATAL_ERROR: strelka EXCEPTION: ERROR: Exemplar alignment contains indel not found in the overlap indel set
        Indel: INDEL pos: 31236730 type: INSERT len: 35 swap_dlen: 0

Reply all
Reply to author
Forward
0 new messages