This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Missing Snp Variants In Two 10Mb Blocks On Chr1 And Chr9 In 1000 Genome Data?

I recently want to use Mach to imputate some SNPs based on 1000 Genome data. Surprisingly, as I followed instruction to chunk ref genome into 10Mb blocks for easy handling, two blocks are of no SNPs at all, one on Chr1 130M-140M, the other on Chr9 50M-60M. So I went back to their allelic.Info file and found there are really big regions without any SNP data. My questions are

  1. Is it normal or something wrong with Mach's 1000G phaseI v3 file?
  2. What's special with these regions? difficult to sequence or too much variations?

Tthanks

snp variant 1000genomes
• 2,886 views • link
updated by Biostar Biostar • written by Zhenyu Zhang Zhenyu Zhang

1 answer

  1. This is normal.
  2. They're centromeres, which are highly repetitive and therefore difficult to sequence cost-effectively. Genome Reference Consortium build 37, which is what the Mach file is based on, simply has large gaps representing the centromeres. Build 38 attempts to model them.
• 0 views • link
written by chrchang523 chrchang523

Thanks a lot. This makes a lot sense.

• 0 views • link
written by Zhenyu Zhang Zhenyu Zhang

Log in to answer this question.

More posts like this

No similar posts found.