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/30/d54995457db7dd4243b60911fe5ff8/.command.sh
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/1d/d9766cfec01d2a4f1d657be37b8350/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
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/1d/d9766cfec01d2a4f1d657be37b8350/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bam
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/30/d54995457db7dd4243b60911fe5ff8/.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 aih-tih-sc-63b0e7-R1_B23WHYVLT4_1
+ all_qc -bam aih-tih-sc-63b0e7-R1_B23WHYVLT4_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 aih-tih-sc-63b0e7-R1_B23WHYVLT4_1 -reference ref_genome.fa -region chr20
mkdir -p aih-tih-sc-63b0e7-R1_B23WHYVLT4_1/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md/
cat Twist_rna_exome_AND_fusion_hg38.bed | cut -f 1-3 | bedtools sort -g ref_genome.fa.fai | bedtools merge > aih-tih-sc-63b0e7-R1_B23WHYVLT4_1/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md/full_bed.bed
seqtool depth -bam aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bam -bed aih-tih-sc-63b0e7-R1_B23WHYVLT4_1/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md/full_bed.bed -threads 4 -bam_threads 2 | bgzip -c > aih-tih-sc-63b0e7-R1_B23WHYVLT4_1/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.depth.gz
saturation_plot -bams aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bam -region chr20 -output aih-tih-sc-63b0e7-R1_B23WHYVLT4_1/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md/saturation_plot.png -threads 1
seqtool saturation_metric -at 1000000 -region chr20 -bam_threads 1 aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bam > aih-tih-sc-63b0e7-R1_B23WHYVLT4_1/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md/saturation_qc.tsv
fragment_size -bams aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bam -region chr20 -output_prefix aih-tih-sc-63b0e7-R1_B23WHYVLT4_1/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md/fragsize -threads 1
samtools mpileup -B --min-MQ 60 --min-BQ 30 -f ref_genome.fa -l dbsnp_hg38_contam_sites.bed aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bam > aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.mpileup && contamination -input aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.mpileup -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz > aih-tih-sc-63b0e7-R1_B23WHYVLT4_1/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md/contamination.tsv
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[mpileup] 1 samples in 1 input files
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_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: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
0it [00:00, ?it/s][W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
136684it [00:00, 1366792.42it/s]
273364it [00:00, 1323285.79it/s]
405788it [00:00, 1288926.13it/s]
569628it [00:00, 1423137.97it/s]
734366it [00:00, 1502584.40it/s]
900067it [00:00, 1554469.58it/s]
1068254it [00:00, 1595780.92it/s]
1228711it [00:00, 1598406.48it/s]
1396475it [00:00, 1622921.62it/s]
1565234it [00:01, 1642691.43it/s]
1729598it [00:01, 1640236.52it/s]
1896890it [00:01, 1650135.86it/s]
2064595it [00:01, 1658108.02it/s]
2235171it [00:01, 1672465.17it/s]
2408421it [00:01, 1690398.77it/s]
2577481it [00:01, 1687423.58it/s]
2746279it [00:01, 1687572.97it/s]
2915047it [00:01, 1682791.54it/s]
3084165it [00:01, 1685284.98it/s]
3252702it [00:02, 1634849.74it/s]
3419161it [00:02, 1643504.84it/s]
3586325it [00:02, 1651800.88it/s]
3756224it [00:02, 1665757.51it/s]
3924213it [00:02, 1669947.35it/s]
4092500it [00:02, 1673787.01it/s]
4262707it [00:02, 1682096.36it/s]
4430974it [00:02, 1682182.07it/s]
4599232it [00:02, 1677616.70it/s]
4771066it [00:02, 1689769.66it/s]
4940068it [00:03, 1676566.42it/s]
5107763it [00:03, 1666758.62it/s]
5278069it [00:03, 1677510.68it/s]
5449465it [00:03, 1688337.19it/s]
5622172it [00:03, 1699869.78it/s]
5796286it [00:03, 1712166.69it/s]
5967525it [00:03, 1711900.50it/s]
6139553it [00:03, 1714372.58it/s]
6312574it [00:03, 1718982.21it/s]
6486241it [00:03, 1724250.96it/s]
6658673it [00:04, 1718056.97it/s]
6830488it [00:04, 1702597.39it/s]
7000782it [00:04, 1666077.97it/s]
7169234it [00:04, 1671439.85it/s]
7339027it [00:04, 1679249.50it/s]
7507565it [00:04, 1681049.73it/s]
7677043it [00:04, 1685121.15it/s]
7850676it [00:04, 1700388.25it/s]
8022115it [00:04, 1704553.42it/s]
8192605it [00:04, 1700717.85it/s]
8362702it [00:05, 1692160.84it/s]
8531943it [00:05, 1690492.79it/s]
8701009it [00:05, 1686533.80it/s]
8869675it [00:05, 1683823.17it/s]
9038066it [00:05, 1674439.74it/s]
9205524it [00:05, 1673814.01it/s]
9376082it [00:05, 1683279.09it/s]
9545115it [00:05, 1685367.71it/s]
9718416it [00:05, 1699607.50it/s]
9891675it [00:05, 1709468.56it/s]
10062631it [00:06, 1682784.40it/s]
10233416it [00:06, 1690198.38it/s]
10402513it [00:06, 1688856.82it/s]
10571452it [00:06, 1671204.18it/s]
10738641it [00:06, 1668748.59it/s]
10905567it [00:06, 1668871.47it/s]
11075764it [00:06, 1678639.69it/s]
11243657it [00:06, 1667314.66it/s]
11410645it [00:06, 1668060.77it/s]
11577475it [00:06, 1662313.65it/s]
11743731it [00:07, 1662366.29it/s]
11910878it [00:07, 1664889.96it/s]
12079652it [00:07, 1671706.54it/s]
12248390it [00:07, 1676378.00it/s]
12416036it [00:07, 1673607.61it/s]
12584622it [00:07, 1677258.63it/s]
12752353it [00:07, 1676906.87it/s]
12821242it [00:07, 1664554.97it/s]
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.bai
tabix -s 1 -b 2 -e 2 aih-tih-sc-63b0e7-R1_B23WHYVLT4_1/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.depth.gz
seqtool depth_hist -nbins 500 aih-tih-sc-63b0e7-R1_B23WHYVLT4_1/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.depth.gz > aih-tih-sc-63b0e7-R1_B23WHYVLT4_1/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.depth_hist.tsv
cat aih-tih-sc-63b0e7-R1_B23WHYVLT4_1/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.depth_hist.tsv | plot -x 0 -xlabel Depth -ylabel "Number of Sites" -xlim 0 2000 -title aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md -output aih-tih-sc-63b0e7-R1_B23WHYVLT4_1/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md/aih-tih-sc-63b0e7-R1_B23WHYVLT4_1.md.depth.png