Griffitts:Phylogenetics: Difference between revisions

From OpenWetWare
Jump to navigationJump to search
(Started page)
 
mNo edit summary
Line 4: Line 4:


==Alignment==
==Alignment==
To start off, you need to perform an alignment of all your sequences. This can be done using MUSCLE. The output file (.txt) then needs to be converted to a .nex file, which can be done using ALTER.
===MUSCLE===
===MUSCLE===
# Go to the [http://www.ebi.ac.uk/Tools/services/web/toolform.ebi?tool=muscle MUSCLE website]
# Go to the [http://www.ebi.ac.uk/Tools/services/web/toolform.ebi?tool=muscle MUSCLE website]
Line 28: Line 30:
# Under '''4. Save converted MSA''' hit '''save…''' and then '''Okay'''
# Under '''4. Save converted MSA''' hit '''save…''' and then '''Okay'''
# Open the .nex file with a text editor
# Open the .nex file with a text editor
# On the sixth line down delete where it says <blockquote><code>match=.</code></blockquote>
# On the sixth line down change <blockquote><code>interleave datatype=NUCLEOTIDE gap=- match=.;</code></blockquote> to <blockquote><code>interleave datatype=NUCLEOTIDE gap=-;</code></blockquote>
# Scroll to the bottom and paste the following<sup>[1]</sup> after <code>end;</code> :<blockquote><code>Begin paup;<br>set autoclose=yes;<br>set criterion=parsimony;<br>set root=outgroup;<br>set storebrlens=yes;<br>hsearch addseq=random nreps=1000 swap=tbr hold=1;<br>showtree;<br>End;</code></blockquote>
# Scroll to the bottom and change <blockquote><code>;<br>end;</code></blockquote>to the following<sup>[1]</sup>:<blockquote><code>;<br>end;<br><br>Begin paup;<br>set autoclose=yes;<br>set criterion=parsimony;<br>set root=outgroup;<br>set storebrlens=yes;<br>hsearch addseq=random nreps=1000 swap=tbr hold=1;<br>showtree;<br>End;</code></blockquote>
# Save<br><br>
# Save<br><br>
NOTE: If you're doing this with multiple alignments, you'll have to refresh the ALTER page before working on a new .txt file.
NOTE: If you're doing this with multiple alignments, you'll have to refresh the ALTER page before working on a new .txt file.

Revision as of 13:25, 16 April 2011

Alignment

To start off, you need to perform an alignment of all your sequences. This can be done using MUSCLE. The output file (.txt) then needs to be converted to a .nex file, which can be done using ALTER.

MUSCLE

  1. Go to the MUSCLE website
  2. Paste your sequences in FASTA format with two returns between each one
    • To generate FASTA format from Strider, select File and then Export as…; when the new window pops up, set the Format bar to FASTA
  3. Set the output to "Pearson/FASTA"
  4. Do NOT change any of the other parameters
  5. Hit Submit
  6. When you have the alignment, click on the Result Summary tab
  7. Right click on the link beneath Alignments in FASTA format and select Save link as…
  8. Give the file the name that you want and then add the extension .txt

ALTER

  1. Go to the ALTER website
  2. Under 1. Select format
    • Choose "MUSCLE" for program
    • Choose "FASTA" for format
  3. Under 2. Paste or Upload MSA hit upload…
  4. Navigate to the .txt file that you created using MUSCLE, select it, and hit open
  5. Under 3. Select output format and convert
    • Choose "PAUP" for program
    • Choose "NEXUS" for format
  6. Hit Convert!
  7. Under 4. Save converted MSA hit save… and then Okay
  8. Open the .nex file with a text editor
  9. On the sixth line down change

    interleave datatype=NUCLEOTIDE gap=- match=.;

    to

    interleave datatype=NUCLEOTIDE gap=-;

  10. Scroll to the bottom and change

    ;
    end;

    to the following[1]:

    ;
    end;

    Begin paup;
    set autoclose=yes;
    set criterion=parsimony;
    set root=outgroup;
    set storebrlens=yes;
    hsearch addseq=random nreps=1000 swap=tbr hold=1;
    showtree;
    End;

  11. Save

NOTE: If you're doing this with multiple alignments, you'll have to refresh the ALTER page before working on a new .txt file.

Trees

Notes

[1] See https://lswiki.byu.edu/index.php/Life_Sciences_Wiki:Community_Portal.