Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/76/4b59989df067a4d7e5679fa8cc525a/aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bam
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/76/7c887275aae62345919d9c994ea3b9/.command.sh
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/76/7c887275aae62345919d9c994ea3b9/.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/76/4b59989df067a4d7e5679fa8cc525a/aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
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-ac6a71-R1_B23WHYVLT4_1
+ all_qc -bam aih-tih-sc-ac6a71-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-ac6a71-R1_B23WHYVLT4_1 -reference ref_genome.fa -region chr20
mkdir -p aih-tih-sc-ac6a71-R1_B23WHYVLT4_1/aih-tih-sc-ac6a71-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-ac6a71-R1_B23WHYVLT4_1/aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md/full_bed.bed
seqtool depth -bam aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bam -bed aih-tih-sc-ac6a71-R1_B23WHYVLT4_1/aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md/full_bed.bed -threads 4 -bam_threads 2 | bgzip -c > aih-tih-sc-ac6a71-R1_B23WHYVLT4_1/aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md/aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.depth.gz
saturation_plot -bams aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bam -region chr20 -output aih-tih-sc-ac6a71-R1_B23WHYVLT4_1/aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md/saturation_plot.png -threads 1
seqtool saturation_metric -at 1000000 -region chr20 -bam_threads 1 aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bam > aih-tih-sc-ac6a71-R1_B23WHYVLT4_1/aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md/saturation_qc.tsv
fragment_size -bams aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bam -region chr20 -output_prefix aih-tih-sc-ac6a71-R1_B23WHYVLT4_1/aih-tih-sc-ac6a71-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-ac6a71-R1_B23WHYVLT4_1.md.bam > aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.mpileup && contamination -input aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.mpileup -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz > aih-tih-sc-ac6a71-R1_B23WHYVLT4_1/aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md/contamination.tsv
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-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-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-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-ac6a71-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-ac6a71-R1_B23WHYVLT4_1.md.bai
146462it [00:00, 1464551.56it/s]
308516it [00:00, 1556190.72it/s]
476266it [00:00, 1611542.35it/s]
646820it [00:00, 1648619.48it/s]
818251it [00:00, 1672143.68it/s]
987440it [00:00, 1678832.88it/s]
1156385it [00:00, 1682284.48it/s]
1324648it [00:00, 1682374.55it/s]
1495722it [00:00, 1691011.47it/s]
1666918it [00:01, 1697455.91it/s]
1836664it [00:01, 1694394.22it/s]
2006105it [00:01, 1681613.62it/s]
2174288it [00:01, 1679339.56it/s]
2344691it [00:01, 1686737.70it/s]
2514342it [00:01, 1689653.56it/s]
2689358it [00:01, 1707795.78it/s]
2860150it [00:01, 1705132.02it/s]
3030673it [00:01, 1702218.94it/s]
3200902it [00:01, 1689049.50it/s]
3369832it [00:02, 1686170.32it/s]
3539996it [00:02, 1690593.24it/s]
3709070it [00:02, 1686132.83it/s]
3880100it [00:02, 1693328.08it/s]
4049444it [00:02, 1691353.40it/s]
4220341it [00:02, 1696490.12it/s]
4389997it [00:02, 1696176.23it/s]
4559620it [00:02, 1694800.00it/s]
4731482it [00:02, 1701912.98it/s]
4901677it [00:02, 1696793.72it/s]
5071362it [00:03, 1684208.65it/s]
5241130it [00:03, 1688205.43it/s]
5411816it [00:03, 1693614.49it/s]
5583790it [00:03, 1701408.16it/s]
5755101it [00:03, 1704894.35it/s]
5927604it [00:03, 1710908.73it/s]
6099035it [00:03, 1711909.96it/s]
6270232it [00:03, 1706912.70it/s]
6442122it [00:03, 1710483.87it/s]
6614112it [00:03, 1713132.65it/s]
6785430it [00:04, 1707917.61it/s]
6956228it [00:04, 1702278.13it/s]
7129785it [00:04, 1712172.49it/s]
7301012it [00:04, 1709932.06it/s]
7472012it [00:04, 1701831.51it/s]
7642207it [00:04, 1698329.05it/s]
7812790it [00:04, 1700533.95it/s]
7982850it [00:04, 1683445.81it/s]
8153049it [00:04, 1688806.15it/s]
8321962it [00:04, 1680057.42it/s]
8491001it [00:05, 1683106.34it/s]
8659841it [00:05, 1684667.49it/s]
8829049it [00:05, 1686863.20it/s]
9000437it [00:05, 1694925.55it/s]
9169941it [00:05, 1691191.63it/s]
9340961it [00:05, 1696838.77it/s]
9511986it [00:05, 1700811.60it/s]
9685044it [00:05, 1709569.25it/s]
9858124it [00:05, 1715909.78it/s]
10029720it [00:05, 1706779.13it/s]
10200412it [00:06, 1694031.92it/s]
10369842it [00:06, 1693637.62it/s]
10542553it [00:06, 1703590.62it/s]
10713356it [00:06, 1704885.74it/s]
10883859it [00:06, 1702050.26it/s]
11054074it [00:06, 1699422.40it/s]
11224024it [00:06, 1687582.30it/s]
11392803it [00:06, 1681626.54it/s]
11562986it [00:06, 1687610.96it/s]
11731899it [00:06, 1688041.91it/s]
11901707it [00:07, 1691018.43it/s]
12055969it [00:07, 1692329.38it/s]
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-ac6a71-R1_B23WHYVLT4_1.md.bai