Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,375,282	A→G	Y126H (TAT→CAT)	epmB ←	EF‑P‑Lys34 lysylation protein, weak lysine 2,3‑aminomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,375,282	0	A	G	100.0%	171.7 / NA	48	Y126H (TAT→CAT)	epmB	EF‑P‑Lys34 lysylation protein, weak lysine 2,3‑aminomutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/24);  total (24/24)

CAGCGGATCGCCGCCGGAGAAAATCATCTCGTCCAGTTCCGGATGCGCAGCAACATACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATAGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAA  >  NC_000913/4375154‑4375411                                                                                                                                 |                                                                                                                                  cAGCGGATCGCCGCCGGAGAAAATCATCTCGTCCAGTTCCGGATGCGCAGCAACATACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAAttt                                                                                                                         >  1:1254313/1‑139 (MQ=255)        tCGCCGCCGGAGAAAATCATCTCGTCCAGTTCCGGATGCGCAGCAACATACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCgcgc                                                                                                                  >  2:1410451/1‑139 (MQ=255)           ccgccgGAGAAAATCATCTCGTCCAGTTCCGGATGCGCAGCAACATACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAgc                                                                                                               <  1:548630/139‑1 (MQ=255)                ggAGAAAATCATCTCGTCCAGTTCCGGATGCGCAGCAACATACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCt                                                                                                          <  2:49936/139‑1 (MQ=255)                           tctcGTCCAGTTCCGGATGCGCAGCAACATACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCaa                                                                                               <  2:119599/139‑1 (MQ=255)                            ctcGTCCAGTTCCGGATGCGCAGCAACATACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAg                                                                                                                <  1:1221587/121‑1 (MQ=255)                            ctcGTCCAGTTCCGGATGCGCAGCAACATACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAg                                                                                                                >  2:1221587/1‑121 (MQ=255)                                tCCAGTTCCGGATGCGCAGCAACATACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGACAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAgcg                                                                                          >  2:309141/1‑139 (MQ=255)                                  cAGTTCCGGATGCGCAGCAACATACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGcc                                                                                        <  1:1040355/139‑1 (MQ=255)                                  cAGTTCCGGATGCGCAGCAACATACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGcc                                                                                        <  2:362082/139‑1 (MQ=255)                                        cGGATGCGCAGCAACATACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTtg                                                                                  <  2:1254313/139‑1 (MQ=255)                                          gATGCGCAGCAACATACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCgg                                                                                     >  1:1199482/1‑134 (MQ=255)                                          gATGCGCAGCAACATACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCgg                                                                                     <  2:1199482/134‑1 (MQ=255)                                             gcgcAGCAACATACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTa                                                                             >  2:153639/1‑139 (MQ=255)                                                       aTACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGcaac                                                                   <  1:808702/139‑1 (MQ=255)                                                        tACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGcaaca                                                                  >  1:1002515/1‑139 (MQ=255)                                                        tACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGcaaca                                                                  >  2:1379244/1‑139 (MQ=255)                                                        tACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGcaaca                                                                  <  1:595993/139‑1 (MQ=255)                                                            cAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAAcc                                                              <  1:230430/139‑1 (MQ=255)                                                             aaGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGcc                                                                                        >  1:496002/1‑112 (MQ=255)                                                             aaGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGcc                                                                                        <  2:496002/112‑1 (MQ=255)                                                               gTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAgg                                                           >  2:1371944/1‑139 (MQ=255)                                                                tGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGc                                                          >  2:34887/1‑139 (MQ=255)                                                                 gCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTg                                                                                                       >  1:1347942/1‑93 (MQ=255)                                                                 gCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTg                                                                                                       <  2:1347942/93‑1 (MQ=255)                                                                  cAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGGAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCATCCGCCTTTGACCAGCAAAAGCGCCCGGGTGTTGTATTTATGCAACAAACCAGGCAc                                                        >  2:855007/1‑139 (MQ=255)                                                                  cAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAAcc                                                              >  1:813893/1‑133 (MQ=255)                                                                  cAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAAcc                                                              <  2:813893/133‑1 (MQ=255)                                                                          cAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACgctgt                                                <  1:1379244/139‑1 (MQ=255)                                                                          cAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACgctgt                                                >  2:237995/1‑139 (MQ=255)                                                                          cAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACg                                                    <  2:1114773/135‑1 (MQ=255)                                                                          cAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACg                                                    >  1:1114773/1‑135 (MQ=255)                                                                                           cTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACgctgtgctg                                            <  2:286334/126‑1 (MQ=255)                                                                                           cTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACgctgtgctg                                            >  1:286334/1‑126 (MQ=255)                                                                                           cTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTggg                               >  2:1239730/1‑139 (MQ=255)                                                                                             gATTTTCGGCATAGGGGAAGTGACGAAGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCGAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTc                             <  1:855007/139‑1 (MQ=255)                                                                                                        tAGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATc                  >  2:896886/1‑139 (MQ=255)                                                                                                         aGGGGAAGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATcc                 >  1:1141926/1‑139 (MQ=255)                                                                                                              aaGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGgcg            <  2:1141926/139‑1 (MQ=255)                                                                                                              aaGTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGgcg            <  1:1410451/139‑1 (MQ=255)                                                                                                                gTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCg                                                                                      <  2:21460/63‑1 (MQ=255)                                                                                                                gTGACGACGGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCg                                                                                      >  1:21460/1‑63 (MQ=255)                                                                                                                   acgacgGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGcc                                                                                        <  2:103491/58‑1 (MQ=255)                                                                                                                   acgacgGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGcc                                                                                        >  1:103491/1‑58 (MQ=255)                                                                                                                   acgacgGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACg                                                    >  2:179079/1‑94 (MQ=255)                                                                                                                   acgacgGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACg                                                    <  1:179079/94‑1 (MQ=255)                                                                                                                      acgGAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACa    >  2:706032/1‑139 (MQ=255)                                                                                                                        ggAAGCAATGGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACaaa  <  1:237995/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                  CAGCGGATCGCCGCCGGAGAAAATCATCTCGTCCAGTTCCGGATGCGCAGCAACATACTCAAGTGCAGTTTGCCAGTTACGCTTGTTGCCCTGATTTTCGGCATAGGGGAAGTGACGACGGAAGCAATAGCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAA  >  NC_000913/4375154‑4375411

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