This is a test version of Biostars. For the public version, visit https://www.biostars.org.
1000G mask with VCFtools

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!

vcftools mask 1000genomes

1 answer

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

Great, thank you!

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

Log in to answer this question.