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_coding_postQC_hg38.bed
Downloading: s3://natera-platform-sandbox/pipeline-resources/beds/altera_subpanels/hg38/altera_v3_targets_postQC_hg38.bed
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/69/e2109afe678e2520271a463d7be2f0/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
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/ee/eb0c2fa8972d8e4e1c745bb8cac313/.command.sh
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/69/e2109afe678e2520271a463d7be2f0/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam
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/ee/eb0c2fa8972d8e4e1c745bb8cac313/.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-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_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/full_bed.bed
seqtool depth -bam 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam -bed 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/full_bed.bed -threads 4 -bam_threads 2 | bgzip -c > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz
saturation_plot -bams 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam -region chr20 -output 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/saturation_plot.png -threads 1
seqtool saturation_metric -at 1000000 -region chr20 -bam_threads 1 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/saturation_qc.tsv
fragment_size -bams 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam -region chr20 -output_prefix 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.mpileup && contamination -input 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.mpileup -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/contamination.tsv
seqtool noise -bam 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam -ref Homo_sapiens_assembly38.fasta -region chr20 -include_bed 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/full_bed.bed -exclude_bed dbsnp138_hg38_common.bed -threads 1 > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.noise_summary.tsv
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
0it [00:00, ?it/s][W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
145422it [00:00, 1454152.05it/s]
306550it [00:00, 1546542.47it/s]
469691it [00:00, 1585265.22it/s]
632468it [00:00, 1601971.94it/s]
795693it [00:00, 1612869.69it/s]
959011it [00:00, 1619751.84it/s]
1124097it [00:00, 1629835.60it/s]
1289596it [00:00, 1637657.14it/s]
1453362it [00:00, 1631820.01it/s]
1618357it [00:01, 1637329.99it/s]
1782789it [00:01, 1639446.81it/s]
1946738it [00:01, 1614091.94it/s]
2108234it [00:01, 1613200.79it/s]
2269614it [00:01, 1612302.20it/s]
2433055it [00:01, 1618881.33it/s]
2594976it [00:01, 1559808.43it/s]
2757756it [00:01, 1579644.04it/s]
2924144it [00:01, 1604419.91it/s]
3089969it [00:01, 1620169.29it/s]
3252934it [00:02, 1622644.54it/s]
3417732it [00:02, 1630175.68it/s]
3580965it [00:02, 1630801.95it/s]
3747088it [00:02, 1639882.81it/s]
3911143it [00:02, 1567406.25it/s]
4070793it [00:02, 1575827.81it/s]
4231044it [00:02, 1583629.21it/s]
4389795it [00:02, 1575286.64it/s]
4547593it [00:02, 1564089.34it/s]
4712072it [00:02, 1587887.44it/s]
4874380it [00:03, 1598137.41it/s]
5038728it [00:03, 1611609.85it/s]
5199994it [00:03, 1571802.19it/s]
5364452it [00:03, 1593149.92it/s]
5527885it [00:03, 1605295.29it/s]
5689280it [00:03, 1607845.06it/s]
5850197it [00:03, 1606512.26it/s]
6010940it [00:03, 1566067.35it/s]
6168559it [00:03, 1569017.85it/s]
6328996it [00:03, 1579284.44it/s]
6487077it [00:04, 1572779.62it/s]
6650212it [00:04, 1590083.66it/s]
6812853it [00:04, 1600879.62it/s]
6973941it [00:04, 1603847.00it/s]
7134379it [00:04, 1517433.86it/s]
7287102it [00:04, 1478855.50it/s]
7435763it [00:04, 1466339.33it/s]
7598806it [00:04, 1513538.75it/s]
7764685it [00:04, 1555834.40it/s]
7932886it [00:04, 1592895.85it/s]
8097782it [00:05, 1609454.76it/s]
8262530it [00:05, 1620729.79it/s]
8426660it [00:05, 1626841.12it/s]
8589906it [00:05, 1628499.10it/s]
8757881it [00:05, 1643751.57it/s]
8925159it [00:05, 1652419.97it/s]
9093075it [00:05, 1660412.95it/s]
9262248it [00:05, 1669780.22it/s]
9432828it [00:05, 1680560.39it/s]
9600910it [00:05, 1667799.08it/s]
9768815it [00:06, 1671137.46it/s]
9935959it [00:06, 1670161.14it/s]
10102996it [00:06, 1666478.21it/s]
10269659it [00:06, 1665782.27it/s]
10436248it [00:06, 1661843.89it/s]
10602441it [00:06, 1657668.70it/s]
10768215it [00:06, 1647331.26it/s]done loading BEDs.
10932964it [00:06, 1646072.54it/s]
11097815it [00:06, 1646784.78it/s]
11262965it [00:06, 1648114.22it/s]
11427783it [00:07, 1642747.39it/s]
11593071it [00:07, 1645756.96it/s]
11758998it [00:07, 1649752.01it/s]
11925628it [00:07, 1654654.28it/s]
12091099it [00:07, 1651490.49it/s]
12256253it [00:07, 1649110.04it/s]
12421168it [00:07, 1647964.57it/s]
12587697it [00:07, 1653133.76it/s]
12753014it [00:07, 1652170.62it/s]
12918233it [00:08, 1649851.14it/s]
13086983it [00:08, 1661107.18it/s]
13256046it [00:08, 1669827.41it/s]
13423033it [00:08, 1659042.09it/s]
13588955it [00:08, 1655131.93it/s]
13755246it [00:08, 1657436.76it/s]
13921112it [00:08, 1657787.46it/s]
14086898it [00:08, 1657521.13it/s]
14252655it [00:08, 1599933.55it/s]
14419978it [00:08, 1621336.50it/s]
14586777it [00:09, 1635060.61it/s]
14752451it [00:09, 1641360.49it/s]
14916788it [00:09, 1639908.08it/s]
15083575it [00:09, 1648220.94it/s]
15248501it [00:09, 1635845.45it/s]
15417918it [00:09, 1653169.98it/s]
15588563it [00:09, 1668931.57it/s]
15758543it [00:09, 1678134.46it/s]
15926746it [00:09, 1679281.97it/s]
16096775it [00:09, 1685554.11it/s]
16265356it [00:10, 1684522.08it/s]
16433827it [00:10, 1674765.41it/s]
16601328it [00:10, 1651621.06it/s]
16766569it [00:10, 1638114.43it/s]
16930444it [00:10, 1630965.00it/s]
17094309it [00:10, 1633221.87it/s]
17258584it [00:10, 1635983.50it/s]
17423174it [00:10, 1638915.27it/s]
17588486it [00:10, 1643060.76it/s]
17754204it [00:10, 1647267.28it/s]
17919391it [00:11, 1648529.98it/s]
18084253it [00:11, 1634090.59it/s]
18249765it [00:11, 1640278.88it/s]
18416093it [00:11, 1647120.38it/s]
18582048it [00:11, 1650819.64it/s]
18747845it [00:11, 1652905.47it/s]
18913149it [00:11, 1648478.57it/s]
19078008it [00:11, 1642845.46it/s]
19245925it [00:11, 1653662.28it/s]
19412482it [00:11, 1657183.27it/s]
19579891it [00:12, 1662219.20it/s]
19746121it [00:12, 1660894.16it/s]
19912216it [00:12, 1654440.01it/s]
20077669it [00:12, 1638728.11it/s]
20243255it [00:12, 1643803.45it/s]
20408235it [00:12, 1645573.44it/s]
20572814it [00:12, 1640022.50it/s]
20736834it [00:12, 1638252.45it/s]
20902635it [00:12, 1644132.22it/s]
21070699it [00:12, 1655030.89it/s]
21238961it [00:13, 1663269.75it/s]
21406709it [00:13, 1667507.23it/s]
21573468it [00:13, 1664234.30it/s]
21739898it [00:13, 1657404.48it/s]
21905648it [00:13, 1650522.56it/s]
22070711it [00:13, 1606931.64it/s]
22237156it [00:13, 1623777.92it/s]
22403237it [00:13, 1634689.68it/s]
22571357it [00:13, 1648465.38it/s]
22736718it [00:13, 1649973.53it/s]
22901807it [00:14, 1574834.16it/s]
23063435it [00:14, 1586799.71it/s]
23230393it [00:14, 1610975.04it/s]
23396374it [00:14, 1625313.60it/s]
23560361it [00:14, 1629607.91it/s]
23725169it [00:14, 1635073.54it/s]
23889762it [00:14, 1638225.37it/s]
24055153it [00:14, 1642835.34it/s]
24220495it [00:14, 1645967.13it/s]
24385158it [00:14, 1633821.64it/s]
24553086it [00:15, 1647329.18it/s]
24720110it [00:15, 1654151.08it/s]
24885567it [00:15, 1640286.50it/s]
25049832it [00:15, 1640859.78it/s]
25213955it [00:15, 1616989.90it/s]
25375746it [00:15, 1616518.26it/s]
25540312it [00:15, 1625152.88it/s]
25705374it [00:15, 1632715.49it/s]
25870476it [00:15, 1638159.22it/s]
26034325it [00:15, 1631360.74it/s][W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
26197489it [00:16, 1631090.65it/s]
26362208it [00:16, 1635881.81it/s]
26527443it [00:16, 1640690.67it/s]
26691992it [00:16, 1642104.59it/s]
26856212it [00:16, 1641847.51it/s]
27021786it [00:16, 1645987.55it/s]
27189029it [00:16, 1653895.56it/s]
27355206it [00:16, 1656236.91it/s]
27520833it [00:16, 1652225.02it/s]
27686152it [00:16, 1652457.79it/s]
27821163it [00:17, 1629763.21it/s]
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
tabix -s 1 -b 2 -e 2 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz
seqtool depth_hist -nbins 500 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.depth_hist.tsv
cat 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.depth_hist.tsv | plot -x 0 -xlabel Depth -ylabel "Number of Sites" -xlim 0 2000 -title 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md -output 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
seqtool read_qc -bed foresight_clarity_targets_hg38.bed -n_threads 1 -bam_threads 2 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
seqtool reads_per_target -bed foresight_clarity_targets_hg38.bed -threads 1 -bam_threads 2 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
seqtool depth_stats -bed foresight_clarity_targets_hg38.bed 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
cat 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "FINAL_hotspot_list_07302025_hg38_sorted.bed" -output "659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "foresight_clarity_targets_hg38.bed" -output "659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias.foresight_clarity_targets_hg38.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
cat 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "altera_v3_targets_coding_postQC_hg38.bed" -output "659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "altera_v3_targets_noncoding_postQC_hg38.bed" -output "659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "altera_v3_targets_postQC_hg38.bed" -output "659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias.altera_v3_targets_postQC_hg38.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
cat 659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed" -output "659_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-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_cfy-T1-TDNA-1_A23MHWWLT4_1/659_cfy-T1-TDNA-1_A23MHWWLT4_1.md/qc_complete.txt