Our staff will be observing the holiday and will be unavailable from 22nd to 25th November. This will cause a delay in reaching out to you and answering your questions immediately. Rest assured we will get back to it on Monday November 26th. We are grateful for your support and patience.
Have a great holiday everyone!!!
Using Picardtools LiftoverVCF between species
I have a VCF file produced by mapping reads from chimpanzee onto the human genome. I want to take this VCF and do a liftover to the chimpanzee genome using picardtools LiftoverVcf. However, when I do this, everything gets dumped in the 'rejected' folder and labeled "MismatchedRefAllele".
I understand why this would be a good feature when you are mapping to the same species, but is there anyway to change this behavior in the program? I tried setting the "LIFTOVER_MIN_MATCH" parameter to 0, but that doesn't seem to do anything.
One possibility is to switch the "REF" and "ALT" bases in the original VCF and ignore all positions that have an allele frequency of 100% (AF=1.000)--i.e. fixed differences between human and chimpanzee. However, I'd prefer not to mess with the file too much myself if possible.
I also tried liftovervariants from and earlier version of GATK (2.9), but it doesn't seem like it supports lifting over between species.
Any thoughts would be much appreciated!