Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • extended cigar in bwa-mem

    Hi all,

    I am aligning pacbio reads to a reference genome using bwa mem. Do you know how to output an extended CIGAR in output sam file from bwa v0.7.13?

    My command is:
    Code:
    bwa mem -x pacbio refgenome.fasta reads.fastq > output.sam
    Thanks in advance

  • #2
    What exactly do you mean by "extended CIGAR"? bwa mem will output SAM files with an appropriate CIGAR string already.

    Comment


    • #3
      He probably means cigar strings with X and = symbols instead of M, which are very handy.

      Comment


      • #4
        Yes, this is what I mean by "extended cigar". By the way, I found a tool that can do that a posteriori, in case someone else is interested: SamFixCigar (http://github.com/lindenb/jvarkit )

        Comment


        • #5
          Are there any plans to extend bwa to write X/= instead of M? At least optional?

          Comment


          • #6
            You can now use BBMap's reformat.sh to add those:

            reformat.sh in=mapped.sam out=extended.sam sam=1.4

            I have not tested it extensively on things like hard-clipping but it should generally work.

            Comment


            • #7
              Brian,

              thanks for the suggestion. I will try :-)

              Nevertheless I'd wish bwa would provide an option to either use M or X=.
              Just to avoid more I/O due to another conversion ...

              Sven

              Comment


              • #8
                After studying some docs I do see that with optional MD:Z I am even as flexible as with X=, with the advantage of not having too "complicated" CIGARs. :-)

                Comment

                Latest Articles

                Collapse

                • 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
                • seqadmin
                  The Impact of AI in Genomic Medicine
                  by seqadmin



                  Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
                  02-26-2024, 02:07 PM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, 03-14-2024, 06:13 AM
                0 responses
                33 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-08-2024, 08:03 AM
                0 responses
                72 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-07-2024, 08:13 AM
                0 responses
                81 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-06-2024, 09:51 AM
                0 responses
                68 views
                0 likes
                Last Post seqadmin  
                Working...
                X