ERROR: Variant physical position must be monotonically increasing
0
0
Entering edit mode
2.3 years ago
reza ▴ 300

I want to calculate XPEHH for each SNP position. When I run the following command

selscan --xpehh --vcf B10_beagle.vcf --vcf-ref D6_beagle.vcf --map MAP.map --threads 8 --out B10vsD6

I get this error

ERROR: Variant physical position must be monotonically increasing
Ch2:66 66 appears after Ch1:124715015 1247150

my solution for this problem is sorting map file based on SNP position by

sort -nk 4 MAP.map > MAP_sorted.map

running selscan with MAP_sorted.map solved this error, now my question is, my way to solve this error is correct?

XPEHH selscan • 791 views
ADD COMMENT
0
Entering edit mode

Should be ok, can verify by writing a quick script to check chromosome & position from vcf and check against the map, then added confirmation by actually confirming using bash or python

ADD REPLY

Login before adding your answer.

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