Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,117,034	T→A	E63D (GAA→GAT)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,034	0	T	A	92.9%	75.9 / ‑3.3	29	E63D (GAA→GAT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (10/16);  total (12/17)
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.52e-01

CGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGAC  >  NC_000913/4116759‑4117324                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                   cGCTCGCGAGAGCCTGCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCAGTACGACGGCACTGCCAGACAATGGCGTGATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTgg                                                                                                                                                                                                                                                                           <  1:71623/301‑1 (MQ=255)     cgcgAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGCTTTAACTGCTCGCAGATTTCTGCGGTACGACGGCATTGCCAGACAATGGCGTTATAGATAGGCGTGCCGGTTTCTTTTTCCCAGGCAATCGTGGTTTCACGCTGTTTTGTAATTCCGATTTCTGCCATTTTGCCGGCACTGAACTCGGCTCTCGCCCGCACGTcccacagtcttgcaagtttggtgtgcc                                                                                                                                                                                                                                                                        >  2:234788/1‑274 (MQ=255)             ccgtgtgctcttccgatcttCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCa                                                                                                                                                                                                                                                               <  1:14644/282‑1 (MQ=255)                cccctacacgacgctcttccgatctcttTGGTGGCAGAAAAGTACGGGTTAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTg                                                                                                                                                                                                                                                            <  2:103369/276‑1 (MQ=255)                          cgtggttCCACTTCACTTTGGTGCCAGAACAGTACGTGTCAATCGCCAGACCGGTATTGCTGCGGATATAATCTTCTACACCGTCACGTTTGAAATGCTCGCAGATTTCTGCGGTACGATGGCCCTTCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGCCAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGTAACTGATATCGGCTTTCGCCGGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTc                                                                                                                                                                                                                                                  <  2:206059/295‑1 (MQ=255)                                tCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGGTTTCCaagatcggaagagggtggtg                                                                                                                                                                                                                                            >  2:14643/1‑282 (MQ=255)                                         cttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTgagatcggaagagcacacgtctgaa                                                                                                                                                                                                                                    >  1:103376/1‑276 (MQ=255)                                             ggTGCCAGAAAAGTACGGGTGAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGATTTCACGCTAGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                                                                                                                                                <  1:98549/300‑1 (MQ=255)                                              gtaccAGAAATATACGGGTGAATCACCCGACCGGGTTGGCTGCGGGTACAAGCTTCTAAATCGTCACGTTTTAAATTCTCGCCGATTTCTTCGGTACGACGGCACTGCCCACCAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTggg                                                                                                                                                                                                                               <  2:171296/297‑1 (MQ=255)                                               tGCCAGAAAAGTACGGGGCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTa                                                                                                                                                                                                                             <  1:35820/301‑1 (MQ=255)                                                  cAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAt                                                                                                                                                                                                                          >  1:189021/1‑301 (MQ=255)                                                               gTCAATCACGAGACCGGTATAGATTAGCATATAAGCTTCTAAACAGTCACGTTTTAAATACTCGCAGATTTCTGCGTTACGACGGTGCGGCCGGACAGTGTCGTTATAGATAGGCTTGCCGGTTACTTTGTCGCAGACAATAGTGAATTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTGCCAGCGTGGAGCTTTGGGTGGCCCAGGTTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTc                                                                                                                                                                                                              <  2:72804/300‑1 (MQ=255)                                                                      aCCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTTGGTCGTGTTCTTCCCAACCTGGTTTTGGGTATATTTGCTCAAATTCGCGCTGc                                                                                                                                                                                                       >  1:153142/1‑300 (MQ=255)                                                                           aCCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACg                                                                                                                                                                                                 <  1:300616/301‑1 (MQ=255)                                                                                               aaTCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTg                                                                                                                                                                             >  1:279730/1‑301 (MQ=255)                                                                                                              cACGTGTTAAATGTTCGCCAGGTTCTTCGGTGAGAAGGCACTGCCAGACAATAGCGTTATAGATATGCTTGCCGGTTTCTTTTACCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTAGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCCGATTTCCATTGGGTCGTGTTCTACCCAACCTAGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACcg                                                                                                                                                               <  2:153132/300‑1 (MQ=255)                                                                                                                      atggagttcagacgtgtgctcttccgatctcACTGCCAGACAATGTCGTTATAGATAGGCTTGCCGATTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGGGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGc                                                                                                                                                      <  1:44547/271‑1 (MQ=255)                                                                                                                           cTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCACATTCGCGCTGAGACACGCGAATGATATTTGCATCGTGATCCATTACGACCGCGTGGGAGCtggtg                                                                                                                                                 >  2:308491/1‑301 (MQ=255)                                                                                                                                                    cACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAATTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGACACGCTCATGATTTTGGCATCGTGATCCATTACTACCGCNCGGGAGccgtagcgaagaggggcgggtagggatagag                                                                                                                        >  2:44544/1‑271 (MQ=255)                                                                                                                                                               aaTGGCGTTATAGATGGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCTCGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCa                                                                                                             <  2:279718/301‑1 (MQ=255)                                                                                                                                                                    cGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCAt                                                                                                        >  1:180070/1‑301 (MQ=255)                                                                                                                                                                           attagTCTTGGCGGTTTCCTGTTCCCGGGTAATCGTGGTTTGACTATTGTTCGGTATGCCGATAGCTGCAATTTGATCGAAACTGTTATCGGCTTTCGCCAGCACATCTGCCGGTGTGGAGCTTTGGGTGGCTCAGATTGCCATGGGGTCGTGTTCTACCCAACCTGGTTTTGAGTAGATTTGCTCAACTTCGCGCTGCTACACGCTAATGATATTGGCATGGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGGTTaa                                                                                                  <  2:202070/298‑1 (MQ=255)                                                                                                                                                                                 ctggCCGGTTTCGTTTGACCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTGGCTCGGAACTGAGATCGGCTTTCGCCAGCACATCTACCCGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTcc                                                                                            <  2:70275/297‑1 (MQ=255)                                                                                                                                                                                    gCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCCTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGAAACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGTAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTa                                                                                        >  1:61258/1‑301 (MQ=255)                                                                                                                                                                                     ccGGTTTTTTTTTCCCAGAAAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACGTATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTa                                                                                        <  2:180060/300‑1 (MQ=255)                                                                                                                                                                                     ccGGTTTCTTTTTCCAAGACATTAGTGGGTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTa                                                                                        <  2:189011/300‑1 (MQ=255)                                                                                                                                                                                            cTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATAtt                                                                                >  1:152349/1‑301 (MQ=255)                                                                                                                                                                                                                          ggTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAACCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGttt                                                  >  1:317602/1‑300 (MQ=255)                                                                                                                                                                                                                                                                          tCGCCAGCACATCTACCAGCGTGGAGCTTGTGGTGGCCCAGATTTCCATTGTGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGAc  <  1:160031/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                   CGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGAC  >  NC_000913/4116759‑4117324

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