VCFtools: how does --diff option consider genotypes
Hello,
I'm trying to compare two vcf files, and I understand that VCFtools output will definitely show vcf1 and vcf2 differ, when they are AA vs TT. But how about the following two cases:
- AA vs AT -- When I see ALT1 column is A, ALT2 column is T, how would I know it is AA vs TT, or AA vs AT?
- What if AT vs AC?
Thank you very much.
• 2,695 views
•
link
1 answer
Well upon a few testing, I found that VCFtools does not distinguish the above two situations. I think it is sometimes better to just write a simple code to process two vcf comparisons.
• 1 views
•
link
Log in to answer this question.