====== 两种情况质控 ======
- 有下机路径,文件命名包含诺禾编号,如FKDO240072420-1A_16SV34,接下来以情况1替代
- 客户直接提供数据,文件命名自定义,如SC120s.1,接下来以情况2替代
根据两种情况的质控方法分别如下:
====== 情况1 ======
=== (1)首先需要配置list,包含七列分别为:文库号、标签、下机路径、诺禾编号、扩增区间、项目批次、测序通道,七列缺一不可 ===
{{ :kappframework-huptgi_1_1_.png?600 |}}
其中第一列、第四列、第五列、第七列都能在路径或文件信息中找到,剩余三列可自定义,第二列要求不同的标签、都三列和第六列填写相同内容
填写完list后使用:
=== (2)生成toml文件 ===
''/TJPROJ7/GB_MICRO/PUBLIC/software/amplicon/python3/bin/python /TJPROJ13/GB_MICRO/USER/yangyuhong/script/qc_toml/tab2list.py list''
----
=== (3)在目录下mkdir log ===
''/TJPROJ7/GB_MICRO/PUBLIC/software/amplicon/python3/bin/python /TJPROJ7/GB_MICRO/PUBLIC/source/amplicon/module/splitQualityControl -PJ qc.toml -OD ./ -SD ./log''
----
=== (4)将rawdata软链接到autoDemultiplex文件夹的最底层目录下 ===
----
=== (5)进入log目录 ===
''for i in $(ls qc*); do qsub -V -cwd -l vf=10G,p=2 -q meta_ass.q $i; done''
一键投递所有qc任务,结果在结果文件夹(项目批次id+扩增区间)下查看
====== 情况2 ======
=== (1)首先需要配置list,包含两列分别为:诺禾编号、扩增区间,两列缺一不可 ===
{{ :ka2.png?200 |}}
----
=== (2)生成toml文件 ===
''/TJPROJ7/GB_MICRO/PUBLIC/software/amplicon/python3/bin/python /TJPROJ7/META_ASS/16s/chenlei/jiaoben/qc/bin/tab2list.py list''
----
=== (3)在目录下mkdir log ===
''/TJPROJ7/GB_MICRO/PUBLIC/software/amplicon/python3/bin/python /TJPROJ7/META_ASS/16s/chenlei/jiaoben/qc/bin/splitQualityControl -PJ qc.toml -OD ./ -SD ./log''
----
=== (4)将rawdata软链接到autoDemultiplex文件夹扩增区间的最底层目录下 ===
----
=== (5)进入log下执行toudi.sh,结果在结果文件夹(扩增区间)下查看 ===