Multiple Sequence Alignment Python
Hi,
Currently I am working on classifying coding/non coding RNAs based on sequence features. I want to try MSA score as a feature for sequence classification.
- Please help me to get MSA score (using Clustlw) for a collection of sequences, using Biopython.
- How can I a get a profile sequence after MSA of training sequences?
- Also suggest other features that can be derived from MSA
• 3,475 views
•
link
updated
by
Ram
•
written
by
Manu Madhavan •
•
written
by
Manu Madhavan •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
long non coding RNAs annotation file
written by annaA •
-
Tool: DASHR 2.0: integrated database of human small non-coding RNA genes and mature products
written by pkouxa •
-
RNA secondary structure - Feature extraction
written by Manu Madhavan •
-
Using BLAST to calculate all-pair sequence similarity score
written by convalescence21 •
-
Get pairwise scores from alignment without re-aligning?
written by mpd •
-
Estimating the distribution of small non-coding RNAs among the reads: Which annotation database / r…
written by f.fassihian •
-
Ideal sequence % identity for profile construction
written by Anand Rao
-
gDNA annotation, exon and intron boundary identification
written by chethusomashekar •
-
Detecting Derived Sites
written by ngsgene
-
Novel Strategy For Improving Sequence Alignment
written by Pals
What have you tried so far?
I have used k-mer (k-length substring) frequencies, GC content, Molecular weight ...,. I would like to include features from MSA (some variance from profile sequence) to this feature set.
There is a Clustal module for python I believe, though I've never personally used it. There should be a way to get the scores back out of that (though I can't help you with the specifics). In the past, I've simply parsed the STOUT output of running a commandline invocation of Clustal. It prints the pairwise alignment scores for all the sequences so you could run with that.
For building a profile sequence, you can pass your MSA through
hmmerto get an HMM of the sequence.As for other features, I can't think of anything specifically that would be valid for all sequences beyond their pairwise agreements. Others may know more...