Clustal format can take numeric names and correction of seq conservation

evolbioinfo · Jul 4, 2018 · c4790bd · c4790bd
1 parent 72d22cf
commit c4790bd
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/align/align.go b/align/align.go
@@ -1295,7 +1295,7 @@ func (a *align) SiteConservation(position int) (conservation int, err error) {
 				}
 			}
 		}
-		if prevchar != ';' && sequence[position] != prevchar {
+		if (prevchar != ';' && sequence[position] != prevchar) || sequence[position] == GAP {
 			same = false
 		}
 		prevchar = sequence[position]

diff --git a/io/clustal/parser.go b/io/clustal/parser.go
@@ -114,7 +114,7 @@ func (p *Parser) Parse() (align.Alignment, error) {
 			currentnbseqs = 0
 		}
 
-		if tok != IDENTIFIER {
+		if tok != IDENTIFIER && tok != NUMERIC {
 			return nil, errors.New("We expect a sequence identifier here")
 		}
 		name = lit