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