Downloading: s3://natera-platform-sandbox/pipeline-resources/ngi-igenomes/igenomes/Homo_sapiens/GATK/GRCh38/Sequence/WholeGenomeFasta/Homo_sapiens_assembly38.fasta
Downloading: s3://natera-platform-sandbox/pipeline-resources/Homo_sapiens_assembly38.dbsnp138.vcf.gz
Downloading: s3://natera-platform-sandbox/pipeline-resources/beds/altera_subpanels/hg38/altera_v3_targets_postQC_hg38.bed
Downloading: s3://natera-platform-sandbox/pipeline-resources/beds/altera_subpanels/hg38/altera_v3_targets_coding_postQC_hg38.bed
Downloading: s3://natera-platform-sandbox/pipeline-resources/Homo_sapiens_assembly38.dbsnp138.vcf.gz.tbi
Downloading: s3://natera-platform-sandbox/pipeline-resources/beds/altera_subpanels/hg38/foresight/foresight_clarity_targets_hg38.bed
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/ea/a359c86709868b3fa49045f91bbe84/.command.sh
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/fb/ef37e6069f2668d49419f98ecb4969/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
Downloading: s3://natera-platform-sandbox/pipeline-resources/ngi-igenomes/igenomes/Homo_sapiens/GATK/GRCh38/Sequence/WholeGenomeFasta/Homo_sapiens_assembly38.fasta.fai
Downloading: s3://natera-platform-sandbox/pipeline-resources/beds/altera_subpanels/hg38/FINAL_hotspot_list_07302025_hg38_sorted.bed
Downloading: s3://natera-platform-sandbox/pipeline-resources/beds/altera_subpanels/hg38/foresight/xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity.bed
Downloading: s3://natera-platform-sandbox/platform-users/ralla/nfcore_rnafusion_resources/GRCh38/gatk4/dbsnp_hg38_contam_sites.bed
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/ea/a359c86709868b3fa49045f91bbe84/.command.run
Downloading: s3://natera-platform-sandbox/pipeline-resources/beds/altera_subpanels/hg38/xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed
Downloading: s3://natera-platform-sandbox/pipeline-resources/beds/dbsnp138_hg38_common.bed
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/fb/ef37e6069f2668d49419f98ecb4969/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam
Downloading: s3://natera-platform-sandbox/pipeline-resources/beds/altera_subpanels/hg38/foresight/xgen-exome-hyb-panel-v2_AND_altera_v3_probes_short_hg38_AND_cdx_spike_custom30off4x_hg38_AND_foresight_clarity_footprint.bed
Downloading: s3://natera-platform-sandbox/pipeline-resources/beds/altera_subpanels/hg38/altera_v3_targets_noncoding_postQC_hg38.bed
==> STAGING COMPLETE (18 inputs)
mkdir -p 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/
cat xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity.bed xgen-exome-hyb-panel-v2_AND_altera_v3_probes_short_hg38_AND_cdx_spike_custom30off4x_hg38_AND_foresight_clarity_footprint.bed altera_v3_targets_postQC_hg38.bed altera_v3_targets_coding_postQC_hg38.bed altera_v3_targets_noncoding_postQC_hg38.bed xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed FINAL_hotspot_list_07302025_hg38_sorted.bed foresight_clarity_targets_hg38.bed | cut -f 1-3 | bedtools sort -g Homo_sapiens_assembly38.fasta.fai | bedtools merge > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/full_bed.bed
seqtool depth -bam 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam -bed 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/full_bed.bed -threads 4 -bam_threads 2 | bgzip -c > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz
saturation_plot -bams 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam -region chr20 -output 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/saturation_plot.png -threads 1
seqtool saturation_metric -at 1000000 -region chr20 -bam_threads 1 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/saturation_qc.tsv
fragment_size -bams 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam -region chr20 -output_prefix 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/fragsize -threads 1
samtools mpileup -B --min-MQ 60 --min-BQ 30 -f Homo_sapiens_assembly38.fasta -l dbsnp_hg38_contam_sites.bed 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.mpileup && contamination -input 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.mpileup -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/contamination.tsv
seqtool noise -bam 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam -ref Homo_sapiens_assembly38.fasta -region chr20 -include_bed 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/full_bed.bed -exclude_bed dbsnp138_hg38_common.bed -threads 1 > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.noise_summary.tsv
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[mpileup] 1 samples in 1 input files
loading BEDs...
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.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: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
0it [00:00, ?it/s]
148391it [00:00, 1483819.43it/s][W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
311408it [00:00, 1569873.70it/s]
477440it [00:00, 1611152.01it/s]
647152it [00:00, 1645072.48it/s]
813167it [00:00, 1650492.96it/s]
978217it [00:00, 1614887.35it/s]
1139849it [00:00, 1588815.82it/s]
1300437it [00:00, 1594116.86it/s]
1469304it [00:00, 1623223.99it/s]
1631747it [00:01, 1608549.53it/s]
1792701it [00:01, 1561761.78it/s]
1949182it [00:01, 1544841.88it/s]
2113505it [00:01, 1573773.94it/s]
2271113it [00:01, 1564111.89it/s]
2434530it [00:01, 1584818.16it/s]
2593161it [00:01, 1544114.31it/s]
2754008it [00:01, 1562919.91it/s]
2914281it [00:01, 1574494.69it/s]
3071932it [00:01, 1539523.02it/s]
3226157it [00:02, 1449789.70it/s]
3392410it [00:02, 1509842.62it/s]
3562306it [00:02, 1564184.50it/s]
3730893it [00:02, 1599501.60it/s]
3901658it [00:02, 1631251.00it/s]
4074823it [00:02, 1660907.96it/s]
4246085it [00:02, 1676245.89it/s]
4417273it [00:02, 1686842.62it/s]
4587263it [00:02, 1690717.27it/s]
4756524it [00:02, 1687939.35it/s]
4925903it [00:03, 1689672.50it/s]
5094964it [00:03, 1686437.34it/s]
5266042it [00:03, 1693702.62it/s]
5435461it [00:03, 1690126.05it/s]
5604509it [00:03, 1685622.31it/s]
5776204it [00:03, 1694968.80it/s]
5945723it [00:03, 1694066.51it/s]
6115524it [00:03, 1695231.42it/s]
6285059it [00:03, 1647389.23it/s]
6457601it [00:03, 1670272.79it/s]
6627014it [00:04, 1677138.28it/s]
6795476it [00:04, 1679346.52it/s]
6963546it [00:04, 1663781.30it/s]
7138182it [00:04, 1688079.15it/s]
7311858it [00:04, 1702504.41it/s]
7485752it [00:04, 1713207.64it/s]
7657142it [00:04, 1676368.30it/s]
7827846it [00:04, 1685390.13it/s]
7996550it [00:04, 1674796.03it/s]
8166324it [00:04, 1681572.03it/s]
8335950it [00:05, 1685727.84it/s]
8507612it [00:05, 1694919.82it/s]
8678033it [00:05, 1697680.78it/s]
8847841it [00:05, 1677990.75it/s]
9017757it [00:05, 1684222.68it/s]
9189892it [00:05, 1695195.72it/s]
9361613it [00:05, 1701629.09it/s]
9532748it [00:05, 1704519.19it/s]
9703228it [00:05, 1702268.47it/s]
9873864it [00:05, 1703474.76it/s]
10044527it [00:06, 1704405.36it/s]
10216462it [00:06, 1708862.30it/s]
10387356it [00:06, 1708682.17it/s]
10558230it [00:06, 1696024.42it/s]
10727857it [00:06, 1695644.75it/s]
10898603it [00:06, 1699159.15it/s]
11068533it [00:06, 1691844.39it/s]
11237732it [00:06, 1684639.19it/s]done loading BEDs.
11410204it [00:06, 1696559.55it/s]
11583369it [00:07, 1707017.86it/s]
11756861it [00:07, 1715338.57it/s]
11928793it [00:07, 1716507.49it/s]
12100455it [00:07, 1708438.41it/s]
12271314it [00:07, 1701026.36it/s]
12441432it [00:07, 1651420.95it/s]
12612377it [00:07, 1668377.19it/s]
12785605it [00:07, 1687160.96it/s]
12956259it [00:07, 1692880.48it/s]
13125702it [00:07, 1614951.67it/s]
13288013it [00:08, 1609076.71it/s]
13459226it [00:08, 1639011.22it/s]
13627561it [00:08, 1651999.63it/s]
13797940it [00:08, 1667269.96it/s]
13966498it [00:08, 1672683.53it/s]
14138039it [00:08, 1685285.43it/s]
14308836it [00:08, 1692031.54it/s]
14479597it [00:08, 1696514.90it/s]
14649327it [00:08, 1691232.43it/s]
14823027it [00:08, 1704884.94it/s]
14995583it [00:09, 1711050.13it/s]
15166722it [00:09, 1690674.90it/s]
15335865it [00:09, 1658347.27it/s]
15502613it [00:09, 1660894.30it/s]
15674129it [00:09, 1676921.72it/s]
15846033it [00:09, 1689301.37it/s]
16015051it [00:09, 1688756.47it/s]
16186638it [00:09, 1696829.20it/s]
16357983it [00:09, 1701778.92it/s]
16528594it [00:09, 1703055.84it/s]
16699486it [00:10, 1704791.60it/s]
16870452it [00:10, 1706236.06it/s]
17042217it [00:10, 1709644.12it/s]
17213191it [00:10, 1705893.36it/s][W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
17383789it [00:10, 1701870.76it/s]
17432438it [00:10, 1663276.12it/s]
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
tabix -s 1 -b 2 -e 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz
seqtool depth_hist -nbins 500 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.depth_hist.tsv
cat 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.depth_hist.tsv | plot -x 0 -xlabel Depth -ylabel "Number of Sites" -xlim 0 2000 -title 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md -output 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.depth.png
[W::hts_idx_load3] The index file is older than the data file: Homo_sapiens_assembly38.dbsnp138.vcf.gz.tbi
seqtool read_qc -bed xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity.bed -n_threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/read_qc.xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity.bed.tsv
seqtool reads_per_target -bed xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity.bed -threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/reads_per_target.xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity.bed.tsv
seqtool depth_stats -bed xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity.bed 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/depth_stats.xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity.bed.tsv
seqtool gc_bias -reference Homo_sapiens_assembly38.fasta -bed xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity.bed -bin_size 0.05 -threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias_qc.xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity.bed.tsv
seqtool read_qc -bed xgen-exome-hyb-panel-v2_AND_altera_v3_probes_short_hg38_AND_cdx_spike_custom30off4x_hg38_AND_foresight_clarity_footprint.bed -n_threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/read_qc.xgen-exome-hyb-panel-v2_AND_altera_v3_probes_short_hg38_AND_cdx_spike_custom30off4x_hg38_AND_foresight_clarity_footprint.bed.tsv
seqtool reads_per_target -bed xgen-exome-hyb-panel-v2_AND_altera_v3_probes_short_hg38_AND_cdx_spike_custom30off4x_hg38_AND_foresight_clarity_footprint.bed -threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/reads_per_target.xgen-exome-hyb-panel-v2_AND_altera_v3_probes_short_hg38_AND_cdx_spike_custom30off4x_hg38_AND_foresight_clarity_footprint.bed.tsv
seqtool depth_stats -bed xgen-exome-hyb-panel-v2_AND_altera_v3_probes_short_hg38_AND_cdx_spike_custom30off4x_hg38_AND_foresight_clarity_footprint.bed 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/depth_stats.xgen-exome-hyb-panel-v2_AND_altera_v3_probes_short_hg38_AND_cdx_spike_custom30off4x_hg38_AND_foresight_clarity_footprint.bed.tsv
seqtool gc_bias -reference Homo_sapiens_assembly38.fasta -bed xgen-exome-hyb-panel-v2_AND_altera_v3_probes_short_hg38_AND_cdx_spike_custom30off4x_hg38_AND_foresight_clarity_footprint.bed -bin_size 0.05 -threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias_qc.xgen-exome-hyb-panel-v2_AND_altera_v3_probes_short_hg38_AND_cdx_spike_custom30off4x_hg38_AND_foresight_clarity_footprint.bed.tsv
seqtool read_qc -bed altera_v3_targets_postQC_hg38.bed -n_threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/read_qc.altera_v3_targets_postQC_hg38.bed.tsv
seqtool reads_per_target -bed altera_v3_targets_postQC_hg38.bed -threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/reads_per_target.altera_v3_targets_postQC_hg38.bed.tsv
seqtool depth_stats -bed altera_v3_targets_postQC_hg38.bed 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/depth_stats.altera_v3_targets_postQC_hg38.bed.tsv
seqtool gc_bias -reference Homo_sapiens_assembly38.fasta -bed altera_v3_targets_postQC_hg38.bed -bin_size 0.05 -threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias_qc.altera_v3_targets_postQC_hg38.bed.tsv
seqtool read_qc -bed altera_v3_targets_coding_postQC_hg38.bed -n_threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/read_qc.altera_v3_targets_coding_postQC_hg38.bed.tsv
seqtool reads_per_target -bed altera_v3_targets_coding_postQC_hg38.bed -threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/reads_per_target.altera_v3_targets_coding_postQC_hg38.bed.tsv
seqtool depth_stats -bed altera_v3_targets_coding_postQC_hg38.bed 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/depth_stats.altera_v3_targets_coding_postQC_hg38.bed.tsv
seqtool gc_bias -reference Homo_sapiens_assembly38.fasta -bed altera_v3_targets_coding_postQC_hg38.bed -bin_size 0.05 -threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias_qc.altera_v3_targets_coding_postQC_hg38.bed.tsv
seqtool read_qc -bed altera_v3_targets_noncoding_postQC_hg38.bed -n_threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/read_qc.altera_v3_targets_noncoding_postQC_hg38.bed.tsv
seqtool reads_per_target -bed altera_v3_targets_noncoding_postQC_hg38.bed -threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/reads_per_target.altera_v3_targets_noncoding_postQC_hg38.bed.tsv
seqtool depth_stats -bed altera_v3_targets_noncoding_postQC_hg38.bed 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/depth_stats.altera_v3_targets_noncoding_postQC_hg38.bed.tsv
seqtool gc_bias -reference Homo_sapiens_assembly38.fasta -bed altera_v3_targets_noncoding_postQC_hg38.bed -bin_size 0.05 -threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias_qc.altera_v3_targets_noncoding_postQC_hg38.bed.tsv
seqtool read_qc -bed xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed -n_threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/read_qc.xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed.tsv
seqtool reads_per_target -bed xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed -threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/reads_per_target.xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed.tsv
seqtool depth_stats -bed xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/depth_stats.xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed.tsv
seqtool gc_bias -reference Homo_sapiens_assembly38.fasta -bed xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed -bin_size 0.05 -threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias_qc.xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed.tsv
seqtool read_qc -bed FINAL_hotspot_list_07302025_hg38_sorted.bed -n_threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/read_qc.FINAL_hotspot_list_07302025_hg38_sorted.bed.tsv
seqtool reads_per_target -bed FINAL_hotspot_list_07302025_hg38_sorted.bed -threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/reads_per_target.FINAL_hotspot_list_07302025_hg38_sorted.bed.tsv
seqtool depth_stats -bed FINAL_hotspot_list_07302025_hg38_sorted.bed 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/depth_stats.FINAL_hotspot_list_07302025_hg38_sorted.bed.tsv
seqtool gc_bias -reference Homo_sapiens_assembly38.fasta -bed FINAL_hotspot_list_07302025_hg38_sorted.bed -bin_size 0.05 -threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias_qc.FINAL_hotspot_list_07302025_hg38_sorted.bed.tsv
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
seqtool read_qc -bed foresight_clarity_targets_hg38.bed -n_threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/read_qc.foresight_clarity_targets_hg38.bed.tsv
seqtool reads_per_target -bed foresight_clarity_targets_hg38.bed -threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/reads_per_target.foresight_clarity_targets_hg38.bed.tsv
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
seqtool depth_stats -bed foresight_clarity_targets_hg38.bed 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/depth_stats.foresight_clarity_targets_hg38.bed.tsv
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
seqtool gc_bias -reference Homo_sapiens_assembly38.fasta -bed foresight_clarity_targets_hg38.bed -bin_size 0.05 -threads 1 -bam_threads 2 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias_qc.foresight_clarity_targets_hg38.bed.tsv
[W::hts_idx_load3] The index file is older than the data file: 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
cat 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias_qc.FINAL_hotspot_list_07302025_hg38_sorted.bed.tsv | plot -x 0 -header infer -ylim 0 2 -title "1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "FINAL_hotspot_list_07302025_hg38_sorted.bed" -output "1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias.FINAL_hotspot_list_07302025_hg38_sorted.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
cat 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias_qc.altera_v3_targets_coding_postQC_hg38.bed.tsv | plot -x 0 -header infer -ylim 0 2 -title "1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "altera_v3_targets_coding_postQC_hg38.bed" -output "1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias.altera_v3_targets_coding_postQC_hg38.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
cat 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias_qc.altera_v3_targets_noncoding_postQC_hg38.bed.tsv | plot -x 0 -header infer -ylim 0 2 -title "1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "altera_v3_targets_noncoding_postQC_hg38.bed" -output "1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias.altera_v3_targets_noncoding_postQC_hg38.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
cat 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias_qc.foresight_clarity_targets_hg38.bed.tsv | plot -x 0 -header infer -ylim 0 2 -title "1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "foresight_clarity_targets_hg38.bed" -output "1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias.foresight_clarity_targets_hg38.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
cat 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias_qc.altera_v3_targets_postQC_hg38.bed.tsv | plot -x 0 -header infer -ylim 0 2 -title "1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "altera_v3_targets_postQC_hg38.bed" -output "1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias.altera_v3_targets_postQC_hg38.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
cat 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias_qc.xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed.tsv | plot -x 0 -header infer -ylim 0 2 -title "1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed" -output "1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias.xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
cat 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias_qc.xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity.bed.tsv | plot -x 0 -header infer -ylim 0 2 -title "1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity.bed" -output "1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias.xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
cat 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias_qc.xgen-exome-hyb-panel-v2_AND_altera_v3_probes_short_hg38_AND_cdx_spike_custom30off4x_hg38_AND_foresight_clarity_footprint.bed.tsv | plot -x 0 -header infer -ylim 0 2 -title "1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "xgen-exome-hyb-panel-v2_AND_altera_v3_probes_short_hg38_AND_cdx_spike_custom30off4x_hg38_AND_foresight_clarity_footprint.bed" -output "1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias.xgen-exome-hyb-panel-v2_AND_altera_v3_probes_short_hg38_AND_cdx_spike_custom30off4x_hg38_AND_foresight_clarity_footprint.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
echo `date` > 1136_4ZB-T1-TDNA-1_A23MHWWLT4_1/1136_4ZB-T1-TDNA-1_A23MHWWLT4_1.md/qc_complete.txt