Tests

src/longsequences/operators.jl

@@ -362,7 +362,7 @@ function _findlast(seq::SeqOrView{<:KNOWN_ALPHABETS}, enc::UInt64)
     while body_i - unroll + 2 > body_stop
         u = zero(UInt64)
         for j in 0:unroll - 1
-            u |= set_zero_encoding(BitsPerSymbol(seq), @inbounds(data[body_i + j]) ⊻ enc)
+            u |= set_zero_encoding(BitsPerSymbol(seq), @inbounds(data[body_i - j]) ⊻ enc)
         end
         iszero(u) || break
         i -= symbols_per_data_element(seq) * unroll

test/longsequences/basics.jl

@@ -341,6 +341,13 @@ end
         ACGTN
         """
         @test a == b
+
+        a = randdnaseq(512)
+        a[111] = DNA_G
+        b = copy(a)
+        @test a == b
+        a[111] = DNA_C
+        @test a != b
     end
 
     @testset "RNA" begin
@@ -373,6 +380,10 @@ end
         X
         """
         @test a == b
+
+        a = randaaseq(131)
+        b = copy(a)
+        @test a == b
     end
 end
 

test/longsequences/find.jl

@@ -64,6 +64,16 @@
     @test findnext(==(DNA_M), seq, 8) == 8
     @test findnext(==(DNA_M), seq, 9) == 21
 
+    # Very large sequence so it hits the inner SIMD loop
+    seq = randdnaseq(541)
+    seq[100] = DNA_Gap
+    @test findfirst(==(DNA_Gap), seq) == findfirst(isgap, seq) == 100
+    seq[100] = DNA_A
+    seq[500] = DNA_Gap
+    @test findlast(==(DNA_Gap), seq) == findlast(isgap, seq) == 500
+    seq[500] = DNA_A
+    @test findfirst(==(DNA_Gap), seq) == findlast(isgap, seq) == nothing
+
     # View with only tail
     #               1234
     seq = view(aa"KWYPAV-L", 3:6)

test/longsequences/seqview.jl

@@ -44,6 +44,20 @@ end
 	@test String(seq) == "ANKYH"
 end
 
+@testset "Equality" begin
+    for size in [41, 504, 7]
+        for offset in [0, 3, 32]
+			seq = randrnaseq(size)
+            seq2 = view(randrnaseq(offset) * seq * randrnaseq(15), offset+1:offset+size)
+            @test seq == seq2
+            seq[4] = RNA_Gap
+            @test seq != seq2
+            seq3 = view(seq, 1:size)
+            @test seq == seq3
+        end
+    end
+end
+
 # Added after issue 260
 @testset "Random construction" begin
 	for i in 1:100