Microbial gene coverage from blast result
0
0
Entering edit mode
12 months ago

Hello all I am new to BLAST and Biostars. I have removed human-mapped reads from RNA-Seq data and did Kraken2 and Bracken analysis using the microbial reads. Using the Kraken Tool, I retrieved the microbial reads. Then I performed BLASTx in order to find the microbial genes present in the data. Now, the main question is to find the coverage of those genes. I used the following formula : (query start - query end)/total gene length in nucleotides [from Uniprot] multiply by 100. But I am getting the coverage more than 100% also for a few genes. Is the formula correct? Or please suggest some tool or formula or code to calculate the gene coverage using BLAST result.

Thank You

gene coverage • 424 views
ADD COMMENT
0
Entering edit mode

What is qcovhsp (high scoring segment pair) in BLAST output. Can we use this parameter to calculate the gene coverage?

ADD REPLY

Login before adding your answer.

Traffic: 1088 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