blastn: how to show all query list in my blastn results, even if they have no match.
2
0
Entering edit mode
6.2 years ago

Hi,

I am comparing a list of sequences from my sequencing results with public database and saving just the best record like this:

blastn -task blastn -query ${i}.fa -db nt -remote -outfmt '6 qseqid sseqid pident evalue staxids stitle' -num_alignments 1 -out ${i}_blastnxnt.csv

But, I would like to have one line per query line, saving all ids, even those not found in the database. Does anyone know how to do it?

Gracias!

alignment • 3.5k views
ADD COMMENT
1
Entering edit mode
ADD COMMENT
0
Entering edit mode

That's it, with a value for the whole query list.

ADD REPLY
1
Entering edit mode
6.1 years ago
Chirag Parsania ★ 2.0k

Before you move ahead to find find one best hit, I would strongly suggest you to read recently published paper talking about finding best hit and its misunderstanding

ADD COMMENT
0
Entering edit mode

+1 for that Chirag..!

ADD REPLY

Login before adding your answer.

Traffic: 2654 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6