From e9cdfd7e03135aac298102da711b691c625ce9eb Mon Sep 17 00:00:00 2001
From: Eric Coissac <eric.coissac@metabarcoding.org>
Date: Tue, 1 Feb 2022 18:49:32 +0100
Subject: [PATCH] Make subseq method dealing with qualities

---
 .gitignore           | 10 ++++++----
 pkg/obiseq/subseq.go |  9 +++++++++
 2 files changed, 15 insertions(+), 4 deletions(-)

diff --git a/.gitignore b/.gitignore
index 2404b8a..9f0a757 100644
--- a/.gitignore
+++ b/.gitignore
@@ -1,11 +1,13 @@
 cpu.pprof
 cpu.trace
-obicount
-obipairing
 test
 bin
 vendor
 *.fastq
 *.fasta
-obiconvert
-obipcr
+
+/obiconvert
+/obicount
+/obimultiplex
+/obipairing
+/obipcr
diff --git a/pkg/obiseq/subseq.go b/pkg/obiseq/subseq.go
index 9571b1d..bdf353d 100644
--- a/pkg/obiseq/subseq.go
+++ b/pkg/obiseq/subseq.go
@@ -27,11 +27,20 @@ func (sequence BioSequence) Subsequence(from, to int, circular bool) (BioSequenc
 		newSeq = MakeEmptyBioSequence()
 		newSeq.Write(sequence.Sequence()[from:to])
 
+		if sequence.HasQualities() {
+			newSeq.WriteQualities(sequence.Qualities()[from:to])
+		}
+
 		newSeq.sequence.id = fmt.Sprintf("%s_sub[%d..%d]", sequence.Id(), from+1, to)
 		newSeq.sequence.definition = sequence.sequence.definition
 	} else {
 		newSeq, _ = sequence.Subsequence(from, sequence.Length(), false)
 		newSeq.Write(sequence.Sequence()[0:to])
+
+		if sequence.HasQualities() {
+			newSeq.WriteQualities(sequence.Qualities()[0:to])
+		}
+
 	}
 
 	if len(sequence.Annotations()) > 0 {