mutation or mismatch rate per base
Hello!
Greetings from New Haven~
I'm struggling, learning, and having fun with my solexa sequencing data
analysis. And I'm learning a lot from this site, thanks to everyone!
Now I'm trying to do some pileups to get the mutation rate for each site.
Using samtools pileup command, I can get the reads number and mutations,
but what I would like to see is the number of mutations for each site.
Anybody knows how to do this?
Thank you in advance!
Greetings from New Haven~
I'm struggling, learning, and having fun with my solexa sequencing data
analysis. And I'm learning a lot from this site, thanks to everyone!
Now I'm trying to do some pileups to get the mutation rate for each site.
Using samtools pileup command, I can get the reads number and mutations,
but what I would like to see is the number of mutations for each site.
Anybody knows how to do this?
Thank you in advance!