Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,069,436 G→A I155I (ATC→ATT fbaA ← fructose‑bisphosphate aldolase, class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,069,4360GA76.2% 34.1 / 4.3 42I155I (ATC→ATTfbaAfructose‑bisphosphate aldolase, class II
Reads supporting (aligned to +/- strand):  ref base G (10/0);  new base A (32/0);  total (42/0)
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.85e-01

CAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGA  >  W3110S.gb/3069385‑3069483
                                                   |                                               
cAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTctct                                        >  1:32624/1‑61 (MQ=255)
cAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTctct                                        >  1:75264/1‑61 (MQ=255)
cAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTctct                                        >  1:54111/1‑61 (MQ=255)
cAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTctct                                        >  1:52816/1‑61 (MQ=255)
cAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTctct                                        >  1:40679/1‑61 (MQ=255)
cAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTctct                                        >  1:7689/1‑61 (MQ=255)
cAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTctct                                        >  1:63399/1‑61 (MQ=255)
cAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTctct                                        >  1:29704/1‑61 (MQ=255)
cAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTctct                                        >  1:25478/1‑61 (MQ=255)
cAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTctct                                        >  1:86547/1‑61 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACa                     >  1:748/1‑38 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTcaa               >  1:60393/1‑42 (MQ=38)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:70030/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:69082/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:71775/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:69757/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:1572/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:74861/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:76863/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:79829/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:82726/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:84714/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:85699/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:87364/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:87666/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:67083/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:13802/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:16963/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:20847/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:2985/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:31247/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:39438/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:40564/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:43901/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:55543/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:55978/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:56378/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:61862/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:64928/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaaga  >  1:67360/1‑57 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaag   >  1:40446/1‑56 (MQ=255)
                                          gCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCAATCATGTgagaag   >  1:9858/1‑56 (MQ=255)
                                                   |                                               
CAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGA  >  W3110S.gb/3069385‑3069483

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