问题描述
我正在对 3 个 RADseq 数据板进行 ipyrad 分析,但条形码是多余的(每个板重复使用同一组 96 个条形码)。我希望一起分析所有 3 个板块,但不确定如何在 ipyrad 中做到这一点。或者是否需要添加额外的条形码以获得288个样本的288个唯一的左+右条形码。
Fastq 序列标识符: @MG01HX02:1030:HCTKYCCX2:8:1101:8410:1432 1:N:0:NTAACGAC+NCTATACC
解决方法
如果有人有类似的问题,我学会了如何做到这一点。 分别为每个板运行 ipyrad 步骤 1,以根据正确的条形码文件分配唯一 ID。确保所有这些输出都在同一个文件夹中。从第一步(到第七步)重新运行 ipyrad,但要引用 fastqs 的这个文件夹,而不是初始的 fastq 文件。