[EMBOSS] seqret again

Derek Gatherer d.gatherer at vir.gla.ac.uk
Wed Jun 15 13:25:36 UTC 2005


Is this a bug?  Compare the following output from seqret when phylip and 
phylip3 are specified.  Shouldn't the first line of the phylip3 output be 
"2 546 YF" and not "1 546" ?

[gath01d at gamma EBV]$ seqret -osformat phylip seqs/balf1.both
Reads and writes (returns) sequences
Output sequence [c-balf1.phylip]: seqs/balf1.phylip
[gath01d at gamma EBV]$ more seqs/balf1.phylip
  2 546
c-balf1.seATGCAGCCAG CCAAGTCTAC CGATTCGGTG TTTGTGAGGA CCCCGGTCGA
balf1.seq ATGAGGCCAG CCAAGTCTAC AGATTCTGTG TTTGTGAGGA CCCCGGTCGA

           GGCGTGGGTC TCACCCTCGC CCCCGGACGA CAAAGTGGCA GAGACCAGCT
           GGCGTGGGTC GCGCCCTCGC CGCCGGACGA CAAGGTGGCT GAGTCCAGCT
[snip]

[gath01d at gamma EBV]$ seqret -osformat phylip3 seqs/balf1.both
Reads and writes (returns) sequences
Output sequence [c-balf1.phylip3]: seqs/balf1.phylip3
[gath01d at gamma EBV]$ more seqs/balf1.phylip3
1 546 YF
c-balf1.seATGCAGCCAG CCAAGTCTAC CGATTCGGTG TTTGTGAGGA CCCCGGTCGA
           GGCGTGGGTC TCACCCTCGC CCCCGGACGA CAAAGTGGCA GAGACCAGCT
           ACCTCCTGTT CAGGGCCCTA TACGCTGTGT TCACCCAGGA CGAGACGGAC
           CTGCCTCTAC CGGCCCTGGT CATGTGCCGG CTCCTGAAGG CCTCCCTGAG

[snip]





More information about the EMBOSS mailing list