1000G mask with VCFtools

Question

[vitti]

Hi all,

I'm trying to apply the 1000 Genomes masks to the phase 3 vcfs using VCFtools with the --mask option. Strangely, my output files are always truncated down to a few hundred KB, whereas they should each span an entire chromosome. I am not sure why this is happening.

I rewrote the mask FASTA files into the format VCFtools requires, simply converting P (= pass) to 0 and everything else to 1. (I also tried using the reverse convention, as the VCFtools documentation for this option is confusing, and tried both --mask and --invert-mask for each of these). I have checked to make sure that all the information in the original FASTA files (number of characters, position of line breaks, etc.) is preserved. Does anyone know what could be causing this problem?

Thanks very much!

Answer 1

[Adam]

It's quite possibly a bug. If you send us a smallish example on the vcftools-help mailing list, we'll look into it.

[vitti]

Great, thank you!

[Adam]

Apparently there was a bug. Has now been fixed in the SVN download version.