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

TTCTTTGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATT  >  NC_000913/3785536‑3785814                                                                                                                                     |                                                                                                                                                   ttctttGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACgcgg                                                                                                                                    <  2:162178/149‑1 (MQ=255)  tctttGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACgcggcg                                                                                                                                  <  2:765873/150‑1 (MQ=255)  tctttGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGgcgc                                                                                                                                 >  2:597428/1‑151 (MQ=255)    tttGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGgcgc                                                                                                                                 >  1:661693/1‑149 (MQ=255)     ttGCTAATCCGGTGCTGACTGGCGCGGTAGTTAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACCCAGCCACGAAGATCTCATCATGGCGATGGTAGAATTGGCAGCTGAACGCGgcgct                                                                                                                                <  1:36929/149‑2 (MQ=255)          aaTCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGaaaaa                                                                                                                          <  1:337219/150‑1 (MQ=255)          aaTCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGaaaaa                                                                                                                          >  1:340914/1‑150 (MQ=255)           aTCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAAt                                                                                                                        <  2:753071/151‑1 (MQ=255)               ggTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTAc                                                                                                                    >  2:372201/1‑151 (MQ=255)                     gACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCAc                                                                                                              >  1:87154/1‑151 (MQ=255)                         ggCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCAt                                                                                                          <  1:517743/151‑1 (MQ=255)                           cgcgGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCACCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTc                                                                                                        >  1:180964/1‑151 (MQ=255)                            gcgGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTcc                                                                                                       <  1:597389/151‑1 (MQ=255)                                 aGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGAc                                                                                                  >  1:199364/1‑151 (MQ=255)                                 aGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGAc                                                                                                  >  2:580871/1‑151 (MQ=255)                                       aGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTcgc                                                                                            <  1:798895/151‑1 (MQ=255)                                       aGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTcgc                                                                                            >  2:320804/1‑151 (MQ=255)                                       aGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTcgc                                                                                            <  1:437743/151‑1 (MQ=255)                                        gCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTcgcg                                                                                           >  2:378957/1‑151 (MQ=255)                                                 gCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGTGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTccgc                                                                                  >  2:134723/1‑151 (MQ=255)                                                     gCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCcgcg                                                                              <  1:203822/151‑1 (MQ=255)                                                       aaaGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCgcgc                                                                             >  1:94956/1‑150 (MQ=255)                                                       aaaGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCa                                                                            <  2:236488/151‑1 (MQ=255)                                                             gTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTg                                                                      <  1:566153/151‑1 (MQ=255)                                                             gTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTg                                                                      <  2:616114/151‑1 (MQ=255)                                                               acacATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGa                                                                     <  2:86971/150‑1 (MQ=255)                                                               acacATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGa                                                                     <  1:237912/150‑1 (MQ=255)                                                                cacaTTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAAt                                                                   <  2:108412/151‑1 (MQ=255)                                                                          gTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGa                                                          >  1:114658/1‑150 (MQ=255)                                                                             tGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGaaaa                                                       >  2:728086/1‑150 (MQ=255)                                                                              gCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGaaa                                                        <  2:241312/148‑1 (MQ=255)                                                                                tctcCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAAtt                                                   >  1:581009/1‑151 (MQ=255)                                                                                    cGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCgaa                                               <  1:335505/151‑1 (MQ=255)                                                                                    cGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCgaa                                               <  2:428950/151‑1 (MQ=255)                                                                                       aGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCgaagaa                                            <  1:783117/151‑1 (MQ=255)                                                                                               tACACAGCCACGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGc                                    >  1:470077/1‑151 (MQ=255)                                                                                                        aCGAAGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTggg                           <  1:496031/151‑1 (MQ=255)                                                                                                           aaGATCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCa                         >  1:423448/1‑150 (MQ=255)                                                                                                              aTCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAgg                     <  1:471796/151‑1 (MQ=255)                                                                                                               tCACATCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGc                    >  1:483014/1‑151 (MQ=255)                                                                                                                cacaTCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCGGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGcc                   >  1:472874/1‑151 (MQ=255)                                                                                                                 acaTCATGGCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCcg                  <  2:214129/151‑1 (MQ=255)                                                                                                                        ggCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCg           >  1:34517/1‑151 (MQ=255)                                                                                                                        ggCGATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCACTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCGGGGCAAAGGCCGCGTAGCg           <  1:601430/151‑1 (MQ=255)                                                                                                                           gATGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTcc        >  1:706482/1‑151 (MQ=255)                                                                                                                            aTGGTAGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCca       <  1:316306/151‑1 (MQ=255)                                                                                                                                 aGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCTGCGCTGGGCAAAGGCCGCGTAGCGTCCATCAtt  >  1:853979/1‑151 (MQ=255)                                                                                                                                 aGAATTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCcatcat   <  2:464525/150‑1 (MQ=255)                                                                                                                                 aGAATTGGCAGCTGAACGCGGCGCAG‑AAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATACGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCAtt  >  1:793909/1‑150 (MQ=255)                                                                                                                                     |                                                                                                                                                   TTCTTTGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATT  >  NC_000913/3785536‑3785814

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