This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Get Query Coverage From Local Blast

Hi all

Does anyone have an idea how I can get the query coverage of my hits from local BLAST (in unix). The -m 8 option gives you a tabular output that does not contain the 'Descriptions' part of the online BLAST. If I blast a gene against a genome with no contigs, I expect to see an overall percentage similarity of the gene against this genome, which is given in the Descriptions section in online BLAST, but cannot get that in local BLAST. One can get a rough idea from -m 8 option but with so many overlaps it can prove to be a bit difficult. Any ideas anyone?

Thanks

blast coverage
• 18,375 views • link
updated by Ram Ram • written by bioinf80 • bioinf80

2 answers

Hi,

Please use the following awk command on the blast tabular output:

awk '{if ($5/$3 > 0.75 && $5/$4 > 0.75 && $6>55 && $7<0.000000000000001) print $0}'  blast_out.tab
  • $5/$3 > 0.75 = Alignment length should be > than 75% of query length;
  • $5/$4 > 0.75 = Alignment length should be > than 75% of Subject length;
  • $6>55 = Percent identity should be > than 55%;
  • $7<0.000000000000001 = e-value less than e-15.

Use the desired combinations and please check the columns in your tab file.

Hope it would help you.

Wishes,
Rahul

• 144 views • link
updated by Ram Ram • written by Rahul Sharma Rahul Sharma

Hi

I am trying to use above awk command to get query coverage but there is no output or any error message

I would like to get 70 % query coverage. My blast output is below.

664     gi|157146441|ref|YP_001453760.1|        39.23   130     79      0       3       392     5       134     6e-24    104
664     gi|157083646|gb|ABV13324.1|     39.23   130     79      0       3       392     5       134     6e-24    104
664     gi|146311054|ref|YP_001176128.1|        38.46   130     80      0       3       392     5       134     1e-23    103
664     gi|145317930|gb|ABP60077.1|     38.46   130     80      0       3       392     5       134     1e-23    103
664     gi|160865276|gb|ABX21899.1|     37.69   130     81      0       3       392     5       134     2e-22   99.8
664     gi|161503929|ref|YP_001571041.1|        37.69   130     81      0       3       392     5       134     2e-22   99.8
664     gi|117623058|ref|YP_851971.1|   37.69   130     81      0       3       392     5       134     3e-22   99.4
664     gi|145566777|sp|A1A9B7.1|MACB1_ECOK1    37.69   130     81      0       3       392     5       134     3e-22   99.4
664     gi|122064321|sp|Q0TJH0.1|MACB_ECOL5     37.69   130     81      0       3       392     5       134     3e-22   99.4
664     gi|75210431|ref|ZP_00710587.1|  37.69   130     81      0       3       392     5       134     3e-22   99.4

Help is appreciated

Thank you

Raju

• 0 views • link
updated by Ram Ram • written by indugun • indugun

Hi Raju,

Please add column headers, I am wondering if you had used the default parameters or specified your own set of columns.

Wishes,

Rahul

• 0 views • link
written by Rahul Sharma Rahul Sharma

Hi R@hul,

I am new to blast tabular format. I have blat output is blast 8 format. I have also provided the headers. Could you guide how to calculate the coverage percentage from the below table. Kindly guide me.

Qry     Tgt    %_Iden.        Align_Length       mismatch        gap_openings    q_start    q_end    s_start    s_end     E-value      bit_score
X_1     Y_4    99.46          1470               8               0               1          1470     1          1470      00           2818
X_3     Y_5    93.75          64                 4               0               1          64       2371       2434      8.10E-24     108
X_4     Y_9    99.94          1574               1               0               3          1576     1          1574      0.00E+00     3033
• 0 views • link
updated by Ram Ram • written by EVR EVR

You can use awk command as suggested by Rahul since you can't get a direct query coverage percentage as you would see in the GUI. You can add the awk to the command line itself

blastn -i input.fasta -m 8 -db db_name -o out | awk '{add appropriate conditions}'

cheers
vijay

• 171 views • link
updated by Ram Ram • written by vijay vijay

Log in to answer this question.

More posts like this