homoplasyFinder
homoplasyFinder copied to clipboard
Output inconsistent sites as VCF file
An important enhancement pointed out by flashton2003. Thanks!
Include option in homoplasyFinder
to export the consistencyIndexReport_DATE.txt
file in a VCF file format.
The VCF file can then easily be run with bcftools csq
to look at the potential consequences of the changes at the inconsistent sites.
Ideally this would require a reference genome file to be provided so the REF column is accurately filled out. It might be possible to avoid this if I assume the allele with the largest count if the reference/ancestral allele.