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/d3/150cf6af1c0ddc3ff8180038cc22eb/.command.sh
  Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/5e/6c81b8bfd3d284f2fab151ad186064/aih-tih-sc-76a805-R1_B23WHYVLT4_1.md.bai
  Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/5e/6c81b8bfd3d284f2fab151ad186064/aih-tih-sc-76a805-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/d3/150cf6af1c0ddc3ff8180038cc22eb/.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-76a805-R1_B23WHYVLT4_1
+ all_qc -bam aih-tih-sc-76a805-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-76a805-R1_B23WHYVLT4_1 -reference ref_genome.fa -region chr20
mkdir -p aih-tih-sc-76a805-R1_B23WHYVLT4_1/aih-tih-sc-76a805-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-76a805-R1_B23WHYVLT4_1/aih-tih-sc-76a805-R1_B23WHYVLT4_1.md/full_bed.bed


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

0it [00:00, ?it/s]
145816it [00:00, 1457817.30it/s]
299397it [00:00, 1503647.48it/s]
461512it [00:00, 1557257.93it/s]
627843it [00:00, 1599094.22it/s]
793880it [00:00, 1621162.33it/s]
964139it [00:00, 1648821.75it/s]
1132716it [00:00, 1660879.04it/s]
1303655it [00:00, 1676228.44it/s]
1471278it [00:00, 1671844.74it/s]
1645741it [00:01, 1694199.23it/s]
1815167it [00:01, 1691756.35it/s]
1986152it [00:01, 1697214.89it/s]
2156849it [00:01, 1700137.10it/s]
2330121it [00:01, 1709948.84it/s]
2503415it [00:01, 1716848.98it/s]
2675103it [00:01, 1704098.83it/s]
2845536it [00:01, 1703289.25it/s]
3016654it [00:01, 1705358.78it/s]
3187812it [00:01, 1707200.04it/s]
3358541it [00:02, 1688646.76it/s]
3527454it [00:02, 1684196.62it/s]
3697427it [00:02, 1688624.82it/s]
3870604it [00:02, 1701453.36it/s]
4043517it [00:02, 1709556.95it/s]
4222466it [00:02, 1733437.95it/s]
4396438it [00:02, 1735289.53it/s]
4569983it [00:02, 1729097.70it/s]
4742909it [00:02, 1722552.54it/s]
4915179it [00:02, 1712510.71it/s]
5086449it [00:03, 1690828.29it/s]
5259896it [00:03, 1703711.36it/s]
5430748it [00:03, 1705120.61it/s]
5603458it [00:03, 1711645.21it/s]
5777551it [00:03, 1720326.28it/s]
5950605it [00:03, 1723361.72it/s]
6122963it [00:03, 1723116.49it/s]
6295290it [00:03, 1714915.85it/s]
6466800it [00:03, 1691900.06it/s]
6641048it [00:03, 1706850.52it/s]
6811804it [00:04, 1704670.56it/s]
6985959it [00:04, 1715630.32it/s]
7161174it [00:04, 1726505.10it/s]
7333860it [00:04, 1716081.30it/s]
7505504it [00:04, 1712738.00it/s]
7676802it [00:04, 1711454.14it/s]
7847964it [00:04, 1706388.43it/s]
8019985it [00:04, 1710490.66it/s]
8191046it [00:04, 1710004.13it/s]
8362136it [00:04, 1710248.20it/s]
8535460it [00:05, 1717102.61it/s]
8707176it [00:05, 1712349.42it/s]
8879291it [00:05, 1714955.35it/s]
9050792it [00:05, 1711374.84it/s]
9221934it [00:05, 1701506.52it/s]
9392099it [00:05, 1697811.81it/s]
9561890it [00:05, 1692301.15it/s]
9731744it [00:05, 1694142.13it/s]
9902708it [00:05, 1698599.46it/s]
10072575it [00:05, 1697735.01it/s]
10244685it [00:06, 1704629.43it/s]
10418073it [00:06, 1713366.70it/s]
10589415it [00:06, 1697870.06it/s]
10762548it [00:06, 1707759.13it/s]
10933355it [00:06, 1685201.40it/s]
11101954it [00:06, 1503596.87it/s]
11256775it [00:06, 1515804.22it/s]
11421008it [00:06, 1551378.69it/s]
11589331it [00:06, 1589093.65it/s]
11749863it [00:07, 1512547.53it/s]
11922627it [00:07, 1573307.12it/s]
12096585it [00:07, 1621100.79it/s]
12270343it [00:07, 1654999.03it/s]
12443669it [00:07, 1677972.29it/s]
12616515it [00:07, 1692801.16it/s]
12788197it [00:07, 1699911.34it/s]
12962935it [00:07, 1714027.29it/s]
13135354it [00:07, 1717007.30it/s]
13307260it [00:07, 1701404.51it/s]
13477573it [00:08, 1692720.06it/s]
13647472it [00:08, 1694528.38it/s]
13817014it [00:08, 1688450.47it/s]
13988852it [00:08, 1697338.24it/s]
14163589it [00:08, 1712122.62it/s]
14334845it [00:08, 1705647.43it/s]
14505705it [00:08, 1706508.28it/s]
14680259it [00:08, 1718150.74it/s]
14852096it [00:08, 1716006.07it/s]
15023712it [00:08, 1712250.66it/s]
15194949it [00:09, 1707979.72it/s]
15365756it [00:09, 1692031.71it/s]
15534994it [00:09, 1686218.72it/s]
15703639it [00:09, 1673353.73it/s]
15871002it [00:09, 1672046.21it/s]
16038392it [00:09, 1672448.52it/s]
16205650it [00:09, 1667314.11it/s]
16372391it [00:09, 1666634.97it/s]
16539061it [00:09, 1636909.22it/s]
16702862it [00:09, 1504874.65it/s]
16865650it [00:10, 1539222.17it/s]
17031483it [00:10, 1573116.42it/s]
17201810it [00:10, 1610780.92it/s][W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-76a805-R1_B23WHYVLT4_1.md.bai

17363969it [00:10, 1612125.43it/s]
17525942it [00:10, 1574089.43it/s]
17684031it [00:10, 1519913.03it/s]
17836776it [00:10, 1500550.00it/s]
17987334it [00:10, 1420386.06it/s][W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-76a805-R1_B23WHYVLT4_1.md.bai

18153236it [00:10, 1487075.02it/s]
18315462it [00:11, 1525528.04it/s]
18479568it [00:11, 1558880.81it/s]
18644383it [00:11, 1584986.23it/s]
18810785it [00:11, 1608267.68it/s]
18981182it [00:11, 1636607.56it/s]
19151316it [00:11, 1655734.25it/s]
19299578it [00:11, 1665465.23it/s]




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


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


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


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


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