Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,226,119 C→G 11.4% intergenic (‑39/‑481) minC ← / → ycgJ inhibitor of FtsZ ring polymerization/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,226,1190CG11.4% 232.0 / 17.6 70intergenic (‑39/‑481)minC/ycgJinhibitor of FtsZ ring polymerization/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (31/31);  new base G (4/4);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

TTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAAT  >  NC_000913/1225976‑1226246
                                                                                                                                               |                                                                                                                               
tttGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCga                                                                                                                             >  2:721052/1‑148 (MQ=255)
    tCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCgatgat                                                                                                                         <  1:721052/148‑1 (MQ=255)
                           cTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGcgcgcg                                                                                                                              <  4:937157/120‑1 (MQ=255)
                           cTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGcgcgcg                                                                                                                              >  3:937157/1‑120 (MQ=255)
                             tAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTta                                                                                                >  4:715176/1‑148 (MQ=255)
                               ggTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatgat                                                                                                                         >  3:367017/1‑121 (MQ=255)
                               ggTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatgat                                                                                                                         <  4:367017/121‑1 (MQ=255)
                                      ccTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGtt                                                                                                                      <  2:96481/117‑1 (MQ=255)
                                      ccTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGtt                                                                                                                      >  1:96481/1‑117 (MQ=255)
                                        tCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatgatt                                                                                                                        <  4:805134/113‑2 (MQ=255)
                                        tCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatgatt                                                                                                                        >  3:805134/1‑112 (MQ=255)
                                          aTGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATtat                                                                                   <  3:715176/148‑1 (MQ=255)
                                                            tAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCt                                                                 >  3:624410/1‑148 (MQ=255)
                                                             aaaGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCga                                                                                                                             >  3:76958/1‑87 (MQ=255)
                                                             aaaGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCga                                                                                                                             <  4:76958/87‑1 (MQ=255)
                                                                 tGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCga                                                                                                                             <  4:777321/83‑1 (MQ=255)
                                                                 tGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCga                                                                                                                             >  3:777321/1‑83 (MQ=255)
                                                                   aaGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCAttt                                                          >  1:781742/1‑148 (MQ=255)
                                                                      cTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATTCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGAGCGCGATGGTGTTCCGa                                                                                                                  <  3:795360/89‑1 (MQ=255)
                                                                      cTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGa                                                                                                                  >  4:795360/1‑89 (MQ=255)
                                                                      cTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTAtt                                                       <  2:781742/148‑1 (MQ=255)
                                                                       tACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTc                                                      <  4:624410/148‑1 (MQ=255)
                                                                          tGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAg                                                                             <  2:775254/122‑1 (MQ=255)
                                                                          tGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAg                                                                             >  1:775254/1‑122 (MQ=255)
                                                                             cTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGatgatg                                                                                                                        <  2:331604/76‑1 (MQ=255)
                                                                             cTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGatgatg                                                                                                                        >  1:331604/1‑76 (MQ=255)
                                                                                        aTTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatga                                                                                                                          >  4:1095758/1‑63 (MQ=37)
                                                                                        aTTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatga                                                                                                                          <  3:1095758/63‑1 (MQ=37)
                                                                                        aTTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACtata                                                                                                          >  1:94895/1‑79 (MQ=255)
                                                                                        aTTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACtata                                                                                                          <  2:94895/79‑1 (MQ=255)
                                                                                         ttGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGa                                                                                                                  >  1:124997/1‑70 (MQ=255)
                                                                                         ttGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGa                                                                                                                  <  2:124997/70‑1 (MQ=255)
                                                                                          tGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGa                                                                                       <  4:615188/96‑1 (MQ=255)
                                                                                          tGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGa                                                                                       >  3:615188/1‑96 (MQ=255)
                                                                                             cGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatgat                                                                                                                         >  3:610747/1‑59 (MQ=25)
                                                                                             cGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatgat                                                                                                                         <  4:610747/59‑1 (MQ=25)
                                                                                             cGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCtg                                   <  3:1113418/145‑1 (MQ=255)
                                                                                             cGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCtg                                   >  4:1113418/1‑145 (MQ=255)
                                                                                                gtTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCa                                                                                                     >  3:561948/1‑76 (MQ=255)
                                                                                                gtTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCa                                                                                                     <  4:561948/76‑1 (MQ=255)
                                                                                                  ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGtt                                                                                                 <  4:97830/78‑1 (MQ=255)
                                                                                                  ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGtt                                                                                                 >  3:97830/1‑78 (MQ=255)
                                                                                                  ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCa                                                                          >  4:929203/1‑101 (MQ=255)
                                                                                                  ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCa                                                                          <  3:929203/101‑1 (MQ=255)
                                                                                                  ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAcc                                                                   <  2:521596/108‑1 (MQ=255)
                                                                                                  ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAcc                                                                   >  1:1029854/1‑108 (MQ=255)
                                                                                                  ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAcc                                                                   >  1:521596/1‑108 (MQ=255)
                                                                                                  ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAcc                                                                   <  2:1029854/108‑1 (MQ=255)
                                                                                                           tGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTc                                                                                           >  4:66432/1‑75 (MQ=255)
                                                                                                           tGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTc                                                                                           >  2:931417/1‑75 (MQ=255)
                                                                                                           tGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTc                                                                                           <  1:931417/75‑1 (MQ=255)
                                                                                                           tGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTc                                                                                           <  3:66432/75‑1 (MQ=255)
                                                                                                             gCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTa                                                         <  2:437302/107‑1 (MQ=255)
                                                                                                             gCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTa                                                         >  1:437302/1‑107 (MQ=255)
                                                                                                                 tACTCAATTAGCTATTAATCATCGCCAGGGCGCGATGATGTTCCGAAGACTATAAg                                                                                                        <  4:628240/56‑1 (MQ=255)
                                                                                                                 tACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAg                                                                                                        >  3:628240/1‑56 (MQ=255)
                                                                                                                 tACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATgg                         <  3:1021190/135‑1 (MQ=255)
                                                                                                                 tACTCAATTAGCTATTAATCATCGCCAGCGCGCGAGGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATgg                         >  4:1021190/1‑135 (MQ=255)
                                                                                                                    tCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTc                                                      >  3:293148/1‑103 (MQ=255)
                                                                                                                    tCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTc                                                      <  4:293148/103‑1 (MQ=255)
                                                                                                                            cTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCtg                                   >  3:955883/1‑114 (MQ=255)
                                                                                                                            cTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCtg                                   <  4:955883/114‑1 (MQ=255)
                                                                                                                               ttAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCa               <  4:920458/131‑1 (MQ=255)
                                                                                                                               ttAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCa               >  3:920458/1‑131 (MQ=255)
                                                                                                                                   tcatcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCaa                                                                         >  3:754403/1‑69 (MQ=255)
                                                                                                                                   tcatcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCaa                                                                         <  4:754403/69‑1 (MQ=255)
                                                                                                                                     atcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCt                                    <  4:1112830/104‑1 (MQ=255)
                                                                                                                                     atcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCt                                    >  3:1112830/1‑104 (MQ=255)
                                                                                                                                      tcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAAt                           <  1:50943/112‑1 (MQ=255)
                                                                                                                                      tcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAAt                           >  2:50943/1‑112 (MQ=255)
                                                                                                                                          cAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCag                                                    <  1:772514/83‑1 (MQ=255)
                                                                                                                                          cAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCag                                                    >  2:772514/1‑83 (MQ=255)
                                                                                                                                            gcgcgcGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAAt  >  2:424557/1‑131 (MQ=255)
                                                                                                                                            gcgcgcGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAAt  <  1:424557/131‑1 (MQ=255)
                                                                                                                                               |                                                                                                                               
TTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAAT  >  NC_000913/1225976‑1226246

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: