Nice site with nice tips.
I’d like to ask and advice.
After having performed the alignment, what is the best way to evaluate the quality of mapping? Any particular tools one can use?
Thanks for the compliment. It depends on what type of alignment you’re talking about and what type of problem.
If you’re talking about multiple sequence alignment, there’s BAliBASE, which can be used as a benchmark for alignment tools. But that’s to evaluate an alignment program. It’s probably not what you’re asking for but I thought I’ll throw it in for good measure.
If you’re talking about local alignment tools like BLAST, there’s the E value.
If you’re talking about short read alignments, have a look at the mapping quality. I wrote about mapping qualities here. There are tools such as SAMStat that can summarise mapping qualities from BAM files.
To get a sense of your read mapping quality scores (in phred scale, I’m assuming), check out FastQC. Once you get a sense of what your reads look like as far as quality is concerned, check out Btrim.
Btrim lets you set up search windows that average out your quality score and if a minimum threshold is reached for quality they get filtered out. The nice thing about Btrim is you can also adjust the minimum length of the reads you’d like. So it will only filter reads if it can also maintain a minimum read length you don’t want to go under.
Good luck,
Alex P
P.S.
Awesome advice, Dave!!! Just found your blog and have benefitted nicely from it. Keep up the good work.
We use cookies on our website to give you the most relevant experience by remembering your preferences and repeat visits. By clicking “Accept All”, you consent to the use of ALL the cookies. However, you may visit "Cookie Settings" to provide a controlled consent.
This website uses cookies to improve your experience while you navigate through the website. Out of these, the cookies that are categorized as necessary are stored on your browser as they are essential for the working of basic functionalities of the website. We also use third-party cookies that help us analyze and understand how you use this website. These cookies will be stored in your browser only with your consent. You also have the option to opt-out of these cookies. But opting out of some of these cookies may affect your browsing experience.
Necessary cookies are absolutely essential for the website to function properly. This category only includes cookies that ensures basic functionalities and security features of the website. These cookies do not store any personal information.
Any cookies that may not be particularly necessary for the website to function and is used specifically to collect user personal data via analytics, ads, other embedded contents are termed as non-necessary cookies. It is mandatory to procure user consent prior to running these cookies on your website.
Nice site with nice tips.
I’d like to ask and advice.
After having performed the alignment, what is the best way to evaluate the quality of mapping? Any particular tools one can use?
Hi Nosa,
Thanks for the compliment. It depends on what type of alignment you’re talking about and what type of problem.
If you’re talking about multiple sequence alignment, there’s BAliBASE, which can be used as a benchmark for alignment tools. But that’s to evaluate an alignment program. It’s probably not what you’re asking for but I thought I’ll throw it in for good measure.
If you’re talking about local alignment tools like BLAST, there’s the E value.
If you’re talking about short read alignments, have a look at the mapping quality. I wrote about mapping qualities here. There are tools such as SAMStat that can summarise mapping qualities from BAM files.
Hope that helps,
Dave
I was looking for a tool on how to assess the quality of the mapping reads. So SAMsat seems useful.
Thanks a lot for the response.
Hi nosa,
To get a sense of your read mapping quality scores (in phred scale, I’m assuming), check out FastQC. Once you get a sense of what your reads look like as far as quality is concerned, check out Btrim.
Btrim lets you set up search windows that average out your quality score and if a minimum threshold is reached for quality they get filtered out. The nice thing about Btrim is you can also adjust the minimum length of the reads you’d like. So it will only filter reads if it can also maintain a minimum read length you don’t want to go under.
Good luck,
Alex P
P.S.
Awesome advice, Dave!!! Just found your blog and have benefitted nicely from it. Keep up the good work.
Hi Alex,
Thanks for the advice and the kind compliment! I’m glad it’s been useful.
Cheers,
Dave
there are two that I know. I had a hard time working with RNAseqc though RSeqc works fine.
http://www.broadinstitute.org/cancer/cga/rna-seqc
https://code.google.com/p/rseqc/