用户工具

站点工具


框架图新流程生成clean

单菌框架图自动化流程重新生成clean

框架图自动化流程分析只交付了rawdata,但是有客户后续找回需要重新交付clean,可以使用框架图新流程手动步骤将首次交付的rawdata链接到当前路径,然后配置rawlist进行重新质控,但是如果遇到以数字开头的样本名,不可作为sjm的任务名进行分析,所以还涉及到样本重新改名的问题,为了简化操作,可以直接跑自动化流程进行质控。

现整理了相关流程步骤如下:

准备工作

1. 切换liuqingqing账号

运行步骤

1. 配置文件:

将自动化分析路径下的mapfile_all,run.sh,project_info,sample_list拷贝到当前路径:

cp mapfile_all run.sh project_info sample_list 当前分析目录/

将mapfile_all中cleandata列中的0替换为1:

sed -i 's/\t0\t/\t1\t/g'  mapfile_all

复制mapfile_all作为mapfile:

cp mapfile_all mapfile

注:首次分析路径下的mapfile_all记录了所有加测路径,可以直接作为重新质控的mapfile

2. 将run.sh中的–auto参数从y改为n;添加参数–ifpushpom n;修改–mapfile和–outdir路径;修改–batchid 分期号:

修改后的run.sh格式如下:

/PUBLIC/software/RESEQ/software/anaconda2/setup/bin/python /TJPROJ2/GB/PUBLIC/source/GB_SEQ/pipline_v5.1/pipline/..//pipline/GB_QCseq_pipline_v1.py \
        --contract H101SC24042805 \
        --batchid X101SC24042805-Z01-F001 \
        --mapfile /TJPROJ7/META_ASS/16s/chenjiawei/X101SC24042805-Z01/X101SC24042805-Z01-F001/kjt-qc-20240506/mapfile \
        --outdir /TJPROJ7/META_ASS/16s/chenjiawei/X101SC24042805-Z01/X101SC24042805-Z01-F001/kjt-qc-20240506 \
        --header n \
        --auto n \
        --type DSW \
        --admin_email 'yanyoudong@novogene.com;zhangxiaoxia@novogene.com;liutiance7240@novogene.com' \
        --zcat_flag y \
        --gn_flag y \
        --fenyangflag N \
        --checkmismatchflag n \
        --showqcresultbylane n \
        --ifpushpom n

3. 运行run.sh:

sh run.sh

注:此步执行后会生成sjm投递的sjm_QC.sh

4. sjm投递任务:

在sjm_QC.sh第一行添加下面环境:

source /TJPROJ5/GB_SEQ/PJ_GB/seq/auto/task/bash_profile_zidong

投递任务

nohup sh sjm_QC.sh &

流程测试路径:

/TJPROJ7/META_ASS/16s/chenjiawei/X101SC24042805-Z01/X101SC24042805-Z01-F001/kjt-qc-20240506
框架图新流程生成clean.txt · 最后更改: 2024/05/06 08:48 由 chenjiawei