Announcement

Collapse

Welcome to the New Seqanswers!

Welcome to the new Seqanswers! We'd love your feedback, please post any you have to this topic: New Seqanswers Feedback.
See more
See less

Stampy: A statistical algorithm for sensitive and fast mapping of Illumina sequence r

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Stampy: A statistical algorithm for sensitive and fast mapping of Illumina sequence r

    Abstract

    High-volume sequencing of DNA and RNA is now within reach of any research laboratory and is quickly becoming established as a key research tool. In many workflows, each of the short sequences (“reads”) resulting from a sequencing run are first “mapped” (aligned) to a reference sequence to infer the read from which the genomic location derived, a challenging task because of the high data volumes and often large genomes. Existing read mapping software excel in either speed (e.g., BWA, Bowtie, ELAND) or sensitivity (e.g., Novoalign), but not in both. In addition, performance often deteriorates in the presence of sequence variation, particularly so for short insertions and deletions (indels). Here, we present a read mapper, Stampy, which uses a hybrid mapping algorithm and a detailed statistical model to achieve both speed and sensitivity, particularly when reads include sequence variation. This results in a higher useable sequence yield and improved accuracy compared to that of existing software.

    http://www.ncbi.nlm.nih.gov/pubmed/20980556
    Last edited by lh3; 10-29-2010, 06:56 AM.

  • #2
    It looks like an interesting aligner and it showed very good performance at indels. I think this is largely due to it's low gap extend penalty, -40 for open, -3 to extend 1 base. Most other aligners would have a higher gap extension penalty relative to gap open and mismatch penalties.
    Originally posted by lh3 View Post
    Abstract

    High-volume sequencing of DNA and RNA is now within reach of any research laboratory and is quickly becoming established as a key research tool. In many workflows, each of the short sequences (“reads”) resulting from a sequencing run are first “mapped” (aligned) to a reference sequence to infer the read from which the genomic location derived, a challenging task because of the high data volumes and often large genomes. Existing read mapping software excel in either speed (e.g., BWA, Bowtie, ELAND) or sensitivity (e.g., Novoalign), but not in both. In addition, performance often deteriorates in the presence of sequence variation, particularly so for short insertions and deletions (indels). Here, we present a read mapper, Stampy, which uses a hybrid mapping algorithm and a detailed statistical model to achieve both speed and sensitivity, particularly when reads include sequence variation. This results in a higher useable sequence yield and improved accuracy compared to that of existing software.

    http://www.ncbi.nlm.nih.gov/pubmed/20980556

    Comment


    • #3
      You are right that we should apply small gap extension penalty. In addition, I believe stampy is the major competitor to novoalign.

      Comment


      • #4
        Agree, looking at 1000 Genomes indels stats a gap extension penalty of 3 looks about correct. Novoalign still showed better sensitivity on SNPs and small indels and lowering gap extend penalty should fix the longer indels but Stampy is getting close.

        Comment

        Working...
        X