Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • dacotahm
    Member
    • Oct 2011
    • 28

    #1

    How can I estimate overall coverage against a reference database?

    I'm doing de novo assembly on RNA-seq data and I do a stand-alone BLAST against a reference transcriptome. I want to know how much of my assembly aligned, preferably as a percentage so I can estimate quality and have simple statistics to present along with the usual average and median contig lengths-

    Velvet-Oases 6.3%
    Trinity 7.9%

    etc. Is there a tool out there that will estimate how much of a reference database I have covered? Essentially I could add up the lengths of each aligned transcript and divide by the length of all the reference transcripts....not sure if this makes sense, but I could use some advice.
  • swbarnes2
    Senior Member
    • May 2008
    • 910

    #2
    BEDTools has a coverage function. It will tell you how many bases are covered at each depth.

    Comment

    Latest Articles

    Collapse

    • SEQadmin2
      How Immunogenomics Decodes Immunity’s Genetic Blueprint
      by SEQadmin2




      The immune system’s power comes from its genetic diversity, allowing myriad threats to be neutralized through first recognizing foreign antigens. That diversity is also what makes the immune system so difficult to study. Recent advances in sequencing technology and computational biology, however, are giving researchers new tools to understand immune responses and immune-related diseases in greater detail.

      This convergence of genetics, immunology, and computation...
      Today, 05:41 AM
    • SEQadmin2
      Beyond CRISPR/Cas9: Understand, Choose, and Use the Right Genome Editing Tool
      by SEQadmin2



      CRISPR/Cas9 sparked the gene editing revolution for both research and therapeutics.1 But this system still showed severe issues that limited its applications. The most prominent were the heavy reliance on PAM sequences, delivery limitations, double-stranded breaks that prompt unintended edits and cell death, and editing inefficiency (both in targeting and in knock-in reliability).

      Despite this, “CRISPR helped turn genome editing from a specialized technique into
      ...
      07-31-2026, 11:01 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, 08-24-2026, 10:32 AM
    0 responses
    42 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 08-20-2026, 11:17 AM
    0 responses
    48 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 08-18-2026, 10:05 AM
    0 responses
    55 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 08-13-2026, 12:22 PM
    0 responses
    50 views
    0 reactions
    Last Post SEQadmin2  
    Working...