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/1e/5e06d055f9176537f30e7bdaf93187/.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/fb/d7d2801b210e9db9056abad315e950/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_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/1e/5e06d055f9176537f30e7bdaf93187/.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/fb/d7d2801b210e9db9056abad315e950/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_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 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1
+ all_qc -bam 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_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 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1 -reference ref_genome.fa -region chr20
mkdir -p 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/
cat Twist_rna_exome_AND_fusion_hg38.bed | cut -f 1-3 | bedtools sort -g ref_genome.fa.fai | bedtools merge > 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/full_bed.bed
saturation_plot -bams 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bam -region chr20 -output 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/saturation_plot.png -threads 1
seqtool depth -bam 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bam -bed 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/full_bed.bed -threads 4 -bam_threads 2 | bgzip -c > 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.depth.gz
seqtool saturation_metric -at 1000000 -region chr20 -bam_threads 1 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bam > 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/saturation_qc.tsv
fragment_size -bams 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bam -region chr20 -output_prefix 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/fragsize -threads 1
samtools mpileup -B --min-MQ 60 --min-BQ 30 -f ref_genome.fa -l dbsnp_hg38_contam_sites.bed 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bam > 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.mpileup && contamination -input 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.mpileup -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz > 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/contamination.tsv
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[mpileup] 1 samples in 1 input files
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_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: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
0it [00:00, ?it/s]
142084it [00:00, 1420492.52it/s]
296064it [00:00, 1490616.99it/s]
458913it [00:00, 1553529.91it/s]
623386it [00:00, 1589397.03it/s]
786574it [00:00, 1604682.73it/s]
952015it [00:00, 1621563.02it/s]
1116515it [00:00, 1629197.23it/s]
1282906it [00:00, 1640218.16it/s]
1448634it [00:00, 1645489.64it/s]
1614053it [00:01, 1648146.83it/s]
1778868it [00:01, 1638517.04it/s]
1943678it [00:01, 1641395.87it/s]
2110296it [00:01, 1648861.34it/s]
2279408it [00:01, 1661578.71it/s]
2447426it [00:01, 1667022.70it/s]
2615100it [00:01, 1669919.86it/s]
2782098it [00:01, 1641146.84it/s]
2950276it [00:01, 1653174.69it/s]
3115688it [00:01, 1648919.51it/s]
3280645it [00:02, 1618849.68it/s]
3442677it [00:02, 1549638.33it/s]
3598271it [00:02, 1539208.11it/s]
3759334it [00:02, 1559828.06it/s]
3915685it [00:02, 1514222.57it/s]
4067561it [00:02, 1474901.89it/s]
4231885it [00:02, 1523045.68it/s]
4386540it [00:02, 1529829.45it/s]
4541705it [00:02, 1536118.13it/s]
4700203it [00:02, 1550460.69it/s]
4862928it [00:03, 1573201.70it/s]
5020432it [00:03, 1531397.70it/s]
5173918it [00:03, 1472597.52it/s]
5324976it [00:03, 1483452.53it/s]
5473802it [00:03, 1453831.59it/s]
5619573it [00:03, 1424258.23it/s]
5762309it [00:03, 1411809.56it/s]
5903683it [00:03, 1391513.65it/s]
6042969it [00:03, 1380047.72it/s]
6200997it [00:04, 1438030.45it/s]
6366527it [00:04, 1501743.49it/s]
6531157it [00:04, 1544421.36it/s]
6685881it [00:04, 1522710.60it/s]
6852559it [00:04, 1565140.61it/s]
7022376it [00:04, 1604475.70it/s]
7191603it [00:04, 1630559.03it/s]
7357325it [00:04, 1638489.32it/s]
7526522it [00:04, 1654441.38it/s]
7692065it [00:04, 1647741.99it/s]
7859389it [00:05, 1655313.86it/s]
8027960it [00:05, 1664361.19it/s]
8199458it [00:05, 1679466.79it/s]
8367438it [00:05, 1667048.28it/s]
8535547it [00:05, 1671209.02it/s]
8705691it [00:05, 1680209.94it/s]
8873739it [00:05, 1677227.54it/s]
9044989it [00:05, 1687744.83it/s]
9213781it [00:05, 1680913.73it/s]
9381890it [00:05, 1645751.73it/s]
9553154it [00:06, 1665367.27it/s]
9719840it [00:06, 1564065.78it/s]
9877521it [00:06, 1438912.07it/s]
10042897it [00:06, 1497158.75it/s]
10209907it [00:06, 1545473.62it/s]
10380879it [00:06, 1592318.10it/s]
10542915it [00:06, 1600411.58it/s]
10704128it [00:06, 1601169.07it/s]
10865067it [00:06, 1489160.90it/s]
11022092it [00:07, 1511916.39it/s]
11192825it [00:07, 1567851.83it/s]
11363441it [00:07, 1607980.89it/s]
11533605it [00:07, 1635415.36it/s]
11705540it [00:07, 1660154.97it/s]
11874308it [00:07, 1668293.44it/s]
12045598it [00:07, 1681540.62it/s]
12215287it [00:07, 1686022.98it/s]
12385174it [00:07, 1689799.67it/s]
12557710it [00:07, 1700404.35it/s]
12727866it [00:08, 1690081.03it/s]
12896968it [00:08, 1677700.97it/s]
13064818it [00:08, 1665520.89it/s]
13233542it [00:08, 1671938.41it/s]
13400790it [00:08, 1669409.45it/s]
13570420it [00:08, 1677397.63it/s]
13738191it [00:08, 1674751.99it/s]
13905688it [00:08, 1648688.15it/s]
14077513it [00:08, 1669235.15it/s]
14249566it [00:08, 1684408.99it/s]
14421084it [00:09, 1693546.00it/s]
14590503it [00:09, 1686053.83it/s]
14760883it [00:09, 1684048.76it/s]
14935834it [00:09, 1703509.29it/s]
15106875it [00:09, 1705543.35it/s]
15277587it [00:09, 1705959.13it/s]
15448204it [00:09, 1685375.32it/s]
15616808it [00:09, 1675934.50it/s]
15784451it [00:09, 1653871.15it/s]
15952251it [00:09, 1660957.83it/s]
16120547it [00:10, 1667451.94it/s]
16290314it [00:10, 1676346.72it/s]
16464189it [00:10, 1694931.47it/s]
16635991it [00:10, 1701806.20it/s]
16806205it [00:10, 1698696.34it/s]
16976099it [00:10, 1681548.73it/s]
17144305it [00:10, 1668028.36it/s]
17311157it [00:10, 1664226.67it/s]
17478584it [00:10, 1667190.15it/s]
17645328it [00:10, 1664697.64it/s]
17812990it [00:11, 1668227.86it/s]
17980815it [00:11, 1671196.95it/s]
18147946it [00:11, 1663321.31it/s]
18315398it [00:11, 1666490.01it/s]
18485634it [00:11, 1677182.11it/s]
18653696it [00:11, 1678184.79it/s]
18821524it [00:11, 1554030.42it/s]
18990519it [00:11, 1592551.27it/s]
19155214it [00:11, 1608077.01it/s]
19324682it [00:11, 1633320.48it/s]
19492468it [00:12, 1646396.17it/s]
19664145it [00:12, 1667191.11it/s]
19834508it [00:12, 1677917.17it/s]
20003532it [00:12, 1681566.70it/s]
20171925it [00:12, 1677996.58it/s]
20345281it [00:12, 1694557.81it/s]
20514862it [00:12, 1691102.42it/s]
20688753it [00:12, 1705357.18it/s]
20861295it [00:12, 1711284.69it/s]
21032473it [00:12, 1705324.14it/s]
21203044it [00:13, 1700318.24it/s]
21377029it [00:13, 1712098.68it/s]
21548265it [00:13, 1648069.20it/s]
21718402it [00:13, 1663575.26it/s]
21869929it [00:13, 1623415.05it/s]
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
tabix -s 1 -b 2 -e 2 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.depth.gz
seqtool depth_hist -nbins 500 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.depth.gz > 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.depth_hist.tsv
cat 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.depth_hist.tsv | plot -x 0 -xlabel Depth -ylabel "Number of Sites" -xlim 0 2000 -title 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md -output 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_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 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bam > 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_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 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bam > 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/reads_per_target.Twist_rna_exome_AND_fusion_hg38.bed.tsv
seqtool depth_stats -bed Twist_rna_exome_AND_fusion_hg38.bed 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.depth.gz> 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_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 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bam > 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_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: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
[W::hts_idx_load3] The index file is older than the data file: 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md.bai
cat 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/gc_bias_qc.Twist_rna_exome_AND_fusion_hg38.bed.tsv | plot -x 0 -header infer -ylim 0 2 -title "1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md" -suptitle "Twist_rna_exome_AND_fusion_hg38.bed" -output "1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/gc_bias.Twist_rna_exome_AND_fusion_hg38.bed.png" -ylabel "Fold enrichment / depletion" -xlabel "GC Quantile"
echo `date` > 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md/qc_complete.txt
+ merge_qc -folders 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1.md -output 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1_full_basic_stats.tsv
+ head -n1 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1_full_basic_stats.tsv
++ basename Twist_rna_exome_AND_fusion_hg38.bed
+ grep -e Twist_rna_exome_AND_fusion_hg38.bed 1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1/1173_CRL_T1_RNA_SLD_01_A23T55JLT4_1_full_basic_stats.tsv