When Maq is used for RNA-seq data analysis, I saw several people set (q) "Minimum mapping quality allowed for a read to be use" to be 30 as the cutoff value in %Maq pileup. Does this mean the reads after this "pileup" with q=30 are those "mapped uniquely"?
But, since "most of the repetitive reads have a mapping score of 0", why just set the q to be a value a little bit > 0 (other than 30) to discard the non-unique reads? Or in another word, what value should we set for q such that we can confidently say that we have filtered out all the non-unique mappings?
Thanks
But, since "most of the repetitive reads have a mapping score of 0", why just set the q to be a value a little bit > 0 (other than 30) to discard the non-unique reads? Or in another word, what value should we set for q such that we can confidently say that we have filtered out all the non-unique mappings?
Thanks