Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • #16
    Look into freebayes. There is samtools/bcftools based option as well.

    Comment


    • #17
      sorted bam and bbmap

      Is there an option for bbmap that will result in a sorted bam? I know adding "out=mapped.bam" this will produce a bam but it appears to be unsorted:

      Code:
      samtools view -H /home/cmccabe/IDP/mapped_unmatched_adapter_removed_quality_trimmed_NA12878-NextSeq-S1_R1_and_R2.bam
      @HD	VN:1.4	SO:unsorted
      Thank you .

      Comment


      • #18
        No. You have to do some work

        Code:
        $ samtools sort your_file.bam your_file_sorted

        Comment


        • #19
          just checking, thank you .

          bbmap is amazingly fast and accurate . Great tool.

          Comment


          • #20
            Originally posted by GenoMax View Post
            No. You have to do some work

            Code:
            $ samtools sort your_file.bam your_file_sorted
            Of course, I do my best to prevent needing people to do work, so...

            Code:
            bbmap.sh in=reads.fq out=mapped.bam bs=bs.sh; sh bs.sh
            That will write a script "bs.sh" that contains the command to make a sorted, indexed bam file, then execute it after BBMap is done. It saves... well, at least a couple of keystrokes

            Comment


            • #21
              Hmm. I should start a new thread for "Easter eggs" in BBMap

              Comment


              • #22
                Couldn't you create the shell script and then execute it within BBMap?

                Comment


                • #23
                  Originally posted by westerman View Post
                  Couldn't you create the shell script and then execute it within BBMap?
                  Unfortunately, not easily, due to memory constraints. BBMap's shell script already grabs all available memory by default. It can launch a samtools process for conversion to bam, which is a low-memory streaming process; but sorting uses a lot of memory (to do it most efficiently), so it might run out of memory. Java cannot reduce its virtual memory allocation after starting, and won't share that memory with other programs.

                  It would work fine on most individual computers since physical memory limits would never be exceeded, but job schedulers often have strict limits on virtual memory and will kill jobs that exceed it. So, it would reduce stability (and typically fail) on clusters unless I did parsing in the shellscript to launch a new process after BBMap terminates, or the user launched BBMap with a custom -Xmx flag. I try to keep the shellscripts as simple as possible. Anyway, there are workarounds, but they all add complexity due to memory issues.

                  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,”...
                    Yesterday, 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, Yesterday, 07:15 AM
                  0 responses
                  12 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 05-23-2024, 10:28 AM
                  0 responses
                  15 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 05-23-2024, 07:35 AM
                  0 responses
                  16 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 05-22-2024, 02:06 PM
                  0 responses
                  8 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X