Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,900,680 C→G 7.3% G302G (GGG→GGC)  yoaE ← putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,900,6800CG 7.3% 182.2 / 4.9 55G302G (GGG→GGC) yoaEputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (28/23);  new base G (2/2);  total (30/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 = 7.08e-01

GCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAG  >  NC_000913/1900537‑1900829
                                                                                                                                               |                                                                                                                                                     
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:209653/1‑150 (MQ=255)
                                                        gtgCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGc                                                                                         >  1:906407/1‑150 (MQ=255)
                                                          gCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCaa                                                                                       >  1:295031/1‑150 (MQ=255)
                                                            ggTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAAt                                                                                     >  1:16957/1‑150 (MQ=255)
                                                            ggTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAAt                                                                                     >  1:1212868/1‑150 (MQ=255)
                                                            ggTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAAt                                                                                     <  2:152198/150‑1 (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)
                                                                                                                                         cATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGttctt                                                                                                >  1:672354/1‑62 (MQ=255)
                                                                                                                                         cATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGttctt                                                                                                <  2:672354/62‑1 (MQ=255)
                                                                                                                                         cATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCaa        <  1:375096/150‑1 (MQ=255)
                                                                                                                                         cATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCGTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGCACATTGGCCTGACGTTTCCCGCCCATCaa        >  1:493220/1‑150 (MQ=255)
                                                                                                                                           tGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAAc      >  2:1219696/1‑150 (MQ=255)
                                                                                                                                             aTCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGc               <  1:402644/139‑1 (MQ=255)
                                                                                                                                             aTCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGc               >  2:402644/1‑139 (MQ=255)
                                                                                                                                               cccGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAg  >  2:48446/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAG  >  NC_000913/1900537‑1900829

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