Header Leaderboard Ad

Collapse

Find large indels from PacBio reads?

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

  • Find large indels from PacBio reads?

    I know Quiver and GATK have been bundled in SMRT pipeline to call variants (the GATK module only finds SNPs though). I'm wondering if there are other tool(s) that can be used for this purpose, especially to find large indels. Thanks.

  • #2
    I have an unsupported version of BLASR that does better alignment across large gaps. I have used this to simply count indel variants (really mobile element insertion) using PacBio reads. PM me if you want access to this; I've been hesitant to post it online since it is on a branch that is not updated with bug fixes, etc.. The unsupported version uses an affine alignment penalty. Without this, bases tend to "bleed" out into gaps and make it difficult to call exactly where a gap is.

    -mark

    Comment

    Latest Articles

    Collapse

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 07:14 AM
    0 responses
    6 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 06-06-2023, 01:08 PM
    0 responses
    6 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 06-01-2023, 08:56 PM
    0 responses
    159 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 06-01-2023, 07:33 AM
    0 responses
    295 views
    0 likes
    Last Post seqadmin  
    Working...
    X