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/7b/e3f71d99425d0892d9612c9501a9e1/.command.sh
Downloading: s3://natera-platform-sandbox/pipeline-resources/AIH/rna/GRCh38/starfusion/ctat_genome_lib_build_dir/ref_genome.dict
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/c9/3770e3d36cea2c0b907cc6affb19ff/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
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/7b/e3f71d99425d0892d9612c9501a9e1/.command.run
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/c9/3770e3d36cea2c0b907cc6affb19ff/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bam
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-platform-sandbox/pipeline-resources/Homo_sapiens_assembly38.dbsnp138.vcf.gz.tbi
==> STAGING COMPLETE (11 inputs)
+ mkdir BVT_FFPE_TRNA_bst_02_50_50ng_50ng
+ all_qc -bam BVT_FFPE_TRNA_bst_02_50_50ng_50ng.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 BVT_FFPE_TRNA_bst_02_50_50ng_50ng -reference ref_genome.fa -region chr20
mkdir -p BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/
cat Twist_rna_exome_AND_fusion_hg38.bed | cut -f 1-3 | bedtools sort -g ref_genome.fa.fai | bedtools merge > BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/full_bed.bed
seqtool depth -bam BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bam -bed BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/full_bed.bed -threads 4 -bam_threads 2 | bgzip -c > BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.depth.gz
saturation_plot -bams BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bam -region chr20 -output BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/saturation_plot.png -threads 1
seqtool saturation_metric -at 1000000 -region chr20 -bam_threads 1 BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bam > BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/saturation_qc.tsv
fragment_size -bams BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bam -region chr20 -output_prefix BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/fragsize -threads 1
samtools mpileup -B --min-MQ 60 --min-BQ 30 -f ref_genome.fa -l dbsnp_hg38_contam_sites.bed BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bam > BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.mpileup && contamination -input BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.mpileup -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz > BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/contamination.tsv
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[mpileup] 1 samples in 1 input files
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.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: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
0it [00:00, ?it/s]
110584it [00:00, 1105727.70it/s]
231204it [00:00, 1161177.90it/s]
352708it [00:00, 1185546.25it/s]
474079it [00:00, 1196623.51it/s]
595365it [00:00, 1202439.44it/s]
717144it [00:00, 1207627.01it/s]
837908it [00:00, 1207076.73it/s]
958888it [00:00, 1207913.33it/s]
1083786it [00:00, 1220722.97it/s]
1205917it [00:01, 1220736.53it/s]
1330217it [00:01, 1227468.05it/s]
1452965it [00:01, 1223455.95it/s]
1575314it [00:01, 1223029.53it/s]
1697620it [00:01, 1218042.08it/s]
1819430it [00:01, 1215290.70it/s]
1971699it [00:01, 1307369.58it/s]
2139501it [00:01, 1418442.50it/s]
2309796it [00:01, 1503704.23it/s]
2480051it [00:01, 1563239.63it/s]
2636430it [00:02, 1416218.41it/s]
2780725it [00:02, 1345225.27it/s]
2917544it [00:02, 1313341.55it/s]
3050425it [00:02, 1288032.45it/s]
3180239it [00:02, 1272577.39it/s]
3308140it [00:02, 1258952.99it/s]
3434439it [00:02, 1251329.60it/s]
3559826it [00:02, 1246672.29it/s]
3684656it [00:02, 1243288.86it/s]
3809090it [00:02, 1231873.84it/s]
3932347it [00:03, 1216516.00it/s]
4054056it [00:03, 1209435.94it/s]
4175031it [00:03, 1201672.17it/s]
4295215it [00:03, 1198653.80it/s]
4415293it [00:03, 1199183.81it/s]
4535218it [00:03, 1197769.06it/s]
4657038it [00:03, 1203821.44it/s]
4777719it [00:03, 1204687.08it/s]
4898194it [00:03, 1202048.31it/s]
4955740it [00:03, 1252905.72it/s]
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
tabix -s 1 -b 2 -e 2 BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.depth.gz
seqtool depth_hist -nbins 500 BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.depth.gz > BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.depth_hist.tsv
cat BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.depth_hist.tsv | plot -x 0 -xlabel Depth -ylabel "Number of Sites" -xlim 0 2000 -title BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md -output BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.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 BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bam > BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.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 BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bam > BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/reads_per_target.Twist_rna_exome_AND_fusion_hg38.bed.tsv
seqtool depth_stats -bed Twist_rna_exome_AND_fusion_hg38.bed BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.depth.gz> BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.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 BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bam > BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.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: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
[W::hts_idx_load3] The index file is older than the data file: BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md.bai
cat BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/gc_bias_qc.Twist_rna_exome_AND_fusion_hg38.bed.tsv | plot -x 0 -header infer -ylim 0 2 -title "BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md" -suptitle "Twist_rna_exome_AND_fusion_hg38.bed" -output "BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/gc_bias.Twist_rna_exome_AND_fusion_hg38.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
echo `date` > BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md/qc_complete.txt
+ merge_qc -folders BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng.md -output BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng_full_basic_stats.tsv
+ head -n1 BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng_full_basic_stats.tsv
++ basename Twist_rna_exome_AND_fusion_hg38.bed
+ grep -e Twist_rna_exome_AND_fusion_hg38.bed BVT_FFPE_TRNA_bst_02_50_50ng_50ng/BVT_FFPE_TRNA_bst_02_50_50ng_50ng_full_basic_stats.tsv