Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Sequence Coverage and Genome Insertion lenght from SAM

    Hi to everyone,

    I'm quite new of this science, and I'm looking for solve a problem.

    I've a SAM file generated by BWA and must to calculate two different things:

    1) Genome insertion lenght ( for plotting mean and standard dev )
    2) Sequence Coverage


    I've calculate the insertion lenght with this operation abs(col[3]-col[7]) made for each row.

    The main problem is that i can understand how to calculate sequence coverage from SAM file, i must implement the algorithm in python or perl, but the prolem is not "the language" is wath operation i need to do?

    Anyone have some ideas?

    thanks in advance to everyone

  • #2
    genomecov (http://bedtools.readthedocs.org/en/l...genomecov.html) or coveragebed (http://bedtools.readthedocs.org/en/l.../coverage.html) tools from Bedtools package can do this operation. You can find graphical representation of what needs to be done (since you want to write your own code) at those links.
    Last edited by GenoMax; 01-14-2016, 04:22 AM.

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Strategies for Sequencing Challenging Samples
      by seqadmin


      Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
      03-22-2024, 06:39 AM
    • seqadmin
      Techniques and Challenges in Conservation Genomics
      by seqadmin



      The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

      Avian Conservation
      Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
      03-08-2024, 10:41 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 06:37 PM
    0 responses
    7 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 06:07 PM
    0 responses
    7 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-22-2024, 10:03 AM
    0 responses
    49 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-21-2024, 07:32 AM
    0 responses
    66 views
    0 likes
    Last Post seqadmin  
    Working...
    X