Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,644,487	A→C	A126A (GCT→GCG)	ndk ←	nucleoside diphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,644,487	0	A	C	100.0%	87.1 / NA	26	A126A (GCT→GCG)	ndk	nucleoside diphosphate kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (15/11);  total (15/11)

GCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGAGCGGCAGATTCGACGGAATCAGAACCGTGGGTACCGTTTTCGGTCAGGCTGTCAGCGTAATCAGCGCGCAGAGTACC  >  NC_000913/2644411‑2644564                                                                             |                                                                              gCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCgg                                                                           >  1:472960/1‑81 (MQ=255)    cgcgGCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCgg                                                                           >  2:472962/1‑78 (MQ=255)    cgcgGCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCgg                                                                           >  1:472964/1‑78 (MQ=255)    cgcgGCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCGGCAg                                                                        <  1:472919/81‑1 (MQ=255)     gcgGCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCGGCAGa                                                                       <  2:472920/81‑1 (MQ=255)       ggCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCGGCAGAtt                                                                     >  2:472968/1‑81 (MQ=255)           tttACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCGGCAGATTCGAc                                                                 <  2:472921/81‑1 (MQ=255)            ttACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCGGCAGATTCGACg                                                                >  2:472976/1‑81 (MQ=255)               cGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCGGCAGATTCGACGGaa                                                             >  1:472974/1‑81 (MQ=255)                        taACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCGGCAGATTCGACGGAATCAGAACCg                                                    >  1:472963/1‑81 (MQ=255)                          aCGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCGGCAGATTCGACGGAATCAGAACCGTg                                                  <  1:472922/81‑1 (MQ=255)                              gTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCGGCAGATTCGACGGAATCAGAACCGTGGGTa                                              <  1:472923/81‑1 (MQ=255)                              gTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCGGCAGATTCGACGGAATCAGAACCGTGGGTa                                              >  2:472972/1‑81 (MQ=255)                                gcgcGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCGGCAGATTCGACGGAATCAGAACCGTGGGTAcc                                            >  2:472981/1‑81 (MQ=255)                                 cgcgGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCGGCAGATTCGACGGAATCAGAACCGTGGGTACCg                                           >  2:472994/1‑81 (MQ=255)                                  gcgGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCGGCAGATTCGACGGAATCAGAACCGTGGGTACCGt                                          <  2:472924/81‑1 (MQ=255)                                    gggCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCGGCAGATTCGACGGAATCAGAACCGTGGGTACCGttt                                        >  1:472978/1‑81 (MQ=255)                                     ggCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCGGCAGATTCGACGGAATCAGAACCGTGGGTACCGtttt                                       <  1:472925/81‑1 (MQ=255)                                      gCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGCGCGGCAGATTCGACGGAATCAGAACCGTGGGTACCGTTTTc                                      <  1:472926/81‑1 (MQ=255)                                                ccttcgccAAAGAAATAAGCGATTTCGCGCGCGGCAGATTCGACGGAATCAGAACCGTGGGTACCGTTTTCGGTCAGGCTg                            <  1:472927/81‑1 (MQ=255)                                                 cttcgccAAAGAAATAAGCGATTTCGCGCGCGGCAGATTCGACGGAATCAGAACCGTGGGTACCGTTTTCGGTCAGGCTGt                           >  2:472977/1‑81 (MQ=255)                                                      ccAAAGAAATAAGCGATTTCGCGCGCGGCAGATTCGACGGAATCAGAACCGTGGGTACCGTTTTCGGTCAGGCTGTCAGCg                      <  1:472928/81‑1 (MQ=255)                                                         aaGAAATAAGCGATTTCGCGCGCGGCAGATTCAACGGAATCAGAACCGTGGGTACCGTTTTCGGTCAGGCTGTCAGCGTaa                   <  2:472929/81‑1 (MQ=255)                                                                     aTTTCGCGCGCGGCAGATTCGACGGAATCAGAACCGTGGGTACCGTTTTCGGTCAGGCTGTCAGCGTAATCAGCGCGCa         >  2:472971/1‑79 (MQ=255)                                                                       ttCGCGCGCGGCAGATTCGACGGAATCAGAACCGTGGGTACCGTTTTCGGTCAGGCTGTCAGCGTAATCAGCGCGCagag      >  1:472965/1‑80 (MQ=255)                                                                       ttCGCGCGCGGCAGATTCGACGGAATCAGAACCGTGGGTACCGTTTTCGGTCAGGCTGTCAGCGTAATCAGCGCGCagag      >  2:472982/1‑80 (MQ=255)                                                                          gcgcgcGGCAGATTCGACGGAATCAGAACCGTGGGTACCGTTTTCGGTCAGGCTGTCAGCGTAATCAGCGCGCAGAGTAcc  <  2:472932/81‑1 (MQ=255)                                                                          gcgcgcGGCAGATTCGACGGAATCAGAACCGTGGGTACCGTTTTCGGTCAGGCTGTCAGCGTAATCAGCGCGCAGAGTAcc  <  2:472930/81‑1 (MQ=255)                                                                          gcgcgcGGCAGATTCGACGGAATCAGAACCGTGGGTACCGTTTTCGGTCAGGCTGTCAGCGTAATCAGCGCGCAGAGTAc   <  1:472931/80‑1 (MQ=255)                                                                             |                                                                              GCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGAGCGGCAGATTCGACGGAATCAGAACCGTGGGTACCGTTTTCGGTCAGGCTGTCAGCGTAATCAGCGCGCAGAGTACC  >  NC_000913/2644411‑2644564

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