How can I pick out prokaryotic sequences from a multiple sequence alignment file?
I have a huge MSA file containing over 1,000 sequences which is in a typical clustal format. If I want to select those prokaryotes sequences from this msa file and create a new msa file. What should I do?
•
written
by
Chong Tang
1 answer
A non sophisticated way of doing it, would be to cut and paste all the accession numbers into the PIR database http://pir.georgetown.edu/pirwww/index.shtml
This will bring up all the details for all the sequences.
Select from drop down box filters, and select TAXA ID. In this field simply enter 2 (for Prokaryotes) and the results will be filtered. Then click select all and then there is an alignment button (in Jalview) on top right of results.
Depending on your alignment naming/ID DAVID might be handy for converting the initial sequences into a format PIR can accept.
Log in to answer this question.
-
Loading multiple sequence alignment fasta file for ggmsa
written by snowpin •
-
Multiple alignment filtering
written by guillaume.rbt
-
Clustal Omega and Biopython
written by oonatoyou •
-
New element for Sequence Alignment
written by virus_n00b •
-
Multiple sequence alignement MSA editing
written by TEman •
-
Get pairwise scores from alignment without re-aligning?
written by mpd •
-
JavaScript access to Clustal Omega REST service
written by zdebruine •
-
Is FastML a good server for reconstructing ancestral protein sequence?
written by lia.elias •
-
Downloading multiple sequence alignment as clustal format file from clustal omega.
written by siddharth.avadhanam •
-
Select Sequence From A Psiblast, And Launch A New Iteration From It
written by mistervanhalen •
What are the names of the sequences? TXIDs? Organism name? The first step would be to get a list of the desired sequences' names.