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/0a/fadb7d9f9f7a1a2279e4b732a6186f/aih-tih-sc-284856-R1_A23WJ53LT4_1.Aligned.sortedByCoord.out.bam
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/48/e7ef2d7b7cd2063be8c26f44d487c6/aih-tih-sc-284856-R1_A23WJ53LT4_1.Aligned.sortedByCoord.out.bam.bai
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/8b/341ef2e12bb90f4143405bc6c78b71/.command.sh
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/0a/fadb7d9f9f7a1a2279e4b732a6186f/aih-tih-sc-284856-R1_A23WJ53LT4_1.Chimeric.out.junction
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/8b/341ef2e12bb90f4143405bc6c78b71/.command.run
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/0a/fadb7d9f9f7a1a2279e4b732a6186f/aih-tih-sc-284856-R1_A23WJ53LT4_1.SJ.out.tab
==> STAGING COMPLETE (7 inputs)
08:04:45 : INFO : -reading targets list: ctat_genome_lib_build_dir/ref_annot.gtf.mini.sortu
-mapping reads to genes
[150000], rate=9000000.00/min
[160000], rate=9600000.00/min
[170000], rate=10200000.00/min
[180000], rate=10800000.00/min
[190000], rate=11400000.00/min
[200000], rate=12000000.00/min
[210000], rate=12600000.00/min
[220000], rate=13200000.00/min
[230000], rate=13800000.00/min
[240000], rate=14400000.00/min
[250000], rate=15000000.00/min
[260000], rate=15600000.00/min
[270000], rate=16200000.00/min
[280000], rate=16800000.00/min
[290000], rate=17400000.00/min
[300000], rate=18000000.00/min
[310000], rate=18600000.00/min
[320000], rate=19200000.00/min
[330000], rate=19800000.00/min
[340000], rate=20400000.00/min
[350000], rate=21000000.00/min
[360000], rate=21600000.00/min
[370000], rate=22200000.00/min
[380000], rate=22800000.00/min
[390000], rate=23400000.00/min
[400000], rate=12000000.00/min
[410000], rate=12300000.00/min
[420000], rate=12600000.00/min
[430000], rate=12900000.00/min
[440000], rate=13200000.00/min
[450000], rate=13500000.00/min
[460000], rate=13800000.00/min
[470000], rate=14100000.00/min
[480000], rate=14400000.00/min
[490000], rate=14700000.00/min
[500000], rate=15000000.00/min
[510000], rate=15300000.00/min
[520000], rate=15600000.00/min
[530000], rate=15900000.00/min
[540000], rate=16200000.00/min
[550000], rate=16500000.00/min
[560000], rate=16800000.00/min
[570000], rate=17100000.00/min
[580000], rate=17400000.00/min
[590000], rate=17700000.00/min
[600000], rate=18000000.00/min
[610000], rate=18300000.00/min
[620000], rate=18600000.00/min
[630000], rate=18900000.00/min
[640000], rate=19200000.00/min
[650000], rate=19500000.00/min
[660000], rate=19800000.00/min
[670000], rate=20100000.00/min
[680000], rate=13600000.00/min
[690000], rate=13800000.00/min
[700000], rate=14000000.00/min
[710000], rate=14200000.00/min
[720000], rate=14400000.00/min
[730000], rate=14600000.00/min
[740000], rate=14800000.00/min
[750000], rate=15000000.00/min
[760000], rate=15200000.00/min
[770000], rate=15400000.00/min
[780000], rate=15600000.00/min
[790000], rate=15800000.00/min
[800000], rate=16000000.00/min
[810000], rate=16200000.00/min
[820000], rate=16400000.00/min
[830000], rate=16600000.00/min
[840000], rate=16800000.00/min
[850000], rate=17000000.00/min
[860000], rate=17200000.00/min
[870000], rate=17400000.00/min
[880000], rate=17600000.00/min
[890000], rate=17800000.00/min
[900000], rate=18000000.00/min
[910000], rate=18200000.00/min
[920000], rate=18400000.00/min
[930000], rate=13950000.00/min
[940000], rate=14100000.00/min
[950000], rate=14250000.00/min
[960000], rate=14400000.00/min
[970000], rate=14550000.00/min
[980000], rate=14700000.00/min
[990000], rate=14850000.00/min
[1000000], rate=15000000.00/min
[1010000], rate=15150000.00/min
[1020000], rate=15300000.00/min
[1030000], rate=15450000.00/min
[1040000], rate=15600000.00/min
[1050000], rate=15750000.00/min
[1060000], rate=15900000.00/min
[1070000], rate=16050000.00/min
[1080000], rate=16200000.00/min
[1090000], rate=16350000.00/min
[1100000], rate=16500000.00/min
[1110000], rate=16650000.00/min
[1120000], rate=16800000.00/min
[1130000], rate=16950000.00/min
[1140000], rate=17100000.00/min
[1150000], rate=17250000.00/min
[1160000], rate=13920000.00/min
[1170000], rate=14040000.00/min
[1180000], rate=14160000.00/min
[1190000], rate=14280000.00/min
[1200000], rate=14400000.00/min
[1210000], rate=14520000.00/min
[1220000], rate=14640000.00/min
[1230000], rate=14760000.00/min
[1240000], rate=14880000.00/min
[1250000], rate=15000000.00/min
[1260000], rate=15120000.00/min
[1270000], rate=15240000.00/min
[1280000], rate=15360000.00/min
[1290000], rate=15480000.00/min
[1300000], rate=15600000.00/min
[1310000], rate=15720000.00/min
[1320000], rate=15840000.00/min
[1330000], rate=15960000.00/min
[1340000], rate=16080000.00/min
[1350000], rate=16200000.00/min
[1360000], rate=16320000.00/min
[1370000], rate=16440000.00/min
[1380000], rate=16560000.00/min
[1390000], rate=16680000.00/min
[1400000], rate=16800000.00/min
[1410000], rate=16920000.00/min
[1420000], rate=14200000.00/min
[1430000], rate=14300000.00/min
[1440000], rate=14400000.00/min
[1450000], rate=14500000.00/min
[1460000], rate=14600000.00/min
[1470000], rate=14700000.00/min
[1480000], rate=14800000.00/min
[1490000], rate=14900000.00/min
[1500000], rate=15000000.00/min
[1510000], rate=15100000.00/min
[1520000], rate=15200000.00/min
[1530000], rate=15300000.00/min
[1540000], rate=15400000.00/min
[1550000], rate=15500000.00/min
[1560000], rate=15600000.00/min
[1570000], rate=15700000.00/min
[1580000], rate=15800000.00/min
[1590000], rate=15900000.00/min
[1600000], rate=16000000.00/min
[1610000], rate=16100000.00/min
[1620000], rate=16200000.00/min
[1630000], rate=16300000.00/min
[1640000], rate=16400000.00/min
[1650000], rate=16500000.00/min
[1660000], rate=16600000.00/min
[1670000], rate=16700000.00/min
[1680000], rate=14400000.00/min
[1690000], rate=14485714.29/min
[1700000], rate=14571428.57/min
[1710000], rate=14657142.86/min
[1720000], rate=14742857.14/min
[1730000], rate=14828571.43/min
[1740000], rate=14914285.71/min
[1750000], rate=15000000.00/min
[1760000], rate=15085714.29/min
[1770000], rate=15171428.57/min
[1780000], rate=15257142.86/min
[1790000], rate=15342857.14/min
[1800000], rate=15428571.43/min
[1810000], rate=15514285.71/min
[1820000], rate=15600000.00/min
[1830000], rate=15685714.29/min
[1840000], rate=15771428.57/min
[1850000], rate=15857142.86/min
[1860000], rate=15942857.14/min
[1870000], rate=16028571.43/min
[1880000], rate=16114285.71/min
[1890000], rate=16200000.00/min
[1900000], rate=16285714.29/min
[1910000], rate=16371428.57/min
[1920000], rate=16457142.86/min
[1930000], rate=16542857.14/min
[1940000], rate=14550000.00/min
[1950000], rate=14625000.00/min
[1960000], rate=14700000.00/min
[1970000], rate=14775000.00/min
[1980000], rate=14850000.00/min
[1990000], rate=14925000.00/min
[2000000], rate=15000000.00/min
[2010000], rate=15075000.00/min
[2020000], rate=15150000.00/min
[2030000], rate=15225000.00/min
[2040000], rate=15300000.00/min
[2050000], rate=15375000.00/min
[2060000], rate=15450000.00/min
[2070000], rate=15525000.00/min
[2080000], rate=15600000.00/min
[2090000], rate=15675000.00/min
[2100000], rate=15750000.00/min
[2110000], rate=15825000.00/min
[2120000], rate=15900000.00/min
[2130000], rate=15975000.00/min
[2140000], rate=16050000.00/min 08:05:02 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/annotate_cancer_introns.pl --introns_file aih-tih-sc-284856-R1_A23WJ53LT4_1.introns --ctat_genome_lib ctat_genome_lib_build_dir --intron_col 0 > aih-tih-sc-284856-R1_A23WJ53LT4_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.
-aih-tih-sc-284856-R1_A23WJ53LT4_1.introns: identified 3989 cancer introns
08:05:02 : INFO : Execution Time = 0.01 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/annotate_cancer_introns.pl --introns_file aih-tih-sc-284856-R1_A23WJ53LT4_1.introns --ctat_genome_lib ctat_genome_lib_build_dir --intron_col 0 > aih-tih-sc-284856-R1_A23WJ53LT4_1.cancer.introns.prelim
08:05:02 : INFO : Running: touch /tmp/nxf.FRRH4IG0ac/aih-tih-sc-284856-R1_A23WJ53LT4_1.chckpts/prelim_introns.ok
08:05:02 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/filter_by_min_total_reads.py --cancer_intron_candidates aih-tih-sc-284856-R1_A23WJ53LT4_1.cancer.introns.prelim --min_total_reads 5 > aih-tih-sc-284856-R1_A23WJ53LT4_1.cancer.introns
08:05:02 : INFO : Execution Time = 0.01 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/filter_by_min_total_reads.py --cancer_intron_candidates aih-tih-sc-284856-R1_A23WJ53LT4_1.cancer.introns.prelim --min_total_reads 5 > aih-tih-sc-284856-R1_A23WJ53LT4_1.cancer.introns
08:05:02 : INFO : Running: touch /tmp/nxf.FRRH4IG0ac/aih-tih-sc-284856-R1_A23WJ53LT4_1.chckpts/introns_filtered.ok
08:05:03 : INFO : -found 651 cancer introns
08:05:03 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/make_igv_splice_bed.py --all_introns aih-tih-sc-284856-R1_A23WJ53LT4_1.introns --cancer_introns aih-tih-sc-284856-R1_A23WJ53LT4_1.cancer.introns --genome_lib_dir ctat_genome_lib_build_dir --output_bed aih-tih-sc-284856-R1_A23WJ53LT4_1.introns.for_IGV.bed
08:05:03 : INFO : Creating the BED File.
08:05:05 : INFO : Saving Bed File as aih-tih-sc-284856-R1_A23WJ53LT4_1.introns.for_IGV.bed
08:05:05 : INFO : Execution Time = 0.04 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/make_igv_splice_bed.py --all_introns aih-tih-sc-284856-R1_A23WJ53LT4_1.introns --cancer_introns aih-tih-sc-284856-R1_A23WJ53LT4_1.cancer.introns --genome_lib_dir ctat_genome_lib_build_dir --output_bed aih-tih-sc-284856-R1_A23WJ53LT4_1.introns.for_IGV.bed
08:05:05 : INFO : Running: touch /tmp/nxf.FRRH4IG0ac/aih-tih-sc-284856-R1_A23WJ53LT4_1.chckpts/intron_igv_bed.ok
08:05:05 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/igv_read_alignment_extractor.py --igv_introns_bed aih-tih-sc-284856-R1_A23WJ53LT4_1.introns.for_IGV.bed --bam aih-tih-sc-284856-R1_A23WJ53LT4_1.Aligned.sortedByCoord.out.bam --output_prefix aih-tih-sc-284856-R1_A23WJ53LT4_1
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-284856-R1_A23WJ53LT4_1.Aligned.sortedByCoord.out.bam.bai
08:08:04 : INFO : Execution Time = 2.99 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/igv_read_alignment_extractor.py --igv_introns_bed aih-tih-sc-284856-R1_A23WJ53LT4_1.introns.for_IGV.bed --bam aih-tih-sc-284856-R1_A23WJ53LT4_1.Aligned.sortedByCoord.out.bam --output_prefix aih-tih-sc-284856-R1_A23WJ53LT4_1
08:08:04 : INFO : Running: touch /tmp/nxf.FRRH4IG0ac/aih-tih-sc-284856-R1_A23WJ53LT4_1.chckpts/reads_alignments_extracted.ok
08:08:04 : INFO : Running: samtools sort -o aih-tih-sc-284856-R1_A23WJ53LT4_1.cancer_intron_reads.sorted.bam aih-tih-sc-284856-R1_A23WJ53LT4_1.cancer_intron_reads.bam
08:08:05 : INFO : Execution Time = 0.00 minutes. CMD: samtools sort -o aih-tih-sc-284856-R1_A23WJ53LT4_1.cancer_intron_reads.sorted.bam aih-tih-sc-284856-R1_A23WJ53LT4_1.cancer_intron_reads.bam
08:08:05 : INFO : Running: touch /tmp/nxf.FRRH4IG0ac/aih-tih-sc-284856-R1_A23WJ53LT4_1.chckpts/sort_cancer_intron_reads.ok
08:08:05 : INFO : Running: samtools sort -o aih-tih-sc-284856-R1_A23WJ53LT4_1.gene_reads.sorted.bam aih-tih-sc-284856-R1_A23WJ53LT4_1.gene_reads.bam
[bam_sort_core] merging from 13 files and 1 in-memory blocks...
08:11:23 : INFO : Execution Time = 3.31 minutes. CMD: samtools sort -o aih-tih-sc-284856-R1_A23WJ53LT4_1.gene_reads.sorted.bam aih-tih-sc-284856-R1_A23WJ53LT4_1.gene_reads.bam
08:11:23 : INFO : Running: touch /tmp/nxf.FRRH4IG0ac/aih-tih-sc-284856-R1_A23WJ53LT4_1.chckpts/sort_gene_reads.ok
08:11:23 : INFO : Running: /usr/local/src/CTAT-SPLICING/bamsifter/bamsifter -c 50 -i 50 -o aih-tih-sc-284856-R1_A23WJ53LT4_1.gene_reads.sorted.sifted.bam.tmp --keep_secondary aih-tih-sc-284856-R1_A23WJ53LT4_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.
08:12:36 : INFO : Execution Time = 1.21 minutes. CMD: /usr/local/src/CTAT-SPLICING/bamsifter/bamsifter -c 50 -i 50 -o aih-tih-sc-284856-R1_A23WJ53LT4_1.gene_reads.sorted.sifted.bam.tmp --keep_secondary aih-tih-sc-284856-R1_A23WJ53LT4_1.gene_reads.sorted.bam
08:12:36 : INFO : Running: touch /tmp/nxf.FRRH4IG0ac/aih-tih-sc-284856-R1_A23WJ53LT4_1.chckpts/aih-tih-sc-284856-R1_A23WJ53LT4_1.gene_reads.sorted.sifted.bam.tmp.ok
08:12:36 : INFO : Running: samtools sort -o aih-tih-sc-284856-R1_A23WJ53LT4_1.gene_reads.sorted.sifted.bam aih-tih-sc-284856-R1_A23WJ53LT4_1.gene_reads.sorted.sifted.bam.tmp
[bam_sort_core] merging from 1 files and 1 in-memory blocks...
08:12:55 : INFO : Execution Time = 0.32 minutes. CMD: samtools sort -o aih-tih-sc-284856-R1_A23WJ53LT4_1.gene_reads.sorted.sifted.bam aih-tih-sc-284856-R1_A23WJ53LT4_1.gene_reads.sorted.sifted.bam.tmp
08:12:55 : INFO : Running: touch /tmp/nxf.FRRH4IG0ac/aih-tih-sc-284856-R1_A23WJ53LT4_1.chckpts/aih-tih-sc-284856-R1_A23WJ53LT4_1.gene_reads.sorted.sifted.bam.tmpsorted.ok
08:12:55 : INFO : Running: samtools index aih-tih-sc-284856-R1_A23WJ53LT4_1.gene_reads.sorted.sifted.bam
08:12:57 : INFO : Execution Time = 0.03 minutes. CMD: samtools index aih-tih-sc-284856-R1_A23WJ53LT4_1.gene_reads.sorted.sifted.bam
08:12:57 : INFO : Running: touch /tmp/nxf.FRRH4IG0ac/aih-tih-sc-284856-R1_A23WJ53LT4_1.chckpts/aih-tih-sc-284856-R1_A23WJ53LT4_1.gene_reads.sorted.sifted.bam.indexed.ok
08:12:57 : INFO : Running: samtools index aih-tih-sc-284856-R1_A23WJ53LT4_1.cancer_intron_reads.sorted.bam
08:12:58 : INFO : Execution Time = 0.00 minutes. CMD: samtools index aih-tih-sc-284856-R1_A23WJ53LT4_1.cancer_intron_reads.sorted.bam
08:12:58 : INFO : Running: touch /tmp/nxf.FRRH4IG0ac/aih-tih-sc-284856-R1_A23WJ53LT4_1.chckpts/aih-tih-sc-284856-R1_A23WJ53LT4_1.cancer_intron_reads.sorted.bam.indexed.ok
08:12:58 : INFO : Running: create_report aih-tih-sc-284856-R1_A23WJ53LT4_1.introns.for_IGV.bed ctat_genome_lib_build_dir/ref_genome.fa --type junction --output aih-tih-sc-284856-R1_A23WJ53LT4_1.ctat-splicing.igv.html --track-config aih-tih-sc-284856-R1_A23WJ53LT4_1.igv.tracks --info-columns gene variant_name uniquely_mapped multi_mapped TCGA GTEx --title 'CTAT_Splicing: aih-tih-sc-284856-R1_A23WJ53LT4_1'
08:16:30 : INFO : Execution Time = 3.53 minutes. CMD: create_report aih-tih-sc-284856-R1_A23WJ53LT4_1.introns.for_IGV.bed ctat_genome_lib_build_dir/ref_genome.fa --type junction --output aih-tih-sc-284856-R1_A23WJ53LT4_1.ctat-splicing.igv.html --track-config aih-tih-sc-284856-R1_A23WJ53LT4_1.igv.tracks --info-columns gene variant_name uniquely_mapped multi_mapped TCGA GTEx --title 'CTAT_Splicing: aih-tih-sc-284856-R1_A23WJ53LT4_1'
08:16:30 : INFO : Running: touch /tmp/nxf.FRRH4IG0ac/aih-tih-sc-284856-R1_A23WJ53LT4_1.chckpts/igv_create_html.ok
08:16:30 : INFO : done.