This is a test version of Biostars. For the public version, visit https://www.biostars.org.
STAR aligned reads minimum maping quality

I aligned RNAseq data with GRCH 37 and wonder what cut off of read quality is generally used when we perform downstream analysis. I know in HISAT2 generally read quality >20 is considered. I could not find any general rule of quality cutoff in STAR. Any suggestion?

rna-seq

when we perform downstream analysis.

I guess you are talking about differential expression analysis, but it doesn't hurt to be specific...

2 answers

Take a look at this blog post by Dr. Simon Andrews and decide for yourself.

Like mentioned in the blog post of genomax, I always go with the cutoff of 255, resulting in only uniquely mapped reads.

Log in to answer this question.