Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,785,668	C→T	L137L (CTG→TTG)	gpmM →	phosphoglycero mutase III, cofactor‑independent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,785,668	0	C	T	100.0%	219.4 / NA	57	L137L (CTG→TTG)	gpmM	phosphoglycero mutase III, cofactor‑independent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/25);  total (32/25)

AATCAAAGATCGCGCTTTCTTTGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATT  >  NC_000913/3785520‑3785814                                                                                                                                                     |                                                                                                                                                   aaTCAAAGATCGCGCTTTCTTTGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTg                                                                                                                                                  <  1:3291/151‑1 (MQ=255)    cAAAGATCGCGCTTTCTTTGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCa                                                                                                                                               <  2:335559/151‑1 (MQ=255)     aaaGATCGCGCTTTCTTTGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGc                                                                                                                                                >  1:896233/1‑149 (MQ=255)              gcTTTCTTTGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAAcgcg                                                                                                                                     <  2:890098/151‑1 (MQ=255)               ctttctttGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACgcgg                                                                                                                                    >  1:62872/1‑151 (MQ=255)                 ttctttGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACgcggcg                                                                                                                                  >  2:291764/1‑151 (MQ=255)                 ttctttGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACgcgg                                                                                                                                    <  2:195642/149‑1 (MQ=255)                   ctttGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGgcgc                                                                                                                                 >  1:554674/1‑150 (MQ=255)                     ttGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGgcgc                                                                                                                                 >  1:271993/1‑148 (MQ=255)                          aaTCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACCTTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCCGaaaaa                                                                                                                          >  2:793636/1‑150 (MQ=255)                             ccGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATc                                                                                                                       <  1:431380/150‑1 (MQ=255)                             ccGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCt                                                                                                                      <  1:407723/151‑1 (MQ=255)                                    tGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGACGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGc                                                                                                                >  1:799853/1‑150 (MQ=255)                                     gACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGc                                                                                                                <  1:325694/149‑1 (MQ=255)                                         ggCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCa                                                                                                           <  1:567732/150‑1 (MQ=255)                                         ggCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCa                                                                                                           <  2:634407/150‑1 (MQ=255)                                            gcgGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTcc                                                                                                       >  2:594919/1‑151 (MQ=255)                                             cgGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCt                                                                                                      <  2:620984/151‑1 (MQ=255)                                               gTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTg                                                                                                    <  2:312636/151‑1 (MQ=255)                                                   aTAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACgg                                                                                                <  1:865565/151‑1 (MQ=255)                                                    tAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGt                                                                                               >  2:662883/1‑151 (MQ=255)                                                     aaaGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGt                                                                                               <  1:898728/150‑1 (MQ=255)                                                     aaaGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTc                                                                                              >  1:411948/1‑151 (MQ=255)                                                     aaaGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTc                                                                                              >  2:345930/1‑151 (MQ=255)                                                     aaaGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTc                                                                                              >  1:342964/1‑151 (MQ=255)                                                         cGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTTGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGCCGGTcgcg                                                                                           >  2:894605/1‑150 (MQ=255)                                                          gAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGATCTCCGCAGGCGGCGTACACAGCCACGAAGATCAGATCAAGGCGATGGTAGAATTGGAAGCTGAACGCGGCGCAGAAAAAATATACCTGCACGCATTGATTGACGGAcgc                                                                                            >  2:139901/1‑148 (MQ=255)                                                              aaCGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTc                                                                                     >  1:529529/1‑151 (MQ=255)                                                               aCGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTcc                                                                                    <  2:260234/151‑1 (MQ=255)                                                                 gCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTcc                                                                                    >  2:803680/1‑149 (MQ=255)                                                                      cAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCgcgc                                                                             >  1:6489/1‑151 (MQ=255)                                                                       aaaGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCcgcg                                                                              <  2:647332/149‑1 (MQ=255)                                                                       aaaGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCgcgc                                                                             >  2:866424/1‑150 (MQ=255)                                                                        aaGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCgcgc                                                                             >  1:319510/1‑149 (MQ=255)                                                                                cacaTTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAAt                                                                   <  1:371630/151‑1 (MQ=255)                                                                                  caTTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATcc                                                                 <  1:89411/151‑1 (MQ=255)                                                                                    ttATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTc                                                               >  1:747996/1‑151 (MQ=255)                                                                                      aTGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGc                                                             <  1:280906/151‑1 (MQ=255)                                                                                             tGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGaaaa                                                       >  1:814686/1‑150 (MQ=255)                                                                                              gCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGaaaaaa                                                     <  1:424894/151‑1 (MQ=255)                                                                                              gCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGaaaaaa                                                     <  2:271906/151‑1 (MQ=255)                                                                                               ctctCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAAt                                                    >  2:472361/1‑151 (MQ=255)                                                                                                    cGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTcaa                                                >  1:586501/1‑148 (MQ=255)                                                                                                     gCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCgaag                                              <  1:816921/151‑1 (MQ=255)                                                                                                     gCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCga                                                >  2:416843/1‑149 (MQ=255)                                                                                                              gTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTg                                     <  2:499589/151‑1 (MQ=255)                                                                                                               tACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGc                                    >  2:722635/1‑151 (MQ=255)                                                                                                                 cacaGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCcg                                  >  2:908518/1‑151 (MQ=255)                                                                                                                   caGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCcgcg                                >  1:204434/1‑151 (MQ=255)                                                                                                                      ccACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGaaa                                           <  2:470025/137‑1 (MQ=255)                                                                                                                      ccACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGaaa                                           >  1:470072/1‑137 (MQ=255)                                                                                                                                 acaTCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCcg                  >  1:67604/1‑151 (MQ=255)                                                                                                                                     catGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGt               <  2:271967/150‑1 (MQ=255)                                                                                                                                      atGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGt               >  1:765098/1‑149 (MQ=255)                                                                                                                                      atGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAg             <  1:481538/151‑1 (MQ=255)                                                                                                                                            aTGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTcc        >  2:413615/1‑150 (MQ=255)                                                                                                                                                 aGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCAtt  >  2:806923/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AATCAAAGATCGCGCTTTCTTTGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATT  >  NC_000913/3785520‑3785814

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