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
      Advanced Methods for the Detection of Infectious Disease
      by seqadmin




      The recent pandemic caused worldwide health, economic, and social disruptions with its reverberations still felt today. A key takeaway from this event is the need for accurate and accessible tools for detecting and tracking infectious diseases. Timely identification is essential for early intervention, managing outbreaks, and preventing their spread. This article reviews several valuable tools employed in the detection and surveillance of infectious diseases.
      ...
      11-27-2023, 01:15 PM
    • seqadmin
      Strategies for Investigating the Microbiome
      by seqadmin




      Microbiome research has led to the discovery of important connections to human and environmental health. Sequencing has become a core investigational tool in microbiome research, a subject that we covered during a recent webinar. Our expert speakers shared a number of advancements including improved experimental workflows, research involving transmission dynamics, and invaluable analysis resources. This article recaps their informative presentations, offering insights...
      11-09-2023, 07:02 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 08:26 AM
    0 responses
    11 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 08:12 AM
    0 responses
    12 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 11-27-2023, 08:12 AM
    0 responses
    19 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 11-22-2023, 09:29 AM
    1 response
    60 views
    0 likes
    Last Post VilliamPast  
    Working...
    X