Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,331,674 A→G 7.4% T110T (ACA→ACG)  PP_3803 → cation ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,331,6740AG 7.4% 150.8 / 2.9 54T110T (ACA→ACG) PP_3803cation ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (28/22);  new base G (2/2);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

GTCGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTATCCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTAC  >  NC_002947/4331525‑4331816
                                                                                                                                                     |                                                                                                                                              
gTCGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTATCCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATaca                                                                                                                                                <  2:1819883/150‑1 (MQ=255)
 tCGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTATCCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCCGTCGCCGGATacac                                                                                                                                               <  1:1036639/150‑1 (MQ=255)
  cGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTATCCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACt                                                                                                                                              <  1:1617948/150‑1 (MQ=255)
      cATGCGGCGCTGACCGGCGTGGCCATCGGTATCCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGAtt                                                                                                                                          <  1:1514928/150‑1 (MQ=255)
         gcggcgCTGACCGGCGTGGCCATCGGTATCCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGc                                                                                                                                       <  1:669889/150‑1 (MQ=255)
         gcggcgCTGACCGGCGTGGCCATCGGTATCCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGc                                                                                                                                       <  1:867304/150‑1 (MQ=255)
                          ggCCATCGGTATCCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTATc                                                                                                                      >  2:9332/1‑150 (MQ=255)
                                   tATCCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCt                                                                                                             <  1:602811/150‑1 (MQ=255)
                                       ctgctgGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGgc                                                                                                         <  2:1063041/150‑1 (MQ=255)
                                           tgGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCa                                                                                                     <  2:1788981/150‑1 (MQ=255)
                                              gTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTACTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGcc                                                                                                  >  1:190901/1‑150 (MQ=255)
                                                   cccTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCctgctg                                                                                             >  2:1766181/1‑150 (MQ=255)
                                                   cccTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCctgctg                                                                                             >  1:414108/1‑150 (MQ=255)
                                                     cTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGCCGCTTGGCGCCAGCctgctgct                                                                                           >  1:1001336/1‑150 (MQ=255)
                                                        taCCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGAt                                                                                        <  1:1426120/150‑1 (MQ=255)
                                                          ccGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCgg                                                                                                                          >  1:1589065/1‑114 (MQ=255)
                                                          ccGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCgg                                                                                                                          <  2:1589065/114‑1 (MQ=255)
                                                                     ggCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAg                                                                           >  1:1029451/1‑150 (MQ=255)
                                                                      gCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATaca                                                                                                                                                <  2:71974/80‑1 (MQ=255)
                                                                      gCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATaca                                                                                                                                                >  1:71974/1‑80 (MQ=255)
                                                                      gCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGa                                                                          <  1:1590777/150‑1 (MQ=255)
                                                                             gTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCa                                                                                                     >  1:402363/1‑116 (MQ=255)
                                                                             gTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCa                                                                                                     <  2:402363/116‑1 (MQ=255)
                                                                              ttCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTg                                                                  >  1:1632788/1‑150 (MQ=255)
                                                                                         cctgcTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCga                                                       >  1:432273/1‑150 (MQ=255)
                                                                                         cctgcTGTTCGGCATCCTGCTCAACTTCCTGCGCAAACGTACCGGGCTGTAGCCGGATACACTGATTGGCCTGTTCCTGTCGGGGTGGCCGGCGCCAGCCACCAGCCTGCTGATGCTGCTGAGATGCAAGGTAAACGGACACAtacaccg                                                       >  2:620046/1‑144 (MQ=255)
                                                                                                ttCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCt                                                                                                                   <  2:551234/83‑1 (MQ=255)
                                                                                                ttCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCt                                                                                                                   >  1:551234/1‑83 (MQ=255)
                                                                                                ttCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTg                                                <  2:190901/150‑1 (MQ=255)
                                                                                                ttCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCCGTCGGTGTCGCTGGCGCTTGGCCCCTGCGTACTACGGAGACGGTCAGGCAAGATCAACGTGCACATCCTCGAGAACGTg                                                >  1:1400394/1‑150 (MQ=255)
                                                                                                ttCGGCATCCTGCTCAACTGCCTGCGCCACCGTCCCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTg                                                <  2:167236/150‑1 (MQ=255)
                                                                                                   ggCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTg                                             >  1:1131447/1‑150 (MQ=255)
                                                                                                     cATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgtt                                           <  1:1766181/150‑1 (MQ=255)
                                                                                                         cTGCTCAACTTCCTGCGCAACCGTACCCGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttcggt                                       >  2:710172/1‑150 (MQ=255)
                                                                                                           gCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACAACCTCGAGAACGTGCTgttcggttc                                     >  2:64396/1‑150 (MQ=255)
                                                                                                              cAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTg                                                                                                                  <  2:1145871/70‑1 (MQ=255)
                                                                                                              cAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTg                                                                                                                  >  1:1145871/1‑70 (MQ=255)
                                                                                                               aaCTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGt                                                                   >  1:1583223/1‑116 (MQ=255)
                                                                                                               aaCTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGt                                                                   <  2:1583223/116‑1 (MQ=255)
                                                                                                               aaCTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTg                                 >  2:863294/1‑150 (MQ=255)
                                                                                                                aCTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGc                                <  1:710172/150‑1 (MQ=255)
                                                                                                                 cTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCt                               <  2:1632788/150‑1 (MQ=255)
                                                                                                                   tCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGa                             >  1:1053558/1‑150 (MQ=255)
                                                                                                                   tCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCATGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCt                                                          <  2:1655998/121‑1 (MQ=255)
                                                                                                                   tCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCATGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCt                                                          >  1:1655998/1‑121 (MQ=255)
                                                                                                                         gcAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTa                       >  1:1448977/1‑150 (MQ=255)
                                                                                                                             ccGTACCGGGCTGTCGCCGGATACACTGATTAGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTg                                                                  <  1:1257438/103‑1 (MQ=255)
                                                                                                                             ccGTACCGGGCTGTCGCCGGATACACTGATTAGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTg                                                                  >  2:1257438/1‑103 (MQ=255)
                                                                                                                                   cGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGa             >  2:333678/1‑150 (MQ=255)
                                                                                                                                    gggCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGAc            >  2:25380/1‑150 (MQ=255)
                                                                                                                                      gCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGTAcca          >  2:77622/1‑148 (MQ=255)
                                                                                                                                        tGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTggcaggc                                                                              >  2:1730629/1‑80 (MQ=255)
                                                                                                                                        tGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTggcaggc                                                                              <  1:1730629/80‑1 (MQ=255)
                                                                                                                                        gatctccGGATACGCTGATTGGCGTGTTCCTGTCGGTGCTGCTGGCGCTTGGCGCCAGCCTACTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGc        <  1:18181/145‑1 (MQ=255)
                                                                                                                                        gatctccGGATACGCTGATTGGCGTGTTCCTGTCGGTATGCTTGGCGCTTGGCGCCAGCCTACTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGc        <  1:18312/145‑1 (MQ=255)
                                                                                                                                             ccGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGAACTTGCTGGTa   >  1:636714/1‑150 (MQ=255)
                                                                                                                                              cGGATACGCTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTACTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTAGAGAACGTGATGTTTGGTTAGGTGGTGGTGGTTAGTGTCGAGTTCTTGCtagatc  >  2:18181/1‑145 (MQ=255)
                                                                                                                                              cGGATACGCTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTACTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTAGAGAACGTACTATTAGGTTCTGTTCATACCTTTAGGGCCCAGtgcttgctagatc  >  2:18312/1‑145 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GTCGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTATCCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTAC  >  NC_002947/4331525‑4331816

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