[BioPython] Performing sequence alignments, etc.
Caitlin
The_Polymorph at rocketmail.com
Mon Oct 15 02:13:24 UTC 2007
Thanks Peter. The genomes are viral. I'll definitely read that
tutorial.
Your help is very appreciated.
~Caitlin
--- Peter <biopython at maubp.freeserve.co.uk> wrote:
> Caitlin wrote:
> > Hi all.
> >
> > I'm relatively new to the field of bioinformatics and I'm trying to
> > perform a multiple sequence alignment on 5-6 sequences (fasta
> format -
> > dna sequences). I'd like the output to be formatted in the
> following
> > manner (clustalw standalone output):
>
> For reading and writing Clustalw alignment files, you could either
> use
> Bio.SeqIO (format name "clustal") or the Bio.Clustalw module.
> http://biopython.org/wiki/SeqIO
>
> > When one more more nucleotides columns are identical, clustalw
> displays
> > an asterisk. If not, a blank space is displayed. Is this a standard
> > feature of BioPython?
>
> There is an example of Clustalw output online here - note there can
> also
> be a column of numbers on the right hand side (not shown here):
> http://www.bioperl.org/wiki/ClustalW_multiple_alignment_format
>
> It sounds like you are describing the simple consensus string which
> clustalw outputs under the alignment (using *:. and space).
>
> Biopython has a SummaryInfo object which can calculate simple
> consensus
> sequences (see the tutorial). Perhaps this would be close to what you
>
> want to do.
>
> > Also, I'm evaluating several sequences but I'd like to obtain the
> most
> > recent complete genomes possible from various countries. Is there a
> > convenient source to use (GenBank?) if I don't know the accession
> > numbers?
>
> What sort of Genomes? Bacteria? Vertebrates? You could start by
> having
> a look at any of the EMBL, NCBI/GenBank or the Japanese DDBJ (these
> three are kept in sync with each other).
>
> Biopython has quite a nice interface for searching and downloading
> sequences from GenBank (again, see the tutorial) so that would be my
> first suggestion.
>
> Peter
>
>
>
>
"Be who you are and say what you feel because those who mind don't
matter and those who matter don't mind."
- Dr. Seuss, "Oh the Places You'll Go"
____________________________________________________________________________________
Don't let your dream ride pass you by. Make it a reality with Yahoo! Autos.
http://autos.yahoo.com/index.html
More information about the Biopython
mailing list