Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • extracting identical reads from bam file

    Hello,
    I have an alignment file of smallRNAs in bam format.
    How can I extact from this file all the identical reads in a single sequence but mantaining the reads count?
    Thank you for the help!

  • #2
    Are you trying to match a specific sequence ?

    samtools view x.bam | grep TTTTCTGCCTGTTGGGCTGGAG | awk '{print $10}' | uniq -c


    All reads with same sequence as another read, try this ...

    samtools view x.bam | awk '{print $10}' | sort --buffer-size=20G | uniq -c | awk '{if ($1!=1) print $0}'


    where x.bam is you bam file.
    fine tune the --buffer-size parameter to sort

    Comment


    • #3
      It would take longer, but swap the second awk in Richard's command line with

      Code:
      sort -nr
      To get the most common reads in order, starting with the most abundant. Tack a

      Code:
      | head -n 100
      to get only the top 100

      Comment


      • #4
        thank you! it worked!
        can I do the same from a fastq file?

        Comment


        • #5
          Yes, you might want to use awk to only print every 4th line.
          Note 1) the "mod" operator and 2) "line count" intrinsic variable in awk.
          Perl/python/C/java if you prefer can address the issue of filtering for only the sequence also.

          Comment


          • #6
            After I extracted the identical reads in a single sequence from the bam file I aligned them again to the genome. When I use igv to visualize the alignment now all the sequences are mapping in sense orientation.. even those sequences that are supposed to be antisense to the genome are shown in sense orientation. Why is that?

            Comment


            • #7
              The sequence for the read in a bam files may be reverse complemented to align to the reference.
              Reads are supposed to be properly noted as reversed in the bitwise flags field in a line/entry of sam/bam.

              You may wish to interrogate this flag for special processing.

              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