obitools/obitools4
2
0
Fork 0
mirror of https://github.com/metabarcoding/obitools4.git synced 2025-06-29 16:20:46 +00:00
Code Issues Packages Projects Releases Wiki Activity

Prevens errors on FASTA formating an empty sequence

Browse Source
This commit is contained in:
Eric Coissac
2022-11-16 09:22:55 +01:00
parent 493c7fddef
commit 7f863fb383
1 changed files with 12 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

18
pkg/obiformats/fastseq_write_fasta.go
View File

@ -30,15 +30,21 @@ func FormatFasta(seq *obiseq.BioSequence, formater FormatHeader) string {
s := seq.Sequence() s := seq.Sequence()
l := len(s) l := len(s)
fragments.Grow(l + int(l/60) + 10) folded := ""
if l == 0 {
log.Println("Writing a BioSequence of length zero")
} else {
fragments.Grow(l + int(l/60)*2 + 100)
for i := 0; i < l; i += 60 { for i := 0; i < l; i += 60 {
to := min(i+60, l) to := min(i+60, l)
fmt.Fprintf(&fragments, "%s\n", string(s[i:to])) fmt.Fprintf(&fragments, "%s\n", string(s[i:to]))
}
folded = fragments.String()
folded = folded[:fragments.Len()-1]
} }
folded := fragments.String()
folded = folded[:fragments.Len()-1]
info := formater(seq) info := formater(seq)
return fmt.Sprintf(">%s %s %s\n%s", return fmt.Sprintf(">%s %s %s\n%s",
seq.Id(), info, seq.Id(), info,