This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Haplotype analysis from phased vcf

Hi,

I want to ask that how can I perform haplotype analysis from phased vcf. I have vcf in this format.

#CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO    FORMAT  S1     S2      S3      S4     S5     S6
10      94870   rs123        A       G       30      PASS    AN=20;AC=11     GT      ./.     ./.     ./.     ./.     ./.     ./.     1|0     ./.     ./.
10      94872   rs456        A       C       30      PASS    AN=20;AC=11     GT      ./.     ./.     ./.     ./.     ./.     ./.     1|0     ./.     ./.
10      95074   rs345        G       A       30      PASS    AN=22;AC=13     GT      0|1     1|0     ./.     ./.     ./.     ./.     1|0     ./.     ./.
10      95096   rs234      A       G       30      PASS    AN=26;AC=19     GT      1|0     1|1     1|1     0|1     1|1     0|1     ./.     1|1     ./.

I want pharmgkb haplotype table this like.

LOCATION    ALLELE  REFERENCE ALLELE
rs2230037   A   G
rs5030868   A   G

Is there any way to do this.

vcf phased haplotype
• 1,656 views • link
written by BAGeno BAGeno

from your input and output, you just need cut -f3-5 and I don't see how it's related to 'phasing'

• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

0 answers

No answers yet.

Log in to answer this question.

More posts like this