Making All Protein Sequence Lengths Same
Is there any code in perl / python to make all protein sequences of same length, otherwise my phylogenetic tool MEGA is not working on them ?
• 4,862 views
•
link
updated
by
Jan Kosinski
•
written
by
Ashru2006 •
•
written
by
Ashru2006 •
1 answer
If you don't mind using interactive programs instead of Python/Perl, you can try with Jalview (cross-platform) or BioEdit (Windows)
Jalview: Menu edit-> Pad Gaps
BioEdit: Menu Alignment -> Flush the alignment
But of course, as others suggested in the comments, ensure first you have your sequences aligned, and different lengths arose for different reasons.
• 1 views
•
link
written
by
Jan Kosinski
Log in to answer this question.
More posts like this
-
Phylogenetic tree for 16s sequences in a metagenomic data set
written by surathsanda •
-
a hash with a pair of sequences as a value
written by natasha.sernova
-
MEGA: Infer ancestral states
written by cv.napping •
-
how to find the ancestor sequence of a molecular phylogenetic tree?
written by crystalrose •
-
how to find conserved region within a protein family for the construction of phylogenetic tree
written by star_143_101 •
-
Distance matrix with different length sequences - novice question
written by mootzk •
-
Code For Generating 'Consensus' Sequence From Pairwise Alignment
written by Woa
-
Building A Ultrametric Phylogenetic Tree With Integer Branch Length
written by Venky •
-
Phylogenetic Tree Of Fragments Of The Same Protein (From A Metagenome)
written by Gimly_Gloin •
-
How Can I Get Protein Sequence In Fasta Format Using Taxon Id?
written by Luke
Have you aligned your sequences? You need to in order to do phylogeny and almost any aligner will pad out sequences to make them the same length.
I would guess that they should be globally-aligned and therefore of the same length. Check your alignment?
David W, I've heard that there are programs that construct alignment and phylogeny in the same time. But I think they are not implemented in MEGA
It happened to me several times to have alignments with sequences of different lengths. These can arise in several situations when you do more sophisticated things, e.g. manual edits, combining several alignments into one using some sequence alignment editors etc. The fact that most of the programs and programming libraries complain about unequal sequence lengths, instead of correcting that automatically, is for me one of the biggest annoyances in bioinformatics.
I liked the way Jan described "manual editing" as one of the more sophisticated things you can do;-)