Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • trouble with bwakit - where are the options?

    Hi all,

    I'm wondering if anyone has some experience with bwakit v0.7.12 and can chime in.

    The command run-bwamem within bwakit seems to have many fewer options than previous stand alone versions of bwa. Most notably, I would like to pass the option -M to bwamem, which marks shorter split hits as secondary. However, when I execute run-bwamem with this option, it throws an error:

    Code:
    /r/bwa-kit/run-bwamem -o test -t8 -MHR"@RG\tID:run1\tSM:${sample}\tLB:LaneX\tPL:Illumina\tPU:none" /r/bwa-kit/hs38DH.fa LP6005206-DNA_A06.1.fastq LP6005206-DNA_A06.2.fastq | sh
    Stout: Unknown option: M

    Is it the case that this option simply doesn't exist for the version of BWA mem that bwakit calls?


  • #2
    I've used bwa mem but not via that kit.

    Does running /r/bwa-kit/run-bwamem give you the available options?

    Comment


    • #3
      Yes, it does list some options, but they are far more limited than what one sees running bwa mem. I suppose that's a big clue, but I assumed those options existed on top of what you can already pass to bwa mem. Seems like that's not the case, though. A big minus for bwa-kit.

      Can anyone else confirm or deny?

      Thanks all!

      Comment


      • #4
        If you look at the run-bwamem script you can see the options it accepts (which is a sub-set of bwa mem): https://github.com/lh3/bwa/blob/mast...kit/run-bwamem

        Looks like Heng Li developed this as an end-to-end solution for HLA typing analysis (https://github.com/lh3/bwa/blob/master/bwakit/README.md).
        Last edited by GenoMax; 11-13-2015, 10:19 AM. Reason: Corrected URL

        Comment


        • #5
          @Genomax is correct.
          bwakit is nothing more than a collection of Perl scripts written to process data with bwa in the context of a specific project.
          It seems simpler and more direct to use bwa mem directly, without using the bwakit scripts, especially if your project is different from theirs.
          Nevertheless, the Perl scripts can easily be modified, if you want to continue working with bwakit for your own project.

          I've forked the bwakit GitHub repository, and added the -M option.
          GitHub is where people build software. More than 100 million people use GitHub to discover, fork, and contribute to over 420 million projects.


          The usual caveat applies that my modifications may always introduce a new bug, but that any bug detected can easily be fixed.

          Comment


          • #6
            Thanks blancha, Genomax, for the clarification. Sounds like the simplest thing for us would be to use the stand alone BWA software for the time being.

            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, 11-01-2024, 06:09 AM
            0 responses
            18 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 10-30-2024, 05:31 AM
            0 responses
            18 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 10-24-2024, 06:58 AM
            0 responses
            24 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 10-23-2024, 08:43 AM
            0 responses
            53 views
            0 likes
            Last Post seqadmin  
            Working...
            X