Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,969,921 G→A 22.7% L99L (CTG→TTG)  cheB ← fused chemotaxis regulator and protein‑glutamate methylesterase in two‑component regulatory system with CheA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,969,9210GA22.7% 36.2 / 7.0 22L99L (CTG→TTG) cheBfused chemotaxis regulator and protein‑glutamate methylesterase in two‑component regulatory system with CheA
Reads supporting (aligned to +/- strand):  ref base G (7/10);  new base A (3/2);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

ATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTGGTGACAAAATCTATCGCCCCCAGCTCCAGCGCGCGCAGCGTGACTTCTGACCCTTTGCCGGTCAGGGAAGAAACCATCACAACGGGCATTG  >  W3110S.gb/1969851‑1969989
                                                                      |                                                                    
aTTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTGGTGACAAAATCTATCGCCCCCAg                                                                      >  1:3646087/1‑71 (MQ=255)
         taCGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTGGTGACAAAATCTATCGCCCCCAGCTCCAgcg                                                              <  1:423672/70‑1 (MQ=255)
                gCATTCCTTCGCGAATACCCAGTTGCGGTTTGGTGACAAAATCTATCGCCCCCAACTCCAgcgc                                                             <  1:3803452/64‑1 (MQ=255)
                gCATTCCTTCGCGAATACCCAGTTGCGGTTTGGTGACAAAATCTATCGCCCCCAACTCCAgcgc                                                             <  1:2704345/64‑1 (MQ=255)
                 cATACCTTCGCGAATACCCAGTTGCGGTTTGGTGACAAAATCTATCGCCCCCAGCTCCAGCGCGCGCAGc                                                      <  1:1650734/70‑1 (MQ=255)
                              aTACCCAGTTGCGGTTTGGTGACAAAATCTATCGCCCCCAGCTCCAGCGCGCGCAGCGTGACTTCTGAccc                                        <  1:382020/71‑1 (MQ=255)
                                   cAGTTGCGGTTTGGTGACAAAATCTATCGCCCCCAGCTCCAGCGCGCGCAGCGt                                                    <  1:1903780/54‑1 (MQ=255)
                                     gTTGCGGTTTGGTGACAAAATCTATCGCCCCCAGCTc                                                                   <  1:2539898/37‑1 (MQ=255)
                                       tGCGGTTTGGTGACAAAATCTATCGCCCCCAGCTCCAGCGCGCGCAGCGTGACTTCTGACCCTTTGCCGGt                               >  1:2939815/1‑71 (MQ=255)
                                       tGCGGTTTGGTGACAAAATCTATCGCCCCCAGCTCCAGCGCGCGCAGCGTGACTTCTGACCCTTTGCCGGt                               >  1:1716371/1‑71 (MQ=255)
                                         cGGTTTGGTGACAAAATCTATCGCCCCCAGCTCCAGCGCGCGCAGCGTGACTTCTGACCCTTTGCCGGt                               >  1:1204002/1‑69 (MQ=255)
                                             ttGGTGACAAAATCTATCGCCCCCAGCTCCAGCGCGCGCAGCGTGACTTCTGACCCTTTGCCGGTCAggg                          <  1:834159/70‑1 (MQ=255)
                                              tGGTGACAAAATCTATCGCCCCCAGCTCCAGCGCGCGCAGCGTGACTTCTGACCCTTTGCCGGTCAggg                          >  1:2865034/1‑69 (MQ=255)
                                              tGGTGACAAAATCTATCGCCCCCAGCTCCAGCGCGCGCAGCGTGACTTCTGACCCTTTGCCGGTCAggg                          >  1:3452688/1‑69 (MQ=255)
                                              tGGTGACAAAATCTATCGCCCCCAACTCCAGCGCGCGCAGCGTGACTTCTGACCCTTTGCCGGTCAggg                          >  1:1834248/1‑69 (MQ=255)
                                                   aCAAAATCTATCGCCCCCAGCTCCAGCGCGCGCAGCGTGACTTCTGACCCTTTGCCGGTCAGGgaa                        >  1:396893/1‑66 (MQ=255)
                                                       aaTCTATCGCCCCCAACTCCAGCGCGCGCAGCGTGACTTCTg                                            >  1:1967688/1‑42 (MQ=255)
                                                       aaTCTATCGCCCCCAACTCCAGCGCGCGCAGCGTGACTTCTg                                            >  1:556579/1‑42 (MQ=255)
                                                             tCGCCCCCAGCTCCAGCGCGCGCAGCGTGACTTCTGACCCTTTGCCGGTCAGGGAAGAAACCATCACAAc          <  1:2737315/70‑1 (MQ=255)
                                                                   ccAGCTCCAGCGCGCGCAGCGTGACTTCTGACCCTTTGCCGGTCAGGGAAGAAACCATCACAACGGGCAtt   <  1:1446992/71‑1 (MQ=255)
                                                                   ccAGCTCCAGCGCGCGCAGCGTGACTTCTGACCCTTTGCCGGTCAGGGAAGAAACCATCACAACGGGCAtt   <  1:1298350/71‑1 (MQ=255)
                                                                    cAGCTCCAGCGCGCGCAGCGTGACTTCTGACCCTTTGCCGGTCAGGGAAGAAACCATCACAACGGGCATTg  >  1:2755072/1‑71 (MQ=255)
                                                                    cAGCTCCAGCGCGCGCAGCGTGACTTCTGACCCTTTGCCGGTCAGGGAAGAAACCATCACAACGGGCATTg  <  1:1451756/71‑1 (MQ=255)
                                                                      |                                                                    
ATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTGGTGACAAAATCTATCGCCCCCAGCTCCAGCGCGCGCAGCGTGACTTCTGACCCTTTGCCGGTCAGGGAAGAAACCATCACAACGGGCATTG  >  W3110S.gb/1969851‑1969989

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