Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,619,627 A→G 8.1% L203P (CTG→CCG)  opdH ← tricarboxylate‑specific outer membrane porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,619,6270AG 8.1% 134.9 / 2.6 49L203P (CTG→CCG) opdHtricarboxylate‑specific outer membrane porin
Reads supporting (aligned to +/- strand):  ref base A (21/24);  new base G (2/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

GATGAAGTGCTGGCGGTAGATATCGTCGAGCTGGGCGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCTTTGGAGGTGAGCAGGCCGCCCTCGAAGGTTTGCGGCAGGATACGCCCGTCGT  >  NC_002947/1619485‑1619762
                                                                                                                                               |                                                                                                                                       
gATGAAGTGCTGGCGGTAGATATCGTCGAGCTGGGCGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCa                                                                                                                                   <  2:323570/149‑1 (MQ=255)
 aTGAAGTGCTGGCGGTAGATATCGTCGAGCTGGGCGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAg                                                                                                                                  <  1:336103/149‑1 (MQ=255)
        gCTGGCGGTAGATATCGTCGAGCTGGGCGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTcc                                                                                                                           <  8:78288/149‑1 (MQ=255)
         cTGGCGGTAGATATCGTCGAGCTGGGCGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCTGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTcc                                                                                                                           <  1:245879/149‑1 (MQ=255)
         cTGGCGGTAGATATCGTCGAGCTGGGCGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTccc                                                                                                                          <  2:305565/149‑1 (MQ=255)
         cTGGCGGTAGATATCGTCGAGCTGGGCGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTcc                                                                                                                           >  7:187881/1‑148 (MQ=255)
            gCGGTAGATATCGTCGAGCTGGGCGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCTCCAGGTCTTCCCAGt                                                                                                                       >  4:31582/1‑149 (MQ=255)
             cGGTAGATATCGTCGAGCTGGGCGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGtt                                                                                                                      <  8:241950/149‑1 (MQ=255)
              ggTAGATATCGTCGAGCTGGGCGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTg                                                                                                                     <  5:165400/149‑1 (MQ=255)
               gTAGATATCGTCGAGCTGGGCGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTgg                                                                                                                    >  5:42522/1‑149 (MQ=255)
                    tatCGTCGAGCTGGGCGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCGGCGCCa                                                                                                                                   <  6:95971/129‑1 (MQ=255)
                    tatCGTCGAGCTGGGCGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCGGCGCCa                                                                                                                                   >  5:95971/1‑129 (MQ=255)
                    tatCGTCGAGCTGGGCGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCg                                                                                                               >  1:85903/1‑149 (MQ=255)
                               tGGGCGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCtt                                                                                                    >  2:210470/1‑149 (MQ=255)
                                 ggCGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCtt                                                                                                    <  4:72576/147‑1 (MQ=255)
                                   cGAAGTTGTAGCTGCAGGTGATGCGGTCGGTGAACTTGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAATGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGTTGTCGTTTCGGTCCtt                                                                                                    <  8:104981/145‑1 (MQ=255)
                                   cGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCt                                                                                                <  4:336808/149‑1 (MQ=255)
                                   agAAGTGGTAGCTGCCGGTGATGCGGTCTGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTATTACAAGTTGGTGTCGTTGCGGTCCtt                                                                                                    >  7:104981/2‑145 (MQ=255)
                                      aGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTg                                                                                             >  5:194143/1‑149 (MQ=255)
                                       gTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGt                                                                                            <  8:187881/149‑1 (MQ=255)
                                                         gCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCg                                                                                                          <  3:104183/117‑1 (MQ=255)
                                                         gCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTCAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCg                                                                                                          >  4:104183/1‑117 (MQ=255)
                                                         gCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGCCCAGGTAGTCAGCGCTGAAGGTGCAACCGAAGCTGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGATGGTGACGTTGCGGTCCTTGGCTTTGACTATGCGGACAACGGTg                                                                           >  4:29447/1‑148 (MQ=255)
                                                                 tGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAgg                                                                  <  1:210470/149‑1 (MQ=255)
                                                                           gTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCGGCGCCAGGTCTTCCCAGTTGGTGTCGTTg                                                                                                            >  1:354556/1‑97 (MQ=255)
                                                                           gTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCGGCGCCAGGTCTTCCCAGTTGGTGTCGTTg                                                                                                            <  2:354556/97‑1 (MQ=255)
                                                                              gAGCCCGCCCCGGGCCAGGTTGGCAGCGCTGAAGGTGCCACCCAAGCGGCCGTCCTTGTTGTTCAGCGCCAGGTCTTCCCAGTGGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCtt                                                     <  6:105173/149‑1 (MQ=255)
                                                                              gAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCtt                                                     >  8:114162/1‑149 (MQ=255)
                                                                                gcccgcccAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCTTTg                                                   <  7:114162/149‑1 (MQ=255)
                                                                                gcccgcccAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCGTTGTCTATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCTTTg                                                   >  2:279685/1‑149 (MQ=255)
                                                                                gcccgcccAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCCCCAGGTCATCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGCCCCACGGTGAAGGTCAGGCCCTTGATCTCTTTg                                                   >  2:308339/1‑149 (MQ=255)
                                                                                 cccgcccAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCTTTgg                                                  <  6:194143/149‑1 (MQ=255)
                                                                                  ccgcccAGGGCCTGGTTGTCAGCGCTGAAGGTGCCTCCGTAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTGGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGTAGGTCAGGCCCTTGATCTCTTTGGa                                                 <  3:29447/149‑1 (MQ=255)
                                                                                   cgcccAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCTTTGGAg                                                >  6:212494/1‑149 (MQ=255)
                                                                                     cccAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCTTTGGAGGt                                              <  2:85903/149‑1 (MQ=255)
                                                                                                    tCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGt                                                                      >  7:245984/1‑110 (MQ=255)
                                                                                                    tCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGt                                                                      <  8:245984/110‑1 (MQ=255)
                                                                                                     cAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCTTTGGAGGTGAGCAGGCCGCCCTc                               >  3:180095/1‑148 (MQ=255)
                                                                                                      aGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCTTTGGAGGTGAGCAGGCCGCCCTCg                              <  6:42522/148‑1 (MQ=255)
                                                                                                           tGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCTTTGGAGGTGAGCAGGCCGCCCTCGAAGGtt                        >  3:100501/1‑149 (MQ=255)
                                                                                                              aGGTGACACCGATGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCTTTGGAGGTGAGCAGGCCGCCCTCGAAGGTTTGc                     <  3:31582/149‑1 (MQ=255)
                                                                                                                gTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCTTTGGAGGTGAGCAGGCCGCCCTCGaa                            >  4:198737/1‑140 (MQ=255)
                                                                                                                gTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCTTTGGAGGTGAGCAGGCCGCCCTCGaa                            <  3:198737/140‑1 (MQ=255)
                                                                                                                gTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCTTTGGAGGTGAGCAGGCCGCCCTCGAAGGTTTGCgg                   <  4:180095/149‑1 (MQ=255)
                                                                                                                     aCCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTATATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCTTTGGAGGTGAGCAGGCCGCCCTCGAAGGTTTGCGGCAGGa              >  3:143412/1‑149 (MQ=255)
                                                                                                                               gCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAg                                                                   >  6:182109/1‑87 (MQ=255)
                                                                                                                               gCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAg                                                                   <  5:182109/87‑1 (MQ=255)
                                                                                                                                 cGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGATCGTGTCGTTGCGGACCTTGGCTTTGTCTAAGCCGCCACACGTGAAGGTCAGGCCCTAGATCGCATTGGAGGTGAGCAGCCCGCCCTCGCAGGAATGCAGCAGGATACgaaggtcg   >  4:245415/1‑142 (MQ=255)
                                                                                                                                  gTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCTTTGGAGGTGAGCAGGCCGCCCTCGAAGGTTTGCGGCAGGATACGCCcgtcgt  >  7:117351/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                       
GATGAAGTGCTGGCGGTAGATATCGTCGAGCTGGGCGAAGTGGTAGCTGCCGGTGATGCGGTCGGTGAACTGGTAGTCGAGCCCGCCCAGGGCCAGGTTGTCAGCGCTGAAGGTGCCACCGAAGCGGCCGTTCTTGTTGTTCAGCGCCAGGTCTTCCCAGTTGGTGTCGTTGCGGTCCTTGGCTTTGTCTATGCGGCCACCGGTGAAGGTCAGGCCCTTGATCTCTTTGGAGGTGAGCAGGCCGCCCTCGAAGGTTTGCGGCAGGATACGCCCGTCGT  >  NC_002947/1619485‑1619762

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