微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

如何使用 varscan 三重奏调用者管理 samtools mileup?

如何解决如何使用 varscan 三重奏调用者管理 samtools mileup?

我已经看到如下使用 varscan 管道 samtools

samtools mpileup -f reference.fasta sample1.bam sample2.bam | java -jar VarScan.jar mpileup2snp

但是如何为 varscan trio 管道以下管道

#Generate a three-sample mpileup
samtools mpileup -B -q 1 \
     -f ref.fasta \
     dad.bam mom.bam child.bam > trio.mpileup

#Run VarScan trio
varscan trio \
     trio.mpileup \ #this is input
     trio.mpileup.output \ #this is output
     --min-coverage 10 --min-var-freq 0.20 --p-value 0.05 \ 
     -adj-var-freq 0.05 -adj-p-value 0.15

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。