This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Blast Parameters For Searching Rna Database

I've a collection of cDNAs (median length ~1500 ) and I wish to filter out those are RNA coding ones thru a blastn search. Can people tell me what typical blast parameters are to be chosen for that kind of search? What will be the typical cutoff values for bit score, evalue, alignment length, low complexity masking(default or soft) etc.?

Thanks in advance

blast rna cdna
• 4,174 views • link
updated by Peter • Peter • written by Woa Woa

Have you considered using Infernal instead of blast?

• 0 views • link
written by Jarretinha Jarretinha

Thanks and yes,I'm looking into Infernal and Rfam

• 0 views • link
written by Woa Woa

2 answers

One easy thing to do is to pick one RNA from the blast database and search for it in the database, then look at the values of evalue and bitscore of the hit with itself. This can give you an idea of the range of the values for an exact match.

Giving an absolute threshold value for the evalue or the bitscore can be a bit risky since they normally depend on the size of the blast database, the size of the query set...

• 0 views • link
written by Marina Manrique Marina Manrique

If you want to search for RNAs with blast, maybe also consider these parameters for the local alignment step: http://stevemount.outfoxing.com/Posting0004.html

• 0 views • link
written by Peter • Peter

Log in to answer this question.

More posts like this