Header Leaderboard Ad

Collapse

Aligner that reports all alignments?

Collapse

Announcement

Collapse

SEQanswers June Challenge Has Begun!

The competition has begun! We're giving away a $50 Amazon gift card to the member who answers the most questions on our site during the month. We want to encourage our community members to share their knowledge and help each other out by answering questions related to sequencing technologies, genomics, and bioinformatics. The competition is open to all members of the site, and the winner will be announced at the beginning of July. Best of luck!

For a list of the official rules, visit (https://www.seqanswers.com/forum/sit...wledge-and-win)
See more
See less
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Aligner that reports all alignments?

    Hi, I'm looking for an aligner that reports multiple/all hits rather than just a single one. I know mrfast or BWA-short can do this for Illumina data, but is there any aligner that can do this that works on non-Illumina reads? My searches haven't turned up anything.

  • #2
    What type of non-Illumina reads? There are now a few categories out there. mrFAST should work for Ion Torrent. If you are using pacbio, the blasr parameter -bestn limits the number of hits you will get, up to 2^32 - 1.

    Originally posted by Rocketknight View Post
    Hi, I'm looking for an aligner that reports multiple/all hits rather than just a single one. I know mrfast or BWA-short can do this for Illumina data, but is there any aligner that can do this that works on non-Illumina reads? My searches haven't turned up anything.

    Comment


    • #3
      Good point. I'm interested in 454/Ion Torrent or PacBio reads, so thanks for the info! I didn't realize mrfast could handle Ion Torrent reads, they're pretty specific about calling it an Illumina read aligner on their site. I'll give it a try though!

      Comment


      • #4
        Have you tried using other aligners with parameters optimized for returning all hits like bowtie2 (-a), bwa mem (-a), novoalign (-r All) and tmap (mapall)? Thanks

        Comment


        • #5
          See Heng Li's announcement about BWA-MEM. http://www.biostars.org/p/67236/

          Comment


          • #6
            Oh god, I'm sorry. I kept looking through that manual and never saw that. Sorry!

            Comment


            • #7
              Glad it helped :-)

              Comment

              Latest Articles

              Collapse

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, 06-01-2023, 08:56 PM
              0 responses
              11 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 06-01-2023, 07:33 AM
              0 responses
              58 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 05-31-2023, 07:50 AM
              0 responses
              98 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 05-26-2023, 09:22 AM
              0 responses
              102 views
              0 likes
              Last Post seqadmin  
              Working...
              X