Question: Does Ncbitblastncommandline Reads The Outfmt (Options 6, 7, And 10 ) Correctly?
2
gravatar for microbeatic
6.6 years ago by
microbeatic80
hanover, NH
microbeatic80 wrote:

I am trying to call

tblastn_hits = NcbitblastnCommandline(cmd='tblastn',
                                         query=query_handle,
                                         db=db_handle,
                                         num_threads=6,
                                         evalue=0.001,
                                         outfmt="6 qstart qend evalue bitscore",
                                         dbsize=10000000,
                                         out="%s%s.txt" % (result_handle,orf_name))

But, when i run the script, it ignores anything after 6 in outfmt, so the blast results only have default columns. Is this a problem with the Biopython version or i am doing something wrong?

biopython blast • 3.2k views
ADD COMMENTlink modified 6.6 years ago by Peter5.8k • written 6.6 years ago by microbeatic80
1

I'm having the same issue when running blast+ executables straight from command line.

ADD REPLYlink written 6.6 years ago by gammyknee200
4
gravatar for Peter
6.6 years ago by
Peter5.8k
Scotland, UK
Peter5.8k wrote:

You need to quote it, the easiest way would be:

tblastn_hits = NcbitblastnCommandline(cmd='tblastn',
                                     query=query_handle,
                                     db=db_handle,
                                     num_threads=6,
                                     evalue=0.001,
                                     outfmt='"6 qstart qend evalue bitscore"',
                                     dbsize=10000000,
                                     out="%s%s.txt" % (result_handle,orf_name))

Here I am using single quotes for the Python string, which means I can use double quotes inside the string.

ADD COMMENTlink written 6.6 years ago by Peter5.8k

Problem with qcovs in blast outfmt 6  ... I'm so grateful to you :D

Seriously , I can't thank you enough . 

ADD REPLYlink modified 5.9 years ago • written 5.9 years ago by Bara'a250

Yes, thank you! I was having this problem on blast+ on the command line when combining it with GNU Parallel.

ADD REPLYlink written 3.3 years ago by johnsenkyle130
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: 1286 users visited in the last hour