Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,623,319	C→T	D92D (GAC→GAT)	ppk →	polyphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,623,319	0	C	T	100.0%	89.3 / NA	26	D92D (GAC→GAT)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/12);  total (14/12)

TATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGCGTCATTATTTTAAGCAGTATCTGCGTCAGCACATTACGC  >  NC_000913/2623191‑2623452                                                                                                                                 |                                                                                                                                      tataAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACaa                                                                                                                            <  2:49408/140‑1 (MQ=255)           gCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAg                                                                                                                        >  1:40410/1‑134 (MQ=255)           gCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAg                                                                                                                        <  2:40410/134‑1 (MQ=255)           gCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTg                                                                                                                  >  1:332747/1‑140 (MQ=255)             ttCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCt                                                                                                                >  1:300741/1‑140 (MQ=255)                       tgaaGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGgca                                                                                                      <  2:332747/140‑2 (MQ=255)                         aaGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGgcgcg                                                                                                    >  1:266457/1‑140 (MQ=255)                            cGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCaa                                                                                                 <  2:70645/140‑1 (MQ=255)                             gACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAAc                                                                                                <  1:272331/140‑1 (MQ=255)                                   tcatTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCaa                                                                                                 >  2:288873/1‑133 (MQ=255)                                   tcatTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCaa                                                                                                 <  1:288873/133‑1 (MQ=255)                                             cGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTa                                                                                >  2:86723/1‑140 (MQ=255)                                               aagaaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTc                                                                                       <  2:259416/131‑1 (MQ=255)                                               aagaaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTc                                                                                       >  1:259416/1‑131 (MQ=255)                                                agaaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATg                                                                             >  2:34745/1‑140 (MQ=255)                                                     aaGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGc                                                                        <  1:325336/140‑1 (MQ=255)                                                     aaGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGc                                                                        <  2:300741/140‑1 (MQ=255)                                                          tCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGct                                                                   <  1:86723/140‑1 (MQ=255)                                                            cAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGctct                                                                 >  2:348968/1‑140 (MQ=255)                                                                          gCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATcaacaa                                                   >  1:351569/1‑140 (MQ=255)                                                                           ccATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACaaa                                                  >  2:61838/1‑140 (MQ=255)                                                                                     ggCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGCg                                        >  1:176788/1‑140 (MQ=255)                                                                                      gCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACCGGCTGCGt                                       >  1:130345/1‑140 (MQ=255)                                                                                         aaaTTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGCGTCa                                     <  2:263481/139‑1 (MQ=255)                                                                                         aaaTTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGCGTCa                                     >  1:263481/1‑139 (MQ=255)                                                                                                                           aTTCGATGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCGCTCCGTCAATCAACAAAACTGGCTGCGTCATTATTTTAAGCAGTATCTGCGTCAGCACATTACGc  <  1:366653/140‑1 (MQ=255)                                                                                                                                 |                                                                                                                                      TATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGCGTCATTATTTTAAGCAGTATCTGCGTCAGCACATTACGC  >  NC_000913/2623191‑2623452

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