.command.log - Workdir - DAQ:BAM_POSTDEDUP_QC:SEQTOOL (659_cAd-T1-TDNA-1_A23MHWWLT4_1)

File Info

Filename: .command.log
Full Path: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/e5/d1c97bf0c509028c5a4a859dbc8c0f/.command.log
Size: 31.7 KB
Task: DAQ:BAM_POSTDEDUP_QC:SEQTOOL (659_cAd-T1-TDNA-1_A23MHWWLT4_1)
Attempt
Content

  Downloading: s3://natera-platform-sandbox/pipeline-resources/Homo_sapiens_assembly38.dbsnp138.vcf.gz
  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/beds/altera_subpanels/hg38/altera_v3_targets_postQC_hg38.bed
  Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/1b/ebf8c638462d79a5888f5503bbb6ae/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
  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/1b/ebf8c638462d79a5888f5503bbb6ae/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam
  Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/e5/d1c97bf0c509028c5a4a859dbc8c0f/.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/e5/d1c97bf0c509028c5a4a859dbc8c0f/.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_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/full_bed.bed


seqtool depth -bam 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam -bed 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/full_bed.bed -threads 4 -bam_threads 2 | bgzip -c > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz
saturation_plot       -bams 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam       -region chr20       -output 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/saturation_plot.png       -threads 1 
seqtool saturation_metric       -at 1000000       -region chr20       -bam_threads 1       659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/saturation_qc.tsv 
fragment_size       -bams 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam       -region chr20       -output_prefix 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.mpileup && contamination -input 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.mpileup -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/contamination.tsv 
seqtool noise   -bam 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam   -ref Homo_sapiens_assembly38.fasta   -region chr20   -include_bed 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/full_bed.bed   -exclude_bed dbsnp138_hg38_common.bed   -threads 1 > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.noise_summary.tsv 
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai

0it [00:00, ?it/s]
150009it [00:00, 1500027.06it/s]
314602it [00:00, 1585827.26it/s]
481415it [00:00, 1623385.50it/s]
647482it [00:00, 1637985.08it/s]
816838it [00:00, 1658005.24it/s]
983361it [00:00, 1660443.24it/s]
1150515it [00:00, 1664050.66it/s]
1316921it [00:00, 1640899.32it/s]
1483890it [00:00, 1649797.71it/s]
1648925it [00:01, 1620507.53it/s]
1812803it [00:01, 1625905.72it/s]
1986511it [00:01, 1659304.68it/s]
2158284it [00:01, 1676838.93it/s]
2327079it [00:01, 1680038.65it/s]
2496756it [00:01, 1684896.25it/s]
2665294it [00:01, 1622855.41it/s]
2829760it [00:01, 1629213.13it/s]
2993056it [00:01, 1625016.86it/s]
3155814it [00:01, 1620711.21it/s]
3318061it [00:02, 1620611.19it/s]
3480245it [00:02, 1609654.77it/s]
3644887it [00:02, 1620445.82it/s]
3812004it [00:02, 1635394.93it/s]
3976421it [00:02, 1637983.77it/s]
4140268it [00:02, 1595399.08it/s]
4300212it [00:02, 1596566.00it/s]
4462514it [00:02, 1604364.92it/s]
4631067it [00:02, 1628440.58it/s]
4794857it [00:02, 1631105.04it/s]
4958058it [00:03, 1626118.22it/s]
5120734it [00:03, 1534139.18it/s]
5275227it [00:03, 1523198.21it/s]
5428277it [00:03, 1518866.21it/s]
5595215it [00:03, 1562660.06it/s]
5766244it [00:03, 1606001.61it/s]
5936213it [00:03, 1633671.29it/s]
6103796it [00:03, 1646168.01it/s]
6271725it [00:03, 1656003.03it/s]
6439029it [00:03, 1661067.03it/s]
6609744it [00:04, 1674761.28it/s]
6780799it [00:04, 1685332.56it/s]
6952707it [00:04, 1695368.30it/s]
7125543it [00:04, 1705216.61it/s]
7299545it [00:04, 1715631.09it/s]
7471135it [00:04, 1702957.99it/s]
7641907it [00:04, 1704283.21it/s]
7812364it [00:04, 1701934.84it/s]
7982578it [00:04, 1695424.94it/s]
8152138it [00:04, 1687845.69it/s]
8320939it [00:05, 1682536.21it/s]
8489205it [00:05, 1674798.20it/s]
8656696it [00:05, 1673696.87it/s]
8824073it [00:05, 1672680.00it/s]
8991346it [00:05, 1669599.46it/s]
9158309it [00:05, 1669266.79it/s]
9326415it [00:05, 1672772.91it/s]
9495180it [00:05, 1677209.07it/s]
9662904it [00:05, 1676026.83it/s]
9830509it [00:05, 1673228.91it/s]
9997834it [00:06, 1670923.76it/s]
10167147it [00:06, 1677343.77it/s]
10334884it [00:06, 1676180.63it/s]
10502894it [00:06, 1677335.61it/s]
10674393it [00:06, 1688599.71it/s]
10843256it [00:06, 1683434.02it/s]
11011605it [00:06, 1673019.38it/s]
11179403it [00:06, 1674386.81it/s]
11347846it [00:06, 1677365.31it/s]done loading BEDs.

