blast alignemnt query coverage

Question

[sharmatina189059]

[Karma]

Why don't you state the error that you are getting while running the command

[RamRS]

The pattern of somehow receiving upvotes for trivial comments continues into your new account, Mask/snijesh. If you're encouraging your colleagues to do this, as you mentioned in your previous comments, please do not do that. If not, I direct this message to the OP sharmatina189059 and ask them to use upvotes with caution and award them to posts that contribute significantly to the subject matter at hand.

[Karma]

I am not sure what you are worried about. I did not tell anyone to upvote my comments.

[RamRS]

Not to debate semantics, but I'm not "worried" about anything.

You remarked on how Mask was your previous account (See: C: RNA-Seq or Exome Seq for Clinical data analysis ). Using the older account, you also said you'd asked your colleagues to up-vote your comments. (See: C: Getting error in running GATK (GenomeAnalysisTK-3.7)? ), which we observed as odd owing to them being up-votes to relatively trivial comments.

Here, the pattern of receiving up-votes to trivial comments continues, so I asked you if your colleagues were still doing this, or if this were a one-off occurrence (a possibility I addressed with the comment directed to /u/sharmatina189059 ).

I hope this clears stuff up for you.

[Karma]

I am not using the account you mentioned and I don't know who upvoted my comments here.

[RamRS]

Thank you for clearing that up. New users tend to upvote any response they get, so that might be the case here.

[Karma]

Just try the same command with quotes The command string needs these quotes, such that the command line looks like:

blastn -outfmt "6 qcovs" ...

[sharmatina189059]

Thanks ..it works but which coloumn gives query coverage

[Karma]

why don't you upload your results.. at-least first 10 lines

[genomax]

Output should be in the same order as what you specified in -outfmt 6 "sseqid qseq qcovs sseq" gives you

gi|672239066|ref|NR_125655.1|   AACTGACAGTTTTAGAAACAGGACT   2   AATTGACAGTTTTAGAAATAGAACT
gi|631253076|ref|NR_114274.1|   AACTGACAGTTTTAGAAACAGGACT   2   AATTGACAGTTTTAGAAATAGAACT
gi|636560310|ref|NR_116370.1|   AACTGACAGTTTTAGAAACAGGACT   2   AATTGACAGTTTTAGAAATAGAACT
gi|343198619|ref|NR_043083.1|   CTAAAAGTTACCCCACCGAC    2   CTAAAGGTTACCCCACCGAC
gi|343201461|ref|NR_042187.1|   AACTGACAGTTTTAGAAACAGGACT   2   AATTGACAGTTTTAGAAATAGAACT

Here qcovs is in column 3.