Selscan: ERROR: Alleles must be coded 0/1 only.
1
0
Entering edit mode
2.5 years ago

I am trying to run nsl statistics using selscan tool on a vcf file (created using freebayes). But i am getting this error

ERROR: Alleles must be coded 0/1 only.

0 2

Kindly give me a solution because genotypes can be 2 3 so how to compute the statistics, in this case.

Many thanks in advance.

VCF nsl Selscan • 1.6k views
ADD COMMENT
0
Entering edit mode

Please edit your post and add your freebayes command.

ADD REPLY
0
Entering edit mode

This was the freebayes command used to generate the vcf file..

  freebayes -f ref.fa aln.bam >var.vcf
ADD REPLY
1
Entering edit mode
2.5 years ago

you could filter out the multiallelic in the VCF using bcftools view -m2 M2 or convert multiallelic to biallelic using bcftools norm

ADD COMMENT
1
Entering edit mode

Side note: Use vt decomp | vt norm instead of bcftools norm to retain information on the original variant entries after the operation. The added INFO/OLD_MULTIALLELIC field serves as a quasi-log of the operation.

ADD REPLY
0
Entering edit mode

Sure, I am grateful to you for helping me out thank you. I will try this out

ADD REPLY

Login before adding your answer.

Traffic: 3558 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