Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,416,763	2 bp→TC	coding (2017‑2018/2145 nt)	pta →	phosphate acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,416,763	0	C	T	100.0%	127.9 / NA	40	L673S (CTG→TCG)	pta	phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/18);  total (22/18)
*	NC_000913	2,416,764	0	T	C	100.0%	145.2 / NA	40	L673S (CTG→TCG)	pta	phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/18);  total (22/18)

TAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCT  >  NC_000913/2416638‑2416894                                                                                                                              ||                                                                                                                                   tAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCgg                                                                                                                        <  1:673021/139‑1 (MQ=255)      gCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGa                                                                                                                   <  2:127181/139‑1 (MQ=255)      gCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGa                                                                                                                   <  2:762717/139‑1 (MQ=255)       cTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGAt                                                                                                                  <  1:46410/139‑1 (MQ=255)       cTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTACAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGAt                                                                                                                  <  1:368083/139‑1 (MQ=255)                 gAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGt                                                                                                        <  2:741492/139‑1 (MQ=255)                                              cAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCg                                                                                                                    >  1:123989/1‑98 (MQ=255)                                              cAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCg                                                                                                                    <  2:123989/98‑1 (MQ=255)                                              cAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAACGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCg                                                                           >  2:499719/1‑139 (MQ=255)                                                ggTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTg                                                                         >  2:393708/1‑139 (MQ=255)                                                     cgcTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCa                                                                    <  1:688926/139‑1 (MQ=255)                                                               ttCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAACGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACg                                                          >  2:189857/1‑139 (MQ=255)                                                                    cTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGAc                                                                                   <  2:624551/109‑1 (MQ=255)                                                                    cTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGAc                                                                                   >  1:624551/1‑109 (MQ=255)                                                                    cTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATc                                                     >  1:170983/1‑139 (MQ=255)                                                                    cTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATc                                                     >  1:719648/1‑139 (MQ=255)                                                                       cccGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCAc                                                                   <  1:625759/122‑1 (MQ=255)                                                                       cccGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCAc                                                                   >  2:625759/1‑122 (MQ=255)                                                                         cGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTa                                                <  1:644429/139‑1 (MQ=255)                                                                                    acCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGgcg                                                                      <  1:380564/106‑1 (MQ=255)                                                                                    acCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGgcg                                                                      >  2:380564/1‑106 (MQ=255)                                                                                      cGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTg                                   >  1:392331/1‑139 (MQ=255)                                                                                      cGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTg                                   >  2:171027/1‑139 (MQ=255)                                                                                        gTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGAt                                                                                                                  <  1:243264/58‑1 (MQ=255)                                                                                        gTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGAt                                                                                                                  >  2:243264/1‑58 (MQ=255)                                                                                        gTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGAc                                 >  1:771747/1‑139 (MQ=255)                                                                                        gTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACAATCGCGCTGAc                                 >  2:520301/1‑139 (MQ=255)                                                                                             accaccTACAAAGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGa                            >  1:358871/1‑139 (MQ=255)                                                                                                      aaaGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTg                   <  2:82999/139‑1 (MQ=255)                                                                                                       aaGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGc                  <  2:719648/139‑1 (MQ=255)                                                                                                       aaGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGc                  >  2:86346/1‑139 (MQ=255)                                                                                                        aGCGGTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGca                 >  1:802552/1‑139 (MQ=255)                                                                                                            gTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGat                                                        >  2:208485/1‑96 (MQ=255)                                                                                                            gTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGat                                                        <  1:208485/96‑1 (MQ=255)                                                                                                            gTACAGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCAcagc             >  1:92879/1‑139 (MQ=255)                                                                                                                aGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCAcagcagct         <  1:499719/139‑2 (MQ=255)                                                                                                                aGCGTTCTGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCAcagcagca         <  1:86346/139‑1 (MQ=255)                                                                                                                       tGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAAtct  >  1:651069/1‑139 (MQ=255)                                                                                                                       tGCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAAtct  >  1:555483/1‑139 (MQ=255)                                                                                                                       cgCCGACTCGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAAtct  >  1:499104/2‑139 (MQ=255)                                                                                                                              ||                                                                                                                                   TAATGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCGGGCCGATGCTGCAGGGTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCT  >  NC_000913/2416638‑2416894

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: ‑