diff --git a/seqkit/cmd/seq.go b/seqkit/cmd/seq.go index 2264fac..35430ea 100644 --- a/seqkit/cmd/seq.go +++ b/seqkit/cmd/seq.go @@ -545,19 +545,23 @@ Filtering records to edit: } if matched && filterMinLen && len(record.Seq.Seq) < minLen { + log.Info("sequence '%s' too short (%d < %d)", record.Seq.Id, len(record.Seq.Seq), minLen) continue } if matched && filterMaxLen && len(record.Seq.Seq) > maxLen { + log.Info("sequence '%s' too long (%d > %d)", record.Seq.Id, len(record.Seq.Seq), maxLen) continue } if matched && (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 } }