Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • LacquerHead
    Member
    • Nov 2015
    • 31

    #1

    BAM file nt offsets

    Hi wondering if someone knows of a quick and easy way to offset a BAM alignment by N nucleotiodes towards the center of a fragment? Thank you.
  • dpryan
    Devon Ryan
    • Jul 2011
    • 3478

    #2
    I guess it depends on what you're trying to achieve. If you're trying to generate sharper coverage tracks then you can use bamCoverage from deepTools with either the --MNase or --centerReads options. Those both operate on fragments.

    Comment

    • LacquerHead
      Member
      • Nov 2015
      • 31

      #3
      Transposase

      I am trying to account for the binding of Transposase by offsetting the read alignments on the +/- strands.


      Originally posted by dpryan View Post
      I guess it depends on what you're trying to achieve. If you're trying to generate sharper coverage tracks then you can use bamCoverage from deepTools with either the --MNase or --centerReads options. Those both operate on fragments.

      Comment

      • dpryan
        Devon Ryan
        • Jul 2011
        • 3478

        #4
        Provided you're happy with a bigWig file then either of the options I presented will work. Otherwise you might have to script something using pysam in python.

        Comment

        • LacquerHead
          Member
          • Nov 2015
          • 31

          #5
          I guess bamCoverage does not give the option to offset reads by specific numbers of bases given the binding event of specific enzymes?

          Comment

          • dpryan
            Devon Ryan
            • Jul 2011
            • 3478

            #6
            The current release doesn't offer this. However I happen to have a test branch that can do this (feature/riboseq_352). This was initially written with RiboSeq in mind, so there's a "--RiboSeq X" that will then use the Xth base from the 5' end of the read. I could change this to allow negative values, which would then support things like GROseq and PROseq. This is all in development, but it'll likely give reasonable results.

            Comment

            Latest Articles

            Collapse

            • SEQadmin2
              Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
              by SEQadmin2



              Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

              Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
              ...
              Today, 09:55 AM
            • SEQadmin2
              New Genomics Technologies Take Aim at Long-Standing Limits
              by SEQadmin2


              Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

              We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
              ...
              09-28-2026, 10:25 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by SEQadmin2, Today, 05:39 AM
            0 responses
            10 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, Yesterday, 09:56 AM
            0 responses
            12 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 10-09-2026, 11:14 AM
            0 responses
            26 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 09-29-2026, 09:51 AM
            0 responses
            77 views
            0 reactions
            Last Post SEQadmin2  
            Working...