Skip to content

Commit 82c98e4

Browse files
committed
skip Viterbi for sequences shorter than 3 nucleotides
1 parent 08124bd commit 82c98e4

File tree

1 file changed

+1
-1
lines changed

1 file changed

+1
-1
lines changed

src/bin/FragGeneScanRs.rs

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -245,7 +245,7 @@ fn run<R: Read + Send, W: WritingBuffer + Send>(
245245
.map(u8::to_ascii_uppercase)
246246
.map(Nuc::from)
247247
.collect();
248-
let read_prediction = if nseq.is_empty() {
248+
let read_prediction = if nseq.len() < 3 {
249249
gene::ReadPrediction::new(head)
250250
} else {
251251
viterbi(&global, &locals, head, nseq, whole_genome)

0 commit comments

Comments
 (0)