I'm using Neanderthal alignment data from here: http://www.eva.mpg.de/neandertal/draft-neandertal-genome/data.html

Specifically, the .bam files. I would like to find alignments for neanderthal sequences to human genomes that state ambiguities (such as M or R). However, the data from that site does not have the ambiguities. Is there a way to find such information? Or are there no ambiguity codes in neanderthal alignments?

My goal is to compare diploid neanderthal sequences to diploid human sequences. How do researchers typically deal with this given the limited data they have for neanderthals?

  • $\begingroup$ Why would you use the draft genome instead of the high coverage Altai genome? $\endgroup$ – AlexDeLarge Jul 27 '16 at 16:45
  • $\begingroup$ Is there more than 1 whole genome sequence (besides Altai) for neanderthals? I'm looking for at least 5. $\endgroup$ – ghgh Jul 27 '16 at 18:03
  • $\begingroup$ Oh, ok. There is the Denisova genome (a close relative of Neandertals), but other than that, you really have to stick to low coverage (it was not clear from the question how many genomes you need). But if you use low coverage, you should avoid using Altai or Denisova genomes for the sake of comparability. $\endgroup$ – AlexDeLarge Jul 28 '16 at 11:25
  • $\begingroup$ Ah ok, thanks. I'm looking for at least 5 or more Neanderthal genomes. I guess I'd have to stick to low coverage then if I want that many. $\endgroup$ – ghgh Jul 28 '16 at 13:37

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.