View Single Post
Old 03-04-2015, 02:24 PM   #11
rahulvrane
Junior Member
 
Location: Melbourne

Join Date: Apr 2012
Posts: 6
Default

Technically you dont need to go that far. Since blastn/blastp scores for latest blast+ packages is already sorted with top hit as the first hit, you can just use the following.

Code:
blastn -query transcripts.fa -db Targets.fa [other options you like ] -outfmt 7 -out blast_output.txt

cat blast_output.txt |awk '/hits found/{getline;print}' | grep -v "#" > top_hits.txt
Rahul

Last edited by rahulvrane; 03-04-2015 at 02:27 PM.
rahulvrane is offline   Reply With Quote