Adding coverage data to your VCF_bam-readcount-优快云博客

本文链接：https://blog.youkuaiyun.com/bio_meimei/article/details/89948165

pVACseq可以从VCF直接解析覆盖率信息。如果VCF缺少覆盖率信息或需要为额外样本或RNA-seq数据添加信息，可以使用vcf-readcount-annotator。该工具结合bam-readcount的输出，在VCF中添加读取计数。在使用vcf-readcount-annotator之前，需要安装vt和bam-readcount。安装过程中可能遇到libcurl和libcrypto的问题，需要手动解决。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

文章目录

Using the vcf-readcount-annotator to add coverage information to your VCF.
Installing vt
Installing bam-readcount
Installing the vcf-readcount-annotator

Adding coverage data to your VCF

pVACseq is able to parse coverage information directly from the VCF. The expected annotation format is outlined below.
pvacseq能够解析你的覆盖信息从你的ＶＣＦ文件中.
有些

Using the vcf-readcount-annotator to add coverage information to your VCF.

Some variant callers will already have added coverage information to your VCF. However, if your VCF doesn’t contain coverage information or if you need to add coverage information for additional samples or for RNA-seq data, you can use the vcf-readcount-annotator to do so.

The vcf-readcount-annotator will take the output from bam-readcount and use it to add readcounts to your VCF.

bam-readcount needs to be run separately for snvs and indels so it is recommended to first split multi-allelic sites by using a tool such as vt decompose.

Installing vt

 #this will create a directory named vt in the directory you cloned the repository
 1. git clone https://github.com/atks/vt.git  
 #change directory to vt
 2. cd vt