在与SLURM进行交互式会话时,如何在我的双消化,双末端Illumina NovaSeq数据中使用Stacks process_radtags?

问题描述

我想使用Stacks process_radtags开始分离合并的DNA样本,以便对Cassin's Sparrow进行亲子鉴定。我是通过哈佛大学使用SLURM的在线集群Cannon进行此操作的。我想在一个互动会议上做这些工作。这是SLURM运行交互式会话的示例:

srun --pty -p test -t 10 --mem 1000 /bin/bash

就process_radtags而言,我已经找到了如何运行该程序的特定示例(我的两种消化酶是SphI-HF和EcoRI-HF):

process_radtags -P -p ./raw  -b ./barcodes/barcodes_lane4 -o ./samples/ \
                  -c -q -r --inline_index --renz_1 nlaIII --renz_2 mluCI

来自该网站:https://catchenlab.life.illinois.edu/stacks/comp/process_radtags.php#hiseq

这是我当前尝试使用的原始数据的样子。我想从mydata的子集开始。也许现在只关注配对的索引文件之一,然后使用process_radtags将它们分开,然后使用Trimmomatic或Stacks进行质量过滤。

1_S1_L002_R1_001.fastq.gz  3_S3_L002_R1_001.fastq.gz  5_S5_L002_R1_001.fastq.gz  7_S7_L002_R1_001.fastq.gz
1_S1_L002_R2_001.fastq.gz  3_S3_L002_R2_001.fastq.gz  5_S5_L002_R2_001.fastq.gz  7_S7_L002_R2_001.fastq.gz
2_S2_L002_R1_001.fastq.gz  4_S4_L002_R1_001.fastq.gz  6_S6_L002_R1_001.fastq.gz  8_S8_L002_R1_001.fastq.gz
2_S2_L002_R2_001.fastq.gz  4_S4_L002_R2_001.fastq.gz  6_S6_L002_R2_001.fastq.gz  8_S8_L002_R2_001.fastq.gz

我是生物信息学和通用计算机科学的新手,对从所有这一切入手感到迷茫。任何帮助将不胜感激!

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)