Downloading: s3://natera-platform-sandbox/pipeline-resources/Homo_sapiens_assembly38.dbsnp138.vcf.gz
Downloading: s3://natera-platform-sandbox/pipeline-resources/AIH/rna/GRCh38/starfusion/ctat_genome_lib_build_dir/ref_genome.fa
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/38/35a54b79f58af46a874aab63d5ee92/.command.sh
Downloading: s3://natera-platform-sandbox/pipeline-resources/AIH/rna/GRCh38/starfusion/ctat_genome_lib_build_dir/ref_genome.dict
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/7a/b743f77ee72f1e406649cc5cbaf4b2/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bam
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/38/35a54b79f58af46a874aab63d5ee92/.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/7a/b743f77ee72f1e406649cc5cbaf4b2/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
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-9114fc-R1_B23WHYVLT4_1
+ all_qc -bam aih-tih-sc-9114fc-R1_B23WHYVLT4_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-9114fc-R1_B23WHYVLT4_1 -reference ref_genome.fa -region chr20
mkdir -p aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_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-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md/full_bed.bed
seqtool depth -bam aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bam -bed aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md/full_bed.bed -threads 4 -bam_threads 2 | bgzip -c > aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.depth.gz
saturation_plot -bams aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bam -region chr20 -output aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md/saturation_plot.png -threads 1
seqtool saturation_metric -at 1000000 -region chr20 -bam_threads 1 aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bam > aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md/saturation_qc.tsv
fragment_size -bams aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bam -region chr20 -output_prefix aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_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-9114fc-R1_B23WHYVLT4_1.md.bam > aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.mpileup && contamination -input aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.mpileup -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz > aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md/contamination.tsv
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_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-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_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-9114fc-R1_B23WHYVLT4_1.md.bai
0it [00:00, ?it/s]
135454it [00:00, 1354324.95it/s][W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
280838it [00:00, 1412624.44it/s]
425679it [00:00, 1428736.38it/s]
573896it [00:00, 1449809.19it/s]
721256it [00:00, 1458268.72it/s]
868315it [00:00, 1462325.14it/s]
1015976it [00:00, 1466976.66it/s]
1162674it [00:00, 1464672.75it/s]
1309618it [00:00, 1466006.80it/s]
1457015it [00:01, 1468441.89it/s]
1603915it [00:01, 1468507.01it/s]
1751758it [00:01, 1471505.47it/s]
1898910it [00:01, 1464573.54it/s]
2047286it [00:01, 1470332.90it/s]
2199355it [00:01, 1485458.93it/s]
2351113it [00:01, 1495094.92it/s]
2500631it [00:01, 1484813.24it/s]
2653073it [00:01, 1496629.26it/s]
2803372it [00:01, 1498517.49it/s]
2953240it [00:02, 1495776.64it/s]
3102829it [00:02, 1491732.37it/s]
3252012it [00:02, 1482836.06it/s]
3401899it [00:02, 1487540.77it/s]
3550667it [00:02, 1480029.37it/s]
3698685it [00:02, 1476067.56it/s]
3847366it [00:02, 1479190.10it/s]
3995295it [00:02, 1475684.55it/s]
4142871it [00:02, 1467055.46it/s]
4289589it [00:02, 1459796.26it/s]
4438018it [00:03, 1467065.39it/s]
4586068it [00:03, 1471055.52it/s]
4733258it [00:03, 1471177.22it/s]
4881126it [00:03, 1473404.65it/s]
5032012it [00:03, 1483954.84it/s]
5180415it [00:03, 1479873.04it/s]
5329152it [00:03, 1482025.90it/s]
5477361it [00:03, 1475425.17it/s]
5625817it [00:03, 1478134.37it/s]
5773639it [00:03, 1473431.32it/s]
5920990it [00:04, 1460424.88it/s]
6067059it [00:04, 1459026.78it/s]
6215365it [00:04, 1466054.54it/s]
6364445it [00:04, 1473423.45it/s]
6511804it [00:04, 1473050.59it/s]
6659120it [00:04, 1456451.35it/s]
6804812it [00:04, 1448396.06it/s]
6949687it [00:04, 1440443.35it/s]
7094204it [00:04, 1441829.14it/s]
7241539it [00:04, 1451191.87it/s]
7392438it [00:05, 1468289.13it/s]
7544776it [00:05, 1484651.50it/s]
7696679it [00:05, 1494918.66it/s]
7846493it [00:05, 1495865.37it/s]
7998059it [00:05, 1501778.54it/s]
8148248it [00:05, 1497613.32it/s]
8302160it [00:05, 1510014.26it/s]
8453172it [00:05, 1503895.40it/s]
8603574it [00:05, 1491032.60it/s]
8754067it [00:05, 1495148.06it/s]
8904503it [00:06, 1497878.76it/s]
9054488it [00:06, 1498450.88it/s]
9204347it [00:06, 1493881.67it/s]
9357718it [00:06, 1505761.78it/s]
9509867it [00:06, 1510446.96it/s]
9660923it [00:06, 1507767.22it/s]
9812248it [00:06, 1509377.81it/s]
9964189it [00:06, 1512365.85it/s]
10115431it [00:06, 1506323.93it/s]
10266072it [00:06, 1496251.28it/s]
10415715it [00:07, 1486647.83it/s]
10564978it [00:07, 1488408.70it/s]
10713835it [00:07, 1469046.00it/s]
10860796it [00:07, 1467411.35it/s]
11008746it [00:07, 1470977.97it/s]
11156772it [00:07, 1473631.30it/s]
11306468it [00:07, 1480575.38it/s]
11455434it [00:07, 1483155.77it/s]
11605329it [00:07, 1487863.76it/s]
11754478it [00:07, 1488934.07it/s]
11903380it [00:08, 1485456.57it/s]
12055552it [00:08, 1496288.81it/s]
12205189it [00:08, 1496119.14it/s]
12358432it [00:08, 1506980.80it/s]
12511654it [00:08, 1514526.24it/s]
12663112it [00:08, 1512757.61it/s]
12815884it [00:08, 1517222.49it/s]
12967610it [00:08, 1516591.50it/s]
13119883it [00:08, 1518373.16it/s]
13271722it [00:08, 1516205.76it/s]
13424485it [00:09, 1519612.11it/s]
13576748it [00:09, 1520499.58it/s]
13728800it [00:09, 1520274.70it/s]
13880829it [00:09, 1514969.11it/s]
14032331it [00:09, 1514021.59it/s]
14183737it [00:09, 1510722.24it/s]
14334813it [00:09, 1504003.49it/s]
14485221it [00:09, 1493053.69it/s]
14634544it [00:09, 1476228.73it/s]
14782209it [00:09, 1456253.11it/s]
14927899it [00:10, 1439427.33it/s]
15071900it [00:10, 1434884.46it/s]
15216466it [00:10, 1437964.43it/s]
15360718it [00:10, 1439295.14it/s]
15504670it [00:10, 1438024.79it/s]
15648488it [00:10, 1428304.45it/s]
15793698it [00:10, 1435357.21it/s]
15938947it [00:10, 1440443.60it/s]
16083971it [00:10, 1443347.52it/s]
16229939it [00:10, 1448136.49it/s]
16377674it [00:11, 1456862.56it/s]
16527074it [00:11, 1467773.95it/s]
16679189it [00:11, 1483748.89it/s]
16830885it [00:11, 1493600.20it/s]
16981216it [00:11, 1496492.35it/s]
17130870it [00:11, 1495629.22it/s]
17280450it [00:11, 1495656.25it/s]
17430018it [00:11, 1483353.50it/s]
17578378it [00:11, 1480619.32it/s]
17726457it [00:11, 1474338.24it/s]
17873906it [00:12, 1463574.48it/s]
18023158it [00:12, 1472151.93it/s]
18170913it [00:12, 1473743.64it/s]
18318787it [00:12, 1475217.72it/s]
18467798it [00:12, 1479652.39it/s]
18615774it [00:12, 1479006.47it/s]
18763682it [00:12, 1478663.71it/s]
18911554it [00:12, 1471882.51it/s]
19058752it [00:12, 1468822.09it/s]
19205642it [00:12, 1465468.09it/s]
19353890it [00:13, 1470529.96it/s]
19500973it [00:13, 1470601.22it/s]
19651070it [00:13, 1479569.40it/s]
19800008it [00:13, 1482485.74it/s]
19948261it [00:13, 1471574.50it/s]
20095439it [00:13, 1467603.45it/s]
20242214it [00:13, 1466759.20it/s]
20392891it [00:13, 1478687.24it/s]
20540772it [00:13, 1441408.72it/s]
20689274it [00:13, 1454210.17it/s]
20834865it [00:14, 1451485.25it/s]
20985469it [00:14, 1467636.85it/s]
21132333it [00:14, 1463150.33it/s]
21280978it [00:14, 1470066.27it/s]
21428040it [00:14, 1461392.60it/s]
21574622it [00:14, 1462599.41it/s]
21723993it [00:14, 1471701.37it/s]
21872907it [00:14, 1476895.00it/s]
22020619it [00:14, 1473079.79it/s][W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
22167944it [00:15, 1469235.87it/s]
22299011it [00:15, 1477871.24it/s]
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
tabix -s 1 -b 2 -e 2 aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.depth.gz
seqtool depth_hist -nbins 500 aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.depth.gz > aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.depth_hist.tsv
cat aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.depth_hist.tsv | plot -x 0 -xlabel Depth -ylabel "Number of Sites" -xlim 0 2000 -title aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md -output aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md/aih-tih-sc-9114fc-R1_B23WHYVLT4_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-9114fc-R1_B23WHYVLT4_1.md.bam > aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_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-9114fc-R1_B23WHYVLT4_1.md.bam > aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_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-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.depth.gz> aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_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-9114fc-R1_B23WHYVLT4_1.md.bam > aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_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-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md.bai
cat aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_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-9114fc-R1_B23WHYVLT4_1.md" -suptitle "Twist_rna_exome_AND_fusion_hg38.bed" -output "aih-tih-sc-9114fc-R1_B23WHYVLT4_1/aih-tih-sc-9114fc-R1_B23WHYVLT4_1.md/gc_bias.Twist_rna_exome_AND_fusion_hg38.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"