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%	154.7 / NA	41	L137L (CTG→TTG)	gpmM	phosphoglycero mutase III, cofactor‑independent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/21);  total (20/21)

AATCAAAGATCGCGCTTTCTTTGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATT  >  NC_000913/3785520‑3785814                                                                                                                                                     |                                                                                                                                                   aaTCAAAGATCGCGCTTTCTTTGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTg                                                                                                                                                  >  1:407160/1‑151 (MQ=255)         aTCGCGCTTTCTTTGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGa                                                                                                                                          <  2:656856/151‑1 (MQ=255)           cgcgCTTTCTTTGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAAc                                                                                                                                        <  1:420602/151‑1 (MQ=255)                tttctttGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACgcggc                                                                                                                                   >  2:315911/1‑151 (MQ=255)                   ctttGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGgcgc                                                                                                                                 >  1:667149/1‑150 (MQ=255)                      tGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGTCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGa                                                                                                                              <  2:199151/150‑1 (MQ=255)                       gCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGaa                                                                                                                             <  1:76035/150‑1 (MQ=255)                             ccGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCt                                                                                                                      <  2:167174/151‑1 (MQ=255)                               ggTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTAc                                                                                                                    >  2:466924/1‑151 (MQ=255)                                         ggCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCAt                                                                                                          <  2:19993/151‑1 (MQ=255)                                          gcgcGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCAtt                                                                                                         <  1:58026/151‑1 (MQ=255)                                                   aTAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACgg                                                                                                <  2:548719/151‑1 (MQ=255)                                                    tAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGt                                                                                               >  1:168025/1‑151 (MQ=255)                                                        gCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTcgcg                                                                                           >  2:446493/1‑151 (MQ=255)                                                                   aGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTccgccg                                                                                >  2:556935/1‑151 (MQ=255)                                                                        aaGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAg                                                                           <  1:741993/151‑1 (MQ=255)                                                                           cGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGc                                                                        >  2:571513/1‑151 (MQ=255)                                                                            ggTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCt                                                                       <  1:220252/151‑1 (MQ=255)                                                                                cacaTTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAAt                                                                   <  2:725511/151‑1 (MQ=255)                                                                                 acaTTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATc                                                                  <  1:581226/151‑1 (MQ=255)                                                                                        gggTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCCGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCt                                                            <  2:663592/150‑1 (MQ=255)                                                                                        ggatCTGCTCCTCGACGGCGGCGTACACAGCCCGGAAGGTCACATGATGGCGATGGTAGAATGGGGAGGTGGAAGCGGCGCAGGGAAAATCTACCTGCACGCATTCCTTGCCGGTCGCGACACTCCGCCGCGCGGTGAGGGATCCt                                                                >  1:354402/4‑146 (MQ=255)                                                                                          gTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGa                                                          >  1:158819/1‑150 (MQ=255)                                                                                                 ctcCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAAt                                                    <  2:133348/149‑1 (MQ=255)                                                                                                     gCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAAt                                                    <  1:155374/145‑1 (MQ=255)                                                                                                         gcggcgTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGaaaa                                          <  1:733023/151‑1 (MQ=255)                                                                                                           ggcgTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAAt                                        >  2:479664/1‑151 (MQ=255)                                                                                                            gcgTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAAtt                                       >  1:774374/1‑151 (MQ=255)                                                                                                              gTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTg                                     <  2:343303/151‑1 (MQ=255)                                                                                                                       cACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTgg                            <  2:568508/151‑1 (MQ=255)                                                                                                                          gAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCa                         <  2:63908/151‑1 (MQ=255)                                                                                                                           aaGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCa                         >  2:386575/1‑150 (MQ=255)                                                                                                                           aaGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCa                         >  1:427940/1‑150 (MQ=255)                                                                                                                             gATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCaa                        >  2:585759/1‑149 (MQ=255)                                                                                                                                    tcatGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGt               >  1:753761/1‑151 (MQ=255)                                                                                                                                      atGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGt               >  1:656547/1‑149 (MQ=255)                                                                                                                                         gCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGt          >  2:338060/1‑151 (MQ=255)                                                                                                                                              ggTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCcat      <  2:417025/150‑1 (MQ=255)                                                                                                                                               gTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCcatc     >  1:670282/1‑150 (MQ=255)                                                                                                                                                tAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCcatcat   >  2:231196/1‑151 (MQ=255)                                                                                                                                                 aGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGTCAAAGGCCGCGTAGCGTCCATCAtt  <  1:303339/151‑1 (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: ‑