Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,468,304 G→C 7.9% intergenic (‑25/+98) PP_2161 ← / ← PP_2162 hypothetical protein/phosphonoacetaldehyde phosphonohydrolase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,468,3040GC 7.9% 152.8 / 4.6 51intergenic (‑25/+98)PP_2161/PP_2162hypothetical protein/phosphonoacetaldehyde phosphonohydrolase‑like protein
Reads supporting (aligned to +/- strand):  ref base G (24/23);  new base C (2/2);  total (26/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 = 8.49e-01

GTTGCGCCTCGATCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTCCACCAGGCAG  >  NC_002947/2468156‑2468450
                                                                                                                                                    |                                                                                                                                                  
gTTGCGCCTCGATCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGc                                                                                                                                                   <  1:999571/150‑1 (MQ=255)
  tGCGCCTCGATCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCaa                                                                                                                                                 <  2:1145768/150‑1 (MQ=255)
     gcCTCGATCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGcc                                                                                                                                              >  1:686176/1‑150 (MQ=255)
           aTCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGAc                                                                                                                                        <  1:155242/150‑1 (MQ=255)
            tCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGCCAAGCCAttgatg                                                                                                                                       <  2:86429/150‑6 (MQ=255)
              tGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGCCAAGCCAttgatgga                                                                                                                                     <  1:153608/150‑8 (MQ=255)
                     ccTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACg                                                                                                                              >  1:93243/1‑150 (MQ=255)
                      ctctTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGt                                                                                                                             >  2:887260/1‑150 (MQ=255)
                                 gggCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTc                                                                                                                  <  1:887260/150‑1 (MQ=255)
                                                  tcCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGt                                                                                                 <  2:65173/150‑1 (MQ=255)
                                                     gAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTg                                                                                                                <  1:1123754/132‑1 (MQ=255)
                                                     gAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTg                                                                                                                >  2:1123754/1‑132 (MQ=255)
                                                         gCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTTCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTAc                                                                                          <  2:686176/150‑1 (MQ=255)
                                                                      gTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGAttt                                                                             <  2:468579/150‑1 (MQ=255)
                                                                      acccAGTACCTCGCTCCAGCCGTGCGGTCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCAGCTTTTGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGAttt                                                                             <  1:379984/148‑1 (MQ=255)
                                                                         cAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCa                                                                          >  1:1082431/1‑150 (MQ=255)
                                                                                tCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCcg                                                                   >  2:664904/1‑150 (MQ=255)
                                                                                   cgcAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGt                                                                >  2:544922/1‑150 (MQ=255)
                                                                                    gcAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAgg                                                                                                   >  2:936140/1‑114 (MQ=255)
                                                                                    gcAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAgg                                                                                                   <  1:936140/114‑1 (MQ=255)
                                                                                            ttCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGc                                                       >  2:957110/1‑150 (MQ=255)
                                                                                               aGGCGCTCTTGCAATTCGCGGGCAGGAATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGt                                                    >  1:1123315/1‑150 (MQ=255)
                                                                                                        tGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGtt                                           <  1:664904/150‑1 (MQ=255)
                                                                                                         gCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCa                                                <  2:1175248/144‑1 (MQ=255)
                                                                                                         gCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCa                                                >  1:1175248/1‑144 (MQ=255)
                                                                                                           aaTTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAggg                                                                                                  <  1:618780/92‑1 (MQ=255)
                                                                                                           aaTTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAggg                                                                                                  <  1:618409/92‑1 (MQ=255)
                                                                                                           aaTTCGCGGGCAGGCATGGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAggg                                                                                                  >  2:618780/1‑92 (MQ=255)
                                                                                                           aaTTCGCGGGCAGGCATGGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAggg                                                                                                  >  2:618409/1‑92 (MQ=255)
                                                                                                               cgcgGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGTATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCAcc                                    <  1:399921/150‑1 (MQ=255)
                                                                                                               cgcgGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCGTGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCAcc                                    <  1:416567/150‑1 (MQ=255)
                                                                                                                        gcgttGGCTATTCTCCATCAATGGCCTGGCAAGCCAGGGACCTCTGCGACGTGGCAAAGTTTCTGGCCTGTCTTAAGCGTAATCTACCTTCTAGCGATTTGCATGATCCCAGTCAAGACGGCTGTATCAGGGTTTTTCACCGTTGCGCCg                           >  2:496876/4‑150 (MQ=255)
                                                                                                                          aTTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGCGTTTTTCACCCTTGCGCCGgc                         >  2:593761/1‑150 (MQ=255)
                                                                                                                              gCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCaa                                                             <  1:1097075/110‑1 (MQ=255)
                                                                                                                              gCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCaa                                                             >  2:1097075/1‑110 (MQ=255)
                                                                                                                                tATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTggg                   >  1:258430/1‑150 (MQ=255)
                                                                                                                                tATTCTCCATCAATGGCTTGGCAAGACATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTggg                   >  1:1186604/1‑150 (MQ=255)
                                                                                                                                  ttCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCcgc                                                                  <  2:114742/101‑1 (MQ=255)
                                                                                                                                  ttCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCcgc                                                                  >  1:114742/1‑101 (MQ=255)
                                                                                                                                   tctcCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGa                <  2:1123315/150‑1 (MQ=255)
                                                                                                                                     tcCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGt                                            <  2:382715/120‑1 (MQ=255)
                                                                                                                                     tcCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGt                                            >  1:382715/1‑120 (MQ=255)
                                                                                                                                     tcCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGgcgc                       <  2:398740/141‑1 (MQ=255)
                                                                                                                                     tcCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGgcgc                       >  1:398740/1‑141 (MQ=255)
                                                                                                                                     tcCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATg              >  2:157601/1‑150 (MQ=255)
                                                                                                                                        gtccaTGGCTTGCCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTcc           >  2:1116069/5‑150 (MQ=255)
                                                                                                                                            aTGGCTTGGCAAGCCATGGACCTCTGCGACGTGGc                                                                                                                          <  2:766519/35‑1 (MQ=255)
                                                                                                                                            aTGGCTTGGCAAGCCATGGACCTCGACGACGTGGc                                                                                                                          >  1:766519/1‑35 (MQ=255)
                                                                                                                                            aTGGCTTGCCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTccacca       >  1:37541/1‑150 (MQ=255)
                                                                                                                                               gCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTCCACcaggc    >  1:201160/1‑150 (MQ=255)
                                                                                                                                                 ttGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTCCACcaggcag  <  1:957110/150‑1 (MQ=255)
                                                                                                                                                    gCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGAt                                                                      >  1:939944/1‑79 (MQ=255)
                                                                                                                                                    gCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGAt                                                                      <  2:939944/79‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTTGCGCCTCGATCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTCCACCAGGCAG  >  NC_002947/2468156‑2468450

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