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/22/31d4b87124f5c7f3e56c4068cf99c0/.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/e3/ccac06bd0f800d8a043625e6dfa095/aih-tih-sc-445ab1-R1_B23WHYVLT4_1.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/22/31d4b87124f5c7f3e56c4068cf99c0/.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-rnd-pltf-dev-nextflow-scratch-01/work/e3/ccac06bd0f800d8a043625e6dfa095/aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bam
  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-445ab1-R1_B23WHYVLT4_1
+ all_qc -bam aih-tih-sc-445ab1-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-445ab1-R1_B23WHYVLT4_1 -reference ref_genome.fa -region chr20
mkdir -p aih-tih-sc-445ab1-R1_B23WHYVLT4_1/aih-tih-sc-445ab1-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-445ab1-R1_B23WHYVLT4_1/aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md/full_bed.bed


seqtool depth -bam aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bam -bed aih-tih-sc-445ab1-R1_B23WHYVLT4_1/aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md/full_bed.bed -threads 4 -bam_threads 2 | bgzip -c > aih-tih-sc-445ab1-R1_B23WHYVLT4_1/aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md/aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.depth.gz
saturation_plot       -bams aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bam       -region chr20       -output aih-tih-sc-445ab1-R1_B23WHYVLT4_1/aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md/saturation_plot.png       -threads 1 
seqtool saturation_metric       -at 1000000       -region chr20       -bam_threads 1       aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bam > aih-tih-sc-445ab1-R1_B23WHYVLT4_1/aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md/saturation_qc.tsv 
fragment_size       -bams aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bam       -region chr20       -output_prefix aih-tih-sc-445ab1-R1_B23WHYVLT4_1/aih-tih-sc-445ab1-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-445ab1-R1_B23WHYVLT4_1.md.bam > aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.mpileup && contamination -input aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.mpileup -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz > aih-tih-sc-445ab1-R1_B23WHYVLT4_1/aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md/contamination.tsv 
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-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-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-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-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai

0it [00:00, ?it/s]
144207it [00:00, 1441971.68it/s]
304578it [00:00, 1537077.72it/s]
469879it [00:00, 1589988.17it/s]
641029it [00:00, 1637935.94it/s]
809989it [00:00, 1656543.69it/s]
976506it [00:00, 1659454.24it/s]
1142452it [00:00, 1646095.37it/s]
1311086it [00:00, 1658818.34it/s]
1483709it [00:00, 1679793.78it/s]
1656400it [00:01, 1694238.14it/s]
1826800it [00:01, 1697074.86it/s]
1996522it [00:01, 1694231.95it/s]
2165956it [00:01, 1684305.03it/s]
2338256it [00:01, 1695783.37it/s]
2507854it [00:01, 1640349.89it/s]
2672275it [00:01, 1603419.56it/s]
2844972it [00:01, 1639278.79it/s]
3014247it [00:01, 1654941.31it/s]
3183703it [00:01, 1666620.67it/s]
3350603it [00:02, 1660195.51it/s]
3517328it [00:02, 1662269.27it/s]
3684233it [00:02, 1664271.84it/s]
3854086it [00:02, 1674476.57it/s]
4024950it [00:02, 1684670.52it/s]
4196718it [00:02, 1694528.20it/s]
4368227it [00:02, 1700668.70it/s]
4538322it [00:02, 1694500.50it/s]
4707796it [00:02, 1688767.25it/s]
4881587it [00:02, 1703424.55it/s]
5051950it [00:03, 1690287.08it/s]
5221013it [00:03, 1668394.38it/s]
5392541it [00:03, 1682237.43it/s]
5566369it [00:03, 1698862.38it/s]
5739224it [00:03, 1707583.61it/s]
5913297it [00:03, 1717438.07it/s]
6086965it [00:03, 1723175.62it/s]
6259315it [00:03, 1721633.41it/s]
6431501it [00:03, 1720486.40it/s]
6603565it [00:03, 1717101.57it/s]
6775287it [00:04, 1710828.42it/s]
6946382it [00:04, 1710247.88it/s]
7118258it [00:04, 1712774.77it/s]
7290383it [00:04, 1715293.22it/s]
7461999it [00:04, 1715523.90it/s]
7633556it [00:04, 1621338.55it/s]
7796754it [00:04, 1607967.98it/s]
7967247it [00:04, 1636014.25it/s]
8137249it [00:04, 1654687.58it/s]
8307148it [00:04, 1667723.47it/s]
8478723it [00:05, 1681829.54it/s]
8650075it [00:05, 1691095.52it/s]
8819379it [00:05, 1678114.60it/s]
8987343it [00:05, 1667468.87it/s]
9159049it [00:05, 1682157.15it/s]
9328220it [00:05, 1684984.53it/s]
9503226it [00:05, 1704377.43it/s]
9678662it [00:05, 1719293.86it/s]
9850640it [00:05, 1718036.75it/s]
10023733it [00:05, 1721751.62it/s]
10195933it [00:06, 1676061.43it/s]
10363816it [00:06, 1598549.76it/s]
10524504it [00:06, 1576378.86it/s]
10684518it [00:06, 1583171.69it/s]
10852742it [00:06, 1611948.14it/s]
11025426it [00:06, 1645637.00it/s]
11191661it [00:06, 1650483.97it/s]
11362034it [00:06, 1666172.31it/s]
11531372it [00:06, 1674250.48it/s]
11702287it [00:06, 1684627.86it/s]
11873038it [00:07, 1691446.19it/s]
12042259it [00:07, 1688887.45it/s]
12215187it [00:07, 1700846.37it/s]
12386518it [00:07, 1704559.55it/s]
12558468it [00:07, 1709014.55it/s]
12734038it [00:07, 1722935.34it/s]
12906348it [00:07, 1720961.91it/s]
13078456it [00:07, 1720095.97it/s]
13250474it [00:07, 1716256.27it/s]
13347848it [00:07, 1677573.45it/s]
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai




[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai


[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-445ab1-R1_B23WHYVLT4_1.md.bai