Question: Receiving low query covers for blast searches on genome files, but not transcript files
0
gravatar for Jacob
22 months ago by
Jacob10
Jacob10 wrote:

I'm performing blast searches with only exon regions for certain genes against various genomes(for different organisms). For example, I will blast regions 1..367, 368..444, 445..828 of the gene ACAN(https://www.ncbi.nlm.nih.gov/nuccore/NM_013227.3) against an organisms genome such as

ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/002/315/GCF_000002315.4_Gallus_gallus-5.0/GCF_000002315.4_Gallus_gallus-5.0_genomic.fna.gz

For a lot of my blast searches, I'm only receiving hits that have a max query cover of about 30-40. If i perform the same blast searches on the transcript file

ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/002/315/GCF_000002315.4_Gallus_gallus-5.0/GCF_000002315.4_Gallus_gallus-5.0_rna.fna.gz

On almost all my blast searches I am receiving hits that have a query cover around 80-100

What could be the reason for this. Exons should be completely connected within the genome shouldn't they? Why would blasting against the rna transcript files return hits with higher query covers than blasting against the genome?

rna-seq transcript blast genome • 544 views
ADD COMMENTlink written 22 months ago by Jacob10
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: 1411 users visited in the last hour