售后

编辑本页后请点击“保存”。请参阅 syntax 了解维基语法。只有在您能改进该页面的前提下才编辑它。如果您想尝试一些东西，请先到 playground 热身。

======简介======
Compare-indel分析时，结果文件会出现空行，抓不到信息的情况。

此流程仅作为compare-indel分析后的修复结果。

======功能======
Compare-indel结果整理。
======数据准备======

**1. 链接基本分析结果：**

在售后操作目录中建立文件夹：

02.varDetect/02.SNPindel/01.detection

02.varDetect/02.SNPindel/02.ANNO

将基本分析中涉及到的indel变异结果vcf文件，注释文件链接到上述目录中；

**2. 生成位点文件：**

将all.vcf中位点的缺失信息提取，将indel.vcf中的变异信息提取；

**3. 合并indel注释文件：**

把sample1和sample2的4个注释结果文件合并；
<code>
cat sample1.avinput.variant_function.indel sample1.avinput.variant_function.indel sample2.avinput.exonic_variant_function.indel sample2.avinput.exonic_variant_function.indel >sample1_sample2.indel
</code>

======数据分析======
**1. 正常配置及运行流程中的compare-indel脚本：**
<code>
qsub run.compare.indel_pipeline.sh
</code>

**2. 对结果进行处理：**
<code>
mkdir deal
qsub cmd.sh
</code>

示例脚本：
<code>
/TJNAS01/AFS_RESEQ/Proj/hanyue/02.PAG/Compare/cmd.sh
</code>

注：根据样本路径及名称修改以下脚本：

tiquweidian_genotype_sample01.pl

tiquweidian_genotype_sample02.pl

tiquweidian.pl

cmd.sh的cut的列情况需要注意

编辑摘要

当您选择开始编辑本页，即寓示你同意将你贡献的内容按下列许可协议发布： CC Attribution-Share Alike 4.0 International

售后

用户工具

站点工具

页面工具