SEQanswers

Go Back   SEQanswers > Introductions



Similar Threads
Thread Thread Starter Forum Replies Last Post
BWA, I don't know what I'm doing asheenlevrai General 0 12-29-2011 11:35 AM
Qiagen SeqTarget system - Don't go there friend !! lek2k Sample Prep / Library Generation 2 08-18-2011 06:32 PM
Samtools don't read the whole data mikael Bioinformatics 6 08-09-2011 12:37 PM
Bam and Sam don't like my fasta file mindlessbrain Bioinformatics 2 12-09-2010 11:47 PM
something wrong with tophat-1.0.14 version, don't use it frankyue50 Bioinformatics 4 10-04-2010 09:33 AM

Reply
 
Thread Tools
Old 12-13-2010, 09:59 AM   #1
anna_
Member
 
Location: germany

Join Date: Dec 2010
Posts: 15
Red face Don't know if I'm right

Hi everyone,

we have just sequenced a whole trancriptome of a conifer de novo and now it is my challenge to find candidate loci in there. I want to perform a stand alone blast search. I'm so new to the subject that I don't even have a concrete question.

Anyway - hello to all!
anna_ is offline   Reply With Quote
Old 12-13-2010, 10:33 AM   #2
westerman
Rick Westerman
 
Location: Purdue University, Indiana, USA

Join Date: Jun 2008
Posts: 1,104
Default

Well, if you are new to sequencing then I suspect that you do not have access to a stand-alone blast server. While it is fairly easy to set one up, you may wish to just try Blasting several of your transcriptomes to a public blast server (e.g., America's http://blast.ncbi.nlm.nih.gov/ or Europe's EMBL http://dove.embl-heidelberg.de/Blast2/ ... and others)

Actually for a one-off project I would suggest just using a public resource instead of your own server unless you have special requirements (e.g., Blasting against a unusual database). The public resource may be slower than your own server but it will not be that much slower and will be much easier to use.
westerman is offline   Reply With Quote
Old 12-13-2010, 11:03 AM   #3
anna_
Member
 
Location: germany

Join Date: Dec 2010
Posts: 15
Default

Thanks for that quick comment.
Our transcriptome has 500 000bp, the LargeContig FASTA-File has 7.5MB. I want to find protein coding sequences in there. The plan is to start with proteins from Arabidopsis and Rice. Therefore I think I need tblastn. I guess that can't be done online. I had the idea to format my FASTA-file to a database and blast the know proteins of the plants against it.

At the moment I'm trying to install the executables available at the NCBI site.
anna_ is offline   Reply With Quote
Old 12-13-2010, 11:06 AM   #4
anna_
Member
 
Location: germany

Join Date: Dec 2010
Posts: 15
Default

Sorry, I meant 500 000 reads (454)
anna_ is offline   Reply With Quote
Old 12-13-2010, 12:32 PM   #5
westerman
Rick Westerman
 
Location: Purdue University, Indiana, USA

Join Date: Jun 2008
Posts: 1,104
Default

tblastn can be run on-line as well as blastx. To do what I suggested you would use blastx -- search the rice and arabidopsis databases with your nucleotide 454 sequences, translated. Your idea of using tblastn to look through your translated 454 sequences with the rice/arabidopsis proteins is also good.

With 500 000 reads it is probably is better to install your own local blast. As I said it is easy enough to do, at least from the command line. Spot checking your results with the results you can get from a public site is recommended. Depending on the power of your computer the blasting can take many hours or days. I suggest reducing your data set to 5000 or so reads for testing purposes before you try all 500 000.

You may also wish to consider a program that can take splicing into account (e.g., splign).


If you do come up with a concrete question then post on SeqAnswers. There are many knowledgeable people here.
westerman is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 07:48 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO