Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,987,439:1 +T 100% coding (284/1047 nt) ansB ← periplasmic L‑asparaginase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,987,4391.T83.9% 89.3 / 14.3 31coding (284/1047 nt)ansBperiplasmic L‑asparaginase II
Reads supporting (aligned to +/- strand):  ref base . (0/5);  new base T (0/26);  total (0/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

TGGGTAATGACGAAGCCGTCGGTCTTATCGCAGTCGGTGTTAATTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGA  >  minE/1987397‑1987490
                                           |                                                    
tGGGTAATGACGAAGCCGTCGGTCTTGTCGCAGTCGGTGTTAATTTTTTTCGCCAGTGTCAGCCAGACAtt                         <  1:101203/71‑1 (MQ=255)
tGGGTAATGACGAAGCCGTCGGTCTTATCGCAGTCGGTGTTAATTTTTTTCGCCAGTGTCAGCCAGACAtt                         <  1:38505/71‑1 (MQ=255)
tGGGTAATGACGAAGCCGTCGGTCTTATCGCAGTCGGTGTTAATTTTTTTCGCCAGTGTCAGCCAGACAtt                         <  1:351237/71‑1 (MQ=255)
cgggTAATGACGAAGCCGTCGGGCTTATCGCAGTCGGTGTTAATTTTTTTCGCCAGTGTCAGCCAGACAtt                         <  1:472354/70‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:3472/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:95348/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:547811/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:545077/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:514088/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:436159/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:431391/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:423630/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:370433/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:357605/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:349590/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:347734/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:300819/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:290000/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:289812/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:287774/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:279957/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:269129/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:261656/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:190428/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:181491/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:178533/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:154959/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:150576/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:112306/69‑1 (MQ=255)
                          aTCGCAGTCGGTGTTAATTTTTTCTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:627058/69‑1 (MQ=255)
                               aGTCGGTGTTAATTTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGa  <  1:236975/64‑1 (MQ=255)
                                           |                                                    
TGGGTAATGACGAAGCCGTCGGTCTTATCGCAGTCGGTGTTAATTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCATGTCCTGGGAGCCGA  >  minE/1987397‑1987490

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: