SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Pacific Biosciences



Similar Threads
Thread Thread Starter Forum Replies Last Post
PacBio Amplicon reads assembly PatrickV Pacific Biosciences 4 07-01-2015 09:30 AM
Filtering PacBio reads ebioman Pacific Biosciences 30 04-20-2015 05:26 PM
Alignment produces BAM file with sorted reads, but cannot see alignment alch Bioinformatics 11 04-14-2015 06:59 PM
BWA_MEM on PacBio Reads amaurizio Bioinformatics 1 09-08-2014 03:06 AM
Would PacBio RII reads be useful at 1X? jwag Pacific Biosciences 7 04-21-2014 07:11 AM

Reply
 
Thread Tools
Old 05-18-2016, 08:42 PM   #1
Rui Guo
Member
 
Location: China

Join Date: Apr 2016
Posts: 18
Default pacbio reads alignment

Which software is proper for pacbio reads alignment to reference genome? I see bwa mem is ok, and there is blasr, but slow. How about MHAP? But from the manual it says the output is blasr m4 format, not sam/bam. Why that? The MHAP paper show better alignment performance to blasr, but it cannot generate sam?
Rui Guo is offline   Reply With Quote
Old 05-19-2016, 09:21 AM   #2
gconcepcion
Member
 
Location: Menlo Park

Join Date: Dec 2010
Posts: 68
Default

Quote:
Originally Posted by Rui Guo View Post
Which software is proper for pacbio reads alignment to reference genome? I see bwa mem is ok, and there is blasr, but slow. How about MHAP? But from the manual it says the output is blasr m4 format, not sam/bam. Why that? The MHAP paper show better alignment performance to blasr, but it cannot generate sam?
Blasr is the recommended option, but as you noted, can be slow. Take a look at some of these tips to speed it up a little bit for nearly identical sequences:
https://github.com/PacificBioscience...ical-sequences

mhap is another option, but as you noted, output is not in SAM format.

The third option would be 'daligner' which is pretty fast, but also uses it's own output format (*.las) which is not what you asked for, but one of our preferred methods for aligning long reads back to a reference. See Gene Meyer's Dazzlerblog for details:
https://dazzlerblog.wordpress.com/about/
gconcepcion is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 05:42 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO