19:07:58 : INFO : -reading targets list: ctat_genome_lib_build_dir/ref_annot.gtf.mini.sortu
-mapping reads to genes
[130000], rate=7800000.00/min
[140000], rate=8400000.00/min
[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=24000000.00/min
[410000], rate=24600000.00/min
[420000], rate=25200000.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=20400000.00/min
[690000], rate=20700000.00/min
[700000], rate=21000000.00/min
[710000], rate=21300000.00/min
[720000], rate=21600000.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=18600000.00/min
[940000], rate=18800000.00/min
[950000], rate=19000000.00/min
[960000], rate=19200000.00/min
[970000], rate=19400000.00/min
[980000], rate=19600000.00/min
[990000], rate=19800000.00/min
[1000000], rate=20000000.00/min
[1010000], rate=20200000.00/min
[1020000], rate=20400000.00/min
[1030000], rate=20600000.00/min
[1040000], rate=20800000.00/min
[1050000], rate=21000000.00/min
[1060000], rate=21200000.00/min
[1070000], rate=21400000.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=17400000.00/min
[1170000], rate=17550000.00/min
[1180000], rate=17700000.00/min
[1190000], rate=17850000.00/min
[1200000], rate=18000000.00/min
[1210000], rate=18150000.00/min
[1220000], rate=18300000.00/min
[1230000], rate=18450000.00/min
[1240000], rate=18600000.00/min
[1250000], rate=18750000.00/min
[1260000], rate=18900000.00/min
[1270000], rate=19050000.00/min
[1280000], rate=19200000.00/min
[1290000], rate=19350000.00/min
[1300000], rate=19500000.00/min
[1310000], rate=19650000.00/min
[1320000], rate=19800000.00/min
[1330000], rate=19950000.00/min
[1340000], rate=20100000.00/min
[1350000], rate=20250000.00/min
[1360000], rate=20400000.00/min
[1370000], rate=20550000.00/min
[1380000], rate=20700000.00/min
[1390000], rate=20850000.00/min
[1400000], rate=21000000.00/min
[1410000], rate=21150000.00/min
[1420000], rate=21300000.00/min
[1430000], rate=21450000.00/min
[1440000], rate=21600000.00/min
[1450000], rate=21750000.00/min
[1460000], rate=21900000.00/min
[1470000], rate=22050000.00/min
[1480000], rate=22200000.00/min
[1490000], rate=17880000.00/min
[1500000], rate=18000000.00/min
[1510000], rate=18120000.00/min
[1520000], rate=18240000.00/min
[1530000], rate=18360000.00/min
[1540000], rate=18480000.00/min
[1550000], rate=18600000.00/min
[1560000], rate=18720000.00/min
[1570000], rate=18840000.00/min
[1580000], rate=18960000.00/min
[1590000], rate=19080000.00/min
[1600000], rate=19200000.00/min
[1610000], rate=19320000.00/min
[1620000], rate=19440000.00/min
[1630000], rate=19560000.00/min
[1640000], rate=19680000.00/min
[1650000], rate=19800000.00/min
[1660000], rate=19920000.00/min
[1670000], rate=20040000.00/min 19:08:07 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/annotate_cancer_introns.pl --introns_file aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.introns --ctat_genome_lib ctat_genome_lib_build_dir --intron_col 0 > aih-tih-sc-d5fda9-R1_A23YTGFLT4_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-d5fda9-R1_A23YTGFLT4_1.introns: identified 370 cancer introns
19:08:07 : INFO : Execution Time = 0.01 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/annotate_cancer_introns.pl --introns_file aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.introns --ctat_genome_lib ctat_genome_lib_build_dir --intron_col 0 > aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.cancer.introns.prelim
19:08:07 : INFO : Running: touch /tmp/nxf.MxGgINGdCU/aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.chckpts/prelim_introns.ok
19:08:07 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/filter_by_min_total_reads.py --cancer_intron_candidates aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.cancer.introns.prelim --min_total_reads 5 > aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.cancer.introns
19:08:07 : 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-d5fda9-R1_A23YTGFLT4_1.cancer.introns.prelim --min_total_reads 5 > aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.cancer.introns
19:08:07 : INFO : Running: touch /tmp/nxf.MxGgINGdCU/aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.chckpts/introns_filtered.ok
19:08:07 : INFO : -found 196 cancer introns
19:08:07 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/make_igv_splice_bed.py --all_introns aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.introns --cancer_introns aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.cancer.introns --genome_lib_dir ctat_genome_lib_build_dir --output_bed aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.introns.for_IGV.bed
19:08:08 : INFO : Creating the BED File.
19:08:09 : INFO : Saving Bed File as aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.introns.for_IGV.bed
19:08:09 : INFO : Execution Time = 0.02 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/make_igv_splice_bed.py --all_introns aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.introns --cancer_introns aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.cancer.introns --genome_lib_dir ctat_genome_lib_build_dir --output_bed aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.introns.for_IGV.bed
19:08:09 : INFO : Running: touch /tmp/nxf.MxGgINGdCU/aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.chckpts/intron_igv_bed.ok
19:08:09 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/igv_read_alignment_extractor.py --igv_introns_bed aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.introns.for_IGV.bed --bam aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.Aligned.sortedByCoord.out.bam --output_prefix aih-tih-sc-d5fda9-R1_A23YTGFLT4_1
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.Aligned.sortedByCoord.out.bam.bai
19:08:24 : INFO : Execution Time = 0.26 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/igv_read_alignment_extractor.py --igv_introns_bed aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.introns.for_IGV.bed --bam aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.Aligned.sortedByCoord.out.bam --output_prefix aih-tih-sc-d5fda9-R1_A23YTGFLT4_1
19:08:24 : INFO : Running: touch /tmp/nxf.MxGgINGdCU/aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.chckpts/reads_alignments_extracted.ok
19:08:24 : INFO : Running: samtools sort -o aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.cancer_intron_reads.sorted.bam aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.cancer_intron_reads.bam
19:08:24 : INFO : Execution Time = 0.00 minutes. CMD: samtools sort -o aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.cancer_intron_reads.sorted.bam aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.cancer_intron_reads.bam
19:08:24 : INFO : Running: touch /tmp/nxf.MxGgINGdCU/aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.chckpts/sort_cancer_intron_reads.ok
19:08:24 : INFO : Running: samtools sort -o aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.gene_reads.sorted.bam aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.gene_reads.bam
[bam_sort_core] merging from 1 files and 1 in-memory blocks...
19:08:42 : INFO : Execution Time = 0.29 minutes. CMD: samtools sort -o aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.gene_reads.sorted.bam aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.gene_reads.bam
19:08:42 : INFO : Running: touch /tmp/nxf.MxGgINGdCU/aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.chckpts/sort_gene_reads.ok
19:08:42 : INFO : Running: /usr/local/src/CTAT-SPLICING/bamsifter/bamsifter -c 50 -i 50 -o aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.gene_reads.sorted.sifted.bam.tmp --keep_secondary aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.gene_reads.sorted.bam
19:08:47 : INFO : Execution Time = 0.08 minutes. CMD: /usr/local/src/CTAT-SPLICING/bamsifter/bamsifter -c 50 -i 50 -o aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.gene_reads.sorted.sifted.bam.tmp --keep_secondary aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.gene_reads.sorted.bam
19:08:47 : INFO : Running: touch /tmp/nxf.MxGgINGdCU/aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.chckpts/aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.gene_reads.sorted.sifted.bam.tmp.ok
19:08:47 : INFO : Running: samtools sort -o aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.gene_reads.sorted.sifted.bam aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.gene_reads.sorted.sifted.bam.tmp
19:08:49 : INFO : Execution Time = 0.04 minutes. CMD: samtools sort -o aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.gene_reads.sorted.sifted.bam aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.gene_reads.sorted.sifted.bam.tmp
19:08:49 : INFO : Running: touch /tmp/nxf.MxGgINGdCU/aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.chckpts/aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.gene_reads.sorted.sifted.bam.tmpsorted.ok
19:08:49 : INFO : Running: samtools index aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.gene_reads.sorted.sifted.bam
19:08:49 : INFO : Execution Time = 0.00 minutes. CMD: samtools index aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.gene_reads.sorted.sifted.bam
19:08:49 : INFO : Running: touch /tmp/nxf.MxGgINGdCU/aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.chckpts/aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.gene_reads.sorted.sifted.bam.indexed.ok
19:08:49 : INFO : Running: samtools index aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.cancer_intron_reads.sorted.bam
19:08:49 : INFO : Execution Time = 0.00 minutes. CMD: samtools index aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.cancer_intron_reads.sorted.bam
19:08:49 : INFO : Running: touch /tmp/nxf.MxGgINGdCU/aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.chckpts/aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.cancer_intron_reads.sorted.bam.indexed.ok
19:08:49 : INFO : Running: create_report aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.introns.for_IGV.bed ctat_genome_lib_build_dir/ref_genome.fa --type junction --output aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.ctat-splicing.igv.html --track-config aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.igv.tracks --info-columns gene variant_name uniquely_mapped multi_mapped TCGA GTEx --title 'CTAT_Splicing: aih-tih-sc-d5fda9-R1_A23YTGFLT4_1'
19:09:12 : INFO : Execution Time = 0.37 minutes. CMD: create_report aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.introns.for_IGV.bed ctat_genome_lib_build_dir/ref_genome.fa --type junction --output aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.ctat-splicing.igv.html --track-config aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.igv.tracks --info-columns gene variant_name uniquely_mapped multi_mapped TCGA GTEx --title 'CTAT_Splicing: aih-tih-sc-d5fda9-R1_A23YTGFLT4_1'
19:09:12 : INFO : Running: touch /tmp/nxf.MxGgINGdCU/aih-tih-sc-d5fda9-R1_A23YTGFLT4_1.chckpts/igv_create_html.ok
19:09:12 : INFO : done.