Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,644,623	T→C	N81S (AAC→AGC)	ndk ←	nucleoside diphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,644,623	0	T	C	100.0%	91.1 / NA	27	N81S (AAC→AGC)	ndk	nucleoside diphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/10);  total (17/10)

GTGGGTACCGTTTTCGGTCAGGCTGTCAGCGTAATCAGCGCGCAGAGTACCAGCCAGTGCGTTTGCCGGATTGGTCGCGCCCAGCAGATCGCGGTGACGCTGAACGGCGTTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCATCGTGTTCAGCATAAAAGCCACGTGCCTGTTCAACGGTCAGGTGC  >  NC_000913/2644514‑2644745                                                                                                              |                                                                                                                           gTGGGTACCGTTTTCGGTCAGGCTGTCAGCGTAATCAGCGCGCAGAGTACCAGCCAGTGCGTTTGCCGGATTGGTCGCGCCCAGCAGATCGCGGTGACGCTGGACGGCGCGTTCGCCGTCCAGCACGGAAAcc                                                                                                     >  1:234329/1‑133 (MQ=255)      tACCGTTTTCGGTCAGGCTGTCAGCGTAATCAGCGCGCAGAGTACCAGCCAGTGCGTTTGCCGGATTGGTCGCGCCCAGCAGATCGCGGTGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGAt                                                                                                <  2:461056/133‑1 (MQ=255)             ttCGGTCAGGCTGTCAGCGTAATCAGCGCGCAGAGTACCAGCCAGTGCGTTTGCCGGATTGGTCGCGCCCAGCAGATCGCGGTGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCa                                                                                         >  1:35373/1‑133 (MQ=255)                      gCTGTCAGCGTAATCAGCGCGCAGAGTACCAGCCAGTGCGTTTGCCGGATTGGTCGCGCCCAGCAGATCGCGGTGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATg                                                                                >  2:460031/1‑133 (MQ=255)                                      gcgcgcAGAGTACCAGCCAGTGCGTTTGCCGGATTGGTCGCGCCCAGCAGATCGCGGTGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCAt                                                                <  2:35373/133‑1 (MQ=255)                                        gcgcAGAGTACCAGCCAGTGCGTTTGCCGGATTGGTCGCGCCCAGCAGATCGCGGTGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCa                                                              >  2:484229/1‑133 (MQ=255)                                                    agccagTGCCGTTGCCGGATTGGTTGGGACCATCAGATCGCGGTGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGtttt                                                  <  2:234329/133‑1 (MQ=255)                                                        agTGCGTTTGCCGGATTGGTCGCGCCCAGCAGATCGCGGTGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCGGACCATCAAAGAACGGTTTTCCAt                                              >  1:18192/1‑133 (MQ=255)                                                        agTGCGTTTGCCGGATTGGTCGCGCCCAGCAGATCGCGGTGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGGCCATCAGAGAACGATTTTCCAt                                              >  1:252520/1‑133 (MQ=255)                                                        agTGCGTTTGCCGGATTGGTCGCGCCCAGCAGATCGCGGTGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCAt                                              >  1:302279/1‑133 (MQ=255)                                                          tGCGTTTGCCGGATTGGTCGCGCCCAGCAGATCGCGGTGACGCTGAACGGCGCTTTCACCTTCCAGTACGGAAACCACGATCGGGCCAGAGCTCATGAATTCAATCAGGACATCAAAGAACGGTTTTCCAtag                                            >  2:241918/1‑131 (MQ=255)                                                          tGCGTTTGCCGGATTGGTCGCGCCCAGCAGATCGCGGTGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCATCg                                            >  1:94864/1‑133 (MQ=255)                                                          tGCGTTTGCCGGATTGGTCGCGCCCAGCAGATCGCGGTGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCATCg                                            >  2:448782/1‑133 (MQ=255)                                                          tGCGTTTGCCGGATTGGTCGCGCCCAGCAGATCGCGGTGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCATCg                                            >  1:323603/1‑133 (MQ=255)                                                          tGCGTTTGCCGGATTGGTCGCGCCCAGCAGATCGCGGTGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCATCg                                            >  1:25994/1‑133 (MQ=255)                                                               ttGCCGGACTGGTCGAGCTCAGAAGATCGCGGTGACGCTGAACGGCGCTTTCACCTTCCAGTACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCATCGTGTTc                                       <  1:241918/133‑1 (MQ=255)                                                                  ccGGATTGGTGGCGCGCAGCAGATAGCGGTGAGGCTGAACGACGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCATCGTGTTCAGc                                    <  2:252520/133‑1 (MQ=255)                                                                             gcgcCCAGCAGATCGCGGGGACCCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCATCGTGTTCAGCATAAAAGCCAc                         <  1:229638/133‑1 (MQ=255)                                                                                         tCGCGGTGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCATCGTGTTCAGCATAAAAGCCACGTGCCTGTTCaa             <  1:484229/133‑1 (MQ=255)                                                                                          cgcgGTGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCATCGTGTTCAGCATAAAAGCCACGTGCCTGTTCAAc            >  1:65026/1‑133 (MQ=255)                                                                                          cgcgGTGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCATCGTGTTCAGCATAAAAGCCACGTGCCTGTTCAAc            >  2:471583/1‑133 (MQ=255)                                                                                             ggTGACGCTGGACGGGGCTTTCACCTTCAAGCACGGAAACCCCGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCATCGTGTTCAGCATAAAAGCCACGTGCCTGTTCAACGGt         <  2:196966/133‑1 (MQ=255)                                                                                              gTGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCATCGTGTTCAGCATAAAAGCCACGTACCTGTTCAACGGTc        <  2:323603/133‑1 (MQ=255)                                                                                              gTGACGCTGAACGGCGCTTCGGGATTCCAGCACTGGCAGCACGAACGGGACAGAGGTCATGAATTCAACCAGACCATCAAAGAACTGCTTTCCACCGTGTTCAGCATAAAAGCCACGTGCATGTGCAAGGGTc        >  2:124031/1‑133 (MQ=255)                                                                                               tGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCATCGTGTTCAGCATAAAAGCCACGTGCCTGTTCAACGGTCa       <  2:302279/133‑1 (MQ=255)                                                                                               tGACGCTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCATCGTGTTCAGCATAAAAGCCACGTGCCTGTTCAACGGTCa       >  2:62447/1‑133 (MQ=255)                                                                                                    cTGAACGGCGCTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCATCGTGTTCAGCATAAAAGCCACGTGCCTGTTCAACGGTCAGGTgc  >  1:390681/1‑133 (MQ=255)                                                                                                              |                                                                                                                           GTGGGTACCGTTTTCGGTCAGGCTGTCAGCGTAATCAGCGCGCAGAGTACCAGCCAGTGCGTTTGCCGGATTGGTCGCGCCCAGCAGATCGCGGTGACGCTGAACGGCGTTTTCACCTTCCAGCACGGAAACCACGATCGGGCCAGAGGTCATGAATTCAACCAGACCATCAAAGAACGGTTTTCCATCGTGTTCAGCATAAAAGCCACGTGCCTGTTCAACGGTCAGGTGC  >  NC_000913/2644514‑2644745

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