浏览代码

summary

master
LUYAO REN 4 年前
父节点
当前提交
b74677e74a
共有 1 个文件被更改,包括 2 次插入0 次删除
  1. +2
    -0
      tasks/region.wdl

+ 2
- 0
tasks/region.wdl 查看文件

@@ -54,6 +54,7 @@ task region {

for i in *vcf.gz
do
cat $i | grep -v '#' >> ${sample}.${tag}.all.count.txt
cat $i | grep -v '#' | grep 'SNV;MC' | wc -l >> ${sample}.${tag}.snv.mc.count.txt
cat $i | grep -v '#' | grep 'SNV;MC' | awk '{ if (($10 == $11) && ($11 == $12) && ($12 == $13)) { print } }' | wc -l >> ${sample}.${tag}.snv.mc.same.genotype.count.txt
cat $i | grep -v '#' | grep 'INDEL;MC' | wc -l >> ${sample}.${tag}.indel.mc.count.txt
@@ -98,5 +99,6 @@ task region {
File indel = "${sample}.${tag}.indel.mc.count.txt"
File indel_mc = "${sample}.${tag}.indel.mc.same.genotype.count.txt"
File file_list = "file.list.txt"
File all = "${sample}.${tag}.all.count.txt"
}
}

正在加载...
取消
保存