sequence similarity in R
I use emboss water on the command line to do pairwise sequence alignment. It worked fine for a few sequences, but I now have a list of 100s of sequences to check pairwise alignment for. Is there a way to automate this in R ?
protein-alignment
homology
protein
similarity
sequence
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Pairwise alignment and Ka/Ks computation
written by maxime.policarpoDear all, I have a multi-fasta file which consists in 1,000 coding sequences. I would like to compute the Ka/Ks between each pair of sequences. …
-
Annotating genome based on Sequence
written by buhbs •Hi all, I was wondering if anyone knows of an R package to annotate genomes based on the sequence of features. For example, I would …
-
how to visualizing residue wise pairwise protein alignment
written by kng •How can I visualize the similarity percentage per residue for results from pairwise sequence alignment (results from emboss water). I would like each residue of …
-
EMBOSS program for calculating dN/dS from MSA
written by shpak.max •Does EMBOSS have an associated program that will calculate dN, dS (non-synonymous/synonymous substitutions) for a multiple sequence alignment? If so, what is the program, and …
-
How to obtain PEG IDs from a list of FASTA Sequences
written by tastafor •Hi, i have a list of FASTA sequences. I have been going to the SEED Viewer (https://pubseed.theseed.org/seedviewer.cgi?page=BlastRun) and running a BLAST of the FASTA sequence …
-
sequence alignment of protein sequence
written by fatimarasool135I want to find the protein sequences of wheat which have IGT motif. I used the scanprosite tool to find the protein sequences which have …
-
Multiple Sequence Alignment
written by Natasha •Hello Everyone, I am performing multiple sequence alignment using Clustal Omega. For pairwise alignment, one can obtain the percentage similarity and percentage identity scores. When …
-
Getting Pairwise Sequence alignment From The Muscle Command-Line
written by bioinformaticsm28 •The **MUSCLE** command line doesn't have an option for returning the **pairwise sequence alignment**. Is there a way to get those alignment(in python)? for example …
-
Aligning several sequences with EMBOSS water program
written by dolevrahat •Hello I am trying to use EMBOSS water program to align several sequences stored in one fasta file to a sequence in another file. I'm …
-
Retrieving Best N Local Alignments
written by Guidobot •<p>Can anyone point me to a program [C/Python] that runs a Smith-Waterman pairwise local alignment but (unlike EMBOSS/water) has options for one or both of …
You can shell-out from R as you can with most languages, but I'd question using R for this really. If you've already done alignments on the commandline, stay on the commandline.
All you need is a shell loop or some
GNU parallelusage to achieve this. It can even be done as a one-line shell command.