Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,595,527 T→C T109A (ACG→GCG)  lgoR ← putative transcriptional activator for L‑galactonate catabolism

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,595,5270TC94.4% 108.4 / ‑2.7 36T109A (ACG→GCG) lgoRputative transcriptional activator for L‑galactonate catabolism
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (18/16);  total (19/17)
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.71e-01

CGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACG  >  NC_000913/4595395‑4595647
                                                                                                                                    |                                                                                                                        
cGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTcccc                                                                                                                    <  1:62814/139‑1 (MQ=255)
          cATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAAc                                                                                                          <  2:897022/139‑1 (MQ=255)
                                         aTAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTcccc                                                                                                                    >  1:834485/1‑98 (MQ=255)
                                         aTAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTcccc                                                                                                                    <  2:834485/98‑1 (MQ=255)
                                              aTTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCa                                                                      >  2:76380/1‑139 (MQ=255)
                                                tataACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCaaa                                                                    >  2:397550/1‑139 (MQ=255)
                                                  taACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAATGCCTGCTCcaa                                                                    <  2:622896/137‑1 (MQ=255)
                                                  taACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCaaa                                                                    >  1:622896/1‑137 (MQ=255)
                                                  taACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTccc                                                                                                                     >  1:582900/1‑88 (MQ=255)
                                                  taACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTccc                                                                                                                     <  2:582900/88‑1 (MQ=255)
                                                  taACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAaca                                                                  >  2:671670/1‑139 (MQ=255)
                                                  taACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAaca                                                                  >  1:144577/1‑139 (MQ=255)
                                                     cGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACt                                                                                                         >  1:1006141/1‑97 (MQ=255)
                                                     cGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACt                                                                                                         <  2:1006141/97‑1 (MQ=255)
                                                                  agcagatagctcgTACTACAACTGGCCCCACACCCGCGGCCCAAGCCAGCTGCAGTTCAGAAAACTCTTTCCCGCGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTATGTTTTGTTCACt                                                  <  2:267261/130‑1 (MQ=255)
                                                                        atTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGa                                            <  1:364809/139‑1 (MQ=255)
                                                                                actacAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGgtag                                    >  2:869675/1‑139 (MQ=255)
                                                                                actacAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGgtag                                    >  1:106335/1‑139 (MQ=255)
                                                                                      aCCGGGCTGACACCCGCGGCCCGCGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGATTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTAAAGCgactgtcactt                              >  1:267261/1‑129 (MQ=255)
                                                                                        cGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCg                            <  1:397550/139‑1 (MQ=255)
                                                                                        cGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCg                            >  2:202901/1‑139 (MQ=255)
                                                                                         gggCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGa                           <  1:76380/139‑1 (MQ=255)
                                                                                          ggCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTcccc                                                                                                                    <  2:962148/49‑1 (MQ=255)
                                                                                          ggCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTcccc                                                                                                                    >  1:962148/1‑49 (MQ=255)
                                                                                          ggCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGa                                            >  2:372262/1‑121 (MQ=255)
                                                                                          ggCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGa                                            <  1:372262/121‑1 (MQ=255)
                                                                                           gCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGcc                                                                            >  2:300679/1‑88 (MQ=255)
                                                                                           gCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGcc                                                                            <  1:300679/88‑1 (MQ=255)
                                                                                               acacCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGt                     >  2:491842/1‑139 (MQ=255)
                                                                                                  cccGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGgc                                                                                                      <  2:326563/55‑1 (MQ=255)
                                                                                                  cccGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGgc                                                                                                      >  1:326563/1‑55 (MQ=255)
                                                                                                   ccGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGtcatc                 <  2:106335/139‑1 (MQ=255)
                                                                                                     gcgGCCCGAGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTcatcat               >  1:262975/1‑139 (MQ=255)
                                                                                                             aGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCaa                                                                     >  1:800012/1‑77 (MQ=255)
                                                                                                             aGCCAGTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCaa                                                                     <  2:800012/77‑1 (MQ=255)
                                                                                                                  gTTGCAGTTCAGAAAACGCTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTAcg  <  2:262975/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                        
CGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACG  >  NC_000913/4595395‑4595647

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