seqret reading from stdin with unkown type

David Mathog mathog at mendel.bio.caltech.edu
Fri May 10 21:35:45 UTC 2002


I don't understand why this one fails:

# ( echo "AGTTT" ; echo "ACGTTTTTGGT" ) \
  | seqret -sequence=stdin  -out=gcg::stdout
Reads and writes (returns) sequences
Error: Unable to read sequence 'stdin'

while these work:


# ( echo "AGTTT" ; echo "ACGTTTTTGGT" ) \
  | seqret -sequence=plain::stdin  -out=gcg::stdout
Reads and writes (returns) sequences
!!NA_SEQUENCE 1.0

  Length: 16  Type: N  Check:  665 ..

  1 AGTTTACGTT TTTGGT

bash-2.03# ( echo ">FOO" ; echo "ACGTTTTTGGT" ) \
  | seqret -sequence=stdin  -out=gcg::stdout
Reads and writes (returns) sequences
!!NA_SEQUENCE 1.0

FOO  Length: 11  Type: N  Check: 5205 ..

  1 ACGTTTTTGG T

According the seqret blurb the program is supposed to try
one format at a time until it finds something that works.  In
the first case it seems not to be trying "plain".  This was
for use with a CGI script, where it's hard to control
what the user enters into the box/uploads.

EMBOSS-2.2.0 on Solaris 5.8.

( 2.3.1 is current but I'm waiting for 2.4 before upgrading).

Thanks,

David Mathog
mathog at caltech.edu
Manager, Sequence Analysis Facility, Biology Division, Caltech




More information about the EMBOSS mailing list