Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • 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.

  • #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


    • #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


      • #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


        • #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


          • #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

            • seqadmin
              Choosing Between NGS and qPCR
              by seqadmin



              Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
              10-18-2024, 07:11 AM
            • 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

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, Yesterday, 06:09 AM
            0 responses
            10 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 10-30-2024, 05:31 AM
            0 responses
            13 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 10-24-2024, 06:58 AM
            0 responses
            22 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 10-23-2024, 08:43 AM
            0 responses
            52 views
            0 likes
            Last Post seqadmin  
            Working...
            X