Multiple sequence alignement MSA editing
I want to remove all rare insertions (when it occurs in less than 5% of the sequences) in a multiple sequence alignment file (clustal .aln) with 699 sequences.
That is, I have a MSA with many columns containing only one or two insertions while the rest of the sequences are blank "-". It is by far too much to do manually.
Any suggestions how to do this?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Remove short sequences with length < 50% from a multi-fasta file
written by b.shrestha •
-
Clustal Omega alignment algorithm - does it account for codons and open reading frames (ORF)?
written by manaswwm
-
How to visualize multiple sequence alignment and read coverage simultaneously
written by callamartyn •
-
Generate nucleotide multiple sequence alignment from protein alignment?
written by molly77 •
-
ClustalW alignment - Unexpected result
written by pentium3-user •
-
Generate HMM profile without an MSA (hmmer)
written by elenatea.russo •
-
Combining different tasks in one line of code in terminal
written by Za
-
How to align new DNA sequences with existing multiple DNA sequence alignment
written by liaoyunshi •
-
Automated sliding window analysis of % identity on a Clustal Omega DNA alignment
written by c.merrick •
-
How can I pick out prokaryotic sequences from a multiple sequence alignment file?
written by Chong Tang
Do you specifically want to do this in R?
If you use BioPython, you can create an ungapped concensus sequence with a threshold for inclusion of a particular residue in a column.