Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,117,320	G→T	R257R (CGC→CGA)	glpF ←	glycerol facilitator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,320	0	G	T	96.6%	199.1 / 0.5	59	R257R (CGC→CGA)	glpF	glycerol facilitator
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (33/24);  total (34/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 = 9.71e-01

CCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCA  >  NC_000913/4117185‑4117452                                                                                                                                        |                                                                                                                                     cccTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGAc                                                                                                                                  >  1:329968/1‑140 (MQ=255)     ggTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAAt                                                                                                                              <  1:773371/140‑1 (MQ=255)     ggTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAAt                                                                                                                              <  1:140932/140‑1 (MQ=255)     ggTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAAt                                                                                                                              <  1:707317/140‑1 (MQ=255)      gTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATc                                                                                                                             >  1:757638/1‑140 (MQ=255)         gAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGt                                                                                                                          >  2:43095/1‑140 (MQ=255)           gcgcAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGttt                                                                                                                        >  2:643041/1‑140 (MQ=255)            cgcAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTg                                                                                                                       <  1:340882/140‑1 (MQ=255)               aaCGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCgg                                                                                                                    <  1:643041/140‑1 (MQ=255)                        ttttttttCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGGAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAAt                                                                                                           <  1:689647/140‑1 (MQ=255)                        ttttttttCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAAt                                                                                                           <  1:296075/140‑1 (MQ=255)                        ttttttttCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAAt                                                                                                           <  2:1018107/140‑1 (MQ=255)                                   tCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACATATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAAt                                                                                                <  1:59326/140‑1 (MQ=255)                                    cATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATc                                                                                               <  1:446252/140‑1 (MQ=255)                                         ttAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGcc                                                                                          <  2:643721/140‑1 (MQ=255)                                                tCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATa                                                                                   >  2:923196/1‑140 (MQ=255)                                                 cccGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAg                                                                                  >  2:1031789/1‑140 (MQ=255)                                                 cccGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAg                                                                                  >  2:505479/1‑140 (MQ=255)                                                  ccGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGcc                                                                                          <  2:37544/131‑1 (MQ=255)                                                  ccGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGcc                                                                                          >  1:37544/1‑131 (MQ=255)                                                  ccGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAgg                                                                                 >  1:593499/1‑140 (MQ=255)                                                          atatTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGaaaa                                                                         >  1:17464/1‑140 (MQ=255)                                                           tatTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAg                                                                        <  1:426000/140‑1 (MQ=255)                                                             ttACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCg                                                                      <  1:758955/140‑1 (MQ=255)                                                                 aGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCAc                                                                  >  2:250932/1‑140 (MQ=255)                                                                  gCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGTACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCAcc                                                                 >  1:9347/1‑140 (MQ=255)                                                                  gCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCAcc                                                                 <  1:923196/140‑1 (MQ=255)                                                                  gCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCAcc                                                                 >  1:1025334/1‑140 (MQ=255)                                                                  gCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCAcc                                                                 <  1:791373/140‑1 (MQ=255)                                                                  gCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCAcc                                                                 >  1:709088/1‑140 (MQ=255)                                                                  gCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCAcc                                                                 >  1:665712/1‑140 (MQ=255)                                                                    gAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAg                                                               <  2:17464/140‑1 (MQ=255)                                                                     aaGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCa                                                                >  2:364028/1‑138 (MQ=255)                                                                     aaGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCa                                                                <  1:364028/138‑1 (MQ=255)                                                                     aaGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAgg                                                              >  2:477237/1‑140 (MQ=255)                                                                     aaGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAgg                                                              >  2:28267/1‑140 (MQ=255)                                                                     aaGCTTTTTGTTCTGAAGGAGCTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAgg                                                              >  2:767496/1‑140 (MQ=255)                                                                      aGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGa                                                             <  1:28267/140‑1 (MQ=255)                                                                      aGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGa                                                             >  1:632691/1‑140 (MQ=255)                                                                        cTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAg                                                           <  1:43095/140‑1 (MQ=255)                                                                                       gAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTg                                            <  2:648801/140‑1 (MQ=255)                                                                                         gTTGTGGTTTCTTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTgcc                                          >  1:759871/1‑140 (MQ=255)                                                                                         gTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTgcc                                          >  1:395383/1‑140 (MQ=255)                                                                                         gTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTgcc                                          >  1:132416/1‑140 (MQ=255)                                                                                                         cttcCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACa                          <  2:759871/140‑1 (MQ=255)                                                                                                         cttcCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACa                          <  1:505479/140‑1 (MQ=255)                                                                                                         cttcCACAACACAGATATCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACa                          <  2:395383/140‑1 (MQ=255)                                                                                                                       atatCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCgc            >  2:383150/1‑140 (MQ=255)                                                                                                                       atatCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCgc            >  2:892880/1‑140 (MQ=255)                                                                                                                          tCGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCgccag         >  1:579413/1‑140 (MQ=255)                                                                                                                           cGCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCTCCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCgccagc        <  1:892880/140‑1 (MQ=255)                                                                                                                            gCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCgccagcc       >  2:571179/1‑140 (MQ=255)                                                                                                                            gCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCgccagcc       >  1:874873/1‑140 (MQ=255)                                                                                                                            gCAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCgccagcc       >  1:828713/1‑140 (MQ=255)                                                                                                                             cAAGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCgccagcca      >  2:220513/1‑140 (MQ=255)                                                                                                                               aGGCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAgg    >  2:111757/1‑140 (MQ=255)                                                                                                                                ggCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATG‑ACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCa  >  1:695651/1‑140 (MQ=255)                                                                                                                                 gCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCa  >  2:397234/1‑140 (MQ=255)                                                                                                                                 gCAAATGTCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCa  >  2:750582/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                     CCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCA  >  NC_000913/4117185‑4117452

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