Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Per base sequence coverage from sam/bam file?

    Hi all,

    Does anyone have recommendations for good progams for calculating the per base sequence coverage for mapped reads in sam or bam format? I'd like for the output format to be human readable (not tdf) so I can input into R for downstream analysis.

    Sorry if this seems basic-- I just haven't been able to lay hands on just the right script and I'm sure there's one out there somewhere that will be faster than something my newbie-self can program!

    Thanks!!!
    Lizzy

  • #2
    "samtools pileup"

    Comment


    • #3
      I think GATK -T DepthOfCoverage also does the job. I haven't tried the per-base coverage in GATK myself (= ran it with -omitBaseOutput), but I think samtools pileup omits bases with no coverage, so might require additional additional scripting after running the pileup. Please correct me if I'm wrong!

      Comment


      • #4
        Hi,

        You can use genomeCoverageBed from the http://code.google.com/p/bedtools/ which can also read in BAM files.

        Cheers

        Comment


        • #5
          Hi lizzy,

          see my older thread at http://seqanswers.com/forums/showthread.php?t=7679

          Hope this helps,
          Boetsie

          Comment


          • #6
            Thanks everyone! I knew you'd all have great suggestions

            Comment


            • #7
              Thanks for the great suggestions. I find it very useful.

              Comment


              • #8
                I find "samtools depth" pretty useful. The fact that you can get coverage for a list of target regions and also look for coverage with different base and mapping quality.

                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
                10 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, Yesterday, 06:07 PM
                0 responses
                9 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-22-2024, 10:03 AM
                0 responses
                51 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-21-2024, 07:32 AM
                0 responses
                67 views
                0 likes
                Last Post seqadmin  
                Working...
                X