Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,798,703 G→C 6.5% intergenic (‑208/‑461) stpA ← / → alaE DNA binding protein, nucleoid‑associated/alanine exporter, alanine‑inducible, stress‑responsive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,798,7030GC 6.5% 211.0 / 3.2 61intergenic (‑208/‑461)stpA/alaEDNA binding protein, nucleoid‑associated/alanine exporter, alanine‑inducible, stress‑responsive
Reads supporting (aligned to +/- strand):  ref base G (29/28);  new base C (2/2);  total (31/30)
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.54e-01

ACAAAATTTAATCTATTTACAGTCCTGCGCGAGATTATTTCAGCTGTTATTTTCTGAAACTCTCAGAATAATTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTCTGCTGTTAATGTACATTGTCTGGCATAAAAAATGCTGAAAATGA  >  NC_000913/2798573‑2798844
                                                                                                                                  |                                                                                                                                             
acaAAATTTAATCTATTTACAGTCCTGCGCGAGATTATTTCAGCTGTTATTTTCTGAAACTCTCAGAATAATTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAAt                                                                                                                             >  1:628909/1‑149 (MQ=255)
 caAAATTTAATCTATTTACAGTCCTGCGCGAGATTATTTCAGCTGTTATTTTCTGAAACTCTCAGAATAATTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAAtt                                                                                                                            <  2:628909/149‑1 (MQ=255)
 caAAATTTAATCTATTTACAGTCCTGCGCGAGATTATTTCAGCTGTTATTTTCTGAAACTCTCAGAATAATTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAAtt                                                                                                                            >  4:1689794/1‑149 (MQ=255)
                      tCCTGCGCGAGATTATTTCAGCTGTTATTTTCTGAAACTCTCAGAATAATTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCaa                                                                                                       >  4:1061251/1‑149 (MQ=255)
                       ccTGCGCGAGATTATTTCAGCTGTTATTTTCTGAAACTCTCAGAATAATTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCaaa                                                                                                      <  3:1689794/149‑1 (MQ=255)
                            gcgAGATTATTTCAGCTGTTATTTTCTGAAACTCTCAGAATAATTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGaa                                                                                                 <  3:1061251/149‑1 (MQ=255)
                            gcgAGATTATTTCAGCTGTTATTTTCTGAAACTCTCAGAATAATTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGaa                                                                                                 >  3:1418919/1‑149 (MQ=255)
                                              ttATTTTCTGAAACTCTCAGAATAATTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACa                                                                               <  4:1418919/149‑1 (MQ=255)
                                                     cTGAAACTCTCAGAATAATTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATaaa                                                                        >  1:64518/1‑149 (MQ=255)
                                                        aaaCTCTCAGAATAATTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGa                                                                     >  2:1666467/1‑149 (MQ=255)
                                                                  aataatTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAAt                                                                                          >  4:825126/1‑118 (MQ=255)
                                                                  aataatTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAAt                                                                                          <  3:825126/118‑1 (MQ=255)
                                                                   ataatTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAAcatccatc                                                          >  2:1748142/1‑149 (MQ=255)
                                                                     aatTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCAc                                                        <  2:64518/149‑1 (MQ=255)
                                                                      atTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGCCGctctgtt                                                                                                                                      <  2:952007/70‑4 (MQ=255)
                                                                      atTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGCCGctctgtt                                                                                                                                      >  1:952007/1‑67 (MQ=255)
                                                                      atTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGCCGctctgtt                                                                                                                                      <  2:1044436/70‑4 (MQ=255)
                                                                      atTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGCCGctctgtt                                                                                                                                      >  1:1044436/1‑67 (MQ=255)
                                                                       tttCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCa                                                                                       >  1:791625/1‑116 (MQ=255)
                                                                       tttCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCa                                                                                       <  2:791625/116‑1 (MQ=255)
                                                                       tttCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTg                                                      >  4:426795/1‑149 (MQ=255)
                                                                        ttCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTgg                                                     >  4:691744/1‑149 (MQ=255)
                                                                         tCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGt                                                    <  3:426795/149‑1 (MQ=255)
                                                                            aGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCt                                                 <  1:1666467/149‑1 (MQ=255)
                                                                             gATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAata                                                                                    >  3:574094/1‑113 (MQ=255)
                                                                             gATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAata                                                                                    <  4:574094/113‑1 (MQ=255)
                                                                              aTTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATTAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTg                                                      <  3:67102/142‑1 (MQ=255)
                                                                              aTTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATTAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTg                                                      >  4:67102/1‑142 (MQ=255)
                                                                               ttATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACa                                                                               >  2:2487/1‑116 (MQ=255)
                                                                               ttATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACa                                                                               <  1:2487/116‑1 (MQ=255)
                                                                                   atatTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTctgct                                          <  1:1748142/149‑1 (MQ=255)
                                                                                      ttAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCaaa                                                                                                      >  1:1345546/1‑86 (MQ=255)
                                                                                      ttAGCCATGTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCCCTAATTCCGCAAGCCTGACAAAATCaaa                                                                                                      <  2:1345546/86‑1 (MQ=255)
                                                                                         gCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCaaa                                                                                                      >  1:1168638/1‑83 (MQ=255)
                                                                                         gCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCaaa                                                                                                      <  2:1168638/83‑1 (MQ=255)
                                                                                          ccATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCaaa                                                                                                      >  3:1047046/1‑82 (MQ=255)
                                                                                          ccATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCaaa                                                                                                      <  4:1047046/82‑1 (MQ=255)
                                                                                             ttttGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTCTGCTGTTAATGTAc                                >  2:192544/1‑149 (MQ=255)
                                                                                                   tgaGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTCTGCTGTTAATGTACATTGTc                          <  1:192544/149‑1 (MQ=255)
                                                                                                   tgaGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTCTGCTGTTAATGTACATTGTc                          >  4:955018/1‑149 (MQ=255)
                                                                                                     aGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGa                                                                                                                <  2:1685073/61‑1 (MQ=255)
                                                                                                     aGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGa                                                                                                                >  1:1685073/1‑61 (MQ=255)
                                                                                                         aTTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTCTGCTGt                                        >  1:1173980/1‑129 (MQ=255)
                                                                                                         aTTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTCTGCTGt                                        <  2:1173980/129‑1 (MQ=255)
                                                                                                         aTTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTCTGCTGt                                        >  3:70277/1‑129 (MQ=255)
                                                                                                         aTTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTCTGCTGt                                        <  4:70277/129‑1 (MQ=255)
                                                                                                              gATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTCTGCTGtt                                       <  3:1037230/125‑1 (MQ=255)
                                                                                                              gATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTCTGCTGtt                                       >  4:1037230/1‑125 (MQ=255)
                                                                                                                ttattGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTCTGCTGTTAATGTACATTGTCTGGCATAAAAAAt             <  3:691744/149‑1 (MQ=255)
                                                                                                                 tattGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGaa                                                                    >  1:40425/1‑93 (MQ=255)
                                                                                                                 tattGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGaa                                                                    <  2:40425/93‑1 (MQ=255)
                                                                                                                  attGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAAcatcc                                                             <  1:428783/99‑1 (MQ=255)
                                                                                                                  attGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAAcatcc                                                             >  2:428783/1‑99 (MQ=255)
                                                                                                                        gAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGaa                                                                            >  3:1195628/1‑78 (MQ=255)
                                                                                                                        gAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGaa                                                                            <  4:1195628/78‑1 (MQ=255)
                                                                                                                           cAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTctgc                                           <  2:404272/108‑1 (MQ=255)
                                                                                                                           cAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTctgc                                           >  1:404272/1‑108 (MQ=255)
                                                                                                                           cAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTCTGCTGTTAATGTACATTGTCTGGCATAAAAAATGCTGAAAATGa  <  3:955018/149‑1 (MQ=255)
                                                                                                                           cAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTCTGCTGTTAATGTACATTGTCTGGCATAAAAAATGCTGAAAATGa  >  4:910406/1‑149 (MQ=255)
                                                                                                                            agagCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTCTGCTGTTAATGTACATTGTCTGGCATAAAAAATGCTGaaa      <  1:1487404/144‑1 (MQ=255)
                                                                                                                            agagCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTCTGCTGTTAATGTACATTGTCTGGCATAAAAAATGCTGaaa      >  2:1487404/1‑144 (MQ=255)
                                                                                                                                 ggcgCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGaaa                                                                   >  3:887838/1‑78 (MQ=255)
                                                                                                                                 ggcgCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGaaa                                                                   <  4:887838/78‑1 (MQ=255)
                                                                                                                                  gcgcTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATaaaa                                                                       <  1:1187591/73‑1 (MQ=255)
                                                                                                                                  gcgcTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATaaaa                                                                       >  2:1187591/1‑73 (MQ=255)
                                                                                                                                  |                                                                                                                                             
ACAAAATTTAATCTATTTACAGTCCTGCGCGAGATTATTTCAGCTGTTATTTTCTGAAACTCTCAGAATAATTTCAAGATTATATATTAGCCATTTTGATGAGTAATTTAGATTATTGGTGAACAGAGCGGCGCTCTCTGCTCGCTAATTCCGCAAGCCTGACAAAATCAAAGTGAATAAAAATGCAATATAACAGAATAAAAGAAAACATCCATCACTGGTGCTTTCTGCTGTTAATGTACATTGTCTGGCATAAAAAATGCTGAAAATGA  >  NC_000913/2798573‑2798844

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