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
      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
    • seqadmin
      The Impact of AI in Genomic Medicine
      by seqadmin



      Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
      02-26-2024, 02:07 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 03-14-2024, 06:13 AM
    0 responses
    33 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-08-2024, 08:03 AM
    0 responses
    72 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-07-2024, 08:13 AM
    0 responses
    81 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-06-2024, 09:51 AM
    0 responses
    68 views
    0 likes
    Last Post seqadmin  
    Working...
    X