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	91.3%	59.7 / ‑3.3	23	Q265Q (CAG→CAA)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (10/11);  total (11/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

AGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTTTGATGCCATTCGCGAACGCGATGTGTTG  >  NC_000913/2623604‑2624069                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                        aGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATg                                                                                                                                                                                                                         >  2:13119/1‑251 (MQ=255)                               tcttctGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTg                                                                                                                                                                                            <  1:128298/250‑1 (MQ=255)                                         aaCATTCTTCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTg                                                                                                                                                                                                                <  2:261378/220‑1 (MQ=255)                                         aaCATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTg                                                                                                                                                                                                                >  1:261378/1‑220 (MQ=255)                                                  cGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTc                                                                                                                                                                        >  2:255568/1‑251 (MQ=255)                                                   gTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCGGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTcc                                                                                                                                                                       <  2:337266/251‑1 (MQ=255)                                                         gCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCTAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGt                                                                                                                                                                  <  1:366370/250‑1 (MQ=255)                                                                           aaGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCAt                                                                                                                                               >  1:131587/1‑251 (MQ=255)                                                                              gCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAAt                                                                                                                                            <  1:255568/251‑1 (MQ=255)                                                                                 tcttTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAAttt                                                                                                                                          <  1:252021/250‑1 (MQ=255)                                                                                             aCGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTttatt                                                                                                                             <  1:398327/251‑1 (MQ=255)                                                                                                cgcTGTATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCAGGAGATGGTAGCCAGCCTGATGGAGTTGATGGCTTCCATTCTCAAGCAGCGTTTAACTGCTGATCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGTCTATTTCCCGTTACGACTCCATCGGCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAAt                                                                                                                          <  2:383416/251‑1 (MQ=255)                                                                                                                             cgcgATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGc                                                                                                      <  2:219517/242‑1 (MQ=255)                                                                                                                             cgcgATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGc                                                                                                      >  1:219517/1‑242 (MQ=255)                                                                                                                                gATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGa                                                                                           >  1:45676/1‑250 (MQ=255)                                                                                                                                                 aGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCACTGTCGGCAAAGCCCATCTGGTGAAAAAACCACTCCCGCGtt                                                                         >  2:365969/1‑251 (MQ=255)                                                                                                                                                                                ttGATGACTTCCAGTCTCAATCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGACCAATGCGCTGGTTGAAGTGTTACGCGCAAAACTGACTATTTACCGCTCCGACTCCATCGTCACCGGCGGTCGTTATCCTAATTTTAAAAACTTTATTAATTTCCCCAATTTCGGCAAAGCAAATCTGTTGACCAAACCAATGCCGCGTTTACGCCAAATTTGGTTTGATAAAGCCCattt                                          >  2:298478/1‑248 (MQ=255)                                                                                                                                                                                   aTGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACTACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCg                                       >  1:369438/1‑251 (MQ=255)                                                                                                                                                                                           cAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGtt                               >  1:440005/1‑251 (MQ=255)                                                                                                                                                                                                 cAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGATTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTTTg                            <  1:413033/248‑1 (MQ=255)                                                                                                                                                                                                                    gAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCcgc                                                                            <  1:365400/181‑1 (MQ=255)                                                                                                                                                                                                                    gAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCcgc                                                                            >  2:365400/1‑181 (MQ=255)                                                                                                                                                                                                                        cGGTGCGTTTGGTTTACAAACGCGGTATGCCCAATGCGCTTGTGTAAGTTTTACGCGAAAAACTTACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTTGTGATCATAATTTTAAAGACTTTATTAATGTCCTCAATGTCGGCTAAGCCAATCTGGTGAACTAACCACTGACGCGTTTTCGCCATATTTTGTTTGTTAAAGCCCAGTTCCGCAATGGTTTTGCTGCTATTCTCGAACGCGATGTTTTg  <  2:285254/251‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                        AGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTTTGATGCCATTCGCGAACGCGATGTGTTG  >  NC_000913/2623604‑2624069

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