Done. -parsing GTF file: /tmp/nxf.XSRTDz90GE/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.gtf * Running CMD: /usr/local/bin/STAR --runThreadN 8 --genomeDir /tmp/nxf.XSRTDz90GE/ctat_genome_lib_build_dir/ref_genome.fa.star.idx --outSAMtype BAM SortedByCoordinate --twopassMode None --alignSJDBoverhangMin 10 --genomeSuffixLengthMax 10000 --limitBAMsortRAM 47271261705 --alignInsertionFlush Right --alignMatesGapMax 100000 --alignIntronMax 100000 --outSAMattributes All --readFilesIn /tmp/nxf.XSRTDz90GE/aih-tih-sc-da06ce-R1_A23YTGFLT4_1_1.fastp.fastq.gz /tmp/nxf.XSRTDz90GE/aih-tih-sc-da06ce-R1_A23YTGFLT4_1_2.fastp.fastq.gz --genomeFastaFiles /tmp/nxf.XSRTDz90GE/fi_workdir/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.fa --outSAMfilter KeepAllAddedReferences --sjdbGTFfile /tmp/nxf.XSRTDz90GE/fi_workdir/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.gtf --alignSJstitchMismatchNmax 5 -1 5 5 --scoreGapNoncan -6 --readFilesCommand 'gunzip -c' * Running CMD: mv Aligned.sortedByCoord.out.bam aih-tih-sc-da06ce-R1_A23YTGFLT4_1.star.sortedByCoord.out.bam * Running CMD: samtools index aih-tih-sc-da06ce-R1_A23YTGFLT4_1.star.sortedByCoord.out.bam 23:18:53 : INFO : Done. -parsing GTF file: /tmp/nxf.XSRTDz90GE/fi_workdir/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.gtf -parsing /tmp/nxf.XSRTDz90GE/fi_workdir/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.star.cSorted.dupsMarked.bam [10000] [20000] [30000] [40000] [50000] [60000] [70000] [80000] [90000] [100000] [110000] [120000] [130000] [140000] [150000] [160000] [170000] [180000] [190000] [200000] [210000] [220000] [230000] [240000] [250000] [260000] [270000] [280000] [290000] [300000] [310000] [320000] [330000] [340000] [350000] [360000] [370000] [380000] -done parsing /tmp/nxf.XSRTDz90GE/fi_workdir/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.star.cSorted.dupsMarked.bam. Extracting junction info. junction read elimination tally: $VAR1 = { 'per_id < 96' => 1796, 'low complexity anchor region' => 1, ' ** passed ** ' => 1465, 'num genes matched < 2' => 1, 'exons hit < 2' => 1, 'num_hits: 2 != num_counted_on_fusion_contigs 1 ' => 101, 'num_hits: 4 != num_counted_on_fusion_contigs 1 ' => 2, 'seq-similar region overlap' => 117, 'small anchor length' => 54, 'excessive soft clipping' => 3759 }; -writing fusion junction support info. -outputting the spanning read info: /tmp/nxf.XSRTDz90GE/fi_workdir/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.star.cSorted.dupsMarked.bam.fusion_spanning_info. - counting read alignments among fusion contigs. [1000] [2000] [3000] [4000] [5000] [6000] [7000] -fusion SPANNING read extraction for scaff: AC092807.3--DDAH1 [1000] [2000] [3000] [4000] [5000] [6000] [7000] [8000] [9000] [10000] [11000] [12000] [13000] [14000] [15000] [16000] [17000] [18000] [19000] [20000] [21000] [22000] [23000] [24000] [25000] [26000] [27000] [28000] [29000] [30000] [31000] [32000] [33000] [34000] [35000] [36000] [37000] [38000] [39000] [40000] [41000] [42000] [43000] [44000] [45000] [46000] [47000] [48000] [49000] [50000] [51000] [52000] [53000] [54000] [55000] [56000] [57000] [58000] [59000] [60000] [61000] [62000] [63000] [64000] [65000] [66000] [67000] [68000] [69000] [70000] [71000] [72000] [73000] [74000] [75000] [76000] [77000] [78000] [79000] [80000] [81000] [82000] [83000] [84000] [85000] [86000] [87000] [88000] [89000] [90000] [91000] [92000] [93000] [94000] [95000] [96000] [97000] [98000] [99000] [100000] [101000] [102000] [103000] [104000] [105000] [106000] [107000] [108000] [109000] [110000] [111000] [112000] [113000] [114000] [115000] [116000] [117000] [118000] [119000] [120000] [121000] [122000] [123000] [124000] [125000] [126000] [127000] [128000] [129000] [130000] [131000] [132000] [133000] [134000] [135000] [136000] [137000] [138000] [139000] [140000] [141000] [142000] [143000] [144000] [145000] [146000] [147000] [148000] [149000] [150000] [151000] [152000] [153000] [154000] [155000] [156000] [157000] [158000] [159000] [160000] [161000] [162000] [163000] [164000] [165000] [166000] [167000] [168000] [169000] [170000] [171000] [172000] [173000] [174000] [175000] [176000] [177000] [178000] [179000] [180000] [181000] [182000] [183000] [184000] [185000] [186000] [187000] [188000] [189000] [190000] [191000] [192000] [193000] [194000] [195000] [196000] [197000] [198000] [199000] [200000] [201000] [202000] [203000] [204000] -fusion SPANNING read extraction for scaff: CCAR2--AC055854.1 [1000] [2000] [3000] [4000] [5000] [6000] [7000] [8000] [9000] [10000] [11000] [12000] [13000] [14000] [15000] [16000] [17000] [18000] [19000] [20000] [21000] [22000] [23000] [24000] [25000] [26000] -fusion SPANNING read extraction for scaff: PAPOLA--RBMXP1 [1000] [2000] [3000] [4000] [5000] [6000] [7000] [8000] [9000] [10000] [11000] [12000] [13000] [14000] [15000] [16000] [17000] [18000] [19000] [20000] [21000] [22000] [23000] [24000] [25000] [26000] -fusion SPANNING read extraction for scaff: PAPOLA--FO393412.1 [1000] [2000] [3000] [4000] [5000] [6000] [7000] [8000] [9000] [10000] [11000] [12000] [13000] [14000] -fusion SPANNING read extraction for scaff: PCDHA10--PCDHAC1 [1000] [2000] [3000] [4000] [5000] [6000] [7000] [8000] [9000] [10000] [11000] [12000] [13000] [14000] -fusion SPANNING read extraction for scaff: PCDHA10--PCDHAC2 [1000] [2000] [3000] [4000] [5000] [6000] [7000] [8000] [9000] [10000] [11000] [12000] [13000] [14000] [15000] [16000] -fusion SPANNING read extraction for scaff: PCDHA10--PCDHA13 [1000] [2000] [3000] [4000] [5000] [6000] [7000] [8000] [9000] [10000] [11000] [12000] [13000] [14000] [15000] [16000] [17000] [18000] [19000] [20000] [21000] [22000] -fusion SPANNING read extraction for scaff: PRSS27--TPSB2 [1000] [2000] [3000] [4000] [5000] [6000] [7000] [8000] [9000] [10000] [11000] [12000] [13000] -fusion SPANNING read extraction for scaff: STON1--GTF2A1L [1000] [2000] [3000] [4000] [5000] [6000] [7000] [8000] [9000] [10000] [11000] [12000] [13000] [14000] [15000] [16000] [17000] [18000] [19000] [20000] -fusion SPANNING read extraction for scaff: TLK2--AC240565.1 [1000] [2000] [3000] [4000] [5000] [6000] [7000] [8000] [9000] [10000] [11000] [12000] [13000] [14000] [15000] [16000] [17000] [18000] -fusion SPANNING read extraction for scaff: TVP23C--CDRT4 -filtered reads reasons: $VAR1 = { 'seq similar region alignment' => 13295, 'lacks exon overlap' => 3691 }; EM: Starting log likelihood: -1603.300294 EM: Round [1] log likelihood: -1603.212216 EM: Round [2] log likelihood: -1603.167964 EM: Round [3] log likelihood: -1603.126749 EM: Round [4] log likelihood: -1603.088181 EM: Round [5] log likelihood: -1603.052008 EM: Round [6] log likelihood: -1603.018010 EM: Round [7] log likelihood: -1602.985996 EM: Round [8] log likelihood: -1602.955798 EM: Round [9] log likelihood: -1602.927269 EM: Round [10] log likelihood: -1602.900276 EM: Round [11] log likelihood: -1602.874703 EM: Round [12] log likelihood: -1602.850445 EM: Round [13] log likelihood: -1602.827408 EM: Round [14] log likelihood: -1602.805508 EM: Round [15] log likelihood: -1602.784666 EM: Round [16] log likelihood: -1602.764814 EM: Round [17] log likelihood: -1602.745887 EM: Round [18] log likelihood: -1602.727828 EM: Round [19] log likelihood: -1602.710583 EM: Round [20] log likelihood: -1602.694104 EM: Round [21] log likelihood: -1602.678346 EM: Round [22] log likelihood: -1602.663266 EM: Round [23] log likelihood: -1602.648827 EM: Round [24] log likelihood: -1602.634994 EM: Round [25] log likelihood: -1602.621732 EM: Round [26] log likelihood: -1602.609013 EM: Round [27] log likelihood: -1602.596806 EM: Round [28] log likelihood: -1602.585086 EM: Round [29] log likelihood: -1602.573828 EM: Round [30] log likelihood: -1602.563008 EM: Round [31] log likelihood: -1602.552605 EM: Round [32] log likelihood: -1602.542598 EM: Round [33] log likelihood: -1602.532969 EM: Round [34] log likelihood: -1602.523699 EM: Round [35] log likelihood: -1602.514771 EM: Round [36] log likelihood: -1602.506170 EM: Round [37] log likelihood: -1602.497881 EM: Round [38] log likelihood: -1602.489890 EM: Round [39] log likelihood: -1602.482182 EM: Round [40] log likelihood: -1602.474746 EM: Round [41] log likelihood: -1602.467570 EM: Round [42] log likelihood: -1602.460642 EM: Round [43] log likelihood: -1602.453951 EM: Round [44] log likelihood: -1602.447489 EM: Round [45] log likelihood: -1602.441244 EM: Round [46] log likelihood: -1602.435208 EM: Round [47] log likelihood: -1602.429372 EM: Round [48] log likelihood: -1602.423728 EM: Round [49] log likelihood: -1602.418269 EM: Round [50] log likelihood: -1602.412986 EM: Round [51] log likelihood: -1602.407873 EM: Round [52] log likelihood: -1602.402923 EM: Round [53] log likelihood: -1602.398129 EM: Round [54] log likelihood: -1602.393486 EM: Round [55] log likelihood: -1602.388988 EM: Round [56] log likelihood: -1602.384628 EM: Round [57] log likelihood: -1602.380403 EM: Round [58] log likelihood: -1602.376306 EM: Round [59] log likelihood: -1602.372333 EM: Round [60] log likelihood: -1602.368479 EM: Round [61] log likelihood: -1602.364740 EM: Round [62] log likelihood: -1602.361112 EM: Round [63] log likelihood: -1602.357590 EM: Round [64] log likelihood: -1602.354171 EM: Round [65] log likelihood: -1602.350851 EM: Round [66] log likelihood: -1602.347626 EM: Round [67] log likelihood: -1602.344494 EM: Round [68] log likelihood: -1602.341450 EM: Round [69] log likelihood: -1602.338492 EM: Round [70] log likelihood: -1602.335617 EM: Round [71] log likelihood: -1602.332822 EM: Round [72] log likelihood: -1602.330104 EM: Round [73] log likelihood: -1602.327460 EM: Round [74] log likelihood: -1602.324889 EM: Round [75] log likelihood: -1602.322387 EM: Round [76] log likelihood: -1602.319952 EM: Round [77] log likelihood: -1602.317583 EM: Round [78] log likelihood: -1602.315276 EM: Round [79] log likelihood: -1602.313030 EM: Round [80] log likelihood: -1602.310843 EM: Round [81] log likelihood: -1602.308713 EM: Round [82] log likelihood: -1602.306638 EM: Round [83] log likelihood: -1602.304617 EM: Round [84] log likelihood: -1602.302647 EM: Round [85] log likelihood: -1602.300727 EM: Round [86] log likelihood: -1602.298856 EM: Round [87] log likelihood: -1602.297031 EM: Round [88] log likelihood: -1602.295252 EM: Round [89] log likelihood: -1602.293518 EM: Round [90] log likelihood: -1602.291826 EM: Round [91] log likelihood: -1602.290175 EM: Round [92] log likelihood: -1602.288565 EM: Round [93] log likelihood: -1602.286993 EM: Round [94] log likelihood: -1602.285460 EM: Round [95] log likelihood: -1602.283963 EM: Round [96] log likelihood: -1602.282502 EM: Round [97] log likelihood: -1602.281076 EM: Round [98] log likelihood: -1602.279683 EM: Round [99] log likelihood: -1602.278323 EM: Round [100] log likelihood: -1602.276995 EM: Round [101] log likelihood: -1602.275697 EM: Round [102] log likelihood: -1602.274430 EM: Round [103] log likelihood: -1602.273191 EM: Round [104] log likelihood: -1602.271981 EM: Round [105] log likelihood: -1602.270798 EM: Round [106] log likelihood: -1602.269642 EM: Round [107] log likelihood: -1602.268512 EM: Round [108] log likelihood: -1602.267408 EM: Round [109] log likelihood: -1602.266328 EM: Round [110] log likelihood: -1602.265272 EM: Round [111] log likelihood: -1602.264239 EM: Round [112] log likelihood: -1602.263229 EM: Round [113] log likelihood: -1602.262241 EM: Round [114] log likelihood: -1602.261274 EM: Round [115] log likelihood: -1602.260329 EM: Round [116] log likelihood: -1602.259404 EM: Round [117] log likelihood: -1602.258498 EM: Round [118] log likelihood: -1602.257613 EM: Round [119] log likelihood: -1602.256746 EM: Round [120] log likelihood: -1602.255897 EM: Round [121] log likelihood: -1602.255066 EM: Round [122] log likelihood: -1602.254253 EM: Round [123] log likelihood: -1602.253457 EM: Round [124] log likelihood: -1602.252678 EM: Round [125] log likelihood: -1602.251915 EM: Round [126] log likelihood: -1602.251167 EM: Round [127] log likelihood: -1602.250435 EM: Round [128] log likelihood: -1602.249718 EM: Round [129] log likelihood: -1602.249016 EM: Round [130] log likelihood: -1602.248328 EM: Round [131] log likelihood: -1602.247654 EM: Round [132] log likelihood: -1602.246994 EM: Round [133] log likelihood: -1602.246347 EM: Round [134] log likelihood: -1602.245713 EM: Round [135] log likelihood: -1602.245092 EM: Round [136] log likelihood: -1602.244484 EM: Round [137] log likelihood: -1602.243887 EM: Round [138] log likelihood: -1602.243303 EM: Round [139] log likelihood: -1602.242730 EM: Round [140] log likelihood: -1602.242168 EM: Round [141] log likelihood: -1602.241617 EM: Round [142] log likelihood: -1602.241078 EM: Round [143] log likelihood: -1602.240548 EM: Round [144] log likelihood: -1602.240030 EM: Round [145] log likelihood: -1602.239521 EM: Round [146] log likelihood: -1602.239022 EM: Round [147] log likelihood: -1602.238533 EM: Round [148] log likelihood: -1602.238053 EM: Round [149] log likelihood: -1602.237583 EM: Round [150] log likelihood: -1602.237122 EM: Round [151] log likelihood: -1602.236669 EM: Round [152] log likelihood: -1602.236225 EM: Round [153] log likelihood: -1602.235790 EM: Round [154] log likelihood: -1602.235363 EM: Round [155] log likelihood: -1602.234944 EM: Round [156] log likelihood: -1602.234533 EM: Round [157] log likelihood: -1602.234130 EM: Round [158] log likelihood: -1602.233734 EM: Round [159] log likelihood: -1602.233346 EM: Round [160] log likelihood: -1602.232966 EM: Round [161] log likelihood: -1602.232592 EM: Round [162] log likelihood: -1602.232226 EM: Round [163] log likelihood: -1602.231866 EM: Round [164] log likelihood: -1602.231513 EM: Round [165] log likelihood: -1602.231166 EM: Round [166] log likelihood: -1602.230827 EM: Round [167] log likelihood: -1602.230493 EM: Round [168] log likelihood: -1602.230165 EM: Round [169] log likelihood: -1602.229844 EM: Round [170] log likelihood: -1602.229529 EM: Round [171] log likelihood: -1602.229219 EM: Round [172] log likelihood: -1602.228915 EM: Round [173] log likelihood: -1602.228617 EM: Round [174] log likelihood: -1602.228324 EM: Round [175] log likelihood: -1602.228037 EM: Round [176] log likelihood: -1602.227755 EM: Round [177] log likelihood: -1602.227478 EM: Round [178] log likelihood: -1602.227206 EM: Round [179] log likelihood: -1602.226939 EM: Round [180] log likelihood: -1602.226677 EM: Round [181] log likelihood: -1602.226419 EM: Round [182] log likelihood: -1602.226167 EM: Round [183] log likelihood: -1602.225919 EM: Round [184] log likelihood: -1602.225675 EM: Round [185] log likelihood: -1602.225436 EM: Round [186] log likelihood: -1602.225201 EM: Round [187] log likelihood: -1602.224971 EM: Round [188] log likelihood: -1602.224744 EM: Round [189] log likelihood: -1602.224522 EM: Round [190] log likelihood: -1602.224304 EM: Round [191] log likelihood: -1602.224089 EM: Round [192] log likelihood: -1602.223879 EM: Round [193] log likelihood: -1602.223672 EM: Round [194] log likelihood: -1602.223469 EM: Round [195] log likelihood: -1602.223270 EM: Round [196] log likelihood: -1602.223074 EM: Round [197] log likelihood: -1602.222882 EM: Round [198] log likelihood: -1602.222693 EM: Round [199] log likelihood: -1602.222507 EM: Round [200] log likelihood: -1602.222325 EM: Round [201] log likelihood: -1602.222146 EM: Round [202] log likelihood: -1602.221970 EM: Round [203] log likelihood: -1602.221798 EM: Round [204] log likelihood: -1602.221628 EM: Round [205] log likelihood: -1602.221462 EM: Round [206] log likelihood: -1602.221298 EM: Round [207] log likelihood: -1602.221138 EM: Round [208] log likelihood: -1602.220980 EM: Round [209] log likelihood: -1602.220825 EM: Round [210] log likelihood: -1602.220673 EM: Round [211] log likelihood: -1602.220523 EM: Round [212] log likelihood: -1602.220376 EM: Round [213] log likelihood: -1602.220232 EM: Round [214] log likelihood: -1602.220090 EM: Round [215] log likelihood: -1602.219951 EM: Round [216] log likelihood: -1602.219814 EM: Round [217] log likelihood: -1602.219680 EM: Round [218] log likelihood: -1602.219548 EM: Round [219] log likelihood: -1602.219418 EM: Round [220] log likelihood: -1602.219291 EM: Round [221] log likelihood: -1602.219165 EM: Round [222] log likelihood: -1602.219043 EM: Round [223] log likelihood: -1602.218922 EM: Round [224] log likelihood: -1602.218803 EM: Round [225] log likelihood: -1602.218687 EM: Round [226] log likelihood: -1602.218572 EM: Round [227] log likelihood: -1602.218460 EM: Round [228] log likelihood: -1602.218349 EM: Round [229] log likelihood: -1602.218240 EM: Round [230] log likelihood: -1602.218134 EM: Round [231] log likelihood: -1602.218029 EM: Round [232] log likelihood: -1602.217926 EM: Round [233] log likelihood: -1602.217825 EM: Round [234] log likelihood: -1602.217725 EM: Stopping iterations at round 234 due to insufficient improvement in likelihood. [bam_sort_core] merging from 0 files and 8 in-memory blocks... -extracting read coordinates from /tmp/nxf.XSRTDz90GE/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.consolidated.bam into /tmp/nxf.XSRTDz90GE/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.consolidated.bam.read_coords CMD: touch /tmp/nxf.XSRTDz90GE/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.consolidated.bam.read_coords.ok CMD: /usr/bin/sort --parallel=1 -S10G -T . -k1,1 -k2,2 -k4,4n /tmp/nxf.XSRTDz90GE/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.consolidated.bam.read_coords > /tmp/nxf.XSRTDz90GE/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.consolidated.bam.read_coords.sort_by_readname CMD: /usr/bin/sort --parallel=1 -S10G -T . -k1,1 -k3,3n /tmp/nxf.XSRTDz90GE/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.consolidated.bam.frag_coords > /tmp/nxf.XSRTDz90GE/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.consolidated.bam.frag_coords.coord_sorted CMD: /usr/local/src/STAR-Fusion/FusionInspector/FusionFilter/util/blast_filter.pl --fusion_preds /tmp/nxf.XSRTDz90GE/fi_workdir/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.fusion_preds.coalesced.summary.EMadj.min_frag_thresh.wSpliceInfo --genome_lib_dir ctat_genome_lib_build_dir CMD: /usr/local/src/STAR-Fusion/FusionInspector/FusionFilter/util/promiscuity_filter.pl --fusion_preds /tmp/nxf.XSRTDz90GE/fi_workdir/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.fusion_preds.coalesced.summary.EMadj.min_frag_thresh.wSpliceInfo.post_blast_filter --genome_lib_dir ctat_genome_lib_build_dir --max_promiscuity 10 --min_pct_dom_promiscuity 50 CMD: cp /tmp/nxf.XSRTDz90GE/fi_workdir/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.fusion_preds.coalesced.summary.EMadj.min_frag_thresh.wSpliceInfo.post_blast_filter.post_promisc_filter /tmp/nxf.XSRTDz90GE/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.post_blast_and_promiscuity_filter -done, see /tmp/nxf.XSRTDz90GE/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.post_blast_and_promiscuity_filter.annot_filter.pass EM: Starting log likelihood: -1560.292373 EM: Round [1] log likelihood: -1560.251741 EM: Round [2] log likelihood: -1560.251647 EM: Stopping iterations at round 2 due to insufficient improvement in likelihood. -total frags in /tmp/nxf.XSRTDz90GE/aih-tih-sc-da06ce-R1_A23YTGFLT4_1_1.fastp.fastq.gz: 296209672 Warning - not locating file: /tmp/nxf.XSRTDz90GE/aih-tih-sc-da06ce-R1_A23YTGFLT4_1.gmap_trinity_GG.fusions.gff3.bed WARNING, no entry stored in dbm for [ENSG00000253125.1] WARNING, no entry stored in dbm for [ENSG00000280136.2] WARNING, no entry stored in dbm for [ENSG00000282057.1] WARNING, no entry stored in dbm for [ENSG00000282057.1] WARNING, no entry stored in dbm for [ENSG00000282057.1] ls: cannot access 'IGV_inputs': No such file or directory