Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Samtools index/sort/rmdup

    I am trying to filter several steps before the mpileup:
    1) samtools sort -n in.bam out.bam # sort the bam file
    2) samtools rmdup -sS in.srt.bam out.bam # remove PCR duplicates
    3) samtools view -b -F 1024 -q 20 *.bam > out.bam # filter duplicates reads and mapping quality
    4) samools index .bam # give a index so that for mpileup

    However, for the index step, it shows
    "[bam_index_core] the alignment is not sorted: reads without coordinates prior to reads with coordinates." Anything wrong here?

  • #2
    samtools sort -n sorts by name, not coordinate.. index requires sort by coordinate..

    Comment


    • #3
      Originally posted by peromhc View Post
      samtools sort -n sorts by name, not coordinate.. index requires sort by coordinate..
      So that I can sort without "-n" again and then index the bam file? Thanks.

      Comment


      • #4
        Originally posted by hanleng View Post
        So that I can sort without "-n" again and then index the bam file? Thanks.
        yep.. you should check out the git version with multithreaded sort -@ flag.. you're doing a lot of sorting!

        Comment


        • #5
          Originally posted by peromhc View Post
          yep.. you should check out the git version with multithreaded sort -@ flag.. you're doing a lot of sorting!
          I do not want to a lot of sorting - I just want to filter some reads before doing the mpileup...

          Comment


          • #6
            Originally posted by hanleng View Post
            I am trying to filter several steps before the mpileup:
            1) samtools sort -n in.bam out.bam # sort the bam file
            2) samtools rmdup -sS in.srt.bam out.bam # remove PCR duplicates
            3) samtools view -b -F 1024 -q 20 *.bam > out.bam # filter duplicates reads and mapping quality
            4) samools index .bam # give a index so that for mpileup
            samtools rmdup doesn't mark duplicates with 1024, it just removes them. So you don't need to filter for that. And why not filter before you sort? It would save a bit of time.

            And mpileup doesn't actually require an index.

            Comment


            • #7
              Originally posted by swbarnes2 View Post
              samtools rmdup doesn't mark duplicates with 1024, it just removes them. So you don't need to filter for that. And why not filter before you sort? It would save a bit of time.

              And mpileup doesn't actually require an index.
              According to my understanding, rmdup remove PCR duplicates, and -F 1024 remove duplicates reads, right? I just want to remove all kinds of duplicates to call SNP, so is this correct?

              Comment


              • #8
                Originally posted by hanleng View Post
                According to my understanding, rmdup remove PCR duplicates, and -F 1024 remove duplicates reads, right? I just want to remove all kinds of duplicates to call SNP, so is this correct?
                -F 1024 will remove all reads that have the 1024 flag set. If you don't run any software which will ever set that flag, then there's nothing with that flag to remove. I don't think anything in the samtools suite will ever set that flag. Picard's MarkDuplicates will, but samtools rmdup won't.

                Comment

                Latest Articles

                Collapse

                • seqadmin
                  Latest Developments in Precision Medicine
                  by seqadmin



                  Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

                  Somatic Genomics
                  “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
                  05-24-2024, 01:16 PM
                • seqadmin
                  Recent Advances in Sequencing Analysis Tools
                  by seqadmin


                  The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
                  05-06-2024, 07:48 AM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, 05-24-2024, 07:15 AM
                0 responses
                198 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 05-23-2024, 10:28 AM
                0 responses
                219 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 05-23-2024, 07:35 AM
                0 responses
                228 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 05-22-2024, 02:06 PM
                0 responses
                12 views
                0 likes
                Last Post seqadmin  
                Working...
                X