Hi everyone, I am having a difficult time in dealing with the SNPs from a homozygous genome. Theoretically speaking, I suppose there should be no SNPs detected in the mitotic homozygous genome if no reads are misassembled in a de novo way. I used CLC and SOAPdenovo to assemble this homozygous genome, but still extracted more than 10,000 SNPs in the high coverage regions. Many frequencies were almost 50/50. Before performing the whole genome sequencing using this assumed homozygous individual, I used two microsatellite markers per chromosome to confirm its homozygosity. The genotyping results came out to be 100% homozygous at those loci. So whether it was due to the really variations on two chromosomes or due to the sequencing error, missassemly or something, I really have no idea. During assembly, there could be more or less mistakes, couldn't there? I think it's impossible to get no SNPs. But the boss wants me to prove its complete homzogyosity, now I am sunk. Could somebody show me some points?
Best,
Lisa
Best,
Lisa
Comment