Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/8f/2e9263d06d21050652c9fc2b0879ff/A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
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/4d/1d30c22abf4bf029871d09917e911c/.command.sh
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/8f/2e9263d06d21050652c9fc2b0879ff/A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bam
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/4d/1d30c22abf4bf029871d09917e911c/.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-platform-sandbox/pipeline-resources/Homo_sapiens_assembly38.dbsnp138.vcf.gz.tbi
==> STAGING COMPLETE (11 inputs)
+ mkdir A673_FFPE_RNA_0001_B23LG7FLT4_1
+ all_qc -bam A673_FFPE_RNA_0001_B23LG7FLT4_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 A673_FFPE_RNA_0001_B23LG7FLT4_1 -reference ref_genome.fa -region chr20
mkdir -p A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/
cat Twist_rna_exome_AND_fusion_hg38.bed | cut -f 1-3 | bedtools sort -g ref_genome.fa.fai | bedtools merge > A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/full_bed.bed
seqtool depth -bam A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bam -bed A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/full_bed.bed -threads 4 -bam_threads 2 | bgzip -c > A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/A673_FFPE_RNA_0001_B23LG7FLT4_1.md.depth.gz
saturation_plot -bams A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bam -region chr20 -output A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/saturation_plot.png -threads 1
seqtool saturation_metric -at 1000000 -region chr20 -bam_threads 1 A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bam > A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/saturation_qc.tsv
fragment_size -bams A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bam -region chr20 -output_prefix A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/fragsize -threads 1
samtools mpileup -B --min-MQ 60 --min-BQ 30 -f ref_genome.fa -l dbsnp_hg38_contam_sites.bed A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bam > A673_FFPE_RNA_0001_B23LG7FLT4_1.md.mpileup && contamination -input A673_FFPE_RNA_0001_B23LG7FLT4_1.md.mpileup -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz > A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/contamination.tsv
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[mpileup] 1 samples in 1 input files
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_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: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
0it [00:00, ?it/s]
83941it [00:00, 839318.75it/s]
171038it [00:00, 857900.69it/s]
258576it [00:00, 865851.01it/s]
345162it [00:00, 857366.51it/s]
436816it [00:00, 878508.76it/s]
524689it [00:00, 867182.54it/s]
611445it [00:00, 801665.39it/s]
695044it [00:00, 811911.41it/s]
784116it [00:00, 835479.84it/s]
874861it [00:01, 857008.87it/s]
964995it [00:01, 870275.67it/s]
1053559it [00:01, 874861.56it/s]
1142948it [00:01, 880526.18it/s]
1232099it [00:01, 883797.41it/s]
1320603it [00:01, 862538.83it/s]
1407049it [00:01, 832783.23it/s]
1490634it [00:01, 832061.83it/s]
1574052it [00:01, 823174.48it/s]
1659077it [00:01, 831056.24it/s]
1747146it [00:02, 845606.88it/s]
1831824it [00:02, 809438.42it/s]
1913140it [00:02, 796646.72it/s]
2004992it [00:02, 831721.15it/s]
2096466it [00:02, 855896.59it/s]
2185046it [00:02, 864666.03it/s]
2274669it [00:02, 873978.24it/s]
2363485it [00:02, 878171.83it/s]
2455376it [00:02, 890289.23it/s]
2544502it [00:02, 887400.54it/s]
2633311it [00:03, 879097.26it/s]
2725829it [00:03, 892750.85it/s]
2815163it [00:03, 886744.67it/s]
2906045it [00:03, 893278.96it/s]
2995412it [00:03, 890815.39it/s]
3084521it [00:03, 888949.15it/s]
3173435it [00:03, 862097.42it/s]
3264759it [00:03, 876981.17it/s]
3354457it [00:03, 882855.77it/s]
3442867it [00:04, 876098.72it/s]
3531396it [00:04, 878794.16it/s]
3619629it [00:04, 879823.01it/s]
3707660it [00:04, 841658.61it/s]
3795703it [00:04, 852856.89it/s]
3881286it [00:04, 837900.65it/s]
3965311it [00:04, 791109.24it/s]
4045023it [00:04, 778267.36it/s]
4123248it [00:04, 776180.25it/s]
4206272it [00:04, 791681.95it/s]
4288342it [00:05, 800104.55it/s]
4369804it [00:05, 804345.70it/s]
4452342it [00:05, 810542.43it/s]
4541262it [00:05, 833872.74it/s]
4629242it [00:05, 847529.19it/s]
4719585it [00:05, 864192.25it/s]
4809733it [00:05, 875308.98it/s]
4897317it [00:05, 861007.04it/s]
4995536it [00:05, 896857.77it/s]
5089404it [00:05, 909273.65it/s]
5180420it [00:06, 894836.35it/s]
5270009it [00:06, 877392.25it/s]
5363822it [00:06, 895135.85it/s]
5453470it [00:06, 891128.38it/s]
5543134it [00:06, 892735.70it/s]
5632474it [00:06, 881993.74it/s]
5720739it [00:06, 854482.87it/s]
5809698it [00:06, 864645.77it/s]
5897532it [00:06, 868639.62it/s]
5984526it [00:06, 854407.92it/s]
6076423it [00:07, 873299.85it/s]
6163885it [00:07, 864580.56it/s]
6250442it [00:07, 855927.31it/s]
6336109it [00:07, 853889.39it/s]
6423841it [00:07, 860784.79it/s]
6509965it [00:07, 856838.63it/s]
6595681it [00:07, 851692.23it/s]
6687070it [00:07, 870107.04it/s]
6778581it [00:07, 883475.41it/s]
6869811it [00:08, 892055.36it/s]
6960357it [00:08, 896041.98it/s]
7053112it [00:08, 905446.23it/s]
7143677it [00:08, 897274.88it/s]
7233434it [00:08, 889628.17it/s]
7322427it [00:08, 851487.90it/s]
7407910it [00:08, 801001.94it/s]
7488723it [00:08, 766320.45it/s]
7565988it [00:08, 758762.28it/s]
7652616it [00:08, 788856.99it/s]
7740864it [00:09, 815592.60it/s]
7829222it [00:09, 835202.55it/s]
7920462it [00:09, 857764.76it/s]
8008949it [00:09, 865739.63it/s]
8099481it [00:09, 877452.73it/s]
8187416it [00:09, 850694.64it/s]
8272780it [00:09, 847353.90it/s]
8359338it [00:09, 852679.91it/s]
8444759it [00:09, 851704.22it/s]
8535894it [00:09, 869345.80it/s]
8622926it [00:10, 862206.86it/s]
8709223it [00:10, 852551.05it/s]
8794545it [00:10, 846566.46it/s]
8879249it [00:10, 819492.89it/s]
8966911it [00:10, 836021.70it/s]
9050701it [00:10, 835475.58it/s]
9137095it [00:10, 843840.97it/s]
9221586it [00:10, 829260.33it/s]
9309522it [00:10, 843954.58it/s]
9398708it [00:11, 858096.84it/s]
9484618it [00:11, 857557.34it/s]
9570444it [00:11, 842449.75it/s]
9654787it [00:11, 788601.60it/s]
9734348it [00:11, 762004.80it/s]
9817376it [00:11, 781107.41it/s]
9896017it [00:11, 782503.94it/s]
9988338it [00:11, 823267.08it/s]
10072869it [00:11, 829703.87it/s]
10156141it [00:11, 809721.66it/s]
10237405it [00:12, 772428.03it/s]
10315134it [00:12, 747120.25it/s]
10390257it [00:12, 733060.12it/s]
10471923it [00:12, 756665.31it/s]
10553976it [00:12, 775002.31it/s]
10640537it [00:12, 801367.48it/s]
10724645it [00:12, 813012.96it/s]
10810478it [00:12, 826394.16it/s]
10900126it [00:12, 847106.33it/s]
10989659it [00:12, 861400.28it/s]
11075915it [00:13, 846156.01it/s]
11166401it [00:13, 863458.81it/s]
11253564it [00:13, 865864.42it/s]
11340240it [00:13, 865556.39it/s]
11430244it [00:13, 875815.05it/s]
11520059it [00:13, 882470.42it/s]
11608344it [00:13, 867022.79it/s]
11695130it [00:13, 833945.36it/s]
11778816it [00:13, 827750.28it/s]
11865934it [00:14, 840344.15it/s]
11950146it [00:14, 837498.47it/s]
12034016it [00:14, 836697.39it/s]
12127017it [00:14, 864138.12it/s]
12217318it [00:14, 875660.24it/s]
12306184it [00:14, 879511.41it/s]
12394194it [00:14, 878464.51it/s]
12482081it [00:14, 877368.26it/s]
12570984it [00:14, 880831.94it/s]
12661160it [00:14, 887070.27it/s]
12749884it [00:15, 881621.94it/s]
12838065it [00:15, 872874.59it/s]
12925381it [00:15, 850668.03it/s]
13014946it [00:15, 863691.60it/s]
13102200it [00:15, 866279.75it/s]
13188917it [00:15, 857374.13it/s]
13276364it [00:15, 862396.50it/s]
13362663it [00:15, 846009.13it/s]
13451815it [00:15, 859353.57it/s]
13537846it [00:15, 855840.35it/s]
13625162it [00:16, 860946.69it/s]
13712618it [00:16, 864976.62it/s]
13799157it [00:16, 860622.99it/s]
13892289it [00:16, 881631.02it/s]
13980492it [00:16, 855919.64it/s]
14071064it [00:16, 870445.38it/s]
14158427it [00:16, 871368.35it/s]
14247493it [00:16, 877067.34it/s]
14337643it [00:16, 884314.32it/s]
14426143it [00:16, 868598.65it/s]
14515243it [00:17, 875188.36it/s]
14550647it [00:17, 850315.20it/s]
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
tabix -s 1 -b 2 -e 2 A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/A673_FFPE_RNA_0001_B23LG7FLT4_1.md.depth.gz
seqtool depth_hist -nbins 500 A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/A673_FFPE_RNA_0001_B23LG7FLT4_1.md.depth.gz > A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/A673_FFPE_RNA_0001_B23LG7FLT4_1.md.depth_hist.tsv
cat A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/A673_FFPE_RNA_0001_B23LG7FLT4_1.md.depth_hist.tsv | plot -x 0 -xlabel Depth -ylabel "Number of Sites" -xlim 0 2000 -title A673_FFPE_RNA_0001_B23LG7FLT4_1.md -output A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/A673_FFPE_RNA_0001_B23LG7FLT4_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 A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bam > A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_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 A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bam > A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/reads_per_target.Twist_rna_exome_AND_fusion_hg38.bed.tsv
seqtool depth_stats -bed Twist_rna_exome_AND_fusion_hg38.bed A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/A673_FFPE_RNA_0001_B23LG7FLT4_1.md.depth.gz> A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_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 A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bam > A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_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: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: A673_FFPE_RNA_0001_B23LG7FLT4_1.md.bai
cat A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/gc_bias_qc.Twist_rna_exome_AND_fusion_hg38.bed.tsv | plot -x 0 -header infer -ylim 0 2 -title "A673_FFPE_RNA_0001_B23LG7FLT4_1.md" -suptitle "Twist_rna_exome_AND_fusion_hg38.bed" -output "A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/gc_bias.Twist_rna_exome_AND_fusion_hg38.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
echo `date` > A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md/qc_complete.txt
+ merge_qc -folders A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1.md -output A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1_full_basic_stats.tsv
+ head -n1 A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1_full_basic_stats.tsv
++ basename Twist_rna_exome_AND_fusion_hg38.bed
+ grep -e Twist_rna_exome_AND_fusion_hg38.bed A673_FFPE_RNA_0001_B23LG7FLT4_1/A673_FFPE_RNA_0001_B23LG7FLT4_1_full_basic_stats.tsv