EddyRivasLab / easel

Sequence analysis library used by Eddy/Rivas lab code
Other
46 stars 26 forks source link

esl-seqstat can't read alignment file from stdin #8

Closed traviswheeler closed 7 years ago

traviswheeler commented 7 years ago

ID e2 TITLE esl-seqstat can't read alignment file from stdin AFFECTS - FIXED_IN - STATUS CLOSED XREF J4/84 REPORTED_BY SRE OPENED_DATE 16 Feb 2009 CLOSED_DATE 16 Feb 2009 DESCRIPTION
sqfile_open tries to guess file format; then when msafile_open is called in GuessAlphabet, stdin is already feof(). Using current design, MSA can't be read sequentially (through sqio interface) unless format is known. SQFILE has a recording buffer mechanism, but it only works in sqio proper, not through interface to esl_msa. Added e2.sh to test for the bug.