Question

phasing genotype

0

Entering edit mode

2.1 years ago

FRANCESCA ▴ 10

How can I "translate" the genotype of my vcf file in alleles? For example: 0|1 vs A|T

genotype • 827 views

ADD COMMENT • link updated 2.1 years ago by Pierre Lindenbaum 164k • written 2.1 years ago by FRANCESCA ▴ 10

score 1 · Answer 1 · 2022-10-20

1

Entering edit mode

2.1 years ago

4galaxy77 2.9k

bcftools query -f'%REF|%ALT\n' in.vcf > output.txt

G|A
T|A
C|T
G|A
C|T
C|T
A|G
T|C
T|C
T|C
G|A

ADD COMMENT • link 2.1 years ago by 4galaxy77 2.9k

score 1 · Answer 2 · 2022-10-20

1

Entering edit mode

2.1 years ago

iraun 6.2k

If you want to reconstruct the haplotypes having a VCF with phased genotypes in place, I can recommend haplosaurus, developed by Ensembl.

ADD COMMENT • link 2.1 years ago by iraun 6.2k

score 1 · Answer 3 · 2022-10-20

1

Entering edit mode

2.1 years ago

Pierre Lindenbaum 164k

bcftools query -f '%CHROM %POS %ID[ %TGT]\n' in.vcf

ADD COMMENT • link 2.1 years ago by Pierre Lindenbaum 164k