Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/29/ed3e66b98810a13421e120c2d3991e/aih-tih-sc-099ac4-R1_B23WHYVLT4_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/66/78561455fa5c37648a5fbd763e75f3/.command.sh
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/29/ed3e66b98810a13421e120c2d3991e/aih-tih-sc-099ac4-R1_B23WHYVLT4_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/66/78561455fa5c37648a5fbd763e75f3/.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-099ac4-R1_B23WHYVLT4_1
+ all_qc -bam aih-tih-sc-099ac4-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-099ac4-R1_B23WHYVLT4_1 -reference ref_genome.fa -region chr20
mkdir -p aih-tih-sc-099ac4-R1_B23WHYVLT4_1/aih-tih-sc-099ac4-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-099ac4-R1_B23WHYVLT4_1/aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md/full_bed.bed
seqtool depth -bam aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md.bam -bed aih-tih-sc-099ac4-R1_B23WHYVLT4_1/aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md/full_bed.bed -threads 4 -bam_threads 2 | bgzip -c > aih-tih-sc-099ac4-R1_B23WHYVLT4_1/aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md/aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md.depth.gz
saturation_plot -bams aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md.bam -region chr20 -output aih-tih-sc-099ac4-R1_B23WHYVLT4_1/aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md/saturation_plot.png -threads 1
seqtool saturation_metric -at 1000000 -region chr20 -bam_threads 1 aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md.bam > aih-tih-sc-099ac4-R1_B23WHYVLT4_1/aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md/saturation_qc.tsv
fragment_size -bams aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md.bam -region chr20 -output_prefix aih-tih-sc-099ac4-R1_B23WHYVLT4_1/aih-tih-sc-099ac4-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-099ac4-R1_B23WHYVLT4_1.md.bam > aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md.mpileup && contamination -input aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md.mpileup -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz > aih-tih-sc-099ac4-R1_B23WHYVLT4_1/aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md/contamination.tsv
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-099ac4-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-099ac4-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-099ac4-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-099ac4-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-099ac4-R1_B23WHYVLT4_1.md.bai
149909it [00:00, 1499005.66it/s]
314955it [00:00, 1588052.96it/s]
480912it [00:00, 1620454.94it/s]
649852it [00:00, 1647542.59it/s]
818116it [00:00, 1660172.32it/s]
988927it [00:00, 1676403.67it/s]
1159740it [00:00, 1686749.09it/s]
1329783it [00:00, 1690977.60it/s]
1499567it [00:00, 1693012.15it/s]
1670995it [00:01, 1699550.21it/s]
1842413it [00:01, 1704002.91it/s]
2013464it [00:01, 1705956.80it/s]
2184060it [00:01, 1703514.90it/s]
2355278it [00:01, 1706106.62it/s]
2533419it [00:01, 1728760.60it/s]
2708254it [00:01, 1734597.17it/s]
2881716it [00:01, 1721544.24it/s]
3053894it [00:01, 1711126.18it/s]
3226210it [00:01, 1714687.12it/s]
3397700it [00:02, 1696405.06it/s]
3567390it [00:02, 1691435.93it/s]
3736566it [00:02, 1689329.99it/s]
3907463it [00:02, 1695038.25it/s]
4077388it [00:02, 1696250.64it/s]
4247028it [00:02, 1691477.63it/s]
4416187it [00:02, 1663254.39it/s]
4582614it [00:02, 1626178.19it/s]
4745438it [00:02, 1619039.07it/s]
4912798it [00:02, 1635016.80it/s]
5077037it [00:03, 1637170.90it/s]
5245023it [00:03, 1649755.07it/s]
5410078it [00:03, 1646489.77it/s]
5580159it [00:03, 1662653.35it/s]
5754982it [00:03, 1688185.14it/s]
5926900it [00:03, 1697434.47it/s]
6098218it [00:03, 1702130.38it/s]
6268458it [00:03, 1699253.42it/s]
6438403it [00:03, 1685817.27it/s]
6609574it [00:03, 1693495.67it/s]
6778954it [00:04, 1652326.36it/s]
6945830it [00:04, 1657132.74it/s]
7113847it [00:04, 1663920.03it/s]
7280366it [00:04, 1591290.03it/s]
7440199it [00:04, 1583579.30it/s]
7599031it [00:04, 1583040.14it/s]
7757666it [00:04, 1569407.06it/s]
7917725it [00:04, 1578450.79it/s]
8075752it [00:04, 1574463.60it/s]
8237312it [00:04, 1586623.96it/s]
8405857it [00:05, 1615977.18it/s]
8574815it [00:05, 1637878.15it/s]
8744715it [00:05, 1656112.01it/s]
8910542it [00:05, 1656738.89it/s]
9080384it [00:05, 1669193.86it/s]
9250495it [00:05, 1678721.50it/s]
9421417it [00:05, 1687730.93it/s]
9590878it [00:05, 1689772.86it/s]
9760451it [00:05, 1691538.85it/s]
9932320it [00:05, 1699652.02it/s]
10105124it [00:06, 1708095.40it/s]
10277541it [00:06, 1712897.48it/s]
10448835it [00:06, 1708514.67it/s]
10621668it [00:06, 1714367.55it/s]
10793110it [00:06, 1711830.74it/s]
10964297it [00:06, 1701296.98it/s]
11134443it [00:06, 1682197.87it/s]
11302713it [00:06, 1665390.08it/s]
11469307it [00:06, 1652383.69it/s]
11634588it [00:06, 1648686.74it/s]
11803858it [00:07, 1661691.14it/s]
11974043it [00:07, 1673608.00it/s]
12144251it [00:07, 1682075.59it/s]
12316567it [00:07, 1694325.96it/s]
12487639it [00:07, 1699210.75it/s]
12659985it [00:07, 1706451.38it/s]
12830646it [00:07, 1702364.79it/s]
13000896it [00:07, 1660395.85it/s]
13168477it [00:07, 1664912.41it/s]
13335166it [00:07, 1665479.27it/s]
13501833it [00:08, 1639749.90it/s]
13665959it [00:08, 1634981.28it/s]
13829558it [00:08, 1619705.13it/s]
13991611it [00:08, 1587304.75it/s]
14150495it [00:08, 1575421.13it/s]
14313443it [00:08, 1591218.48it/s]
14472672it [00:08, 1579537.95it/s]
14630702it [00:08, 1573153.48it/s]
14788066it [00:08, 1570408.63it/s]
14946956it [00:09, 1575863.41it/s]
15108351it [00:09, 1587066.99it/s]
15276444it [00:09, 1615039.49it/s]
15443548it [00:09, 1631754.64it/s]
15612526it [00:09, 1649046.48it/s]
15782418it [00:09, 1663960.45it/s]
15949678it [00:09, 1666531.19it/s]
16116348it [00:09, 1645734.37it/s][W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md.bai
16281955it [00:09, 1648787.38it/s]
16446885it [00:09, 1643382.87it/s]
16612108it [00:10, 1646001.26it/s]
16776737it [00:10, 1646071.52it/s]
16946651it [00:10, 1661917.11it/s]
17116398it [00:10, 1672477.45it/s]
17283662it [00:10, 1671981.04it/s]
17450871it [00:10, 1669205.39it/s]
17620737it [00:10, 1678003.73it/s]
17788864it [00:10, 1678934.62it/s]
17959110it [00:10, 1685895.82it/s]
18082130it [00:10, 1660872.49it/s]
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-099ac4-R1_B23WHYVLT4_1.md.bai