homoplasyFinder icon indicating copy to clipboard operation
homoplasyFinder copied to clipboard

Output inconsistent sites as VCF file

Open JosephCrispell opened this issue 5 years ago • 3 comments

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.

JosephCrispell avatar Nov 05 '19 16:11 JosephCrispell