Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/eb/3f9b8fcca2aee8f8e9e55406e975c2/aih-tih-sc-3e91d2-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/7f/d7c88e002b98064e5561f615e17014/.command.sh
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/eb/3f9b8fcca2aee8f8e9e55406e975c2/aih-tih-sc-3e91d2-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/7f/d7c88e002b98064e5561f615e17014/.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-3e91d2-R1_B23WHYVLT4_1
+ all_qc -bam aih-tih-sc-3e91d2-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-3e91d2-R1_B23WHYVLT4_1 -reference ref_genome.fa -region chr20
mkdir -p aih-tih-sc-3e91d2-R1_B23WHYVLT4_1/aih-tih-sc-3e91d2-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-3e91d2-R1_B23WHYVLT4_1/aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md/full_bed.bed
seqtool depth -bam aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bam -bed aih-tih-sc-3e91d2-R1_B23WHYVLT4_1/aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md/full_bed.bed -threads 4 -bam_threads 2 | bgzip -c > aih-tih-sc-3e91d2-R1_B23WHYVLT4_1/aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md/aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.depth.gz
saturation_plot -bams aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bam -region chr20 -output aih-tih-sc-3e91d2-R1_B23WHYVLT4_1/aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md/saturation_plot.png -threads 1
seqtool saturation_metric -at 1000000 -region chr20 -bam_threads 1 aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bam > aih-tih-sc-3e91d2-R1_B23WHYVLT4_1/aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md/saturation_qc.tsv
fragment_size -bams aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bam -region chr20 -output_prefix aih-tih-sc-3e91d2-R1_B23WHYVLT4_1/aih-tih-sc-3e91d2-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-3e91d2-R1_B23WHYVLT4_1.md.bam > aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.mpileup && contamination -input aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.mpileup -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz > aih-tih-sc-3e91d2-R1_B23WHYVLT4_1/aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md/contamination.tsv
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-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-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-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-3e91d2-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-3e91d2-R1_B23WHYVLT4_1.md.bai
152190it [00:00, 1521854.28it/s]
317258it [00:00, 1597595.99it/s]
481766it [00:00, 1619207.96it/s]
643687it [00:00, 1531265.76it/s]
812678it [00:00, 1586392.90it/s]
971952it [00:00, 1562971.97it/s]
1134713it [00:00, 1583442.90it/s]
1299536it [00:00, 1603673.68it/s]
1465775it [00:00, 1621774.54it/s]
1631261it [00:01, 1631805.50it/s]
1797498it [00:01, 1641086.58it/s]
1965182it [00:01, 1651900.84it/s]
2130457it [00:01, 1648443.73it/s]
2301791it [00:01, 1667958.11it/s]
2470759it [00:01, 1674474.89it/s]
2639262it [00:01, 1677625.06it/s]
2807051it [00:01, 1671105.09it/s]
2974185it [00:01, 1666136.83it/s]
3140817it [00:01, 1661462.45it/s]
3306977it [00:02, 1555359.25it/s]
3475365it [00:02, 1592027.79it/s]
3635699it [00:02, 1567800.37it/s]
3800083it [00:02, 1589771.45it/s]
3971497it [00:02, 1626106.47it/s]
4143553it [00:02, 1653910.38it/s]
4309381it [00:02, 1645627.96it/s]
4477178it [00:02, 1655194.06it/s]
4646865it [00:02, 1667570.93it/s]
4813793it [00:02, 1662816.43it/s]
4980195it [00:03, 1658495.43it/s]
5146129it [00:03, 1650638.06it/s]
5313831it [00:03, 1658371.33it/s]
5482189it [00:03, 1665874.45it/s]
5650427it [00:03, 1670793.55it/s]
5820509it [00:03, 1679760.91it/s]
5989402it [00:03, 1682494.64it/s]
6157669it [00:03, 1673486.10it/s]
6325040it [00:03, 1668496.07it/s]
6492106it [00:03, 1669129.13it/s]
6659031it [00:04, 1667447.74it/s]
6825784it [00:04, 1660096.93it/s]
6991805it [00:04, 1657624.31it/s]
7157575it [00:04, 1653456.59it/s]
7327663it [00:04, 1667598.26it/s]
7494434it [00:04, 1651584.99it/s]
7659630it [00:04, 1634503.26it/s]
7823132it [00:04, 1615780.33it/s]
7984771it [00:04, 1610516.99it/s]
8145860it [00:04, 1605680.01it/s]
8306451it [00:05, 1582356.05it/s]
8464752it [00:05, 1565737.21it/s]
8622153it [00:05, 1568154.69it/s]
8779013it [00:05, 1561969.08it/s]
8935238it [00:05, 1560782.41it/s]
9091335it [00:05, 1556421.45it/s]
9249335it [00:05, 1563421.84it/s]
9405692it [00:05, 1560562.39it/s]
9562523it [00:05, 1562861.75it/s]
9719861it [00:05, 1565991.76it/s]
9881280it [00:06, 1580397.85it/s]
10039329it [00:06, 1573633.02it/s]
10197364it [00:06, 1575627.39it/s]
10355351it [00:06, 1576885.46it/s]
10514844it [00:06, 1582275.16it/s]
10673078it [00:06, 1575258.15it/s]
10830615it [00:06, 1534900.46it/s]
10984326it [00:06, 1523352.94it/s]
11136813it [00:06, 1520140.25it/s]
11291576it [00:07, 1528245.14it/s]
11446090it [00:07, 1533246.23it/s]
11604941it [00:07, 1549644.27it/s]
11760814it [00:07, 1552345.33it/s]
11917693it [00:07, 1557247.70it/s]
12073450it [00:07, 1543094.06it/s]
12229475it [00:07, 1548114.92it/s]
12387001it [00:07, 1556194.34it/s]
12542652it [00:07, 1554843.79it/s]
12700477it [00:07, 1561781.08it/s]
12857387it [00:08, 1563960.52it/s]
13015844it [00:08, 1570118.33it/s]
13172866it [00:08, 1559286.46it/s]
13333654it [00:08, 1573770.38it/s]
13500040it [00:08, 1600677.78it/s]
13668018it [00:08, 1624332.74it/s]
13833058it [00:08, 1632024.90it/s]
14004097it [00:08, 1655374.68it/s]
14172872it [00:08, 1665065.58it/s]
14343286it [00:08, 1676731.71it/s]
14512834it [00:09, 1682341.60it/s]
14682827it [00:09, 1687579.18it/s]
14851591it [00:09, 1683836.50it/s]
15019981it [00:09, 1668006.81it/s]
15186817it [00:09, 1666710.81it/s]
15353512it [00:09, 1658257.07it/s]
15519360it [00:09, 1657042.15it/s]
15685655it [00:09, 1658792.47it/s]
15851546it [00:09, 1658294.21it/s]
16017384it [00:09, 1655801.08it/s]
16182970it [00:10, 1640919.42it/s][W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
16349770it [00:10, 1648873.46it/s]
16514685it [00:10, 1638983.87it/s]
16678611it [00:10, 1634287.94it/s]
16842059it [00:10, 1631555.65it/s]
17005227it [00:10, 1620612.42it/s]
17168832it [00:10, 1625126.56it/s]
17331362it [00:10, 1618264.21it/s]
17494052it [00:10, 1620824.05it/s]
17659787it [00:10, 1631710.55it/s]
17822973it [00:11, 1622297.51it/s]
17986241it [00:11, 1625377.97it/s]
18148795it [00:11, 1621692.13it/s]
18314772it [00:11, 1633047.11it/s]
18481193it [00:11, 1642219.08it/s]
18645428it [00:11, 1641888.55it/s]
18713543it [00:11, 1617295.08it/s]
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3e91d2-R1_B23WHYVLT4_1.md.bai