Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,051,158 A→G D2074G (GAC→GGC)  yeeJ → putative adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,051,1580AG88.7% 170.5 / 9.0 62D2074G (GAC→GGC) yeeJputative adhesin
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base G (26/29);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

GAGCGATGACAGGTAGCGTCACGGTAAGCGCAGTCACGACCGCTGGTGGAATGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGACTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATCAGGTACTAACGTGCCCTATATAAAAATTAGCGCAATTGATTACA  >  NC_000913/2051020‑2051292
                                                                                                                                          |                                                                                                                                      
gAGCGATGACAGGTAGCGTCACGGTAAGCGCAGTCACGACCGCTGGTGGAATGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGa                                                                                                                                        <  1:628806/139‑1 (MQ=255)
    gATGACAGGTAGCGTCACGGTAAGCGCAGTCACGACCGCTGGTGGAATGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCtat                                                                                                                                    <  1:713860/139‑1 (MQ=255)
         cAGGTAGCGTCACGGTAAGCGCAGTCACGACCGCTGGTGGAATGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGa                                                                                                                               <  2:382356/139‑1 (MQ=255)
                cGTCACGGTAAGCGCAGTCACGACCGCTGGTGGAATGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCt                                                                                                                        >  1:722380/1‑139 (MQ=255)
                    aCGGTAAGCGCAGTCACGACCGCTGGTGGAATGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCA‑‑‑‑GTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACg                                                                                                                >  2:652890/1‑139 (MQ=255)
                      ggTAAGCGCAGTCACGACCGCTGGTGGAATGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTa                                                                                                                  >  2:805094/1‑139 (MQ=255)
                           gcgcAGTCACGACCGCTGGTGGAATGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAAtt                                                                                                                   >  2:422347/1‑133 (MQ=255)
                           gcgcAGTCACGACCGCTGGTGGAATGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAAtt                                                                                                                   <  1:422347/133‑1 (MQ=255)
                                    cGACCGCTGGTGGAATGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCa                                                                                                    <  1:805094/139‑1 (MQ=255)
                                     gACCGCTGGTGGAATGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCAc                                                                                                   <  1:1166305/139‑1 (MQ=255)
                                      aCCGCTGGTGGAATGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGACTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACg                                                                                                  >  1:499528/1‑139 (MQ=255)
                                           tggtggAATGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACg                                                                                                  >  2:1270664/1‑134 (MQ=255)
                                           tggtggAATGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACg                                                                                                  <  1:1270664/134‑1 (MQ=255)
                                            ggtggAATGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGatatat                                                                                            <  2:1371828/139‑1 (MQ=255)
                                                 aaTGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGc                                                                                       <  2:1263962/139‑1 (MQ=255)
                                                 aaTGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGACTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGc                                                                                       <  2:499528/139‑1 (MQ=255)
                                                     cAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCCGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATAACGAGAGTGCTGAATTACGTCTTTTTCTGCACGATATATCAGGCAAta                                                                                   >  2:1107913/1‑138 (MQ=255)
                                                           gTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCa                                                                             >  2:762262/1‑139 (MQ=255)
                                                              gATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACg                                                                                                  >  1:1070249/1‑115 (MQ=255)
                                                              gATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACg                                                                                                  <  2:1070249/115‑1 (MQ=255)
                                                                      gCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCtat                                                                                                                                    >  1:177087/1‑73 (MQ=255)
                                                                      gCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCtat                                                                                                                                    <  2:177087/73‑1 (MQ=255)
                                                                         ggtggCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATCGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAggg                                                               >  1:986017/1‑139 (MQ=255)
                                                                          gtggCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAAt                                                                                    <  1:1369777/117‑1 (MQ=255)
                                                                          gtggCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAAt                                                                                    >  2:1369777/1‑117 (MQ=255)
                                                                          gtggCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGa                                                              >  2:865576/1‑139 (MQ=255)
                                                                          gtggCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGACTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGa                                                              >  1:804464/1‑139 (MQ=255)
                                                                          gtggCTGGCCCGGCAGACACCTCCCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGATTTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGa                                                              >  1:920439/1‑139 (MQ=255)
                                                                           tggctggcCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGAt                                                             <  1:1107913/139‑1 (MQ=255)
                                                                               tggcCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTg                                                                                                      >  1:1145357/1‑94 (MQ=255)
                                                                               tggcCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTg                                                                                                      <  2:1145357/94‑1 (MQ=255)
                                                                               tggcCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGaa                                                         <  2:722380/139‑1 (MQ=255)
                                                                                ggcCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAAt                                                        <  2:539601/139‑1 (MQ=255)
                                                                                 gcCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAAtt                                                       >  1:1166991/1‑139 (MQ=255)
                                                                                        aGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGACTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCaa                                                <  2:804464/139‑1 (MQ=255)
                                                                                          acacCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATCGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATc                                              <  2:986017/139‑1 (MQ=255)
                                                                                          acacCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATc                                              <  1:865576/139‑1 (MQ=255)
                                                                                               tggCAGTCCAGCCTTAAGAGGAATCGGGCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATGACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATCAGGTa                                         <  2:920439/137‑1 (MQ=255)
                                                                                                cGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATCAGGTAc                                        >  2:394579/1‑139 (MQ=255)
                                                                                                     tccgtccTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGa                                                                   <  2:481761/107‑1 (MQ=255)
                                                                                                     tccgtccTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGa                                                                   >  1:481761/1‑107 (MQ=255)
                                                                                                      ccgtccTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGtt                                                                                                         <  1:232732/68‑1 (MQ=255)
                                                                                                      ccgtccTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGtt                                                                                                         >  2:232732/1‑68 (MQ=255)
                                                                                                       agtccTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATCAGGTACTAACGTg                                 <  1:652890/138‑1 (MQ=255)
                                                                                                         tccTTAAGAGCAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATCAGGTACTAACGTGcc                               >  1:1336802/1‑139 (MQ=255)
                                                                                                                 agCAATCGGTCAACACTGAAAGGGGACTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATCAGGTACTAACGTGCCCTATATaa                       >  2:300040/1‑139 (MQ=255)
                                                                                                                   cAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGaa                                                                                                                     <  1:861499/43‑1 (MQ=255)
                                                                                                                   cAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGaa                                                                                                                     >  2:861499/1‑43 (MQ=255)
                                                                                                                   cAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGAt                                                                               <  2:1345532/81‑1 (MQ=255)
                                                                                                                   cAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGAt                                                                               >  1:1345532/1‑81 (MQ=255)
                                                                                                                   cAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAg                                                                 <  2:1394466/95‑1 (MQ=255)
                                                                                                                   cAATCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAg                                                                 >  1:1394466/1‑95 (MQ=255)
                                                                                                                      tCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCa                                                 >  2:1432356/1‑108 (MQ=255)
                                                                                                                      tCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCa                                                 <  1:1432356/108‑1 (MQ=255)
                                                                                                                      tCGGTCATCACTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATCAGGTACTAACGTGCCCTATATAAAAAtt                  >  1:1406029/1‑139 (MQ=255)
                                                                                                                           catcaCTGAAAGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATCAGGTACTAACGTGCCCTATATAAAAATTAgcgc             >  1:794693/1‑139 (MQ=255)
                                                                                                                                   aaaGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATCAGGTACTAAc                                    <  2:1595015/108‑1 (MQ=255)
                                                                                                                                   aaaGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATCAGGTACTAAc                                    >  1:1595015/1‑108 (MQ=255)
                                                                                                                                    aaGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATCAGGTACTAACGTGCCCTATATAAAAATTAGCGCAATTGATTa    <  2:1336802/139‑1 (MQ=255)
                                                                                                                                    aaGGGGGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATCAGGTACTAACGTGCCCTATATAAAAATTAGCGCAATTGATTa    <  2:794693/139‑1 (MQ=255)
                                                                                                                                      ggggGCTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATCAGGTACTAACGTGCCCTATATAAAAATTAGCGCAATTGATTACa  <  1:394579/139‑1 (MQ=255)
                                                                                                                                        ggACTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATc                                              <  2:1106950/93‑1 (MQ=255)
                                                                                                                                        ggACTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATc                                              >  1:1106950/1‑93 (MQ=255)
                                                                                                                                          |                                                                                                                                      
GAGCGATGACAGGTAGCGTCACGGTAAGCGCAGTCACGACCGCTGGTGGAATGCAAACAGTAGATATAACGCTGGTGGCTGGCCCGGCAGACACCTCGCAGTCCGTCCTTAAGAGCAATCGGTCATCACTGAAAGGGGACTATACCGATAGTGCTGAATTACGTCTTGTTCTGCACGATATATCAGGCAATCCGATCAAAGTTTCTGAAGGGATGGAATTTGTGCAATCAGGTACTAACGTGCCCTATATAAAAATTAGCGCAATTGATTACA  >  NC_000913/2051020‑2051292

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