Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,623,838	G→A	Q265Q (CAG→CAA)	ppk →	polyphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,623,838	0	G	A	100.0%	36.6 / NA	13	Q265Q (CAG→CAA)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/6);  total (7/6)

GAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTT  >  NC_000913/2623592‑2624041                                                                                                                                                                                                                                                       |                                                                                                                                                                                                            gAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAacgc                                                                                                                                                                                                          <  2:48577/250‑1 (MQ=255)               aaGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATgcg                                                                                                                                                                                           >  1:42972/1‑251 (MQ=255)                       cgTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTTGTTGa                                                                                                                                                                                   <  1:61766/251‑1 (MQ=255)                                        gatTCTTCTGGATAACATTCTGCGTTTCTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGTGTTTAACTGCTGAGCCGGTGCGTTTTGGTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAAc                                                                                                                                                                  <  1:169300/251‑1 (MQ=255)                                           tcttctGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTg                                                                                                                                                                <  2:84455/250‑1 (MQ=255)                                                  gATAACATTCTGCGTTACTGCATTGATGATATTTTCAAAGGCTTCTTTGCTTATGACGCGCTGAAGGCATATTCAATGAAGATGAACCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGAGCGTTTTGTTTATCAACGAGAACTGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTc                                                                                                                                                        >  1:220720/1‑251 (MQ=255)                                                            tGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGatat                                                                                                                                                                                                     >  1:136719/1‑196 (MQ=255)                                                            tGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGGTTAACTGCTGAGCCGGTGCTTTTTGTTTATCAACGCGatat                                                                                                                                                                                                     <  2:136719/196‑1 (MQ=255)                                                                                             tcttTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGGCCCCGGCGGTCGTTATCATAAtt                                                                                                               <  2:17622/249‑1 (MQ=255)                                                                                                                      cTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATg                                                                                    >  1:85536/1‑251 (MQ=255)                                                                                                                                      aCCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAAACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCt                                                                    >  2:149138/1‑251 (MQ=255)                                                                                                                                         cgcgATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTGTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGGGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGt                                                                 <  1:240653/251‑1 (MQ=255)                                                                                                                                                                              ccagccTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGttt                            >  1:25978/1‑251 (MQ=255)                                                                                                                                                                                                        aGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGttt  >  2:261335/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                            GAATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTT  >  NC_000913/2623592‑2624041

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