将 vcf 转换为 hap 文件，折叠基因型

Question

我用照片问问题

我有一个vcf file.
还有一个R-code.
这是我的输出结果：output

我希望我的输出文件只有 number/line/sample。

如何删除新数据...C...和“”

有什么建议吗？

注意：我是第一次在这里提问。我试着遵守规则。我还在学习。

Answer 1

试试这个例子：

library(data.table)

# example vcf
hap <- fread('
##fileformat=VCFv4.0
##fileDate=20090805
##source=myImputationProgramV3.1
##reference=1000GenomesPilot-NCBI36
##phasing=partial
##INFO=<ID=NS,Number=1,Type=Integer,Description="Number of Samples With Data">
##INFO=<ID=AN,Number=1,Type=Integer,Description="Total number of alleles in called genotypes">
##FILTER=<ID=q10,Description="Quality below 10">
##FILTER=<ID=s50,Description="Less than 50% of samples have data">
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
#CHROM  POS ID  REF ALT QUAL    FILTER  INFO    FORMAT  NA00001 NA00002 NA00003
19  111 .   A   C   9.6 .   .   GT  0|0 0|0 0|1
19  112 .   A   G   10  .   .   GT  0|0 1|0 1|1
19  112 .   A   G   4   .   .   GT  1|0 1|0 1|1
')

data.table(gsub("|", "", do.call(paste0, hap[, -c(1:9)]), fixed = TRUE))
#        V1
# 1: 000001
# 2: 001011
# 3: 101011

将 vcf 转换为 hap 文件，折叠基因型

Convert vcf to hap file, collapse genotypes

r

bioinformatics

vcf-variant-call-format