Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,917,857	G→A	S116S (TCG→TCA)	yebT →	MCE domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,917,857	0	G	A	97.0%	222.0 / ‑2.9	67	S116S (TCG→TCA)	yebT	MCE domain protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (31/34);  total (32/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.99e-01

ATCTTCGTAAGATTGAAGTCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCGTTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACT  >  NC_000913/1917760‑1917952                                                                                                  |                                                                                                aTCTTCGTAAGATTGAAGTCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCAtt                                                                                               <  2:1180918/100‑1 (MQ=255)  tCTTCGTAAGATTGAAGTCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTg                                                                                              <  2:346359/100‑1 (MQ=255)  tCTTCGTAAGATTGAAGTCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTg                                                                                              <  1:1635501/100‑1 (MQ=255)     tCGTAAGATTGAAGTCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCa                                                                                           <  1:1560108/100‑1 (MQ=255)      cGTAAGATTGAAGTCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAg                                                                                          <  2:1560079/100‑1 (MQ=255)       gTAAGATTGAAGTCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAgg                                                                                         >  1:1065009/1‑100 (MQ=255)       gTAAGATTGAAGTCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAgg                                                                                         <  1:110512/100‑1 (MQ=255)       gTAAGATTGAAGTCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAgg                                                                                         <  2:1360518/100‑1 (MQ=255)        tAAGATTGAAGTCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGgt                                                                                        >  1:1199781/1‑100 (MQ=255)        tAAGATTGAAGTCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGgt                                                                                        <  1:1604468/100‑1 (MQ=255)         aaGATTGAAGTCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGgtg                                                                                       <  2:1079521/100‑1 (MQ=255)         aaGATTGAAGTCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGgtg                                                                                       >  1:226361/1‑100 (MQ=255)           gATTGAAGTCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGtc                                                                                     >  1:925501/1‑100 (MQ=255)               gAAGTCCAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCg                                                                                 >  1:403840/1‑100 (MQ=255)                 aGTCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCggg                                                                               >  1:369340/1‑100 (MQ=255)                  gTCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGc                                                                              >  1:931094/1‑100 (MQ=255)                   tCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCt                                                                             >  1:1812528/1‑100 (MQ=255)                   tCAAGGGAAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCt                                                                             <  2:1661325/100‑1 (MQ=255)                    cAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTg                                                                            <  1:1773330/100‑1 (MQ=255)                    cAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTg                                                                            <  2:1028801/100‑1 (MQ=255)                    cAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTg                                                                            >  2:621353/1‑100 (MQ=255)                     aaGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTgg                                                                           <  1:752051/100‑1 (MQ=255)                        gTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACg                                                                        <  2:1647649/100‑1 (MQ=255)                         tCAGAATTAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGa                                                                          <  1:203070/97‑1 (MQ=255)                          cAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGcc                                                                      <  2:696091/100‑1 (MQ=255)                            gCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCt                                                                    >  1:434609/1‑100 (MQ=255)                            gCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCt                                                                    <  1:44805/100‑1 (MQ=255)                               tCAAGTCCGATATGAAAGATGCGCTGTGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCtcg                                                                  <  2:740334/99‑1 (MQ=255)                                 aaGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCtcgacg                                                               >  2:1388169/1‑100 (MQ=255)                                  aGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCgg                                                              <  2:494415/100‑1 (MQ=255)                                     ccGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTgg                                                           >  1:631055/1‑100 (MQ=255)                                      cGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCGTTGGCAGATGTCTCCGGGCTGGACGCCCTCGTCGGAgg                                                           >  2:1441128/1‑99 (MQ=255)                                       gATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGa                                                         >  2:1820253/1‑100 (MQ=255)                                        atatGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTATCCGGGCTGGACGCCCTCGTCGGTGGGaa                                                        >  2:1647122/1‑100 (MQ=255)                                         tatGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAAc                                                       <  1:246404/100‑1 (MQ=255)                                            gAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGaa                                                        >  1:1144147/1‑96 (MQ=255)                                             aaaGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACtata                                                   <  1:312007/100‑1 (MQ=255)                                              aaGCTGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTatat                                                  >  1:1628962/1‑100 (MQ=255)                                              aaGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTatat                                                  >  1:405801/1‑100 (MQ=255)                                                    cgctgcgcGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCat                                            >  1:1717779/1‑100 (MQ=255)                                                           cgAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGTCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCg                                     >  2:1302824/1‑100 (MQ=255)                                                                 gaCTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGggtaaa                               >  1:1447920/1‑100 (MQ=255)                                                                   cTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGggtaaagg                             >  2:599301/1‑100 (MQ=255)                                                                   cTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGggtaaagg                             >  1:1298570/1‑100 (MQ=255)                                                                   cTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGggtaaagg                             <  1:1790852/100‑1 (MQ=255)                                                                   cTCAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGggtaaagg                             <  2:671987/100‑1 (MQ=255)                                                                     cAGTTCTGGCTGGTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGggtaaaggta                           <  2:1667075/100‑1 (MQ=255)                                                                          ctggctggTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAagag                      <  1:1820240/100‑1 (MQ=255)                                                                                ggTGACGCCAAAAGCATCATTGGCAGGTGTTTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGa                 <  1:1548059/99‑1 (MQ=255)                                                                                ggTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGAt                >  1:823302/1‑100 (MQ=255)                                                                                 gTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGATGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATc               <  2:771538/100‑1 (MQ=255)                                                                                 gTGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATc               <  2:772059/100‑1 (MQ=255)                                                                                  tGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCa              >  2:661824/1‑100 (MQ=255)                                                                                  tGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCa              >  1:729495/1‑100 (MQ=255)                                                                                  tGACGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCa              >  2:173847/1‑100 (MQ=255)                                                                                   gACGCCAAAAGCATCATTGGCAGGTGTCTACGGGCTGGACGCCCTCGTCGGTGGGAACTTTATAGGCATGATGACGGGTAAAGGTAAAGAGCAgg                  >  2:1517780/1‑95 (MQ=255)                                                                                    aCGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACt            <  1:111819/100‑1 (MQ=255)                                                                                    aCGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATACCGGGTAAAGGAAAAAAACAGGAACACt            >  2:533121/1‑100 (MQ=255)                                                                                     cGCCAAAAGCATCGTTGGCAGGTGTCCCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACtt           <  1:1351533/100‑1 (MQ=255)                                                                                     cGCCAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACtt           >  1:557800/1‑100 (MQ=255)                                                                                        cAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGt        <  1:101033/100‑1 (MQ=255)                                                                                        cAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGt        >  1:1600646/1‑100 (MQ=255)                                                                                        cAAAAGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGt        >  1:1281069/1‑100 (MQ=255)                                                                                         aaaaGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTc       <  2:1199792/100‑1 (MQ=255)                                                                                          aaaGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCg      <  2:686706/100‑1 (MQ=255)                                                                                          aaaGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCg      <  1:1149990/100‑1 (MQ=255)                                                                                          aaaGCATCATTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCg      <  2:824789/100‑1 (MQ=255)                                                                                              catcatTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACt  >  1:500082/1‑100 (MQ=255)                                                                                                  |                                                                                                ATCTTCGTAAGATTGAAGTCAAGGTCAGCATCAAGTCCGATATGAAAGATGCGCTGCGCGAAGAGACTCAGTTCTGGCTGGTGACGCCAAAAGCATCGTTGGCAGGTGTCTCCGGGCTGGACGCCCTCGTCGGTGGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACT  >  NC_000913/1917760‑1917952

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