Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Need to compare multiple mpileups for coverage overlap

    I have 11 mpileups generated by novoalign from whole genome sequencing of Daphnia, which has a genome size of ~200MB. Each is the result of alignment of about 130M paired-end illumina reads. I find that all of the mpileups have about 110M lines, implying that we have info for about 55% of the reference. I want to quantify the degree of overlap in coverage among the mpileups. Is there an easy way to do this with samtools or something similar? Also, is there a way to construct an mpileup file that has a line for every site in the reference, rather than just those with read coverage?

  • #2
    Look into the bedtools package for easy ways to do things like this:

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Non-Coding RNA Research and Technologies
      by seqadmin




      Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.

      Nobel Prize for MicroRNA Discovery
      This week,...
      10-07-2024, 08:07 AM
    • seqadmin
      Recent Developments in Metagenomics
      by seqadmin





      Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...
      09-23-2024, 06:35 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Today, 06:55 AM
    0 responses
    7 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 10-02-2024, 04:51 AM
    0 responses
    104 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 10-01-2024, 07:10 AM
    0 responses
    112 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 09-30-2024, 08:33 AM
    1 response
    117 views
    0 likes
    Last Post EmiTom
    by EmiTom
     
    Working...
    X