SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
blastall output: NCBI vs command line ElMichael Bioinformatics 8 11-17-2011 12:52 PM
BLAST+ vs BLASTALL (legacy BLAST) Symphysodon Bioinformatics 4 10-25-2011 03:52 PM

Reply
 
Thread Tools
Old 03-09-2011, 08:21 AM   #1
empyrean
Member
 
Location: EU

Join Date: Sep 2010
Posts: 52
Default Blastall help !!!

Hello All,

I need some help in optimizing my blastall query. I have a file of 10000 genes, and i wanted to blast them against themselves. So, i created blast db to perform blastall. But i dont want the genes that map to themselves. I would like to get second best hit from the results. Can i do this in blastall using some parameters or do i need to write a perl script for this?

Code:
blastall -p blastn -m 9 -d testgenesdb -i test.fasta -o outputtest.QUERY
That is the query i uses and it brings out blast results in tabular format..
please help me if u know optimizing blast queries.

Thank you all
empyrean is offline   Reply With Quote
Old 03-09-2011, 11:01 AM   #2
klaas.vandepoele
Junior Member
 
Location: Ghent

Join Date: Jan 2011
Posts: 6
Default example

code exclude_self.pl

while (<STDIN>)
{
my @l= split(/\t/,$_);
if ($l[0] ne $l[1])
{
print join("\t",@l)."\n";
}
}


--

perl exclude_self.pl < outputtest.QUERY > outputtest.QUERY.noself

Something like this should work.

best,
Klaas
klaas.vandepoele is offline   Reply With Quote
Old 03-09-2011, 12:45 PM   #3
empyrean
Member
 
Location: EU

Join Date: Sep 2010
Posts: 52
Default

Thank you for the code.. this is what i wanted but i need only second best hit.. i dont want all the results...how can i achieve it?

Last edited by empyrean; 03-10-2011 at 08:12 PM.
empyrean is offline   Reply With Quote
Old 03-10-2011, 08:13 PM   #4
empyrean
Member
 
Location: EU

Join Date: Sep 2010
Posts: 52
Default

How to get the second best hit or after running above code, i need the first row with best evalue for all the combinations..
empyrean 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 06:08 AM.


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