Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,623,776	G→A	E245K (GAG→AAG)	ppk →	polyphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,623,776	0	G	A	100.0%	75.7 / NA	25	E245K (GAG→AAG)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/10);  total (15/10)

ACCATCCGTTACGCGCTGCTGGAGATCCCATCAGATAAAGTGCCGCGCTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAG  >  NC_000913/2623539‑2624022                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       aCCATCCGTTACGCGCTGCTGGAGATCCCATCAGATAAAGTGCCGCGCTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTg                                                                                                                                                                                                                                                   >  2:338520/1‑243 (MQ=255) aCCATCCGTTACGCGCTGCTGGAGATCCCATCAGATAAAGTGCCGCGCTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTg                                                                                                                                                                                                                                                   <  1:338520/243‑1 (MQ=255)                            ccATCAGATAAAGTGCCGCGCTTTGTGAATTTACCGCCAGCAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAActgctg                                                                                                                                                                                                                 >  2:94027/1‑250 (MQ=255)                                   aTAAAGTGCCGCGTTTTGGAAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTg                                                                                                                                                                                                         <  2:392552/251‑1 (MQ=255)                                                              cGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGccc                                                                                                                                                                              <  2:114889/251‑1 (MQ=255)                                                                 cAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAAt                                                                                                                                                                           >  2:296401/1‑251 (MQ=255)                                                                                  cGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGtt                                                                                                                                                          <  1:112344/251‑1 (MQ=255)                                                                                   gCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTa                                                                                                                                                         <  2:395622/251‑1 (MQ=255)                                                                                           atgatTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGaaa                                                                                                                                                  >  2:226267/1‑250 (MQ=255)                                                                                                  ttctGGATAACATTCGGAGTACCTGCCTTGATGATACTTTCCAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATAAAGATGACCAGCGATGCAGAATACGCTTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGTGATATGCCCATTGCGCTGGTTGATGTGTTACGCGAAAAACTGACt                                                                                                                                          >  2:179459/1‑251 (MQ=255)                                                                                                   tctGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACt                                                                                                                                          <  1:296401/250‑1 (MQ=255)                                                                                                       gATAACATTCTGCGTTACTGCCTTGATGATATTTGCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTTCGTGTTGTTTATCAGCGCGATATGCCCAAAGCGCTGGTTAAAGTGTTACTCAAAAAACGGACTATTTc                                                                                                                                     <  1:333175/251‑1 (MQ=255)                                                                                                              ttCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTAc                                                                                                                              >  1:79982/1‑251 (MQ=255)                                                                                                                   cGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAACACCTGACTATTTACCGCTACGACTc                                                                                                                         >  1:388718/1‑251 (MQ=255)                                                                                                                       aCTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATc                                                                                                                     >  2:389944/1‑251 (MQ=255)                                                                                                                               gatgatATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCcgg                                                                                                             >  2:27811/1‑251 (MQ=255)                                                                                                                                   atatTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGt                                                                                                         >  1:350442/1‑251 (MQ=255)                                                                                                                                            aaGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGCCCCGCGATGCCGAATACGATTTAGTACATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCCAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGCAAACCTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCa                                                                                                 >  2:340509/1‑250 (MQ=255)                                                                                                                                                   cttTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAAtttt                                                                                          >  1:102838/1‑250 (MQ=255)                                                                                                                                                                      aaTGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTccc                                                                      <  2:102838/251‑1 (MQ=255)                                                                                                                                                                               tCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTTACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCACTGTCgg                                                             >  1:304034/1‑251 (MQ=255)                                                                                                                                                                                      aGATGACCCGCCATGCCGACTACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAACAACTGACTATTTCACGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTACAGACTTTATTAATTTCCCCAATGTCGGCAAAGcc                                                      >  2:91935/1‑251 (MQ=255)                                                                                                                                                                                                                                aaGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGt                                                                <  1:285859/199‑1 (MQ=255)                                                                                                                                                                                                                                aaGCCAGCCTGATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGt                                                                >  2:285859/1‑199 (MQ=255)                                                                                                                                                                                                                                          gATGAAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAg  <  2:364633/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       ACCATCCGTTACGCGCTGCTGGAGATCCCATCAGATAAAGTGCCGCGCTTTGTGAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAG  >  NC_000913/2623539‑2624022

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