Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,900,671 C→G 8.0% T305T (ACG→ACC) ‡ yoaE ← putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,900,6710CG 8.0% 164.2 / 5.1 50T305T (ACG→ACC) ‡yoaEputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (26/20);  new base G (2/2);  total (28/22)
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.00e-01

GCAGTTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCC  >  NC_000913/1900526‑1900818
                                                                                                                                                 |                                                                                                                                                   
gcagTTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCa                                                                                                                                                 <  1:746174/150‑1 (MQ=255)
gcagTTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCa                                                                                                                                                 >  1:1336222/1‑150 (MQ=255)
  agTTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATg                                                                                                                                               <  1:724630/150‑1 (MQ=255)
   gTTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGa                                                                                                                                              >  2:573213/1‑150 (MQ=255)
           gCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGGATCATGACGCCgcgc                                                                                                                                      >  2:859181/1‑150 (MQ=255)
           gCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCgcgc                                                                                                                                      >  1:152198/1‑150 (MQ=255)
              cGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCaga                                                                                                                                   <  2:81880/150‑1 (MQ=255)
                tACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGGATCATGACGCCgcgcggtga                                                                                                                                 <  2:187309/150‑6 (MQ=255)
                      aCCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGGATCATGACGCCGCGCGGTGaaataag                                                                                                                           <  1:859181/150‑6 (MQ=255)
                       ccGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCCCCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGAc                                                                                                                          >  2:375096/1‑150 (MQ=255)
                              tttCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCg                                                                                                                   >  1:253622/1‑150 (MQ=255)
                                     cAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTAc                                                                                                            >  1:1387826/1‑150 (MQ=255)
                                            ccGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTa                                                                                                     >  2:1011186/1‑150 (MQ=255)
                                            ccGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTa                                                                                                     >  1:34715/1‑150 (MQ=255)
                                             cgcgACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTaa                                                                                                    <  2:977604/150‑1 (MQ=255)
                                                gACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCa                                                                                                 <  2:173390/150‑1 (MQ=255)
                                                   aTACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGt                                                                                              >  1:97798/1‑150 (MQ=255)
                                                             cAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGttcttc                                                                                    >  1:556170/1‑150 (MQ=255)
                                                                aCTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCTCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGttcttcttc                                                                                 <  1:667211/150‑1 (MQ=255)
                                                                aCTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGttcttcttc                                                                                 <  2:411710/150‑1 (MQ=255)
                                                                   gtgCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGc                                                                              >  1:906407/1‑150 (MQ=255)
                                                                   gtgCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGc                                                                              >  1:209653/1‑150 (MQ=255)
                                                                     gCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCaa                                                                            >  1:295031/1‑150 (MQ=255)
                                                                       ggTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAAt                                                                          <  2:152198/150‑1 (MQ=255)
                                                                       ggTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAAt                                                                          >  1:16957/1‑150 (MQ=255)
                                                                       ggTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAAt                                                                          >  1:1212868/1‑150 (MQ=255)
                                                                          gAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCa                                                                       >  2:1201233/1‑150 (MQ=255)
                                                                          gAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCa                                                                       >  1:1156112/1‑150 (MQ=255)
                                                                                  cAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCgg                                                               <  2:271962/150‑1 (MQ=255)
                                                                                          cgcgGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGatcggca                                                       >  1:67790/1‑150 (MQ=255)
                                                                                            cgGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGatcggcatc                                                     <  2:16957/150‑1 (MQ=255)
                                                                                             ggATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGatcggcatcg                                                    <  2:1387826/150‑1 (MQ=255)
                                                                                              gATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGatcggcatcgg                                                   <  2:1336222/150‑1 (MQ=255)
                                                                                                       gACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTATGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGtt                                          >  1:1128300/1‑150 (MQ=255)
                                                                                                       gACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGtt                                          >  1:1045500/1‑150 (MQ=255)
                                                                                                          cccGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATc                                       <  2:295031/150‑1 (MQ=255)
                                                                                                          cccGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATc                                       <  1:1340735/150‑1 (MQ=255)
                                                                                                          cccGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGACATCGGCATCGGGTTATc                                       >  1:385559/1‑150 (MQ=255)
                                                                                                                   aGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTg                              >  2:790258/1‑150 (MQ=255)
                                                                                                                     cGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCt                            >  2:347049/1‑150 (MQ=255)
                                                                                                                       tCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGa                          <  2:97798/150‑1 (MQ=255)
                                                                                                                         aaCCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAAc                        <  1:1011186/150‑1 (MQ=255)
                                                                                                                          aCCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACa                       >  1:1075191/1‑150 (MQ=255)
                                                                                                                                cTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGcc                 <  2:1321168/150‑1 (MQ=255)
                                                                                                                                    tttCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGAc             >  2:109747/1‑150 (MQ=255)
                                                                                                                                     ttCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACg            <  2:209653/150‑1 (MQ=255)
                                                                                                                                           gcgcgGGATCATGACGCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTccc      >  1:976510/1‑150 (MQ=255)
                                                                                                                                           gcgcgGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTccc      <  2:556170/150‑1 (MQ=255)
                                                                                                                                           gcgcgGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTccc      >  2:653380/1‑150 (MQ=255)
                                                                                                                                            cgcgGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCg     <  1:1323145/150‑1 (MQ=255)
                                                                                                                                              cggcgTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGcc   <  2:1212868/150‑1 (MQ=255)
                                                                                                                                               ggcgTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGccc  <  2:906407/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCAGTTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCC  >  NC_000913/1900526‑1900818

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