Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 325,128 G→A 8.9% V138V (GTG→GTA)  oprQ → outer‑membrane porin D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947325,1280GA 8.9% 157.2 / 4.8 57V138V (GTG→GTA) oprQouter‑membrane porin D
Reads supporting (aligned to +/- strand):  ref base G (25/26);  new base A (3/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

GATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAA  >  NC_002947/324985‑325272
                                                                                                                                               |                                                                                                                                                
gATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTg                                                                                                                                               >  2:243380/1‑147 (MQ=255)
gATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGaa                                                                                                                                             <  2:41520/149‑1 (MQ=255)
gATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGaa                                                                                                                                             <  2:306733/149‑1 (MQ=255)
       gcAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGt                                                                                                                                      <  1:136371/149‑1 (MQ=255)
       gcAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGt                                                                                                                                      >  7:55628/1‑149 (MQ=255)
        cAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTg                                                                                                                                     <  1:243380/149‑1 (MQ=255)
          aaGGCCACAACGGTGGTGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGAc                                                                                                                                   <  7:245411/149‑1 (MQ=255)
               cacaACGGTGGTGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGAt                                                                                                                              >  4:48154/1‑149 (MQ=255)
                acaACGGTGGTGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATg                                                                                                                             >  7:326669/1‑149 (MQ=255)
                 caACGGTGGTGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTACTGAAGTACGGTGACCAGATGc                                                                                                                            <  8:294104/149‑1 (MQ=255)
                      gtggtgGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCa                                                                                                                       <  5:161906/149‑1 (MQ=255)
                          tggCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGc                                                                                                                   <  7:308822/149‑1 (MQ=255)
                              ggcgTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGt                                                                                                               <  5:354034/149‑1 (MQ=255)
                               gcgTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGtg                                                                                                              <  8:55628/149‑1 (MQ=255)
                                     aCTTCTTCAAGCCGCAGCAAACCACCAATGCCCAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGc                                                                                                                   <  2:193883/138‑1 (MQ=255)
                                     aCTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGc                                                                                                                   >  1:193883/1‑138 (MQ=255)
                                     aCTTCTTCAACCCGCACGAAACAACCAATGCCGAAGGCAACGCCAACTCGCCGCACAACATGTCGCGTGGTGGTGATGCTGTGAAGATCCCCATTTCCACCACCGTGCTGAAGTAAGGGGACAAGATGCCAGCACTGCCAGtactg                                                                                                           >  3:139616/1‑146 (MQ=255)
                                               gCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGaa                                                                                                                                             >  6:246078/1‑102 (MQ=255)
                                               gCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGaa                                                                                                                                             <  5:246078/102‑1 (MQ=255)
                                                  gCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCAc                                                                                           >  8:121352/1‑149 (MQ=255)
                                                       aaaCCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTctg                                                                                      >  7:93808/1‑149 (MQ=255)
                                                        aaCCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTctgc                                                                                     >  2:173216/1‑149 (MQ=255)
                                                         accaccAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTctgct                                                                                    <  8:326669/149‑1 (MQ=255)
                                                              caATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCtgctgc                                                                                  >  6:341846/1‑146 (MQ=255)
                                                               aaTGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAg                                                                               >  8:12657/1‑148 (MQ=255)
                                                                aTGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGa                                                                              >  5:12864/1‑148 (MQ=255)
                                                                 tGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGaaa                                                                            <  8:93808/149‑1 (MQ=255)
                                                                     gAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCtt                                                                        <  1:173216/149‑1 (MQ=255)
                                                                      aaGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTc                                                                       <  6:12864/149‑1 (MQ=255)
                                                                           aaCGCCTGCTCGCCGCACAACCTGGCGGGTGTTGTTGCTGCTGTGAAGTTCCGCATTTCCAACAACGTTCTTAAGTAGGTTGATGAGATGAGAGCACTGTCTGTGCTGCAGTACGACGAGGCAAGTCTTGTGGCAGAAATCTTGAcggt                                                                  <  4:139616/149‑4 (MQ=255)
                                                                           aaCGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATCCAAGTCTGCTGCCAGAAAGCTTCACCgg                                                                  >  6:240134/1‑149 (MQ=255)
                                                                                                cTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCAACGGCACCCTGATCACCTCCAAaga                                             <  7:121352/149‑1 (MQ=255)
                                                                                                          gtggtgCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCa                                                       <  1:59417/129‑1 (MQ=255)
                                                                                                          gtggtgCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCa                                                       >  2:59417/1‑129 (MQ=255)
                                                                                                          gtggtgCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGc                                   >  6:176264/1‑149 (MQ=255)
                                                                                                               gctgctGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTg                                <  4:92041/147‑1 (MQ=255)
                                                                                                               gctgctGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATCCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGa                              >  3:92041/1‑149 (MQ=255)
                                                                                                                ctgctgTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGaa                             <  3:48154/149‑1 (MQ=255)
                                                                                                                  gctgTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCAcc                                                     <  7:162037/123‑1 (MQ=255)
                                                                                                                  gctgTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCAcc                                                     >  8:162037/1‑123 (MQ=255)
                                                                                                                   ctgTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACt                           >  6:332213/1‑148 (MQ=255)
                                                                                                                    tgtgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCtt                                                                        <  1:185489/102‑1 (MQ=255)
                                                                                                                    tgtgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCtt                                                                        >  2:185489/1‑102 (MQ=255)
                                                                                                                      tgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAgaga                                           <  7:244672/129‑1 (MQ=255)
                                                                                                                      tgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAgaga                                           >  8:244672/1‑129 (MQ=255)
                                                                                                                      tgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAAc                       <  2:229911/149‑1 (MQ=255)
                                                                                                                      tgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAAc                       >  1:229911/1‑149 (MQ=255)
                                                                                                                               gCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAg                                      >  7:190493/1‑125 (MQ=255)
                                                                                                                               gCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAg                                      <  8:190493/125‑1 (MQ=255)
                                                                                                                                      ggtcacCGTACTTCAGCACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACt                           >  8:144043/4‑129 (MQ=255)
                                                                                                                                      ggtcacCGTACTTCAGCACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACt                           <  7:144043/126‑1 (MQ=255)
                                                                                                                                       gtcacCGTACTTCAGCACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCAcc        >  8:14343/3‑147 (MQ=255)
                                                                                                                                       aaCACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCAcc        >  3:85813/1‑147 (MQ=255)
                                                                                                                                        tcacCGTACTTCAGCACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAg     >  2:215142/2‑149 (MQ=255)
                                                                                                                                        acacCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCAc                                                                     <  8:107344/85‑1 (MQ=255)
                                                                                                                                        acacCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCAc                                                                     >  7:107344/1‑85 (MQ=255)
                                                                                                                                           ccGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGaa  <  4:356164/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAA  >  NC_002947/324985‑325272

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