View Single Post
Old 08-03-2011, 01:15 AM   #1
vebaev
Senior Res.
 
Location: Plovdiv, Bulgaria

Join Date: Oct 2008
Posts: 110
Default mapping reads to known miRNAs precursors

Dear all,
I'm trying to estimate the all reads count that are mapped to a specific precursor. I have done it my using bowtie with no mistmaches -v 0
But in some papers and tools I see that thay allowed 1 mistmache?! I compared the results in both way. Most of the counts are similar but for some miRNAs thay are doubled like in mir143 bellow:

with 1 mistmach mapping:
miR-143 - 4,697,060
miR-10b - 520,833

with 0 mistmach mapping:
mir-143 - 1,976,574
mir-10b - 501,099

So which is the right way? I thought it should be with 0 mistmaches since non specific tags can align when mistmaches occure but some people do it? And if I want to estimate diff.expression I should know which method is most relevant?

Thanks
vebaev is offline   Reply With Quote