Downloading: s3://natera-platform-sandbox/pipeline-resources/AIH/rna/GRCh38/starfusion/ctat_genome_lib_build_dir/ref_genome.fa
Downloading: s3://natera-platform-sandbox/pipeline-resources/Homo_sapiens_assembly38.dbsnp138.vcf.gz
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/7a/bb2491d34a95ade35b70e747e454d4/.command.sh
Downloading: s3://natera-platform-sandbox/pipeline-resources/AIH/rna/GRCh38/starfusion/ctat_genome_lib_build_dir/ref_genome.dict
Downloading: s3://natera-platform-sandbox/pipeline-resources/AIH/rna/GRCh38/gatk4/dbsnp_hg38_contam_sites.bed
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/9d/3e2fd006bfc5d1781c96aadaed5795/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/7a/bb2491d34a95ade35b70e747e454d4/.command.run
Downloading: s3://natera-platform-sandbox/pipeline-resources/AIH/rna/GRCh38/starfusion/ctat_genome_lib_build_dir/ref_genome.fa.fai
Downloading: s3://natera-platform-sandbox/pipeline-resources/beds/twist/Twist_rna_exome_AND_fusion_hg38.bed
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/9d/3e2fd006bfc5d1781c96aadaed5795/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bam
Downloading: s3://natera-platform-sandbox/pipeline-resources/Homo_sapiens_assembly38.dbsnp138.vcf.gz.tbi
==> STAGING COMPLETE (11 inputs)
+ mkdir PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1
+ all_qc -bam PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bam -beds Twist_rna_exome_AND_fusion_hg38.bed -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz -contamination_bed dbsnp_hg38_contam_sites.bed -output_folder PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1 -reference ref_genome.fa -region chr20
mkdir -p PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/
cat Twist_rna_exome_AND_fusion_hg38.bed | cut -f 1-3 | bedtools sort -g ref_genome.fa.fai | bedtools merge > PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/full_bed.bed
seqtool depth -bam PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bam -bed PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/full_bed.bed -threads 4 -bam_threads 2 | bgzip -c > PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.depth.gz
saturation_plot -bams PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bam -region chr20 -output PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/saturation_plot.png -threads 1
seqtool saturation_metric -at 1000000 -region chr20 -bam_threads 1 PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bam > PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/saturation_qc.tsv
fragment_size -bams PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bam -region chr20 -output_prefix PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/fragsize -threads 1
samtools mpileup -B --min-MQ 60 --min-BQ 30 -f ref_genome.fa -l dbsnp_hg38_contam_sites.bed PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bam > PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.mpileup && contamination -input PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.mpileup -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz > PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/contamination.tsv
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[mpileup] 1 samples in 1 input files
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
/opt/conda/lib/python3.9/site-packages/pyranges/__init__.py:8: UserWarning: pkg_resources is deprecated as an API. See https://setuptools.pypa.io/en/latest/pkg_resources.html. The pkg_resources package is slated for removal as early as 2025-11-30. Refrain from using this package or pin to Setuptools<81.
import pkg_resources
/opt/conda/lib/python3.9/site-packages/pyranges/__init__.py:8: UserWarning: pkg_resources is deprecated as an API. See https://setuptools.pypa.io/en/latest/pkg_resources.html. The pkg_resources package is slated for removal as early as 2025-11-30. Refrain from using this package or pin to Setuptools<81.
import pkg_resources
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
0it [00:00, ?it/s][W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
102439it [00:00, 1024285.97it/s]
210100it [00:00, 1055024.17it/s]
321877it [00:00, 1083642.62it/s]
436699it [00:00, 1109036.15it/s]
551758it [00:00, 1123997.65it/s]
667800it [00:00, 1136358.68it/s]
781436it [00:00, 1132977.41it/s]
895459it [00:00, 1135204.47it/s]
1008982it [00:00, 1134716.21it/s]
1122851it [00:01, 1135922.51it/s]
1236564it [00:01, 1136272.17it/s]
1353285it [00:01, 1145655.51it/s]
1468644it [00:01, 1148035.97it/s]
1583449it [00:01, 1137938.23it/s]
1697264it [00:01, 1127833.65it/s]
1810542it [00:01, 1129290.81it/s]
1923493it [00:01, 1100542.37it/s]
2033705it [00:01, 1095588.35it/s]
2148355it [00:01, 1110557.00it/s]
2265073it [00:02, 1127297.89it/s]
2383180it [00:02, 1143263.62it/s]
2503372it [00:02, 1160728.28it/s]
2626353it [00:02, 1181341.56it/s]
2748343it [00:02, 1192855.24it/s]
2869592it [00:02, 1198711.30it/s]
2989497it [00:02, 1185570.91it/s]
3108109it [00:02, 1183650.18it/s]
3226511it [00:02, 1180148.80it/s]
3344552it [00:02, 1088703.36it/s]
3454797it [00:03, 1022834.31it/s]
3558599it [00:03, 962546.47it/s]
3656281it [00:03, 914901.30it/s]
3748882it [00:03, 911610.21it/s]
3840786it [00:03, 909867.38it/s]
3933370it [00:03, 914287.49it/s]
4041471it [00:03, 961941.99it/s]
4157226it [00:03, 1018681.23it/s]
4272933it [00:03, 1059218.30it/s]
4390669it [00:04, 1094061.56it/s]
4507334it [00:04, 1115566.49it/s]
4624036it [00:04, 1130861.44it/s]
4739928it [00:04, 1139211.83it/s]
4855419it [00:04, 1143883.96it/s]
4970388it [00:04, 1145600.75it/s]
5085029it [00:04, 1101376.83it/s]
5195583it [00:04, 1077600.02it/s]
5303692it [00:04, 1046812.88it/s]
5408718it [00:04, 1045092.88it/s]
5513462it [00:05, 1042710.52it/s]
5618521it [00:05, 1044988.92it/s]
5732358it [00:05, 1072462.82it/s]
5850323it [00:05, 1104178.68it/s]
5965783it [00:05, 1119045.75it/s]
6080774it [00:05, 1128222.09it/s]
6194146it [00:05, 1129799.01it/s]
6307184it [00:05, 1129127.84it/s]
6422703it [00:05, 1136902.87it/s]
6539593it [00:05, 1146448.06it/s]
6654262it [00:06, 1084385.31it/s]
6763407it [00:06, 1005958.73it/s]
6865438it [00:06, 976721.17it/s]
6964113it [00:06, 977921.29it/s]
7062617it [00:06, 929643.93it/s]
7157196it [00:06, 934098.80it/s]
7272784it [00:06, 996839.09it/s]
7389626it [00:06, 1046213.43it/s]
7504789it [00:06, 1076899.64it/s]
7620758it [00:07, 1101133.52it/s]
7737479it [00:07, 1120658.55it/s]
7849906it [00:07, 1118735.17it/s]
7962032it [00:07, 1093548.78it/s]
8075972it [00:07, 1106891.31it/s]
8193312it [00:07, 1126504.40it/s]
8306156it [00:07, 1122334.49it/s]
8419298it [00:07, 1124965.78it/s][W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
8536545it [00:07, 1139089.39it/s]
8653059it [00:07, 1146845.08it/s]
8768287it [00:08, 1148447.65it/s]
8883174it [00:08, 1147151.27it/s]
8997919it [00:08, 1137094.95it/s]
9113363it [00:08, 1142235.52it/s]
9227618it [00:08, 1137698.34it/s]
9341412it [00:08, 1129911.09it/s]
9454886it [00:08, 1131326.90it/s]
9570069it [00:08, 1137384.81it/s]
9683825it [00:08, 1136613.79it/s]
9797498it [00:08, 1131019.18it/s]
9912909it [00:09, 1137881.37it/s]
10027666it [00:09, 1140648.95it/s]
10143718it [00:09, 1146494.71it/s]
10262075it [00:09, 1157574.49it/s]
10377842it [00:09, 1148174.14it/s]
10498151it [00:09, 1164514.52it/s]
10616932it [00:09, 1171453.29it/s]
10734100it [00:09, 1170741.94it/s]
10851426it [00:09, 1171477.17it/s]
10968585it [00:09, 1171201.12it/s]
11085713it [00:10, 1169608.04it/s]
11202680it [00:10, 1164188.25it/s]
11319849it [00:10, 1166410.24it/s]
11438505it [00:10, 1172417.93it/s]
11557582it [00:10, 1177888.37it/s]
11675722it [00:10, 1178920.25it/s]
11793619it [00:10, 1168221.68it/s]
11910466it [00:10, 1155684.79it/s]
12026074it [00:10, 1148653.66it/s]
12140969it [00:10, 1148515.54it/s]
12144340it [00:10, 1106565.03it/s]
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
tabix -s 1 -b 2 -e 2 PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.depth.gz
seqtool depth_hist -nbins 500 PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.depth.gz > PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.depth_hist.tsv
cat PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.depth_hist.tsv | plot -x 0 -xlabel Depth -ylabel "Number of Sites" -xlim 0 2000 -title PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md -output PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.depth.png
[W::hts_idx_load3] The index file is older than the data file: Homo_sapiens_assembly38.dbsnp138.vcf.gz.tbi
seqtool read_qc -bed Twist_rna_exome_AND_fusion_hg38.bed -n_threads 1 -bam_threads 2 PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bam > PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/read_qc.Twist_rna_exome_AND_fusion_hg38.bed.tsv
seqtool reads_per_target -bed Twist_rna_exome_AND_fusion_hg38.bed -threads 1 -bam_threads 2 PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bam > PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/reads_per_target.Twist_rna_exome_AND_fusion_hg38.bed.tsv
seqtool depth_stats -bed Twist_rna_exome_AND_fusion_hg38.bed PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.depth.gz> PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/depth_stats.Twist_rna_exome_AND_fusion_hg38.bed.tsv
seqtool gc_bias -reference ref_genome.fa -bed Twist_rna_exome_AND_fusion_hg38.bed -bin_size 0.05 -threads 1 -bam_threads 2 PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bam > PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/gc_bias_qc.Twist_rna_exome_AND_fusion_hg38.bed.tsv
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md.bai
cat PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/gc_bias_qc.Twist_rna_exome_AND_fusion_hg38.bed.tsv | plot -x 0 -header infer -ylim 0 2 -title "PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md" -suptitle "Twist_rna_exome_AND_fusion_hg38.bed" -output "PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/gc_bias.Twist_rna_exome_AND_fusion_hg38.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
echo `date` > PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md/qc_complete.txt
+ merge_qc -folders PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1.md -output PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1_full_basic_stats.tsv
+ head -n1 PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1_full_basic_stats.tsv
++ basename Twist_rna_exome_AND_fusion_hg38.bed
+ grep -e Twist_rna_exome_AND_fusion_hg38.bed PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1/PX_FFPE_pan_TRNA_nat_03_A23MVY2LT4_1_full_basic_stats.tsv