[Bioperl-l] How to write a Bio::DB:Fasta sequence into a file?

Mark A. Jensen maj at fortinbras.us
Sat Nov 29 14:34:33 UTC 2008


## sorry about the multiple posts

What you are doing is a little unclear. If you have obtained a
sequence object, from Bio::DB::Fasta or anywhere, you can write
it to a FASTA file by using the write_sequence call. Maybe you
have done

$seq = read_sequence("theirseqs.fas", 'fasta');

or

($seq, @more_seqs) = read_all_sequences("theirseqs.fas", 'fasta');

or even

$db = new Bio::DB::Fasta("theirseqs.fas");
$seq = $db->get_Seq_by_id('THX1138');

then write it with

write_sequence(">myseqs.fas", 'fasta', $seq);

or append it to your file with

write_sequence(">>myseqs.fas", 'fasta', $seq);

If your problem is, you have a sequence string, and you want to
write it to a file in FASTA format, then you create a sequence
object directly:

$seq = new_sequence( 'atcgtgcaat', 'THX1138' )

and write it using write_sequence() as above. 

Hope this helps-
Mark




More information about the Bioperl-l mailing list