Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,475,315	G→T	V58L (GTG→TTG)	ydbD →	PF10971 family putative periplasmic methylglyoxal resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,475,315	0	G	T	97.0%	218.5 / ‑2.4	67	V58L (GTG→TTG)	ydbD	PF10971 family putative periplasmic methylglyoxal resistance protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (35/30);  total (36/31)
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.77e-01

ACTTCGATGTGTTGCAAATGCAGTTGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCGTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAATATAGAGAGTAATCAAAAATTGAT  >  NC_000913/1475218‑1475412                                                                                                  |                                                                                                  aCTTCGATGTGTTGCAAATGCAGTTGATATTGAATCATTTTTTTCTGCGGCTACCGGTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTg                                                                                                 <  1:650437/100‑1 (MQ=255) aCTTCGATGTGTTGCAAATGCAGTTGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTg                                                                                                 <  2:564513/100‑1 (MQ=255)  cTTCGATGTGTTGCAAATGCAGTTGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGa                                                                                                >  2:460891/1‑100 (MQ=255)  cTTCGATGTGTTGCAAATGCAGTTGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGa                                                                                                >  2:893225/1‑100 (MQ=255)    tCGATGTGTTGCAAATGCAGTTGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAAt                                                                                              >  1:1050489/1‑100 (MQ=255)     cGATGTGTTGCAAATGTAGTTGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCGTGAATc                                                                                             <  1:63829/100‑1 (MQ=255)      ccgatcttGCAAATGCAGTTGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCt                                                                                            <  1:686839/94‑1 (MQ=255)           gtTGCAAATGCAGTTGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTcc                                                                                       >  1:68390/1‑100 (MQ=255)             tGCAAATGCAGTTGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCtagatcgg                                                                                     >  2:686835/1‑93 (MQ=255)             tGCAAATGCAGTTGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTcccc                                                                                     <  1:1212675/100‑1 (MQ=255)              gCAAATGCAGTTGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTccc                                                                                      >  2:58598/1‑98 (MQ=255)                  aTGCAGTTGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCgg                                                                                <  2:403481/100‑1 (MQ=255)                   tGCAGTTGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGt                                                                               <  2:43183/100‑1 (MQ=255)                   tGCAGTTGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGt                                                                               <  1:164014/100‑1 (MQ=255)                    gCAGTTGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGtt                                                                              <  2:1052240/100‑1 (MQ=255)                      aGTTGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTa                                                                            >  1:1278666/1‑100 (MQ=255)                         tGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCt                                                                         >  1:1104822/1‑100 (MQ=255)                         tGATATTGAATCATTTTTTTCCGCGGCTACCGCTGAAGCTAAACAACATGTTGAACAGGCTACCAACACTAGTTAGAATCTTTTCCCCCTCGCTTTAACt                                                                         >  2:1277318/1‑100 (MQ=255)                          gATATTGAAGCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTg                                                                        >  2:894564/1‑100 (MQ=255)                                gAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCg                                                                   >  2:596672/1‑99 (MQ=255)                                gAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGa                                                                  <  1:270952/100‑1 (MQ=255)                                 aaTCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGaa                                                                 >  1:32828/1‑100 (MQ=255)                                    cATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTg                                                              >  2:778091/1‑100 (MQ=255)                                     aTTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTgg                                                             >  2:35361/1‑100 (MQ=255)                                      tttttttCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGa                                                            <  1:1253931/100‑1 (MQ=255)                                      tttttttCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGa                                                            >  2:1110087/1‑100 (MQ=255)                                      tttttttCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGa                                                            >  1:228396/1‑100 (MQ=255)                                       ttttttCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGaa                                                           >  1:1023969/1‑100 (MQ=255)                                         ttttCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAg                                                         <  2:393336/100‑1 (MQ=255)                                             cTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCa                                                     >  1:1194875/1‑100 (MQ=255)                                                 ggCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGt                                                 >  2:629446/1‑100 (MQ=255)                                                   cTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTgg                                               <  2:1003702/100‑1 (MQ=255)                                                    tACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGc                                              <  2:1023960/100‑1 (MQ=255)                                                      ccGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGa                                            <  2:36946/100‑1 (MQ=255)                                                      ccGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGa                                            <  2:1126879/100‑1 (MQ=255)                                                        gCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGAtt                                          >  2:693321/1‑100 (MQ=255)                                                         cTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGAttt                                         >  2:63624/1‑100 (MQ=255)                                                            aaGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGt                                      <  2:32827/100‑1 (MQ=255)                                                              gATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTgg                                    <  2:1104812/100‑1 (MQ=255)                                                              gATAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTgg                                    <  2:302046/100‑1 (MQ=255)                                                               aTAAACAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGt                                   <  2:1291126/100‑1 (MQ=255)                                                                 aaaCAACAAGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAg                                 <  1:778096/100‑1 (MQ=255)                                                                     aacaaGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAgg                             <  1:893232/100‑1 (MQ=255)                                                                      acaaGTTGAACAAGCTATCAACAGTCGCTTGAATCTTGTGGCCAGCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATGTAGTGGTAGAAGGt                            >  2:301200/1‑100 (MQ=255)                                                                      acaaGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGt                            >  1:508790/1‑100 (MQ=255)                                                                      acaaGTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGt                            >  1:1084642/1‑100 (MQ=255)                                                                          gTTGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAata                        <  1:681376/100‑1 (MQ=255)                                                                            tGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAAtata                      >  2:773506/1‑100 (MQ=255)                                                                            tGAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAAtata                      >  2:773365/1‑100 (MQ=255)                                                                             gAACAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAATATag                     >  2:1116547/1‑100 (MQ=255)                                                                              aaCAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAAGTGGAAAGTGCATCGTGGCGATTTAGTAGTAGAAGGTAAtata                      <  2:841872/98‑1 (MQ=255)                                                                              aaCAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAATATaga                    >  2:1113924/1‑100 (MQ=255)                                                                               aCAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAATATagag                   <  2:1084041/100‑1 (MQ=255)                                                                                cAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAATATagaga                  >  2:551875/1‑100 (MQ=255)                                                                                cAAGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAATATagaga                  >  2:827991/1‑100 (MQ=255)                                                                                 aaGCTATCAACTGTAGCTTGAATCTTGTCCCCTTCGGGTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAATATagagag                 <  2:265520/100‑1 (MQ=255)                                                                                  aGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAATATagaga                  >  2:134794/1‑98 (MQ=255)                                                                                  aGCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAATATAGAGAGt                <  1:1048349/100‑1 (MQ=255)                                                                                   gCTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAATATAGAGAGTa               <  1:1002323/100‑1 (MQ=255)                                                                                    cTATCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAATATAGAGAGTaa              >  2:903555/1‑100 (MQ=255)                                                                                      aTCAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCTATTTAGTGGTAGAAGGTAATATAGAGAGTAATc            >  2:546521/1‑100 (MQ=255)                                                                                        cAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAATATAGAGAGTAATCaa          >  1:1006558/1‑100 (MQ=255)                                                                                        cAACAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATCTAGTGGTAGAAGGTAATATAGAGAGTAATaaa          >  2:103254/1‑100 (MQ=255)                                                                                          aCAGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAATATAGAGAGTAATCaaaa        <  1:460895/100‑1 (MQ=255)                                                                                            aGTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAATATAGAGAGTAATCAAAAat      <  2:68389/100‑1 (MQ=255)                                                                                             gTAGCTTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAATATAGAGAGTAATCAAAAatt     <  2:678975/100‑1 (MQ=255)                                                                                                gCGTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAATATAGAGAGTAATCAAAAattgat  >  1:834797/1‑100 (MQ=255)                                                                                                  |                                                                                                  ACTTCGATGTGTTGCAAATGCAGTTGATATTGAATCATTTTTTTCTGCGGCTACCGCTGAAGATAAACAACAAGTTGAACAAGCTATCAACAGTAGCGTGAATCTTGTCCCCTTCGGTTTATCTGCATCGAACTGGAAAGTGCATCGTGGCGATTTAGTGGTAGAAGGTAATATAGAGAGTAATCAAAAATTGAT  >  NC_000913/1475218‑1475412

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