Haplotype analysis from phased vcf
0
0
Entering edit mode
6.2 years ago
BAGeno ▴ 190

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.4k views
ADD COMMENT
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 2111 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6