Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/96/d7b7ccc2a38ec3ac0fea5d427e55bd/aih-tih-sc-3b3c1e-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/64/e4adf2fa8df87c2b5e8f9f42560576/.command.sh
  Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/96/d7b7ccc2a38ec3ac0fea5d427e55bd/aih-tih-sc-3b3c1e-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/64/e4adf2fa8df87c2b5e8f9f42560576/.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-3b3c1e-R1_B23WHYVLT4_1
+ all_qc -bam aih-tih-sc-3b3c1e-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-3b3c1e-R1_B23WHYVLT4_1 -reference ref_genome.fa -region chr20
mkdir -p aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1/aih-tih-sc-3b3c1e-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-3b3c1e-R1_B23WHYVLT4_1/aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1.md/full_bed.bed


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

0it [00:00, ?it/s]
157056it [00:00, 1570494.10it/s]
326298it [00:00, 1642187.26it/s]
495075it [00:00, 1662979.92it/s]
667648it [00:00, 1687633.99it/s]
839575it [00:00, 1699023.11it/s]
1011757it [00:00, 1706611.04it/s]
1186878it [00:00, 1721173.21it/s]
1358996it [00:00, 1719076.62it/s]
1530905it [00:00, 1718831.72it/s]
1702789it [00:01, 1717984.37it/s]
1874588it [00:01, 1711862.30it/s]
2046213it [00:01, 1713177.39it/s]
2223127it [00:01, 1729967.39it/s]
2397343it [00:01, 1733628.74it/s]
2570711it [00:01, 1733440.59it/s]
2744059it [00:01, 1721699.08it/s]
2917626it [00:01, 1725865.88it/s]
3091244it [00:01, 1728939.24it/s]
3264151it [00:01, 1723885.51it/s]
3436551it [00:02, 1706627.12it/s]
3607255it [00:02, 1703968.14it/s]
3779572it [00:02, 1709653.84it/s]
3950561it [00:02, 1704256.59it/s]
4121004it [00:02, 1699052.84it/s]
4291374it [00:02, 1700413.24it/s]
4461426it [00:02, 1692846.69it/s]
4630722it [00:02, 1688691.94it/s]
4799599it [00:02, 1682168.57it/s]
4967824it [00:02, 1672907.67it/s]
5135170it [00:03, 1672971.57it/s]
5302475it [00:03, 1664700.43it/s]
5470923it [00:03, 1670558.78it/s]
5638065it [00:03, 1670760.34it/s]
5805150it [00:03, 1656379.24it/s]
5970817it [00:03, 1647220.18it/s]
6135582it [00:03, 1647329.64it/s]
6302219it [00:03, 1652985.44it/s]
6467534it [00:03, 1649938.47it/s]
6633434it [00:03, 1652625.57it/s]
6798706it [00:04, 1646385.93it/s]
6966522it [00:04, 1655852.00it/s]
7137371it [00:04, 1671509.65it/s]
7307561it [00:04, 1680585.26it/s]
7477389it [00:04, 1685870.38it/s]
7645986it [00:04, 1682060.35it/s]
7814200it [00:04, 1679163.76it/s]
7983173it [00:04, 1682263.61it/s]
8154602it [00:04, 1691821.79it/s]
8324317it [00:04, 1693401.45it/s]
8494656it [00:05, 1696243.89it/s]
8667173it [00:05, 1704894.65it/s]
8838586it [00:05, 1707650.76it/s]
9009354it [00:05, 1691520.73it/s]
9181468it [00:05, 1700331.08it/s]
9356433it [00:05, 1715035.89it/s]
9527966it [00:05, 1705858.67it/s]
9698581it [00:05, 1696996.51it/s]
9869668it [00:05, 1701116.10it/s]
10039801it [00:05, 1699871.03it/s]
10209803it [00:06, 1696368.38it/s]
10379450it [00:06, 1695594.00it/s]
10549571it [00:06, 1697259.99it/s]
10721239it [00:06, 1703059.79it/s]
10892942it [00:06, 1707117.23it/s]
11064069it [00:06, 1708347.31it/s]
11234907it [00:06, 1697745.50it/s]
11405562it [00:06, 1700359.64it/s]
11576870it [00:06, 1704150.73it/s]
11747296it [00:06, 1694898.86it/s]
11917181it [00:07, 1696066.06it/s]
12089942it [00:07, 1705477.14it/s]
12229518it [00:07, 1693977.05it/s]
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1.md.bai




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


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


tabix -s 1 -b 2 -e 2 aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1/aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1.md/aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1.md.depth.gz


seqtool depth_hist     -nbins 500     aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1/aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1.md/aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1.md.depth.gz     > aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1/aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1.md/aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1.md.depth_hist.tsv


cat aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1/aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1.md/aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1.md.depth_hist.tsv | plot         -x 0         -xlabel Depth -ylabel "Number of Sites"         -xlim 0 2000         -title aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1.md         -output aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1/aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1.md/aih-tih-sc-3b3c1e-R1_B23WHYVLT4_1.md.depth.png