Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Trimmoatic: Command for cutting bases from the start and end of the read

    Hi everyone,

    I am working with small non coding RNA data. I removed my adapter sequences from my data. However, there are random miRS of 4 bases attached to both the start and end of the reads. Hence, i wanted to remove 4 bases both from the start and end of the read. I used the command HEADCROP 4 to remove bases from the start of the read. However, i am not quite sure which command to use, to remove bases from the end of the read. the CROP command doesn't cut the bases at the end of the read, however it cuts the read to specified length by cutting the bases from the end of the read. The TRAILING command could be used, but i am not quite sure what would be threshold quality value that could trim the 4 random mirs. Because i need all the 4 bases to be uniformly trimmed across all the reads. Any thoughts? Opinion? Suggestions? I know cutadapt could be used in this case. But since i have already used trimmomatic to cut adapters i want to use the same for trimming the 4 bases from the start and end of the read.

    Thanks

  • #2
    ^^^^ or in other words i am wondering if there is a TRAILCROP option in Trimmomatic

    Comment


    • #3
      The problem with removing bases from the 3' end of the reads is that, even after adapter trimming, there may be some reads that had only a few bases of adapter sequence, so too short to be removed by trimmomatic, and you would be removing those instead of the miRS sequences.

      As far as I know there isn't an option that does what you want in Trimmomatic, unless you start with reads all the same length without any adapter contamination.

      Comment


      • #4
        Originally posted by newfie View Post
        ^^^^ or in other words i am wondering if there is a TRAILCROP option in Trimmomatic
        CROP: Cut the read to a specified length

        Comment


        • #5
          Originally posted by Zaag View Post
          CROP: Cut the read to a specified length
          Yes, but that won't help if your reads are all different lengths because you have removed variable lengths of adapters from the 3' ends.

          Comment


          • #6
            True, but you can leave the adapters on and just cut adapter + 4 bp

            I know PrinSeq has this option:
            -trim_right

            Trim sequence at the 3'-end by trim_right positions


            but I never used it.

            Cheers.

            Comment


            • #7
              @newfie: You could use the CROP command first followed by HEADCROP/trimming if you want to stick with trimmomatic. (If you want to remove an additional 4 bp from what was left after the trimming, then this won't work).

              bbduk.sh will allow you to forcetrimright=n past a certain position.

              Comment

              Latest Articles

              Collapse

              • 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
              • seqadmin
                Recent Developments in Metagenomics
                by seqadmin





                Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...
                09-23-2024, 06:35 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, Yesterday, 02:44 PM
              0 responses
              7 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 10-11-2024, 06:55 AM
              0 responses
              14 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 10-02-2024, 04:51 AM
              0 responses
              110 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 10-01-2024, 07:10 AM
              0 responses
              117 views
              0 likes
              Last Post seqadmin  
              Working...
              X