Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,687,693 C→G 10.2% I223M (ATC→ATG)  PP_4148 → ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,687,6930CG10.2% 172.7 / 9.8 59I223M (ATC→ATG) PP_4148ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (28/25);  new base G (3/3);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

TGTTCGCCATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACGTGT  >  NC_002947/4687548‑4687832
                                                                                                                                                 |                                                                                                                                            
tGTTCGCCATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGAc                                                                                                                                            >  7:123103/1‑148 (MQ=255)
     gCCATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCt                                                                                                                                     >  3:37039/1‑149 (MQ=255)
      ccATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTg                                                                                                                                    >  1:129543/1‑149 (MQ=255)
              ctgATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCg                                                                                                                             <  1:245948/148‑1 (MQ=255)
              ctgATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCg                                                                                                                             <  4:37039/148‑1 (MQ=255)
                 aTCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGACCAGCGC‑GGTGATCGGCag                                                                                                                         <  4:303593/149‑2 (MQ=39)
                 aTCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGACCAGCGC‑GGTGATCGGCag                                                                                                                         <  3:68139/149‑2 (MQ=39)
                         gTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTg                                                                                                                                    <  8:262/130‑1 (MQ=255)
                         gTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTg                                                                                                                                    >  7:262/1‑130 (MQ=255)
                         gTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGcca                                                                                                                 >  5:55140/1‑149 (MQ=255)
                             gCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGccacca                                                                                                              <  3:217743/148‑1 (MQ=255)
                                 gTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGAcc                                                                                                         <  6:55140/149‑1 (MQ=255)
                                           cTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGAcc                                                                                                >  4:241096/1‑148 (MQ=255)
                                           cTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGAc                                                                                                 <  6:150875/147‑1 (MQ=255)
                                           cTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCa                                                                                               <  3:314986/149‑1 (MQ=255)
                                            ttCGAATGGTTCCAGCTGCGTGGCCTCACCTCCAACAAATTCGACGAACTGAGCACCACCGGCAAGTTGCTGGACTACTTATGGAACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGACAACTTCGCACACATGACCATGCTCACCaa                                                                                              >  8:250754/1‑149 (MQ=255)
                                             tCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGAcc                                                                                                >  3:157445/1‑146 (MQ=255)
                                               gACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGa                                                                                            <  4:157445/148‑1 (MQ=255)
                                                aCTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAAc                                                                                          <  2:91156/149‑1 (MQ=255)
                                                  aggTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGCTGCTGCCGATCAcc                                                                                                                                           >  4:275106/2‑99 (MQ=255)
                                                    gTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACATCGCCACCAGGACCCTGCTGCCCAAGAACAGCt                                                                                      >  8:90196/1‑149 (MQ=255)
                                                         cGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTc                                                                                 >  4:277871/1‑149 (MQ=255)
                                                                    cTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATcaa                                                                      <  2:129543/149‑1 (MQ=255)
                                                                     tCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATcaa                                                                      >  8:188499/1‑148 (MQ=255)
                                                                       aCCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATcaaca                                                                    >  4:13716/1‑148 (MQ=255)
                                                                        ccTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCaacaaa                                                                  <  7:188499/149‑1 (MQ=255)
                                                                         cTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCaacaaac                                                                 <  3:277871/149‑1 (MQ=255)
                                                                          tCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACa                                                                                         <  7:268011/124‑1 (MQ=255)
                                                                          tCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACa                                                                                         >  8:268011/1‑124 (MQ=255)
                                                                                  cTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTgg                                                        >  4:118690/1‑149 (MQ=255)
                                                                                      gacgaACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCg                                                                                                                             <  2:115388/76‑1 (MQ=255)
                                                                                      gacgaACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCg                                                                                                                             >  1:115388/1‑76 (MQ=255)
                                                                                         gaACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGt                                                          >  6:178716/1‑140 (MQ=255)
                                                                                         gaACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGt                                                          <  5:178716/140‑1 (MQ=255)
                                                                                          aaCTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgc                                                 >  4:142214/1‑148 (MQ=255)
                                                                                            cTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCt                                                                                      >  7:208388/1‑109 (MQ=255)
                                                                                            cTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCt                                                                                      <  8:208388/109‑1 (MQ=255)
                                                                                            cTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgccaa                                              <  3:13716/149‑1 (MQ=255)
                                                                                               aGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgccaaggc                                           >  6:48409/1‑149 (MQ=255)
                                                                                                   ccaccGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgccaaggccaag                                       >  1:312377/1‑149 (MQ=255)
                                                                                                        ggCAAGGTGCTGGACTACTTCTGGGACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCt                                  <  5:195438/149‑1 (MQ=255)
                                                                                                                 cTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGAt                          >  7:534/1‑148 (MQ=255)
                                                                                                                  tGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCa                        >  3:245126/1‑149 (MQ=255)
                                                                                                                   ggACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGc                              <  7:137259/142‑1 (MQ=255)
                                                                                                                   ggACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGc                              >  8:137259/1‑142 (MQ=255)
                                                                                                                    gACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCAc                       <  4:245126/148‑1 (MQ=255)
                                                                                                                     actactTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACgg                     <  3:241096/149‑1 (MQ=255)
                                                                                                                            cTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACg                      >  3:298031/1‑141 (MQ=255)
                                                                                                                            cTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACg                      <  4:298031/141‑1 (MQ=255)
                                                                                                                             tGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCa                                                                                               <  3:153424/67‑1 (MQ=255)
                                                                                                                             tGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCa                                                                                               >  4:153424/1‑67 (MQ=255)
                                                                                                                                  ccTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGcc        <  5:48409/149‑1 (MQ=255)
                                                                                                                                     gaagttgcCGATGACCAGCGC‑GGTGATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgccaaggccaag                                       <  5:179630/110‑1 (MQ=38)
                                                                                                                                     gaagttgcCGATGACCAGCGC‑GGTGATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgccaaggccaag                                       >  6:179630/6‑115 (MQ=38)
                                                                                                                                     gaagttgcCGATGACCAGCGC‑GGTGATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACg     >  8:293811/6‑149 (MQ=255)
                                                                                                                                      gTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACgt    >  1:136951/1‑149 (MQ=255)
                                                                                                                                        gctgcCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACgtgt  >  7:19599/1‑149 (MQ=255)
                                                                                                                                            ccGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCa                                                                                               >  7:277266/1‑52 (MQ=255)
                                                                                                                                            ccGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCa                                                                                               <  8:277266/52‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
TGTTCGCCATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACGTGT  >  NC_002947/4687548‑4687832

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: ‑