This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How does plink guess the reference base in "--recode vcf-iid"

Here: http://oryzasnp-atcg-irri-org.s3-website-ap-southeast-1.amazonaws.com/3krg-3k_filt_snp-v1/README-3kRG-filtered-SNP-v1.txt

The recommended SOP for recoding PED / MAP to VCF is:

plink --file <ped file name without extension> --recode vcf-iid --out <output file name without extension>

However, without specifying the reference sequence, how is the reference base in the VCF guessed?

Does plink have an extra option for passing the reference sequence when calling --recode?

I'm reading about the associated file formats here:

and I don't see the reference base stored anywhere.

vcf plink
• 4,019 views • link
updated by Ram Ram • written by Dan Dan

1 answer

From plink --help | less

  --recode ...

    The A2 allele is saved as the reference and normally flagged as not based
    on a real reference genome ('PR' INFO field value).  When it is important
    for reference alleles to be correct, you'll also want to include
    --a2-allele and --real-ref-alleles in your command.

  --a2-allele [filename] {a2col} {IDcol} {skip} :
    Force alleles in the file to A2.  ("--a2-allele [VCF filename] 4 3 '#'",
    which scrapes reference allele assignments from a VCF file, is especially
    useful.

Sigh

• 0 views • link
updated by Ram Ram • written by Dan Dan

Log in to answer this question.

More posts like this