Question: BLAST qcov equals to 0
1
gravatar for Ming
6 months ago by
Ming30
Ming30 wrote:

Hello All,

I just performed a blastn, and I have some alignment returning qcovs = 0. I am not sure how this is possible because the alignment reports a pident and a bitscore for that alignment. For instance:

    qseqid sseqid pident bitscore evalue qcovs qcovhsp length
k127_1554010 ref|NZ_AONQ01000017.1| 77.007 409 2.85e-110 0 0

From my understanding, qcovs = query coverage. A qcovs = 0 means that 0% of the query is aligned to the subject. But why is the pident and bitscore reported as positive values?

Thank you in advanced!

blast • 203 views
ADD COMMENTlink modified 4 months ago by h.mon27k • written 6 months ago by Ming30
1

Can you provide the parameters that you have defined for your blast?

ADD REPLYlink written 4 months ago by mgdias.jose10
1

that output is not corresponding to the header you provide, length has no value (from the top of my head, I even think 'length' is not a valid field)

can you provide your actual blast command (as also mgdias.jose requested)

ADD REPLYlink modified 4 months ago • written 4 months ago by lieven.sterck6.0k
1
gravatar for h.mon
4 months ago by
h.mon27k
Brazil
h.mon27k wrote:

As lieven.sterck noted above, the output does not correspond to the header - which may indicate you did something wrong / the output fields are not exactly what you think they are.

I will have an attempt at another (rather wild) guess:

The name of your query and subject indicates you are blasting MEGAHIT contigs against bacterial genomes - potentially very long sequences. If you have a strong short similarity, blast will output the hit, but qcov may be set to zero because it is a tiny fraction of a percentage. Try to use qstart, qend, sstart and send to see how long are the alignments.

ADD COMMENTlink written 4 months ago by h.mon27k
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: 1962 users visited in the last hour