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-rnd-pltf-dev-nextflow-scratch-01/work/f7/93124cd57bf207d807bf685f8c8d91/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam
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/4b/dd6156a7e37329bf0eeb565d4c83e4/.command.sh
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/4b/dd6156a7e37329bf0eeb565d4c83e4/.command.run
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/f7/93124cd57bf207d807bf685f8c8d91/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
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-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/full_bed.bed
seqtool depth -bam 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam -bed 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/full_bed.bed -threads 4 -bam_threads 2 | bgzip -c > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz
saturation_plot -bams 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam -region chr20 -output 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/saturation_plot.png -threads 1
seqtool saturation_metric -at 1000000 -region chr20 -bam_threads 1 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/saturation_qc.tsv
fragment_size -bams 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam -region chr20 -output_prefix 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.mpileup && contamination -input 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.mpileup -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/contamination.tsv
seqtool noise -bam 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam -ref Homo_sapiens_assembly38.fasta -region chr20 -include_bed 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/full_bed.bed -exclude_bed dbsnp138_hg38_common.bed -threads 1 > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.noise_summary.tsv
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-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: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-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: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
0it [00:00, ?it/s]
152237it [00:00, 1522255.31it/s]
318972it [00:00, 1607562.80it/s]
490314it [00:00, 1655869.84it/s]
661661it [00:00, 1678587.60it/s]
833974it [00:00, 1694453.41it/s]
1006638it [00:00, 1705375.66it/s]
1182363it [00:00, 1722313.80it/s]
1354595it [00:00, 1705424.84it/s]
1529021it [00:00, 1717442.57it/s]
1700798it [00:01, 1678516.52it/s]
1868844it [00:01, 1659831.42it/s]
2039053it [00:01, 1672397.60it/s]
2206424it [00:01, 1620533.21it/s]
2368850it [00:01, 1612179.56it/s]
2537914it [00:01, 1635176.81it/s]
2715007it [00:01, 1675251.04it/s]
2891202it [00:01, 1700958.83it/s]
3061501it [00:01, 1692509.32it/s]
3230897it [00:01, 1678491.45it/s]
3398861it [00:02, 1595297.13it/s]
3562242it [00:02, 1606248.83it/s]
3723520it [00:02, 1565728.79it/s]
3888324it [00:02, 1589406.56it/s]
4048056it [00:02, 1591696.48it/s]
4218191it [00:02, 1623782.10it/s]
4388201it [00:02, 1646334.18it/s]
4553090it [00:02, 1607845.66it/s]
4722971it [00:02, 1634467.78it/s]
4894842it [00:02, 1659298.05it/s]
5061466it [00:03, 1661283.69it/s]
5227778it [00:03, 1588404.40it/s]
5387380it [00:03, 1587965.74it/s]
5549963it [00:03, 1598995.38it/s]
5710996it [00:03, 1602312.46it/s]
5880260it [00:03, 1629022.39it/s]
6051295it [00:03, 1653168.52it/s]
6216803it [00:03, 1540767.39it/s]
6372506it [00:03, 1494447.74it/s]
6523203it [00:04, 1473375.05it/s]
6695119it [00:04, 1543072.69it/s]
6872019it [00:04, 1608204.88it/s]
7047009it [00:04, 1649430.41it/s]
7220949it [00:04, 1675879.72it/s]
7395206it [00:04, 1695519.54it/s]
7573721it [00:04, 1722131.80it/s]
7752685it [00:04, 1742232.13it/s]
7932656it [00:04, 1759378.59it/s]
8108778it [00:04, 1754475.41it/s]
8287948it [00:05, 1765555.11it/s]
8464600it [00:05, 1763383.71it/s]
8641006it [00:05, 1763411.69it/s]
8817395it [00:05, 1753778.65it/s]
8992817it [00:05, 1749232.40it/s]
9167771it [00:05, 1741046.27it/s]
9343026it [00:05, 1744457.75it/s]
9519300it [00:05, 1749898.53it/s]
9694308it [00:05, 1736946.93it/s]
9870067it [00:05, 1743075.11it/s]
10046885it [00:06, 1750549.94it/s]
10221962it [00:06, 1742423.41it/s]
10396226it [00:06, 1739889.43it/s]
10571895it [00:06, 1744886.43it/s]
10746397it [00:06, 1742955.65it/s]
10923583it [00:06, 1751434.80it/s]
11099320it [00:06, 1753096.47it/s]done loading BEDs.
11274637it [00:06, 1748459.88it/s]
11449514it [00:06, 1748453.39it/s]
11625316it [00:06, 1751298.82it/s]
11800799it [00:07, 1752337.02it/s]
11976036it [00:07, 1665569.52it/s]
12154497it [00:07, 1699937.87it/s]
12330327it [00:07, 1716924.74it/s]
12502591it [00:07, 1718361.54it/s]
12674830it [00:07, 1698519.27it/s]
12854419it [00:07, 1727202.36it/s]
13034775it [00:07, 1749811.59it/s]
13213571it [00:07, 1761098.69it/s]
13394140it [00:07, 1774378.48it/s]
13571700it [00:08, 1768934.87it/s]
13748681it [00:08, 1737290.96it/s]
13922584it [00:08, 1717158.33it/s]
14094445it [00:08, 1700356.33it/s]
14267948it [00:08, 1710483.63it/s]
14442843it [00:08, 1721710.40it/s]
14617058it [00:08, 1727762.83it/s]
14791576it [00:08, 1732931.67it/s]
14964919it [00:08, 1706320.74it/s]
15140690it [00:08, 1721497.01it/s]
15316680it [00:09, 1732886.54it/s]
15492481it [00:09, 1740351.69it/s]
15666580it [00:09, 1726552.41it/s]
15839735it [00:09, 1728022.54it/s]
16018375it [00:09, 1745406.36it/s]
16194241it [00:09, 1749347.56it/s]
16369208it [00:09, 1738245.73it/s]
16543068it [00:09, 1726662.34it/s]
16717080it [00:09, 1730642.53it/s]
16892630it [00:09, 1738033.68it/s][W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
17066459it [00:10, 1729058.11it/s]
17241794it [00:10, 1736141.74it/s]
17421126it [00:10, 1753194.17it/s]
17598167it [00:10, 1758223.83it/s]
17774008it [00:10, 1746345.74it/s]
17948671it [00:10, 1733457.30it/s]
18122051it [00:10, 1671269.98it/s]
18299311it [00:10, 1700672.88it/s]
18474205it [00:10, 1714800.23it/s]
18646009it [00:11, 1632292.26it/s]
18810204it [00:11, 1612270.12it/s]
18986453it [00:11, 1655424.37it/s]
19158210it [00:11, 1673496.13it/s]
19330719it [00:11, 1688489.62it/s]
19505127it [00:11, 1704902.22it/s]
19682648it [00:11, 1725705.81it/s]
19857110it [00:11, 1731316.56it/s]
20035593it [00:11, 1747122.49it/s]
20210428it [00:11, 1741096.45it/s]
20384627it [00:12, 1718719.99it/s]
20556617it [00:12, 1677052.80it/s]
20728794it [00:12, 1690097.21it/s]
20905291it [00:12, 1712127.81it/s]
21080984it [00:12, 1725376.35it/s]
21255063it [00:12, 1729942.81it/s]
21429279it [00:12, 1733543.08it/s]
21605018it [00:12, 1740645.55it/s]
21781719it [00:12, 1748511.61it/s]
21959106it [00:12, 1756083.87it/s]
22135223it [00:13, 1757588.99it/s]
22311591it [00:13, 1759398.44it/s]
22456597it [00:13, 1697587.55it/s]
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
tabix -s 1 -b 2 -e 2 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz
seqtool depth_hist -nbins 500 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.depth_hist.tsv
cat 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.depth_hist.tsv | plot -x 0 -xlabel Depth -ylabel "Number of Sites" -xlim 0 2000 -title 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md -output 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
seqtool read_qc -bed foresight_clarity_targets_hg38.bed -n_threads 1 -bam_threads 2 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/read_qc.foresight_clarity_targets_hg38.bed.tsv
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
seqtool reads_per_target -bed foresight_clarity_targets_hg38.bed -threads 1 -bam_threads 2 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
seqtool depth_stats -bed foresight_clarity_targets_hg38.bed 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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: 659_dm4-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
cat 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 "659_dm4-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "FINAL_hotspot_list_07302025_hg38_sorted.bed" -output "659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias.FINAL_hotspot_list_07302025_hg38_sorted.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
cat 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias_qc.foresight_clarity_targets_hg38.bed.tsv | plot -x 0 -header infer -ylim 0 2 -title "659_dm4-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "foresight_clarity_targets_hg38.bed" -output "659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias.foresight_clarity_targets_hg38.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
cat 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 "659_dm4-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "altera_v3_targets_coding_postQC_hg38.bed" -output "659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias.altera_v3_targets_coding_postQC_hg38.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
cat 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 "659_dm4-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "altera_v3_targets_noncoding_postQC_hg38.bed" -output "659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias.altera_v3_targets_noncoding_postQC_hg38.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
cat 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 "659_dm4-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "altera_v3_targets_postQC_hg38.bed" -output "659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias.altera_v3_targets_postQC_hg38.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
cat 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 "659_dm4-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed" -output "659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 "659_dm4-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 "659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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 "659_dm4-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 "659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-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` > 659_dm4-T1-TDNA-1_A23MHWWLT4_1/659_dm4-T1-TDNA-1_A23MHWWLT4_1.md/qc_complete.txt