11515878it [00:06, 1678116.61it/s]
11683697it [00:07, 1620249.68it/s]
11854560it [00:07, 1646056.59it/s]
12026195it [00:07, 1666740.06it/s]
12196299it [00:07, 1676876.45it/s]
12365127it [00:07, 1680183.58it/s]
12533303it [00:07, 1675038.80it/s]
12700918it [00:07, 1659735.22it/s]
12871488it [00:07, 1673348.90it/s]
13044210it [00:07, 1689362.89it/s]
13217460it [00:07, 1702117.12it/s]
13387730it [00:08, 1701700.16it/s]
13560428it [00:08, 1709243.38it/s]
13731611it [00:08, 1709999.90it/s]
13902633it [00:08, 1679164.43it/s]
14070686it [00:08, 1666841.91it/s]
14237471it [00:08, 1657490.77it/s]
14404996it [00:08, 1662657.11it/s]
14572645it [00:08, 1666747.15it/s]
14740165it [00:08, 1669198.80it/s]
14908633it [00:09, 1673794.99it/s]
15077181it [00:09, 1677262.13it/s]
15244926it [00:09, 1659758.87it/s][W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai

15410953it [00:09, 1659407.92it/s]
15579300it [00:09, 1666559.21it/s]
15747194it [00:09, 1670212.69it/s]
15915005it [00:09, 1672556.73it/s]
16082277it [00:09, 1672552.36it/s]
16249544it [00:09, 1670843.72it/s]
16420654it [00:09, 1682880.44it/s]
16589679it [00:10, 1685070.82it/s]
16760809it [00:10, 1692838.37it/s]
16930684it [00:10, 1694437.93it/s]
17100155it [00:10, 1694503.03it/s]
17269608it [00:10, 1676762.77it/s]
17437346it [00:10, 1676861.31it/s]
17605063it [00:10, 1676032.22it/s]
17773395it [00:10, 1678190.81it/s]
17941230it [00:10, 1672278.25it/s]
18109109it [00:10, 1674205.51it/s]
18280114it [00:11, 1684908.26it/s]
18450996it [00:11, 1691938.21it/s]
18621295it [00:11, 1695233.12it/s]
18790826it [00:11, 1688512.00it/s]
18959688it [00:11, 1678517.68it/s]
19127558it [00:11, 1652789.96it/s]
19292923it [00:11, 1636914.41it/s]
19462202it [00:11, 1653330.18it/s]
19630273it [00:11, 1661427.01it/s]
19801706it [00:11, 1677074.93it/s]
19970552it [00:12, 1680449.50it/s]
20138641it [00:12, 1623701.19it/s]
20301437it [00:12, 1605211.60it/s]
20471141it [00:12, 1631989.70it/s]
20639288it [00:12, 1646530.43it/s]
20805522it [00:12, 1651086.71it/s]
20970798it [00:12, 1650602.19it/s]
21137655it [00:12, 1655938.98it/s]
21305847it [00:12, 1663548.80it/s]
21472266it [00:12, 1661765.51it/s]
21640968it [00:13, 1669302.66it/s]
21812487it [00:13, 1683022.51it/s]
21980816it [00:13, 1681638.12it/s]
22148999it [00:13, 1669362.18it/s]
22317267it [00:13, 1673187.44it/s]
22484611it [00:13, 1642984.74it/s]
22650622it [00:13, 1647961.10it/s]
22818125it [00:13, 1655975.98it/s]
22983846it [00:13, 1656249.71it/s]
23151793it [00:13, 1663125.43it/s]
23319346it [00:14, 1666812.48it/s]
23486595it [00:14, 1668498.42it/s]
23654639it [00:14, 1671992.07it/s]
23825145it [00:14, 1681882.30it/s]
23993346it [00:14, 1678177.31it/s]
24161193it [00:14, 1678249.83it/s]
24260568it [00:14, 1660053.19it/s]


[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai




[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai


[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai


tabix -s 1 -b 2 -e 2 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz


seqtool depth_hist     -nbins 500     659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz     > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.depth_hist.tsv


cat 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.depth_hist.tsv | plot         -x 0         -xlabel Depth -ylabel "Number of Sites"         -xlim 0 2000         -title 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md         -output 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam  > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam  > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam  > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam  > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam  > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam  > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam  > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
seqtool read_qc    -bed foresight_clarity_targets_hg38.bed     -n_threads 1     -bam_threads 2     659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
seqtool reads_per_target     -bed foresight_clarity_targets_hg38.bed     -threads 1     -bam_threads 2     659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
seqtool depth_stats  -bed foresight_clarity_targets_hg38.bed  659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.depth.gz> 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai
[W::hts_idx_load3] The index file is older than the data file: 659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam  > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md.bam.bai


















cat 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "FINAL_hotspot_list_07302025_hg38_sorted.bed" -output "659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "foresight_clarity_targets_hg38.bed" -output "659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias.foresight_clarity_targets_hg38.bed.png"  -ylabel "Fold enrichment / depletion"  -xlabel "GC Quantile" 








cat 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "altera_v3_targets_coding_postQC_hg38.bed" -output "659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "altera_v3_targets_postQC_hg38.bed" -output "659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/gc_bias.altera_v3_targets_postQC_hg38.bed.png"  -ylabel "Fold enrichment / depletion"  -xlabel "GC Quantile" 




cat 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "altera_v3_targets_noncoding_postQC_hg38.bed" -output "659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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" 




cat 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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_cAd-T1-TDNA-1_A23MHWWLT4_1.md" -suptitle "xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed" -output "659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-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" 


















echo `date` > 659_cAd-T1-TDNA-1_A23MHWWLT4_1/659_cAd-T1-TDNA-1_A23MHWWLT4_1.md/qc_complete.txt