Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/bf/e91e45954716242b1c070b429fcf94/659_cAd-T1-TRNA-1_B23WHTKLT4_1.Chimeric.out.junction
Downloading: s3://natera-platform-sandbox/pipeline-resources/AIH/rna/GRCh38/starfusion/ctat_genome_lib_build_dir
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/b2/83cc8d0a0a2e7c01c8e232a6dcd078/.command.sh
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/2e/226b73e72c865dbbc9d24e127018f0/659_cAd-T1-TRNA-1_B23WHTKLT4_1.Aligned.sortedByCoord.out.bam.bai
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/bf/e91e45954716242b1c070b429fcf94/659_cAd-T1-TRNA-1_B23WHTKLT4_1.Aligned.sortedByCoord.out.bam
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/b2/83cc8d0a0a2e7c01c8e232a6dcd078/.command.run
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/bf/e91e45954716242b1c070b429fcf94/659_cAd-T1-TRNA-1_B23WHTKLT4_1.SJ.out.tab
==> STAGING COMPLETE (7 inputs)
00:05:19 : INFO : -reading targets list: ctat_genome_lib_build_dir/ref_annot.gtf.mini.sortu
-mapping reads to genes
[360000], rate=21600000.00/min
[370000], rate=22200000.00/min
[380000], rate=22800000.00/min
[390000], rate=23400000.00/min
[400000], rate=24000000.00/min
[410000], rate=24600000.00/min
[420000], rate=25200000.00/min
[430000], rate=25800000.00/min
[440000], rate=26400000.00/min
[450000], rate=27000000.00/min
[460000], rate=27600000.00/min
[470000], rate=28200000.00/min
[480000], rate=28800000.00/min
[490000], rate=29400000.00/min
[500000], rate=30000000.00/min
[510000], rate=30600000.00/min
[520000], rate=31200000.00/min
[530000], rate=31800000.00/min
[540000], rate=32400000.00/min
[550000], rate=33000000.00/min
[560000], rate=33600000.00/min
[570000], rate=34200000.00/min
[580000], rate=34800000.00/min
[590000], rate=35400000.00/min
[600000], rate=36000000.00/min
[610000], rate=36600000.00/min
[620000], rate=37200000.00/min
[630000], rate=37800000.00/min
[640000], rate=38400000.00/min
[650000], rate=39000000.00/min
[660000], rate=39600000.00/min
[670000], rate=40200000.00/min
[680000], rate=40800000.00/min
[690000], rate=41400000.00/min
[700000], rate=42000000.00/min
[710000], rate=42600000.00/min
[720000], rate=43200000.00/min
[730000], rate=21900000.00/min
[740000], rate=22200000.00/min
[750000], rate=22500000.00/min
[760000], rate=22800000.00/min
[770000], rate=23100000.00/min
[780000], rate=23400000.00/min
[790000], rate=23700000.00/min
[800000], rate=24000000.00/min
[810000], rate=24300000.00/min
[820000], rate=24600000.00/min
[830000], rate=24900000.00/min
[840000], rate=25200000.00/min
[850000], rate=25500000.00/min
[860000], rate=25800000.00/min
[870000], rate=26100000.00/min
[880000], rate=26400000.00/min
[890000], rate=26700000.00/min
[900000], rate=27000000.00/min
[910000], rate=27300000.00/min
[920000], rate=27600000.00/min
[930000], rate=27900000.00/min
[940000], rate=28200000.00/min
[950000], rate=28500000.00/min
[960000], rate=28800000.00/min
[970000], rate=29100000.00/min
[980000], rate=29400000.00/min
[990000], rate=29700000.00/min
[1000000], rate=30000000.00/min
[1010000], rate=30300000.00/min
[1020000], rate=30600000.00/min
[1030000], rate=30900000.00/min
[1040000], rate=31200000.00/min
[1050000], rate=31500000.00/min
[1060000], rate=31800000.00/min
[1070000], rate=32100000.00/min
[1080000], rate=32400000.00/min
[1090000], rate=32700000.00/min
[1100000], rate=22000000.00/min
[1110000], rate=22200000.00/min
[1120000], rate=22400000.00/min
[1130000], rate=22600000.00/min
[1140000], rate=22800000.00/min
[1150000], rate=23000000.00/min
[1160000], rate=23200000.00/min
[1170000], rate=23400000.00/min
[1180000], rate=23600000.00/min
[1190000], rate=23800000.00/min
[1200000], rate=24000000.00/min
[1210000], rate=24200000.00/min
[1220000], rate=24400000.00/min
[1230000], rate=24600000.00/min
[1240000], rate=24800000.00/min
[1250000], rate=25000000.00/min
[1260000], rate=25200000.00/min
[1270000], rate=25400000.00/min
[1280000], rate=25600000.00/min
[1290000], rate=25800000.00/min
[1300000], rate=26000000.00/min
[1310000], rate=26200000.00/min
[1320000], rate=26400000.00/min
[1330000], rate=26600000.00/min
[1340000], rate=26800000.00/min
[1350000], rate=27000000.00/min
[1360000], rate=27200000.00/min
[1370000], rate=27400000.00/min
[1380000], rate=27600000.00/min
[1390000], rate=27800000.00/min
[1400000], rate=28000000.00/min
[1410000], rate=28200000.00/min
[1420000], rate=28400000.00/min
[1430000], rate=28600000.00/min
[1440000], rate=28800000.00/min
[1450000], rate=29000000.00/min
[1460000], rate=29200000.00/min
[1470000], rate=22050000.00/min
[1480000], rate=22200000.00/min
[1490000], rate=22350000.00/min
[1500000], rate=22500000.00/min
[1510000], rate=22650000.00/min
[1520000], rate=22800000.00/min
[1530000], rate=22950000.00/min
[1540000], rate=23100000.00/min
[1550000], rate=23250000.00/min
[1560000], rate=23400000.00/min
[1570000], rate=23550000.00/min
[1580000], rate=23700000.00/min
[1590000], rate=23850000.00/min
[1600000], rate=24000000.00/min
[1610000], rate=24150000.00/min
[1620000], rate=24300000.00/min
[1630000], rate=24450000.00/min
[1640000], rate=24600000.00/min
[1650000], rate=24750000.00/min
[1660000], rate=24900000.00/min
[1670000], rate=25050000.00/min
[1680000], rate=25200000.00/min
[1690000], rate=25350000.00/min
[1700000], rate=25500000.00/min
[1710000], rate=25650000.00/min
[1720000], rate=25800000.00/min
[1730000], rate=25950000.00/min
[1740000], rate=26100000.00/min
[1750000], rate=26250000.00/min
[1760000], rate=26400000.00/min
[1770000], rate=26550000.00/min
[1780000], rate=26700000.00/min
[1790000], rate=26850000.00/min
[1800000], rate=27000000.00/min
[1810000], rate=27150000.00/min
[1820000], rate=27300000.00/min
[1830000], rate=21960000.00/min
[1840000], rate=22080000.00/min
[1850000], rate=22200000.00/min
[1860000], rate=22320000.00/min
[1870000], rate=22440000.00/min
[1880000], rate=22560000.00/min
[1890000], rate=22680000.00/min
[1900000], rate=22800000.00/min
[1910000], rate=22920000.00/min
[1920000], rate=23040000.00/min
[1930000], rate=23160000.00/min
[1940000], rate=23280000.00/min
[1950000], rate=23400000.00/min
[1960000], rate=23520000.00/min
[1970000], rate=23640000.00/min
[1980000], rate=23760000.00/min
[1990000], rate=23880000.00/min
[2000000], rate=24000000.00/min
[2010000], rate=24120000.00/min
[2020000], rate=24240000.00/min
[2030000], rate=24360000.00/min
[2040000], rate=24480000.00/min
[2050000], rate=24600000.00/min
[2060000], rate=24720000.00/min
[2070000], rate=24840000.00/min
[2080000], rate=24960000.00/min
[2090000], rate=25080000.00/min
[2100000], rate=25200000.00/min
[2110000], rate=25320000.00/min
[2120000], rate=21200000.00/min
[2130000], rate=21300000.00/min
[2140000], rate=21400000.00/min
[2150000], rate=21500000.00/min
[2160000], rate=21600000.00/min
[2170000], rate=21700000.00/min
[2180000], rate=21800000.00/min
[2190000], rate=21900000.00/min
[2200000], rate=22000000.00/min
[2210000], rate=22100000.00/min
[2220000], rate=22200000.00/min
[2230000], rate=22300000.00/min
[2240000], rate=22400000.00/min
[2250000], rate=22500000.00/min
[2260000], rate=22600000.00/min
[2270000], rate=22700000.00/min
[2280000], rate=22800000.00/min
[2290000], rate=22900000.00/min
[2300000], rate=23000000.00/min
[2310000], rate=23100000.00/min
[2320000], rate=23200000.00/min
[2330000], rate=23300000.00/min
[2340000], rate=23400000.00/min
[2350000], rate=23500000.00/min
[2360000], rate=23600000.00/min
[2370000], rate=23700000.00/min
[2380000], rate=23800000.00/min
[2390000], rate=23900000.00/min
[2400000], rate=24000000.00/min
[2410000], rate=24100000.00/min
[2420000], rate=24200000.00/min
[2430000], rate=24300000.00/min
[2440000], rate=24400000.00/min
[2450000], rate=24500000.00/min
[2460000], rate=24600000.00/min
[2470000], rate=24700000.00/min
[2480000], rate=21257142.86/min
[2490000], rate=21342857.14/min
[2500000], rate=21428571.43/min
[2510000], rate=21514285.71/min
[2520000], rate=21600000.00/min
[2530000], rate=21685714.29/min
[2540000], rate=21771428.57/min
[2550000], rate=21857142.86/min
[2560000], rate=21942857.14/min
[2570000], rate=22028571.43/min
[2580000], rate=22114285.71/min
[2590000], rate=22200000.00/min
[2600000], rate=22285714.29/min
[2610000], rate=22371428.57/min
[2620000], rate=22457142.86/min
[2630000], rate=22542857.14/min
[2640000], rate=22628571.43/min
[2650000], rate=22714285.71/min
[2660000], rate=22800000.00/min
[2670000], rate=22885714.29/min
[2680000], rate=22971428.57/min
[2690000], rate=23057142.86/min
[2700000], rate=23142857.14/min
[2710000], rate=23228571.43/min
[2720000], rate=23314285.71/min
[2730000], rate=23400000.00/min
[2740000], rate=23485714.29/min
[2750000], rate=23571428.57/min
[2760000], rate=23657142.86/min
[2770000], rate=23742857.14/min
[2780000], rate=23828571.43/min
[2790000], rate=23914285.71/min
[2800000], rate=24000000.00/min
[2810000], rate=21075000.00/min
[2820000], rate=21150000.00/min
[2830000], rate=21225000.00/min
[2840000], rate=21300000.00/min
[2850000], rate=21375000.00/min
[2860000], rate=21450000.00/min
[2870000], rate=21525000.00/min
[2880000], rate=21600000.00/min
[2890000], rate=21675000.00/min
[2900000], rate=21750000.00/min
[2910000], rate=21825000.00/min
[2920000], rate=21900000.00/min
[2930000], rate=21975000.00/min
[2940000], rate=22050000.00/min
[2950000], rate=22125000.00/min
[2960000], rate=22200000.00/min
[2970000], rate=22275000.00/min
[2980000], rate=22350000.00/min
[2990000], rate=22425000.00/min
[3000000], rate=22500000.00/min
[3010000], rate=22575000.00/min
[3020000], rate=22650000.00/min
[3030000], rate=22725000.00/min
[3040000], rate=22800000.00/min
[3050000], rate=22875000.00/min
[3060000], rate=22950000.00/min
[3070000], rate=23025000.00/min
[3080000], rate=23100000.00/min
[3090000], rate=23175000.00/min
[3100000], rate=23250000.00/min
[3110000], rate=23325000.00/min
[3120000], rate=23400000.00/min
[3130000], rate=23475000.00/min
[3140000], rate=23550000.00/min
[3150000], rate=23625000.00/min
[3160000], rate=23700000.00/min
[3170000], rate=23775000.00/min
[3180000], rate=21200000.00/min
[3190000], rate=21266666.67/min
[3200000], rate=21333333.33/min
[3210000], rate=21400000.00/min
[3220000], rate=21466666.67/min
[3230000], rate=21533333.33/min
[3240000], rate=21600000.00/min
[3250000], rate=21666666.67/min
[3260000], rate=21733333.33/min
[3270000], rate=21800000.00/min 00:05:34 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/annotate_cancer_introns.pl --introns_file 659_cAd-T1-TRNA-1_B23WHTKLT4_1.introns --ctat_genome_lib ctat_genome_lib_build_dir --intron_col 0 > 659_cAd-T1-TRNA-1_B23WHTKLT4_1.cancer.introns.prelim
Use of uninitialized value $ctat_genome_lib in concatenation (.) or string at /usr/local/src/CTAT-SPLICING/util/annotate_cancer_introns.pl line 18.
-659_cAd-T1-TRNA-1_B23WHTKLT4_1.introns: identified 3076 cancer introns
00:05:34 : INFO : Execution Time = 0.01 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/annotate_cancer_introns.pl --introns_file 659_cAd-T1-TRNA-1_B23WHTKLT4_1.introns --ctat_genome_lib ctat_genome_lib_build_dir --intron_col 0 > 659_cAd-T1-TRNA-1_B23WHTKLT4_1.cancer.introns.prelim
00:05:34 : INFO : Running: touch /tmp/nxf.0U1CdND7ED/659_cAd-T1-TRNA-1_B23WHTKLT4_1.chckpts/prelim_introns.ok
00:05:34 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/filter_by_min_total_reads.py --cancer_intron_candidates 659_cAd-T1-TRNA-1_B23WHTKLT4_1.cancer.introns.prelim --min_total_reads 5 > 659_cAd-T1-TRNA-1_B23WHTKLT4_1.cancer.introns
00:05:34 : INFO : Execution Time = 0.00 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/filter_by_min_total_reads.py --cancer_intron_candidates 659_cAd-T1-TRNA-1_B23WHTKLT4_1.cancer.introns.prelim --min_total_reads 5 > 659_cAd-T1-TRNA-1_B23WHTKLT4_1.cancer.introns
00:05:34 : INFO : Running: touch /tmp/nxf.0U1CdND7ED/659_cAd-T1-TRNA-1_B23WHTKLT4_1.chckpts/introns_filtered.ok
00:05:34 : INFO : -found 1162 cancer introns
00:05:34 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/make_igv_splice_bed.py --all_introns 659_cAd-T1-TRNA-1_B23WHTKLT4_1.introns --cancer_introns 659_cAd-T1-TRNA-1_B23WHTKLT4_1.cancer.introns --genome_lib_dir ctat_genome_lib_build_dir --output_bed 659_cAd-T1-TRNA-1_B23WHTKLT4_1.introns.for_IGV.bed
00:05:35 : INFO : Creating the BED File.
00:05:37 : INFO : Saving Bed File as 659_cAd-T1-TRNA-1_B23WHTKLT4_1.introns.for_IGV.bed
00:05:37 : INFO : Execution Time = 0.05 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/make_igv_splice_bed.py --all_introns 659_cAd-T1-TRNA-1_B23WHTKLT4_1.introns --cancer_introns 659_cAd-T1-TRNA-1_B23WHTKLT4_1.cancer.introns --genome_lib_dir ctat_genome_lib_build_dir --output_bed 659_cAd-T1-TRNA-1_B23WHTKLT4_1.introns.for_IGV.bed
00:05:37 : INFO : Running: touch /tmp/nxf.0U1CdND7ED/659_cAd-T1-TRNA-1_B23WHTKLT4_1.chckpts/intron_igv_bed.ok
00:05:37 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/igv_read_alignment_extractor.py --igv_introns_bed 659_cAd-T1-TRNA-1_B23WHTKLT4_1.introns.for_IGV.bed --bam 659_cAd-T1-TRNA-1_B23WHTKLT4_1.Aligned.sortedByCoord.out.bam --output_prefix 659_cAd-T1-TRNA-1_B23WHTKLT4_1
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TRNA-1_B23WHTKLT4_1.Aligned.sortedByCoord.out.bam.bai
00:08:16 : INFO : Execution Time = 2.65 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/igv_read_alignment_extractor.py --igv_introns_bed 659_cAd-T1-TRNA-1_B23WHTKLT4_1.introns.for_IGV.bed --bam 659_cAd-T1-TRNA-1_B23WHTKLT4_1.Aligned.sortedByCoord.out.bam --output_prefix 659_cAd-T1-TRNA-1_B23WHTKLT4_1
00:08:16 : INFO : Running: touch /tmp/nxf.0U1CdND7ED/659_cAd-T1-TRNA-1_B23WHTKLT4_1.chckpts/reads_alignments_extracted.ok
00:08:16 : INFO : Running: samtools sort -o 659_cAd-T1-TRNA-1_B23WHTKLT4_1.cancer_intron_reads.sorted.bam 659_cAd-T1-TRNA-1_B23WHTKLT4_1.cancer_intron_reads.bam
00:08:16 : INFO : Execution Time = 0.00 minutes. CMD: samtools sort -o 659_cAd-T1-TRNA-1_B23WHTKLT4_1.cancer_intron_reads.sorted.bam 659_cAd-T1-TRNA-1_B23WHTKLT4_1.cancer_intron_reads.bam
00:08:16 : INFO : Running: touch /tmp/nxf.0U1CdND7ED/659_cAd-T1-TRNA-1_B23WHTKLT4_1.chckpts/sort_cancer_intron_reads.ok
00:08:16 : INFO : Running: samtools sort -o 659_cAd-T1-TRNA-1_B23WHTKLT4_1.gene_reads.sorted.bam 659_cAd-T1-TRNA-1_B23WHTKLT4_1.gene_reads.bam
[bam_sort_core] merging from 15 files and 1 in-memory blocks...
00:11:52 : INFO : Execution Time = 3.60 minutes. CMD: samtools sort -o 659_cAd-T1-TRNA-1_B23WHTKLT4_1.gene_reads.sorted.bam 659_cAd-T1-TRNA-1_B23WHTKLT4_1.gene_reads.bam
00:11:52 : INFO : Running: touch /tmp/nxf.0U1CdND7ED/659_cAd-T1-TRNA-1_B23WHTKLT4_1.chckpts/sort_gene_reads.ok
00:11:52 : INFO : Running: /usr/local/src/CTAT-SPLICING/bamsifter/bamsifter -c 50 -i 50 -o 659_cAd-T1-TRNA-1_B23WHTKLT4_1.gene_reads.sorted.sifted.bam.tmp --keep_secondary 659_cAd-T1-TRNA-1_B23WHTKLT4_1.gene_reads.sorted.bam
Done with chr chr1.
Done with chr chr2.
Done with chr chr3.
Done with chr chr4.
Done with chr chr5.
Done with chr chr6.
Done with chr chr7.
Done with chr chr8.
Done with chr chr9.
Done with chr chr10.
Done with chr chr11.
Done with chr chr12.
Done with chr chr13.
Done with chr chr14.
Done with chr chr15.
Done with chr chr16.
Done with chr chr17.
Done with chr chr18.
Done with chr chr19.
Done with chr chr20.
Done with chr chr21.
Done with chr chr22.
00:12:58 : INFO : Execution Time = 1.11 minutes. CMD: /usr/local/src/CTAT-SPLICING/bamsifter/bamsifter -c 50 -i 50 -o 659_cAd-T1-TRNA-1_B23WHTKLT4_1.gene_reads.sorted.sifted.bam.tmp --keep_secondary 659_cAd-T1-TRNA-1_B23WHTKLT4_1.gene_reads.sorted.bam
00:12:58 : INFO : Running: touch /tmp/nxf.0U1CdND7ED/659_cAd-T1-TRNA-1_B23WHTKLT4_1.chckpts/659_cAd-T1-TRNA-1_B23WHTKLT4_1.gene_reads.sorted.sifted.bam.tmp.ok
00:12:58 : INFO : Running: samtools sort -o 659_cAd-T1-TRNA-1_B23WHTKLT4_1.gene_reads.sorted.sifted.bam 659_cAd-T1-TRNA-1_B23WHTKLT4_1.gene_reads.sorted.sifted.bam.tmp
00:13:09 : INFO : Execution Time = 0.18 minutes. CMD: samtools sort -o 659_cAd-T1-TRNA-1_B23WHTKLT4_1.gene_reads.sorted.sifted.bam 659_cAd-T1-TRNA-1_B23WHTKLT4_1.gene_reads.sorted.sifted.bam.tmp
00:13:09 : INFO : Running: touch /tmp/nxf.0U1CdND7ED/659_cAd-T1-TRNA-1_B23WHTKLT4_1.chckpts/659_cAd-T1-TRNA-1_B23WHTKLT4_1.gene_reads.sorted.sifted.bam.tmpsorted.ok
00:13:09 : INFO : Running: samtools index 659_cAd-T1-TRNA-1_B23WHTKLT4_1.gene_reads.sorted.sifted.bam
00:13:11 : INFO : Execution Time = 0.03 minutes. CMD: samtools index 659_cAd-T1-TRNA-1_B23WHTKLT4_1.gene_reads.sorted.sifted.bam
00:13:11 : INFO : Running: touch /tmp/nxf.0U1CdND7ED/659_cAd-T1-TRNA-1_B23WHTKLT4_1.chckpts/659_cAd-T1-TRNA-1_B23WHTKLT4_1.gene_reads.sorted.sifted.bam.indexed.ok
00:13:11 : INFO : Running: samtools index 659_cAd-T1-TRNA-1_B23WHTKLT4_1.cancer_intron_reads.sorted.bam
00:13:11 : INFO : Execution Time = 0.00 minutes. CMD: samtools index 659_cAd-T1-TRNA-1_B23WHTKLT4_1.cancer_intron_reads.sorted.bam
00:13:11 : INFO : Running: touch /tmp/nxf.0U1CdND7ED/659_cAd-T1-TRNA-1_B23WHTKLT4_1.chckpts/659_cAd-T1-TRNA-1_B23WHTKLT4_1.cancer_intron_reads.sorted.bam.indexed.ok
00:13:11 : INFO : Running: create_report 659_cAd-T1-TRNA-1_B23WHTKLT4_1.introns.for_IGV.bed ctat_genome_lib_build_dir/ref_genome.fa --type junction --output 659_cAd-T1-TRNA-1_B23WHTKLT4_1.ctat-splicing.igv.html --track-config 659_cAd-T1-TRNA-1_B23WHTKLT4_1.igv.tracks --info-columns gene variant_name uniquely_mapped multi_mapped TCGA GTEx --title 'CTAT_Splicing: 659_cAd-T1-TRNA-1_B23WHTKLT4_1'
00:14:36 : INFO : Execution Time = 1.42 minutes. CMD: create_report 659_cAd-T1-TRNA-1_B23WHTKLT4_1.introns.for_IGV.bed ctat_genome_lib_build_dir/ref_genome.fa --type junction --output 659_cAd-T1-TRNA-1_B23WHTKLT4_1.ctat-splicing.igv.html --track-config 659_cAd-T1-TRNA-1_B23WHTKLT4_1.igv.tracks --info-columns gene variant_name uniquely_mapped multi_mapped TCGA GTEx --title 'CTAT_Splicing: 659_cAd-T1-TRNA-1_B23WHTKLT4_1'
00:14:36 : INFO : Running: touch /tmp/nxf.0U1CdND7ED/659_cAd-T1-TRNA-1_B23WHTKLT4_1.chckpts/igv_create_html.ok
00:14:36 : INFO : done.