From 1d6b319ee1fafbaac462815c594ff7d038192efb Mon Sep 17 00:00:00 2001 From: "Filipe G. Vieira" <1151762+fgvieira@users.noreply.github.com> Date: Mon, 21 Oct 2024 13:36:26 +0200 Subject: [PATCH] Add logs to filtering --- seqkit/cmd/seq.go | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/seqkit/cmd/seq.go b/seqkit/cmd/seq.go index 7b38b53a..351e60d1 100644 --- a/seqkit/cmd/seq.go +++ b/seqkit/cmd/seq.go @@ -269,19 +269,23 @@ var seqCmd = &cobra.Command{ } if filterMinLen && len(record.Seq.Seq) < minLen { + log.Info("sequence '%s' too short (%d < %d)", record.Seq.Id, len(record.Seq.Seq), minLen) continue } if filterMaxLen && len(record.Seq.Seq) > maxLen { + log.Info("sequence '%s' too long (%d > %d)", record.Seq.Id, len(record.Seq.Seq), maxLen) continue } if filterMinQual || filterMaxQual { avgQual := record.Seq.AvgQual(qBase) if filterMinQual && avgQual < minQual { + log.Info("average quality of sequence '%s' too low (%d < %d)", record.Seq.Id, avgQual, minQual) continue } if filterMaxQual && avgQual >= maxQual { + log.Info("average quality of sequence '%s' too high (%d > %d)", record.Seq.Id, avgQual, maxQual) continue } }