Question: blasting genome contigs against local SILVA 16S RNA database
0
gravatar for robjohn70000
2.9 years ago by
robjohn7000090
Nigeria
robjohn7000090 wrote:

Hi All,

I'm new to bioinformatics and interested in determining the bacteria strains in a collection of bacteria genome contigs. I wonder if is logical and appropriate to create a local SILVA 16S RNA database and then blast (using blastn) the genome contigs against the 16S RNA  database. My main concern is that, will it be logical to blast DNA sequences against RNA sequences in the 16S RNA database.

I welcome some advice. Thanks


 

sequencing blast genome • 1.5k views
ADD COMMENTlink modified 2.9 years ago by marina.v.yurieva380 • written 2.9 years ago by robjohn7000090
0
gravatar for Istvan Albert
2.9 years ago by
Istvan Albert ♦♦ 73k
University Park, USA
Istvan Albert ♦♦ 73k wrote:

if your RNA sequences are represented in the DNA alphabet and then there shouldn't be any problem - that's standard protocol

ADD COMMENTlink written 2.9 years ago by Istvan Albert ♦♦ 73k

Thanks @Istvan. I noticed that the RNA sequences have 'U's in place of 'T's. 
 

ADD REPLYlink written 2.9 years ago by robjohn7000090

NCBI provides taxonomically annotated 16S databases that can be used directly.

ADD REPLYlink written 2.9 years ago by Istvan Albert ♦♦ 73k
0
gravatar for marina.v.yurieva
2.9 years ago by
Farmington, CT
marina.v.yurieva380 wrote:

If you haven't done this type of analysis before you might look into developed pipelines: mothur and QIIME. For standard 16S analysis blasting is not very efficient.

ADD COMMENTlink written 2.9 years ago by marina.v.yurieva380

Indeed, for 454 type of data it may work, but for large number of reads such as from MiSeq the process requires a lot of resources

ADD REPLYlink written 2.9 years ago by Istvan Albert ♦♦ 73k

Thank you all very much for your contributions and advice. I will be looking at the the suggested tools.
 

ADD REPLYlink written 2.9 years ago by robjohn7000090
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1394 users visited in the last hour