Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/17/22687b8c29ceaae0e87261a73c5910/659_Kh-T1-TRNA-1_B23MHV2LT4_1.SJ.out.tab 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/17/22687b8c29ceaae0e87261a73c5910/659_Kh-T1-TRNA-1_B23MHV2LT4_1.Aligned.sortedByCoord.out.bam Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/b6/2888402fd5e42177eb04267a8665fc/.command.sh Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/bf/f44acd1442ba55f7900514ebeca01e/659_Kh-T1-TRNA-1_B23MHV2LT4_1.Aligned.sortedByCoord.out.bam.bai Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/17/22687b8c29ceaae0e87261a73c5910/659_Kh-T1-TRNA-1_B23MHV2LT4_1.Chimeric.out.junction Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/b6/2888402fd5e42177eb04267a8665fc/.command.run ==> STAGING COMPLETE (7 inputs) 07:20:31 : INFO : -reading targets list: ctat_genome_lib_build_dir/ref_annot.gtf.mini.sortu -mapping reads to genes [120000], rate=7200000.00/min [130000], rate=7800000.00/min [140000], rate=8400000.00/min [150000], rate=9000000.00/min [160000], rate=9600000.00/min 07:20:36 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/annotate_cancer_introns.pl --introns_file 659_Kh-T1-TRNA-1_B23MHV2LT4_1.introns --ctat_genome_lib ctat_genome_lib_build_dir --intron_col 0 > 659_Kh-T1-TRNA-1_B23MHV2LT4_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_Kh-T1-TRNA-1_B23MHV2LT4_1.introns: identified 81 cancer introns 07:20:36 : INFO : Execution Time = 0.00 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/annotate_cancer_introns.pl --introns_file 659_Kh-T1-TRNA-1_B23MHV2LT4_1.introns --ctat_genome_lib ctat_genome_lib_build_dir --intron_col 0 > 659_Kh-T1-TRNA-1_B23MHV2LT4_1.cancer.introns.prelim 07:20:36 : INFO : Running: touch /tmp/nxf.MOs3xSBgLI/659_Kh-T1-TRNA-1_B23MHV2LT4_1.chckpts/prelim_introns.ok 07:20:36 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/filter_by_min_total_reads.py --cancer_intron_candidates 659_Kh-T1-TRNA-1_B23MHV2LT4_1.cancer.introns.prelim --min_total_reads 5 > 659_Kh-T1-TRNA-1_B23MHV2LT4_1.cancer.introns 07:20:36 : INFO : Execution Time = 0.01 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/filter_by_min_total_reads.py --cancer_intron_candidates 659_Kh-T1-TRNA-1_B23MHV2LT4_1.cancer.introns.prelim --min_total_reads 5 > 659_Kh-T1-TRNA-1_B23MHV2LT4_1.cancer.introns 07:20:36 : INFO : Running: touch /tmp/nxf.MOs3xSBgLI/659_Kh-T1-TRNA-1_B23MHV2LT4_1.chckpts/introns_filtered.ok 07:20:36 : INFO : -found 42 cancer introns 07:20:36 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/make_igv_splice_bed.py --all_introns 659_Kh-T1-TRNA-1_B23MHV2LT4_1.introns --cancer_introns 659_Kh-T1-TRNA-1_B23MHV2LT4_1.cancer.introns --genome_lib_dir ctat_genome_lib_build_dir --output_bed 659_Kh-T1-TRNA-1_B23MHV2LT4_1.introns.for_IGV.bed 07:20:37 : INFO : Creating the BED File. 07:20:37 : INFO : Saving Bed File as 659_Kh-T1-TRNA-1_B23MHV2LT4_1.introns.for_IGV.bed 07:20:37 : INFO : Execution Time = 0.01 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/make_igv_splice_bed.py --all_introns 659_Kh-T1-TRNA-1_B23MHV2LT4_1.introns --cancer_introns 659_Kh-T1-TRNA-1_B23MHV2LT4_1.cancer.introns --genome_lib_dir ctat_genome_lib_build_dir --output_bed 659_Kh-T1-TRNA-1_B23MHV2LT4_1.introns.for_IGV.bed 07:20:37 : INFO : Running: touch /tmp/nxf.MOs3xSBgLI/659_Kh-T1-TRNA-1_B23MHV2LT4_1.chckpts/intron_igv_bed.ok 07:20:37 : INFO : Running: /usr/local/src/CTAT-SPLICING/util/igv_read_alignment_extractor.py --igv_introns_bed 659_Kh-T1-TRNA-1_B23MHV2LT4_1.introns.for_IGV.bed --bam 659_Kh-T1-TRNA-1_B23MHV2LT4_1.Aligned.sortedByCoord.out.bam --output_prefix 659_Kh-T1-TRNA-1_B23MHV2LT4_1 [W::hts_idx_load3] The index file is older than the data file: 659_Kh-T1-TRNA-1_B23MHV2LT4_1.Aligned.sortedByCoord.out.bam.bai 07:20:38 : INFO : Execution Time = 0.01 minutes. CMD: /usr/local/src/CTAT-SPLICING/util/igv_read_alignment_extractor.py --igv_introns_bed 659_Kh-T1-TRNA-1_B23MHV2LT4_1.introns.for_IGV.bed --bam 659_Kh-T1-TRNA-1_B23MHV2LT4_1.Aligned.sortedByCoord.out.bam --output_prefix 659_Kh-T1-TRNA-1_B23MHV2LT4_1 07:20:38 : INFO : Running: touch /tmp/nxf.MOs3xSBgLI/659_Kh-T1-TRNA-1_B23MHV2LT4_1.chckpts/reads_alignments_extracted.ok 07:20:38 : INFO : Running: samtools sort -o 659_Kh-T1-TRNA-1_B23MHV2LT4_1.cancer_intron_reads.sorted.bam 659_Kh-T1-TRNA-1_B23MHV2LT4_1.cancer_intron_reads.bam 07:20:38 : INFO : Execution Time = 0.00 minutes. CMD: samtools sort -o 659_Kh-T1-TRNA-1_B23MHV2LT4_1.cancer_intron_reads.sorted.bam 659_Kh-T1-TRNA-1_B23MHV2LT4_1.cancer_intron_reads.bam 07:20:38 : INFO : Running: touch /tmp/nxf.MOs3xSBgLI/659_Kh-T1-TRNA-1_B23MHV2LT4_1.chckpts/sort_cancer_intron_reads.ok 07:20:38 : INFO : Running: samtools sort -o 659_Kh-T1-TRNA-1_B23MHV2LT4_1.gene_reads.sorted.bam 659_Kh-T1-TRNA-1_B23MHV2LT4_1.gene_reads.bam 07:20:39 : INFO : Execution Time = 0.01 minutes. CMD: samtools sort -o 659_Kh-T1-TRNA-1_B23MHV2LT4_1.gene_reads.sorted.bam 659_Kh-T1-TRNA-1_B23MHV2LT4_1.gene_reads.bam 07:20:39 : INFO : Running: touch /tmp/nxf.MOs3xSBgLI/659_Kh-T1-TRNA-1_B23MHV2LT4_1.chckpts/sort_gene_reads.ok 07:20:39 : INFO : Running: /usr/local/src/CTAT-SPLICING/bamsifter/bamsifter -c 50 -i 50 -o 659_Kh-T1-TRNA-1_B23MHV2LT4_1.gene_reads.sorted.sifted.bam.tmp --keep_secondary 659_Kh-T1-TRNA-1_B23MHV2LT4_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 chr11. Done with chr chr14. Done with chr chr16. Done with chr chr17. Done with chr chr19. Done with chr chr21. Done with chr chr22. 07:20:39 : INFO : Execution Time = 0.00 minutes. CMD: /usr/local/src/CTAT-SPLICING/bamsifter/bamsifter -c 50 -i 50 -o 659_Kh-T1-TRNA-1_B23MHV2LT4_1.gene_reads.sorted.sifted.bam.tmp --keep_secondary 659_Kh-T1-TRNA-1_B23MHV2LT4_1.gene_reads.sorted.bam 07:20:39 : INFO : Running: touch /tmp/nxf.MOs3xSBgLI/659_Kh-T1-TRNA-1_B23MHV2LT4_1.chckpts/659_Kh-T1-TRNA-1_B23MHV2LT4_1.gene_reads.sorted.sifted.bam.tmp.ok 07:20:39 : INFO : Running: samtools sort -o 659_Kh-T1-TRNA-1_B23MHV2LT4_1.gene_reads.sorted.sifted.bam 659_Kh-T1-TRNA-1_B23MHV2LT4_1.gene_reads.sorted.sifted.bam.tmp 07:20:39 : INFO : Execution Time = 0.00 minutes. CMD: samtools sort -o 659_Kh-T1-TRNA-1_B23MHV2LT4_1.gene_reads.sorted.sifted.bam 659_Kh-T1-TRNA-1_B23MHV2LT4_1.gene_reads.sorted.sifted.bam.tmp 07:20:39 : INFO : Running: touch /tmp/nxf.MOs3xSBgLI/659_Kh-T1-TRNA-1_B23MHV2LT4_1.chckpts/659_Kh-T1-TRNA-1_B23MHV2LT4_1.gene_reads.sorted.sifted.bam.tmpsorted.ok 07:20:39 : INFO : Running: samtools index 659_Kh-T1-TRNA-1_B23MHV2LT4_1.gene_reads.sorted.sifted.bam 07:20:39 : INFO : Execution Time = 0.00 minutes. CMD: samtools index 659_Kh-T1-TRNA-1_B23MHV2LT4_1.gene_reads.sorted.sifted.bam 07:20:39 : INFO : Running: touch /tmp/nxf.MOs3xSBgLI/659_Kh-T1-TRNA-1_B23MHV2LT4_1.chckpts/659_Kh-T1-TRNA-1_B23MHV2LT4_1.gene_reads.sorted.sifted.bam.indexed.ok 07:20:39 : INFO : Running: samtools index 659_Kh-T1-TRNA-1_B23MHV2LT4_1.cancer_intron_reads.sorted.bam 07:20:39 : INFO : Execution Time = 0.00 minutes. CMD: samtools index 659_Kh-T1-TRNA-1_B23MHV2LT4_1.cancer_intron_reads.sorted.bam 07:20:39 : INFO : Running: touch /tmp/nxf.MOs3xSBgLI/659_Kh-T1-TRNA-1_B23MHV2LT4_1.chckpts/659_Kh-T1-TRNA-1_B23MHV2LT4_1.cancer_intron_reads.sorted.bam.indexed.ok 07:20:39 : INFO : Running: create_report 659_Kh-T1-TRNA-1_B23MHV2LT4_1.introns.for_IGV.bed ctat_genome_lib_build_dir/ref_genome.fa --type junction --output 659_Kh-T1-TRNA-1_B23MHV2LT4_1.ctat-splicing.igv.html --track-config 659_Kh-T1-TRNA-1_B23MHV2LT4_1.igv.tracks --info-columns gene variant_name uniquely_mapped multi_mapped TCGA GTEx --title 'CTAT_Splicing: 659_Kh-T1-TRNA-1_B23MHV2LT4_1' 07:20:47 : INFO : Execution Time = 0.13 minutes. CMD: create_report 659_Kh-T1-TRNA-1_B23MHV2LT4_1.introns.for_IGV.bed ctat_genome_lib_build_dir/ref_genome.fa --type junction --output 659_Kh-T1-TRNA-1_B23MHV2LT4_1.ctat-splicing.igv.html --track-config 659_Kh-T1-TRNA-1_B23MHV2LT4_1.igv.tracks --info-columns gene variant_name uniquely_mapped multi_mapped TCGA GTEx --title 'CTAT_Splicing: 659_Kh-T1-TRNA-1_B23MHV2LT4_1' 07:20:47 : INFO : Running: touch /tmp/nxf.MOs3xSBgLI/659_Kh-T1-TRNA-1_B23MHV2LT4_1.chckpts/igv_create_html.ok 07:20:47 : INFO : done.