15:56:52 : INFO : -reading targets list: ctat_genome_lib_build_dir/ref_annot.gtf.mini.sortu
-mapping reads to genes
[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 15:56:56 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/annotate_cancer_introns.pl --introns_file rna_s1221_S42.introns --ctat_genome_lib ctat_genome_lib_build_dir --intron_col 0 > rna_s1221_S42.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.
-rna_s1221_S42.introns: identified 1037 cancer introns
15:56:56 : INFO : Execution Time = 0.00 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/annotate_cancer_introns.pl --introns_file rna_s1221_S42.introns --ctat_genome_lib ctat_genome_lib_build_dir --intron_col 0 > rna_s1221_S42.cancer.introns.prelim
15:56:56 : INFO : Running: touch /tmp/nxf.7R9CaDrOLm/rna_s1221_S42.chckpts/prelim_introns.ok
15:56:56 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/filter_by_min_total_reads.py --cancer_intron_candidates rna_s1221_S42.cancer.introns.prelim --min_total_reads 5 > rna_s1221_S42.cancer.introns
15:56:57 : INFO : Execution Time = 0.00 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/filter_by_min_total_reads.py --cancer_intron_candidates rna_s1221_S42.cancer.introns.prelim --min_total_reads 5 > rna_s1221_S42.cancer.introns
15:56:57 : INFO : Running: touch /tmp/nxf.7R9CaDrOLm/rna_s1221_S42.chckpts/introns_filtered.ok
15:56:57 : INFO : -found 19 cancer introns
15:56:57 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/make_igv_splice_bed.py --all_introns rna_s1221_S42.introns --cancer_introns rna_s1221_S42.cancer.introns --genome_lib_dir ctat_genome_lib_build_dir --output_bed rna_s1221_S42.introns.for_IGV.bed
15:56:57 : INFO : Creating the BED File.
15:56:57 : INFO : Saving Bed File as rna_s1221_S42.introns.for_IGV.bed
15:56:57 : INFO : Execution Time = 0.01 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/make_igv_splice_bed.py --all_introns rna_s1221_S42.introns --cancer_introns rna_s1221_S42.cancer.introns --genome_lib_dir ctat_genome_lib_build_dir --output_bed rna_s1221_S42.introns.for_IGV.bed
15:56:57 : INFO : Running: touch /tmp/nxf.7R9CaDrOLm/rna_s1221_S42.chckpts/intron_igv_bed.ok
15:56:57 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/igv_read_alignment_extractor.py --igv_introns_bed rna_s1221_S42.introns.for_IGV.bed --bam rna_s1221_S42.Aligned.sortedByCoord.out.bam --output_prefix rna_s1221_S42
[W::hts_idx_load3] The index file is older than the data file: rna_s1221_S42.Aligned.sortedByCoord.out.bam.bai
15:56:58 : INFO : Execution Time = 0.01 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/igv_read_alignment_extractor.py --igv_introns_bed rna_s1221_S42.introns.for_IGV.bed --bam rna_s1221_S42.Aligned.sortedByCoord.out.bam --output_prefix rna_s1221_S42
15:56:58 : INFO : Running: touch /tmp/nxf.7R9CaDrOLm/rna_s1221_S42.chckpts/reads_alignments_extracted.ok
15:56:58 : INFO : Running: samtools sort -o rna_s1221_S42.cancer_intron_reads.sorted.bam rna_s1221_S42.cancer_intron_reads.bam
15:56:58 : INFO : Execution Time = 0.00 minutes. CMD: samtools sort -o rna_s1221_S42.cancer_intron_reads.sorted.bam rna_s1221_S42.cancer_intron_reads.bam
15:56:58 : INFO : Running: touch /tmp/nxf.7R9CaDrOLm/rna_s1221_S42.chckpts/sort_cancer_intron_reads.ok
15:56:58 : INFO : Running: samtools sort -o rna_s1221_S42.gene_reads.sorted.bam rna_s1221_S42.gene_reads.bam
15:56:58 : INFO : Execution Time = 0.01 minutes. CMD: samtools sort -o rna_s1221_S42.gene_reads.sorted.bam rna_s1221_S42.gene_reads.bam
15:56:58 : INFO : Running: touch /tmp/nxf.7R9CaDrOLm/rna_s1221_S42.chckpts/sort_gene_reads.ok
15:56:58 : INFO : Running: /usr/local/src/CTAT-SPLICING/bamsifter/bamsifter -c 50 -i 50 -o rna_s1221_S42.gene_reads.sorted.sifted.bam.tmp --keep_secondary rna_s1221_S42.gene_reads.sorted.bam
15:56:58 : INFO : Execution Time = 0.00 minutes. CMD: /usr/local/src/CTAT-SPLICING/bamsifter/bamsifter -c 50 -i 50 -o rna_s1221_S42.gene_reads.sorted.sifted.bam.tmp --keep_secondary rna_s1221_S42.gene_reads.sorted.bam
15:56:58 : INFO : Running: touch /tmp/nxf.7R9CaDrOLm/rna_s1221_S42.chckpts/rna_s1221_S42.gene_reads.sorted.sifted.bam.tmp.ok
15:56:58 : INFO : Running: samtools sort -o rna_s1221_S42.gene_reads.sorted.sifted.bam rna_s1221_S42.gene_reads.sorted.sifted.bam.tmp
15:56:58 : INFO : Execution Time = 0.00 minutes. CMD: samtools sort -o rna_s1221_S42.gene_reads.sorted.sifted.bam rna_s1221_S42.gene_reads.sorted.sifted.bam.tmp
15:56:58 : INFO : Running: touch /tmp/nxf.7R9CaDrOLm/rna_s1221_S42.chckpts/rna_s1221_S42.gene_reads.sorted.sifted.bam.tmpsorted.ok
15:56:58 : INFO : Running: samtools index rna_s1221_S42.gene_reads.sorted.sifted.bam
15:56:58 : INFO : Execution Time = 0.00 minutes. CMD: samtools index rna_s1221_S42.gene_reads.sorted.sifted.bam
15:56:58 : INFO : Running: touch /tmp/nxf.7R9CaDrOLm/rna_s1221_S42.chckpts/rna_s1221_S42.gene_reads.sorted.sifted.bam.indexed.ok
15:56:58 : INFO : Running: samtools index rna_s1221_S42.cancer_intron_reads.sorted.bam
15:56:58 : INFO : Execution Time = 0.00 minutes. CMD: samtools index rna_s1221_S42.cancer_intron_reads.sorted.bam
15:56:58 : INFO : Running: touch /tmp/nxf.7R9CaDrOLm/rna_s1221_S42.chckpts/rna_s1221_S42.cancer_intron_reads.sorted.bam.indexed.ok
15:56:58 : INFO : Running: create_report rna_s1221_S42.introns.for_IGV.bed ctat_genome_lib_build_dir/ref_genome.fa --type junction --output rna_s1221_S42.ctat-splicing.igv.html --track-config rna_s1221_S42.igv.tracks --info-columns gene variant_name uniquely_mapped multi_mapped TCGA GTEx --title 'CTAT_Splicing: rna_s1221_S42'
15:57:00 : INFO : Execution Time = 0.03 minutes. CMD: create_report rna_s1221_S42.introns.for_IGV.bed ctat_genome_lib_build_dir/ref_genome.fa --type junction --output rna_s1221_S42.ctat-splicing.igv.html --track-config rna_s1221_S42.igv.tracks --info-columns gene variant_name uniquely_mapped multi_mapped TCGA GTEx --title 'CTAT_Splicing: rna_s1221_S42'
15:57:00 : INFO : Running: touch /tmp/nxf.7R9CaDrOLm/rna_s1221_S42.chckpts/igv_create_html.ok
15:57:00 : INFO : done.