Question: Strands are reversed
0
gravatar for richyanicky
12 months ago by
richyanicky10
richyanicky10 wrote:

Hello,

I have converted a Ukbiobank bgen file to vcf for processing against a reference panel using beagle.

to make the pgen file from bgen I used the following.

plink2 --bgen ukb_imp_chr17_v3.bgen --sample ukb####_imp_chr17_v3_s#####.sample --make-pgen --keep-allele-order

There is a message that --keep-allele-order has no effect ..

This is the step for the vcf:

plink2 --pgen plink2.pgen --pvar plink2.pvar --psam plink2.psam  --export vcf-4.2 vcf-dosage=DS bgz --from-bp 50 --to-bp 10000 --chr 17

The issue I see between the vcf I create and the reference panel is that at each location the Ref and Alt are opposite. This causes Beagle to fail and state there are no matches. Which is correct.

How can I switch or maintain the order of the Ref & Alt so beagle will match them?

Thanks,

Richard

beagle plink • 322 views
ADD COMMENTlink modified 12 months ago by genomax85k • written 12 months ago by richyanicky10

Use the 'ref-first' --bgen modifier.

"plink2 --bgen <bgen filename> ref-first --sample ..."

ADD REPLYlink written 12 months ago by chrchang5237.0k

That cleared up the issue Thank You!

ADD REPLYlink written 12 months ago by richyanicky10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1667 users visited in the last hour