Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,226,114	C→G	11.4%	intergenic (‑34/‑486)	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_000913	1,226,114	0	C	G	11.4%	235.5 / 13.0	70	intergenic (‑34/‑486)	minC/ycgJ	inhibitor 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 = 1.00e+00

TTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCA  >  NC_000913/1225976‑1226233                                                                                                                                           |                                                                                                                        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                                                                                                            <  4:367017/121‑1 (MQ=255)                                ggTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatgat                                                                                                            >  3:367017/1‑121 (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                                                                                                                >  3:777321/1‑83 (MQ=255)                                                                  tGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCga                                                                                                                <  4:777321/83‑1 (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                                                                                                           >  1:331604/1‑76 (MQ=255)                                                                              cTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGatgatg                                                                                                           <  2:331604/76‑1 (MQ=255)                                                                                         aTTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatga                                                                                                             <  3:1095758/63‑1 (MQ=37)                                                                                         aTTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatga                                                                                                             >  4:1095758/1‑63 (MQ=37)                                                                                         aTTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACtata                                                                                             <  2:94895/79‑1 (MQ=255)                                                                                         aTTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACtata                                                                                             >  1:94895/1‑79 (MQ=255)                                                                                          ttGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGa                                                                                                     >  1:124997/1‑70 (MQ=255)                                                                                          ttGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGa                                                                                                     <  2:124997/70‑1 (MQ=255)                                                                                           tGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGa                                                                          >  3:615188/1‑96 (MQ=255)                                                                                           tGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGa                                                                          <  4:615188/96‑1 (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                                                                                    >  3:97830/1‑78 (MQ=255)                                                                                                   ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGtt                                                                                    <  4:97830/78‑1 (MQ=255)                                                                                                   ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCa                                                             >  4:929203/1‑101 (MQ=255)                                                                                                   ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCa                                                             <  3:929203/101‑1 (MQ=255)                                                                                                   ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAcc                                                      >  1:1029854/1‑108 (MQ=255)                                                                                                   ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAcc                                                      >  1:521596/1‑108 (MQ=255)                                                                                                   ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAcc                                                      <  2:1029854/108‑1 (MQ=255)                                                                                                   ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAcc                                                      <  2:521596/108‑1 (MQ=255)                                                                                                            tGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTc                                                                              <  3:66432/75‑1 (MQ=255)                                                                                                            tGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTc                                                                              <  1:931417/75‑1 (MQ=255)                                                                                                            tGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTc                                                                              >  4:66432/1‑75 (MQ=255)                                                                                                            tGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTc                                                                              >  2:931417/1‑75 (MQ=255)                                                                                                              gCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTa                                            >  1:437302/1‑107 (MQ=255)                                                                                                              gCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTa                                            <  2:437302/107‑1 (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                                         <  4:293148/103‑1 (MQ=255)                                                                                                                     tCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTc                                         >  3:293148/1‑103 (MQ=255)                                                                                                                             cTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCtg                      <  4:955883/114‑1 (MQ=255)                                                                                                                             cTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCtg                      >  3:955883/1‑114 (MQ=255)                                                                                                                                ttAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCa  >  3:920458/1‑131 (MQ=255)                                                                                                                                ttAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCa  <  4:920458/131‑1 (MQ=255)                                                                                                                                    tcatcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCaa                                                            >  3:754403/1‑69 (MQ=255)                                                                                                                                    tcatcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCaa                                                            <  4:754403/69‑1 (MQ=255)                                                                                                                                      atcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCt                       >  3:1112830/1‑104 (MQ=255)                                                                                                                                      atcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCt                       <  4:1112830/104‑1 (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)                                                                                                                                           |                                                                                                                        TTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCA  >  NC_000913/1225976‑1226233

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: ‑