Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Alignment based on 1 mismatch.

    Hi
    I have a requirement where on aligning reads to reference the mapping is based only on one mismatch with the reference. I need a tool which can map the reads with the reference genome with max of one mismatch. Can some one suggest one.
    Thanks in advance.

    Regards
    Vishwesh

  • #2
    Bowtie2 will do this fine in end-to-end mode. Just set the minimum score to -1 (--score-min C,-1,0) and maximum/minimum mismatch penalties to both -1 (--mp -1,-1).

    Comment


    • #3
      Thanks David Eccles.. Is it possible to do this in bwa aln... if so what parameters should be given...

      Comment


      • #4
        From the bwa manual it looks like the edit distance option is what you want (-n 1), coupled with zero extensions allowed (-o 0). I have no experience with BWA, so can only do my best to interpret the specifications.

        Comment


        • #5
          Thanks David Eccles.. I am new to bowtie.. the options are bit confusing.. can you please help me the command.. The files I use are not paired ends..

          Comment


          • #6
            What have you tried? Give me an example, together with the error message, and we'll be better able to solve this problem.

            Comment


            • #7
              I tried this..
              bowtie -S all_other_rna -score-min C,-1,0 -mp -1,-1 /home/vishwesh/soft/TLL28_GTAGAG_L004_R1_001_c2.fastaq test.sam

              Comment


              • #8
                please use bowtie2 for doing this:



                Options for bowtie are quite different from bowtie2.

                Comment


                • #9
                  I tried using bowtie2..
                  here is the error i get..

                  Error: Encountered internal Bowtie 2 exception (#1)
                  Command: /data/apps/bowtie2-2.1.0/bin/bowtie2-align --wrapper basic-0 -S all_other_rna -score-min C,-1,0 -mp /home/vishwesh/soft/TLL28_GTAGAG_L004_R1_001_c2.fastaq test.sam
                  bowtie2-align exited with value 1

                  command is
                  bowtie2 -S all_other_rna -score-min C,-1,0 -mp -1,-1 /home/vishwesh/soft/TLL28_GTAGAG_L004_R1_001_c2.fastaq test.sam

                  let me know if anything is wrong in this command..

                  Comment


                  • #10
                    your bowtie2 command should look something like this:
                    Code:
                    bowtie2 [options]* -x <bt2-idx> {-1 <m1> -2 <m2> | -U <r>} -S [<hit>]
                    [taken verbatim from the manual]

                    It looks like you haven't specified an index file on the command line, or the index is attached to the wrong arguments. Perhaps something like this would work better (assuming all_other_rna is your bowtie2 index generated with bowtie2-build):

                    Code:
                    bowtie2 --score-min C,-1,0 --mp -1,-1 -x all_other_rna -S /home/vishwesh/soft/TLL28_GTAGAG_L004_R1_001_c2.fastaq > test.sam

                    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
                    195 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 05-23-2024, 10:28 AM
                    0 responses
                    217 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 05-23-2024, 07:35 AM
                    0 responses
                    219 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