Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 289,771 C→G 7.3% A397G (GCC→GGC)  oprE → outer‑membrane porin E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947289,7710CG 7.3% 171.4 / 3.7 55A397G (GCC→GGC) oprEouter‑membrane porin E
Reads supporting (aligned to +/- strand):  ref base C (28/23);  new base G (2/2);  total (30/25)
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.94e-01

ATACAACCTTTGGTCAGTACTCCTATGACTTTGCGGCTCTGGGTGTTCCAGGCCTGAAGGCTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGC  >  NC_002947/289633‑289918
                                                                                                                                          |                                                                                                                                                   
aTACAACCTTTGGTCAGTACTCCTATGACTTTGCGGCTCTGGGTGTTCCAGGCCTGAAGGCTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGAc                                                                                                                                           <  5:61298/149‑1 (MQ=255)
  aCAACCTTTGGTCAGTACTCCTATGACTTTGCGGCTCTGGGTGTTCCAGGCCTGAAGGCTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTa                                                                                                                                         >  4:330010/1‑149 (MQ=255)
      ccTTTGGTCAGTACTCCTATGACTTTGCGGCTCTGGGTGTTCCAGGCCTGAAGGCTTCGGTGTCCTACCTGCGTGGTGATTATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGAAGCCCGCGTCGACTACAcc                                                                                                                                     >  2:161656/1‑149 (MQ=255)
               aGTACTCCTATGACTTTGCGGCTCTGGGTGTTCCAGGCCTGAAGGCTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAgg                                                                                                                              >  8:212065/1‑147 (MQ=255)
                gTACTCCTATGACTTTGCGGCTCTGGGTGTTCCAGGCCTGAAGGCTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAg                                                                                                                           <  2:267320/149‑1 (MQ=255)
                 tACTCCTATGACTTTGCGGCTCTGGGTGTTCCAGGCCTGAAGGCTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAgg                                                                                                                          <  7:212065/149‑1 (MQ=255)
                    tCCTATGACTTTGCGGCTCTGGGTGTTCCAGGCCTGAAGGCTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCAc                                                                                                                       >  8:47518/1‑149 (MQ=255)
                       tATGACTTTGCGGCTCTGGGTGTTCCAGGCCTGAAGGCTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTaca                                                                                                                                       <  6:22878/130‑1 (MQ=255)
                       tATGACTTTGCGGCTCTGGGTGTTCCAGGCCTGAAGGCTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTaca                                                                                                                                       <  6:17826/130‑1 (MQ=255)
                       tATGACTTTGCGGCTCTGGGTGTTCCAGGCCTGAAGGCTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTaca                                                                                                                                       >  5:22878/1‑130 (MQ=255)
                       tATGACTTTGCGGCTCTGGGTGTTCCAGGCCTGAAGGCTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTaca                                                                                                                                       >  5:17826/1‑130 (MQ=255)
                                  ggCTCTGGGTGTTCCAGGCCTGAAGGCTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACTCCCTTCAGGAAGGCACTCTCAAAGGCCTgg                                                                                                         <  1:30426/149‑1 (MQ=255)
                                              tCCAGGCCTGAAGGCTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGCCGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGc                                                                                             <  2:244575/149‑1 (MQ=255)
                                              tCCAGGCCTGAAGGCTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGc                                                                                             >  7:164489/1‑149 (MQ=255)
                                                cAGGCCTGAAGGCTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGGTTCAGCCTGCGt                                                                                           >  6:119691/1‑149 (MQ=255)
                                                        aaGGCTTCGGTGTCCTACCTGCGTGGTGTTGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGgtgt                                                                                   >  1:203197/1‑149 (MQ=255)
                                                         aGGCTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTc                                                                                  <  5:119691/149‑1 (MQ=255)
                                                            cTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCCGCGTCAAGGTGTCTAc                                                                               <  7:47518/149‑1 (MQ=255)
                                                             ttCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCtt                                                                                                     <  1:161344/126‑1 (MQ=255)
                                                             ttCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCtt                                                                                                     >  2:161344/1‑126 (MQ=255)
                                                                      cTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAg                                                                                                                           <  8:291183/95‑1 (MQ=255)
                                                                      cTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAg                                                                                                                           >  7:291183/1‑95 (MQ=255)
                                                                      cTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAg                                                                                       >  8:311153/1‑131 (MQ=255)
                                                                      cTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAg                                                                                       <  7:311153/131‑1 (MQ=255)
                                                                          cTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAgcggcg                                                                  >  8:181554/1‑148 (MQ=255)
                                                                                   gatgatGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCt                                                        >  6:173007/1‑149 (MQ=255)
                                                                                    atgatgTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTg                                                       >  6:260734/1‑149 (MQ=255)
                                                                                         gTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACa                                                   >  3:53156/1‑148 (MQ=255)
                                                                                             gCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGa                                               <  2:124165/148‑1 (MQ=255)
                                                                                               aaCAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGAc                                              >  3:218783/1‑147 (MQ=255)
                                                                                                     tCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCg                                      >  4:79480/1‑149 (MQ=255)
                                                                                                       gggTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGc                                     <  5:173007/148‑1 (MQ=255)
                                                                                                       gggTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGc                                     <  5:260734/148‑1 (MQ=255)
                                                                                                            aCGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAgcgg                                                                    <  1:302445/112‑1 (MQ=255)
                                                                                                            aCGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAgcgg                                                                    >  2:302445/1‑112 (MQ=255)
                                                                                                            aCGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATc                               <  2:203197/149‑1 (MQ=255)
                                                                                                              ggCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGt                             >  7:186867/1‑149 (MQ=255)
                                                                                                                      cAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTaca                     >  6:218463/1‑149 (MQ=255)
                                                                                                                         cGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTc                                  >  5:347003/1‑133 (MQ=255)
                                                                                                                         cGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTc                                  <  6:347003/133‑1 (MQ=255)
                                                                                                                         cGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACtt                  <  4:53156/149‑1 (MQ=255)
                                                                                                                                 aaCGTGACGGCCGCGTCGACTACACCCTTCAGGAAGGCCCTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAgg                                                                                      >  5:366006/1‑73 (MQ=255)
                                                                                                                                 aaCGTGACGGCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAgg                                                                                      <  6:361059/73‑1 (MQ=255)
                                                                                                                                 aaCGTGACGGCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAgg                                                                                      <  6:366006/73‑1 (MQ=255)
                                                                                                                                 aaCGTGACGGCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAgg                                                                                      >  5:361059/1‑73 (MQ=255)
                                                                                                                                 aaCGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTc          <  4:218783/149‑1 (MQ=255)
                                                                                                                                    gTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGc                                     >  1:264749/1‑119 (MQ=255)
                                                                                                                                    gTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGc                                     <  2:264749/119‑1 (MQ=255)
                                                                                                                                    gTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGc                                     <  2:259686/119‑1 (MQ=255)
                                                                                                                                    gTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCATGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGc                                     >  1:259686/1‑119 (MQ=255)
                                                                                                                                    gTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCAATCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGAAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACATCCTg       >  6:228485/1‑149 (MQ=255)
                                                                                                                                       aCGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGt      >  2:208305/1‑147 (MQ=255)
                                                                                                                                       aCGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTaa    >  7:341100/1‑149 (MQ=255)
                                                                                                                                       aCGCCAGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGCCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGt      >  8:47465/1‑147 (MQ=255)
                                                                                                                                         gCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGGCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGc  <  8:164489/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
ATACAACCTTTGGTCAGTACTCCTATGACTTTGCGGCTCTGGGTGTTCCAGGCCTGAAGGCTTCGGTGTCCTACCTGCGTGGTGATGATGTTCGCAACAGGTCGGGTGACGGCACCTACAACGAATGGGAACGTGACGCCCGCGTCGACTACACCCTTCAGGAAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGC  >  NC_002947/289633‑289918

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