Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Very Short Read aligner

    Hi All,

    I am in process of building my own short read aligner for a lab working with cancer-genome sequencing. And have following questions:

    1. As part for a our problem - we have with us read length ranging from 18-22 bp.

    I am aware of several aligners currently available, and I was wondering if they were any issues at all, working with read length of above range. From all the research papers of aligners that I have gone through; usually the range has varied from 30 bp and upwards.

    It would be really helpful if other members , who might have worked on similar read length, could advice me.

    2. Also, are there any species-specific statistical heuristics that are known when employing an approach for short read alignment. To be more specific for example lets say for non-mammalian sequence source , I would like to set different set of parameters when doing a short read alignment as compared to mammalian source.

    3. As I mentioned, I am in process of developing my own aligner, any piece of advise or suggestions would be very valuable from members who have embarked on the same . I am still in brain-storming phase, and trying the scope my problem range that this short read aligner would address. For now formally this aligner should be able to:
    a. align read's length ranging 18-22 bp to reference genome
    b. ungapped alignment
    c. Applies a species - specific statistical scoring heuristic (If it all it makes sense to use one in first place.)

    I am looking forward to any advice or suggestion or perhaps even a general discussion on process of developing a short read aligner from scratch.

    thank you

    regards
    Rupinder

  • #2
    Rupinder,

    Originally posted by Rupinder View Post
    I am looking forward to any advice or suggestion or perhaps even a general discussion on process of developing a short read aligner from scratch.
    My advice is not to write (yet another) short read aligner, unless you are doing it to "learn" or as part of a student project.

    There are so many already: MAQ, Shrimp, SOAP, Bowtie, ELAND, Novocraft etc. They are all doing essentially the same thing, and most of them are very efficient indeed, and you are unlikely to better them. I'm sure you could find the appropriate parameters to suit your data, and 20-24 bp reads are not a problem. The reason most of them are 30+ bp is that most Illumina data is around that length, but SOLID3 data is often 24 bp and the aligners work well. Shrimp and others have post-processing scripts to help correct for bias with sequences with non-uniform statistics.

    I think you would be better off using the existing tools with appropriate settings, rather than write your own, and get on with the science further down the line.

    (But if you have the time, and want to learn more about alignment, optimization, programming etc, go for it! Especially if you want to write a GPU enabled version)

    --Torsten

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Strategies for Sequencing Challenging Samples
      by seqadmin


      Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
      03-22-2024, 06:39 AM
    • seqadmin
      Techniques and Challenges in Conservation Genomics
      by seqadmin



      The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

      Avian Conservation
      Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
      03-08-2024, 10:41 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 06:37 PM
    0 responses
    8 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 06:07 PM
    0 responses
    8 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-22-2024, 10:03 AM
    0 responses
    49 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-21-2024, 07:32 AM
    0 responses
    66 views
    0 likes
    Last Post seqadmin  
    Working...
    X