Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Sequence Alignment

  1. Dec 31, 2004 #1
    I tryed to align two sequences
    and CATAGA
    using Needleman-Wunsch algorithm but it only gives me 2 alignment resutls

    while I thought it should be 4 results because the first C in the second sequence is also permutable, right ?

    So can you tell me how to take out the other two alignments ?

    Thank you very much inadvance
  2. jcsd
  3. Dec 31, 2004 #2


    User Avatar
    Staff Emeritus
    Science Advisor
    Gold Member

    You should be able to change some parameter in the software option. You can change the gap penalty. Some alignment software also have problem when it comes to the first base pair permutation.
  4. Jan 1, 2005 #3
    Thank iansmith, I solve that problem with a pen and a piece of paper only..
    But i am still really wanting to know how a software parameter change would be able to handle that problem, what does such a software do to get 4 alignments ?
    This is what I took from my biology class, talking with me about software thingies I think is like pouring water on a duck's head because I actually don't understand a thing at all.

    Thank you so very much in advance and a "Happy New Year" from me to you too...:wink:
  5. Jan 2, 2005 #4
    I never get any help from anyone except iansmith !!!!!
  6. Jan 2, 2005 #5
    During class, we were pointed to ClustalW (http://www.ebi.ac.uk/clustalw/) to do sequence alignments. Click the FAQ, and walla! It gives you help on exactly what all those insane options do...

    And don't dispair, bioinformatics is a tough subject - and gets worse. Once I waltzed into my uni's computer labs to check my email and saw my friends (who were doing biomedical science) crying during their bioinformatics prac :surprised
    Last edited: Jan 2, 2005
  7. Jan 2, 2005 #6


    User Avatar
    Staff Emeritus
    Science Advisor
    Gold Member

    For 2 sequence aligment you can use Blast 2 sequence and you easily change so parameter.
    http://www.ncbi.nlm.nih.gov/blast/bl2seq/wblast2.cgi [Broken]

    For Clustal, I prefer to use the version available at the Biology Workbeanch of UCSD

    You need to register and it is free. Also you get a bunch of other bioinformatic tools.

    I also took a grad level class. You need to understand the principle and theory behind bioinformatics. It make more sense once you understand the basics. My class was mostly about sequence alignment, BLAST and making threes.
    Last edited by a moderator: May 1, 2017
  8. Jan 3, 2005 #7
    Thank you jikx
    iansmith, I know that was my mistake, I am really sorry. Can you tell me how to do now ??? I don't know what i should do. I am really stuck!
    By the way, would you also tell me what lies behind bioinformatics ?
  9. Jan 6, 2005 #8
    Just paste this in the big white box (from the site i mentioned above):

    >FOSB_Wunsch1 Protein fosA

    >FOSB_Wunsch2 Protein fosB

    and you get an output like this:
    http://www.ebi.ac.uk/cgi-bin/clustalw/result?tool=clustalw&jobid=clustalw-20050106-11352860&poll=yes [Broken]

    which is an output before you muck around with all the settings like window length, gap extensions etc (the FAQ can help you there) .. however, these programs are more for sequences which are hundreds of of basepairs long!
    Last edited by a moderator: May 1, 2017
  10. Jan 6, 2005 #9


    User Avatar
    Staff Emeritus
    Science Advisor
    Gold Member

    Wunsch, do you have any theoritical part for class.

    If you don't:
    First, you need to understand the mathematical models that are behind bioinformatics. You should look the Wright-Fisher model, Hardy-Weinberg Equilibrium, and Neutral model.
    Second, you need to understand probability and how it applies to bioinformatics.
    Third, try to understand what is behind the scoring methods such as the blossom scoring matrix.
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook