大橙子网站建设,新征程启航

为企业提供网站建设、域名注册、服务器等服务

freebayes怎么安装使用

本篇内容介绍了“freebayes怎么安装使用”的有关知识,在实际案例的操作过程中,不少人都会遇到这样的困境,接下来就让小编带领大家学习一下如何处理这些情况吧!希望大家仔细阅读,能够学有所成!

大悟网站制作公司哪家好,找创新互联!从网页设计、网站建设、微信开发、APP开发、自适应网站建设等网站项目制作,到程序开发,运营维护。创新互联公司2013年成立到现在10年的时间,我们拥有了丰富的建站经验和运维经验,来保证我们的工作的顺利进行。专注于网站建设就选创新互联。

freebayes 是一款snp calling 软件,其灵敏度高,用法简便,所以广受欢迎。

软件的安装过程如下

git clone --recursive git://github.com/ekg/freebayes.git
cd freebayes/
make

编译成功之后,在bin目录就是可执行文件。为了使用方便,可以将bin目录添加到PATH环境变量中。

对于freebayes而言,只需要两个输入文件,一个是参考基因组的fasta文件,另外一个是比对产生的bam文件。基本用法如下

freebayes -f ref.fasta align.bam >var.vcf

参考基因组的fasta文件需要有后缀为.fai的索引文件,可以通过samtools来构建,命令如下

samtools faidx ref.fasta

如果你提供的fasta文件没有对应的索引,程序会自动去构建。对于大型参考基因组而言,建议是先构建好索引。比对的bam文件可以按照GATK官方推荐的预处理流程得到。

输出结果是VCF格式的,示例如下

头部

##fileformat=VCFv4.2
##fileDate=20180626
##source=freeBayes v1.2.0
##reference=ref.fasta
##contig=
##phasing=none
##commandline="./freebayes --fasta-reference ref.fasta align.bam"
##INFO=

正文

#CHROM    POS    ID    REF    ALT    QUAL    FILTER    INFO    FORMAT    GW1
NC_023084.1    107    .    G    T    331.137    .    AB=0;ABP=0;AC=2;AF=1;AN=2;AO=13;CIGAR=1X;DP=13;DPB=13;DPRA=0;EPP=3.17734;EPPR=0;GTI=0;LEN=1;MEANALT=1;MQM=60;MQMR=0;NS=1;NUMALT=1;ODDS=22.627;PAIRED=0.615385;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=423;QR=0;RO=0;RPL=0;RPP=31.2394;RPPR=0;RPR=13;RUN=1;SAF=6;SAP=3.17734;SAR=7;SRF=0;SRP=0;SRR=0;TYPE=snp;technology.illumina=1    GT:DP:AD:RO:QR:AO:QA:GL    1/1:13:0,13:0:0:13:423:-38.3849,-3.91339,0

VCF格式之前的文章中已经详细介绍过,每个字段的含义可以参考头部的注释信息。

“freebayes怎么安装使用”的内容就介绍到这里了,感谢大家的阅读。如果想了解更多行业相关的知识可以关注创新互联网站,小编将为大家输出更多高质量的实用文章!


当前名称:freebayes怎么安装使用
标题网址:http://dzwzjz.com/article/gesgei.html
在线咨询
服务热线
服务热线:028-86922220
TOP