Downloading: s3://natera-platform-sandbox/pipeline-resources/AIH/rna/GRCh38/starfusion/ctat_genome_lib_build_dir/ref_genome.fa
Downloading: s3://natera-platform-sandbox/pipeline-resources/Homo_sapiens_assembly38.dbsnp138.vcf.gz
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/e5/547fa35c3aaeeee82b3d462e0d21be/.command.sh
Downloading: s3://natera-platform-sandbox/pipeline-resources/AIH/rna/GRCh38/starfusion/ctat_genome_lib_build_dir/ref_genome.dict
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/7e/a2b361272ec8e9bd366f3c859e0021/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
Downloading: s3://natera-platform-sandbox/pipeline-resources/AIH/rna/GRCh38/gatk4/dbsnp_hg38_contam_sites.bed
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/e5/547fa35c3aaeeee82b3d462e0d21be/.command.run
Downloading: s3://natera-platform-sandbox/pipeline-resources/AIH/rna/GRCh38/starfusion/ctat_genome_lib_build_dir/ref_genome.fa.fai
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/7e/a2b361272ec8e9bd366f3c859e0021/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bam
Downloading: s3://natera-platform-sandbox/pipeline-resources/beds/twist/Twist_rna_exome_AND_fusion_hg38.bed
Downloading: s3://natera-platform-sandbox/pipeline-resources/Homo_sapiens_assembly38.dbsnp138.vcf.gz.tbi
==> STAGING COMPLETE (11 inputs)
+ mkdir aih-tih-sc-150e73-R1_A23YTGFLT4_1
+ all_qc -bam aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bam -beds Twist_rna_exome_AND_fusion_hg38.bed -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz -contamination_bed dbsnp_hg38_contam_sites.bed -output_folder aih-tih-sc-150e73-R1_A23YTGFLT4_1 -reference ref_genome.fa -region chr20
mkdir -p aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/
cat Twist_rna_exome_AND_fusion_hg38.bed | cut -f 1-3 | bedtools sort -g ref_genome.fa.fai | bedtools merge > aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/full_bed.bed
seqtool depth -bam aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bam -bed aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/full_bed.bed -threads 4 -bam_threads 2 | bgzip -c > aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.depth.gz
saturation_plot -bams aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bam -region chr20 -output aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/saturation_plot.png -threads 1
seqtool saturation_metric -at 1000000 -region chr20 -bam_threads 1 aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bam > aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/saturation_qc.tsv
fragment_size -bams aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bam -region chr20 -output_prefix aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/fragsize -threads 1
samtools mpileup -B --min-MQ 60 --min-BQ 30 -f ref_genome.fa -l dbsnp_hg38_contam_sites.bed aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bam > aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.mpileup && contamination -input aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.mpileup -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz > aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/contamination.tsv
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[mpileup] 1 samples in 1 input files
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
/opt/conda/lib/python3.9/site-packages/pyranges/__init__.py:8: UserWarning: pkg_resources is deprecated as an API. See https://setuptools.pypa.io/en/latest/pkg_resources.html. The pkg_resources package is slated for removal as early as 2025-11-30. Refrain from using this package or pin to Setuptools<81.
import pkg_resources
/opt/conda/lib/python3.9/site-packages/pyranges/__init__.py:8: UserWarning: pkg_resources is deprecated as an API. See https://setuptools.pypa.io/en/latest/pkg_resources.html. The pkg_resources package is slated for removal as early as 2025-11-30. Refrain from using this package or pin to Setuptools<81.
import pkg_resources
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
0it [00:00, ?it/s]
58994it [00:00, 589873.06it/s]
122903it [00:00, 618793.69it/s]
186713it [00:00, 627589.37it/s]
249472it [00:00, 620259.20it/s]
313256it [00:00, 626534.32it/s]
377978it [00:00, 633517.17it/s]
441577it [00:00, 634308.76it/s]
505835it [00:00, 636923.46it/s]
570156it [00:00, 638871.48it/s]
635296it [00:01, 642723.29it/s]
700244it [00:01, 644777.04it/s]
765142it [00:01, 646042.44it/s]
829749it [00:01, 644385.95it/s]
894575it [00:01, 645540.39it/s]
959131it [00:01, 643674.37it/s]
1023816it [00:01, 644614.12it/s]
1088477it [00:01, 645198.87it/s]
1155627it [00:01, 653077.15it/s]
1222440it [00:01, 657578.41it/s]
1289943it [00:02, 662799.04it/s]
1356225it [00:02, 658213.41it/s]
1423206it [00:02, 661663.43it/s]
1490597it [00:02, 665310.33it/s]
1558265it [00:02, 668698.41it/s]
1625678it [00:02, 670310.94it/s]
1692714it [00:02, 658326.98it/s]
1760677it [00:02, 664619.54it/s]
1827963it [00:02, 667055.72it/s]
1895411it [00:02, 669256.49it/s]
1962361it [00:03, 669099.32it/s]
2030014it [00:03, 671307.48it/s]
2098487it [00:03, 675310.47it/s]
2166889it [00:03, 677904.91it/s]
2234687it [00:03, 677434.47it/s]
2302436it [00:03, 677260.70it/s]
2370166it [00:03, 676137.02it/s]
2438056it [00:03, 676950.45it/s]
2505754it [00:03, 675386.40it/s]
2573295it [00:03, 650590.88it/s]
2641107it [00:04, 658604.68it/s]
2708434it [00:04, 662903.20it/s]
2775809it [00:04, 666096.88it/s]
2842507it [00:04, 664926.77it/s]
2909317it [00:04, 665857.29it/s]
2976866it [00:04, 668690.87it/s]
3044658it [00:04, 671434.43it/s]
3113087it [00:04, 675165.89it/s]
3182336it [00:04, 680336.02it/s]
3252290it [00:04, 686025.85it/s]
3321391it [00:05, 687387.98it/s]
3390137it [00:05, 684936.09it/s]
3458638it [00:05, 683024.07it/s]
3526946it [00:05, 681320.20it/s]
3595083it [00:05, 680064.42it/s]
3663093it [00:05, 678947.98it/s]
3730990it [00:05, 678244.75it/s]
3799790it [00:05, 681144.37it/s]
3868993it [00:05, 684386.36it/s]
3937650it [00:05, 685027.69it/s]
4006616it [00:06, 686401.01it/s]
4075258it [00:06, 680676.21it/s]
4143337it [00:06, 678930.11it/s]
4211238it [00:06, 677303.59it/s]
4278974it [00:06, 676081.19it/s]
4346827it [00:06, 676743.07it/s]
4415126it [00:06, 678518.47it/s]
4483093it [00:06, 678845.20it/s]
4551646it [00:06, 680831.66it/s]
4620394it [00:06, 682808.34it/s]
4688677it [00:07, 676836.14it/s]
4756374it [00:07, 672634.29it/s]
4824120it [00:07, 674055.02it/s]
4892546it [00:07, 677017.03it/s]
4961954it [00:07, 682096.92it/s]
5030563it [00:07, 683276.53it/s]
5098897it [00:07, 682821.61it/s]
5167184it [00:07, 682499.19it/s]
5235437it [00:07, 655889.74it/s]
5301247it [00:07, 640338.88it/s]
5369173it [00:08, 651563.21it/s]
5436079it [00:08, 656660.37it/s]
5503751it [00:08, 662524.43it/s]
5571556it [00:08, 667106.81it/s]
5638854it [00:08, 668834.49it/s]
5705798it [00:08, 665345.08it/s]
5773501it [00:08, 668807.85it/s]
5841731it [00:08, 672723.54it/s]
5909531it [00:08, 674171.69it/s]
5977473it [00:08, 675727.96it/s]
6045410it [00:09, 676802.89it/s]
6113543it [00:09, 678142.15it/s]
6181365it [00:09, 678146.77it/s]
6249550it [00:09, 679242.00it/s]
6318481it [00:09, 682243.31it/s]
6386954it [00:09, 682975.04it/s]
6455254it [00:09, 681771.79it/s]
6523936it [00:09, 683267.16it/s]
6592265it [00:09, 680175.98it/s]
6660677it [00:09, 681341.17it/s]
6728815it [00:10, 663931.57it/s]
6795305it [00:10, 653364.11it/s]
6860728it [00:10, 643022.78it/s]
6925104it [00:10, 634977.46it/s]
6988655it [00:10, 631065.98it/s]
7051793it [00:10, 626515.05it/s]
7114831it [00:10, 627634.05it/s]
7177610it [00:10, 623698.32it/s]
7239990it [00:10, 620484.06it/s]
7302045it [00:11, 618049.43it/s]
7364273it [00:11, 619288.57it/s]
7426206it [00:11, 618911.59it/s]
7488396it [00:11, 619787.56it/s]
7550377it [00:11, 619368.13it/s]
7613009it [00:11, 621432.90it/s]
7675665it [00:11, 622952.33it/s]
7738806it [00:11, 625472.38it/s]
7802004it [00:11, 627408.73it/s]
7865787it [00:11, 630518.72it/s]
7929520it [00:12, 632543.68it/s]
7993610it [00:12, 635035.87it/s]
8059300it [00:12, 641580.02it/s]
8125889it [00:12, 648857.75it/s]
8194035it [00:12, 658622.95it/s]
8260062it [00:12, 659103.98it/s]
8326679it [00:12, 661210.00it/s]
8392936it [00:12, 661600.16it/s]
8459097it [00:12, 661224.98it/s]
8525453it [00:12, 661912.10it/s]
8591645it [00:13, 661356.19it/s]
8657781it [00:13, 659373.64it/s]
8723857it [00:13, 659772.62it/s]
8790720it [00:13, 662407.02it/s]
8858295it [00:13, 666391.87it/s]
8924936it [00:13, 663619.13it/s]
8996575it [00:13, 679377.78it/s]
9080544it [00:13, 727329.65it/s]
9153295it [00:13, 602549.69it/s]
9217121it [00:14, 573658.83it/s]
9278519it [00:14, 584110.38it/s]
9339154it [00:14, 590116.08it/s]
9399606it [00:14, 567053.28it/s]
9457406it [00:14, 569310.01it/s]
9516056it [00:14, 574145.66it/s]
9574051it [00:14, 575789.12it/s]
9632449it [00:14, 578161.66it/s]
9690563it [00:14, 567889.46it/s]
9747584it [00:14, 548750.07it/s]
9802719it [00:15, 543043.02it/s]
9857196it [00:15, 541229.47it/s]
9914420it [00:15, 550231.88it/s]
9972133it [00:15, 558112.70it/s]
10028948it [00:15, 561067.18it/s]
10086148it [00:15, 564265.94it/s]
10142628it [00:15, 552035.48it/s]
10197920it [00:15, 539463.62it/s]
10251970it [00:15, 537785.32it/s]
10308446it [00:15, 545669.20it/s]
10365196it [00:16, 552093.51it/s]
10423201it [00:16, 560368.28it/s]
10480520it [00:16, 564071.35it/s]
10540200it [00:16, 573744.35it/s]
10597608it [00:16, 565139.54it/s]
10654171it [00:16, 565227.00it/s]
10711808it [00:16, 568525.48it/s]
10769583it [00:16, 571261.49it/s]
10827034it [00:16, 572220.11it/s]
10884272it [00:16, 561786.84it/s]
10940501it [00:17, 543571.17it/s]
10994999it [00:17, 534664.68it/s]
11048570it [00:17, 531573.99it/s]
11103186it [00:17, 535808.71it/s]
11159679it [00:17, 544350.41it/s]
11217452it [00:17, 554204.41it/s]
11274858it [00:17, 560086.60it/s]
11330913it [00:17, 554464.81it/s]
11386401it [00:17, 540629.22it/s]
11440555it [00:18, 538680.15it/s]
11494484it [00:18, 534741.44it/s]
11551364it [00:18, 544738.49it/s]
11609278it [00:18, 554814.03it/s]
11666819it [00:18, 560917.83it/s]
11731567it [00:18, 586690.08it/s]
11819757it [00:18, 674792.42it/s]
11909931it [00:18, 742538.08it/s]
12005235it [00:18, 805509.46it/s]
12099795it [00:18, 847438.96it/s]
12195136it [00:19, 879173.00it/s]
12283097it [00:19, 840783.15it/s]
12371404it [00:19, 853050.66it/s]
12457018it [00:19, 849206.08it/s]
12542151it [00:19, 825781.55it/s]
12633498it [00:19, 851182.13it/s]
12724987it [00:19, 869817.57it/s]
12812201it [00:19, 796851.67it/s]
12893172it [00:19, 732141.23it/s]
12968047it [00:20, 694685.97it/s]
13051142it [00:20, 730695.65it/s]
13137789it [00:20, 767902.56it/s]
13236338it [00:20, 829116.00it/s]
13329868it [00:20, 859439.60it/s]
13419427it [00:20, 869928.52it/s]
13510827it [00:20, 882840.52it/s]
13599686it [00:20, 881862.64it/s]
13688274it [00:20, 876835.23it/s]
13776241it [00:20, 876947.91it/s]
13869635it [00:21, 893850.06it/s]
13962002it [00:21, 902707.01it/s][W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
14054775it [00:21, 910158.22it/s]
14146991it [00:21, 913729.54it/s]
14238425it [00:21, 911022.21it/s]
14329572it [00:21, 907758.67it/s]
14420381it [00:21, 881238.77it/s]
14517405it [00:21, 907270.32it/s]
14616858it [00:21, 932980.14it/s]
14711970it [00:21, 938338.29it/s]
14810515it [00:22, 952293.74it/s]
14908904it [00:22, 961699.85it/s]
15005153it [00:22, 956510.76it/s]
15100945it [00:22, 956911.16it/s]
15197038it [00:22, 958092.80it/s]
15292878it [00:22, 929636.09it/s]
15386040it [00:22, 906043.59it/s]
15476869it [00:22, 887124.22it/s]
15565770it [00:22, 840013.72it/s]
15650284it [00:23, 754389.43it/s]
15727404it [00:23, 709041.85it/s]
15799666it [00:23, 682352.75it/s]
15868780it [00:23, 662278.82it/s]
15935532it [00:23, 647401.44it/s]
16000558it [00:23, 638240.29it/s]
16064524it [00:23, 627361.48it/s]
16127316it [00:23, 623102.57it/s]
16189642it [00:23, 619374.05it/s]
16251575it [00:24, 617497.25it/s]
16313524it [00:24, 617994.51it/s]
16375605it [00:24, 618806.81it/s]
16437484it [00:24, 617764.96it/s]
16499258it [00:24, 614960.47it/s]
16560752it [00:24, 610251.60it/s]
16621779it [00:24, 604904.81it/s]
16682274it [00:24, 603556.63it/s]
16743168it [00:24, 605138.91it/s]
16828258it [00:24, 678196.21it/s]
16922179it [00:25, 755995.10it/s]
17015777it [00:25, 809734.90it/s]
17110248it [00:25, 850085.29it/s]
17202205it [00:25, 870864.71it/s]
17289902it [00:25, 872606.46it/s]
17380211it [00:25, 881725.44it/s]
17447427it [00:25, 681518.28it/s]
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
tabix -s 1 -b 2 -e 2 aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.depth.gz
seqtool depth_hist -nbins 500 aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.depth.gz > aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.depth_hist.tsv
cat aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.depth_hist.tsv | plot -x 0 -xlabel Depth -ylabel "Number of Sites" -xlim 0 2000 -title aih-tih-sc-150e73-R1_A23YTGFLT4_1.md -output aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/aih-tih-sc-150e73-R1_A23YTGFLT4_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 Twist_rna_exome_AND_fusion_hg38.bed -n_threads 1 -bam_threads 2 aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bam > aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/read_qc.Twist_rna_exome_AND_fusion_hg38.bed.tsv
seqtool reads_per_target -bed Twist_rna_exome_AND_fusion_hg38.bed -threads 1 -bam_threads 2 aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bam > aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/reads_per_target.Twist_rna_exome_AND_fusion_hg38.bed.tsv
seqtool depth_stats -bed Twist_rna_exome_AND_fusion_hg38.bed aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.depth.gz> aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/depth_stats.Twist_rna_exome_AND_fusion_hg38.bed.tsv
seqtool gc_bias -reference ref_genome.fa -bed Twist_rna_exome_AND_fusion_hg38.bed -bin_size 0.05 -threads 1 -bam_threads 2 aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bam > aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/gc_bias_qc.Twist_rna_exome_AND_fusion_hg38.bed.tsv
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-150e73-R1_A23YTGFLT4_1.md.bai
cat aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/gc_bias_qc.Twist_rna_exome_AND_fusion_hg38.bed.tsv | plot -x 0 -header infer -ylim 0 2 -title "aih-tih-sc-150e73-R1_A23YTGFLT4_1.md" -suptitle "Twist_rna_exome_AND_fusion_hg38.bed" -output "aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/gc_bias.Twist_rna_exome_AND_fusion_hg38.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
echo `date` > aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md/qc_complete.txt
+ merge_qc -folders aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1.md -output aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1_full_basic_stats.tsv
+ head -n1 aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1_full_basic_stats.tsv
++ basename Twist_rna_exome_AND_fusion_hg38.bed
+ grep -e Twist_rna_exome_AND_fusion_hg38.bed aih-tih-sc-150e73-R1_A23YTGFLT4_1/aih-tih-sc-150e73-R1_A23YTGFLT4_1_full_basic_stats.tsv