Atkinson-Lab / Tractor

Scripts for implementing the Tractor pipeline
MIT License
44 stars 5 forks source link

fix msp scan #10

Closed mike-w-wilson closed 2 years ago

mike-w-wilson commented 3 years ago

RFMix does not necessary start the ancestry windows where the input VCF does. Currently, the ExtractTracts.py script will read through the entire msp file if the first sites in the input VCF is before the msp windows start. It will then use the last windows ancestry calls as the ancestry calls for all sites. This fixes the bug by skipping the VCF site if it is before the first windows start.