Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,617,928	G→A	P202P (CCG→CCA)	deoC →	deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,617,928	0	G	A	100.0%	79.4 / NA	26	P202P (CCG→CCA)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/15);  total (11/15)

CCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTG‑AAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGG  >  NC_000913/4617682‑4618165                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               ccTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACcagcgg                                                                                                                                                                                                                                            <  2:170452/251‑1 (MQ=255)                       gCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGa                                                                                                                                                                                                                         >  1:252778/1‑248 (MQ=255)                       gCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGa                                                                                                                                                                                                                      <  1:269486/251‑1 (MQ=255)                          tGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGGAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTGCCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTGCAAACGAGCGGGCGGCGTGCGAACTGCGGAAGAt                                                                                                                                                                                                                     <  2:124166/249‑1 (MQ=255)                              cgGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCAGGACATCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATgcgc                                                                                                                                                                                                                 >  1:12707/1‑249 (MQ=255)                                              gctgAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCAt                                                                                                                                                                                                 <  1:222595/249‑1 (MQ=255)                                                     gTGATCATAGAAAACGACGGACTGCAAGACGACGCGCTGATCAGTAAAGCGTCTGAAATATCCATAAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGGGCACGCGACGCCGGTGAGCGCGCGCATTATGATGGTAGTGATACGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCCGTGGGCGGGAGTAATGCAGAAGATGCGCAGAAATATCTCGCCATCGCAGaat                                                                                                                                                                                         >  2:229893/1‑248 (MQ=255)                                                      tGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGa                                                                                                                                                                                        >  1:235873/1‑250 (MQ=255)                                                                            gAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTGCGGTGCTGACTGGGCAg                                                                                                                                                                 <  2:226463/251‑1 (MQ=255)                                                                                 aCGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATgcg                                                                                                                                                            >  2:207987/1‑251 (MQ=255)                                                                                           gATCCGTAAAGCGTCTGAAATCTCCATCACAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTAcc                                                                                                                                                  >  1:88928/1‑251 (MQ=255)                                                                                                            aaaTCTCCATCAAAGCGGGTGCGGACTTCATCAAATCCTCTTCCGGTAAAGTGGCTGTTATCGCGAGGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACGGCGGAAGATGCGCATAAATATCTCGCCATCGCAGAGGAGCTGTTCGGTGCTGACTGGGCAGATGCGCGGCTCTACCGCTTTGGTGCTTCCAGc                                                                                                                                 <  2:130476/251‑1 (MQ=255)                                                                                                                                  ggACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGAGATGGGCGTAGAAAAAGCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTAcc                                                                                                                                                  <  1:156227/212‑1 (MQ=255)                                                                                                                                  ggACTACATCAAAACCTATAACGGTAACGTGGATGTGAACGCGACGCCGCAAAGCGCGCGCTTCATGATGGAAGTGATCCGTGCTATGGGCGTAGAAAAAACCGTTGGTTTCACACCAGCGGGCGGCGTGCGTACTGCGGAAGCTGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGTTGCGCGTCACTAcc                                                                                                                                                  >  2:156227/1‑212 (MQ=255)                                                                                                                                                                     tgAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTG‑AAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCagcagc                                                                       >  2:205030/1‑251 (MQ=255)                                                                                                                                                                                  ggAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGACAAAACCGTTGGTTTCAAACAAGCGGGCGGCGTGCGTCCTGCGGAAGATGCGCAGAAATACCTCGCCATCGCAGCTGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCCAGCCTGTTG‑AAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACtaagta                                                              >  1:129171/1‑247 (MQ=255)                                                                                                                                                                                   gAAAGCGCGCGCATCATGATGGTAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGGCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTG‑AAAGCGCTGGGCCACGTTGACGGTAAGAGCGCCAGCAGATACTAAGTAAGATg                                                         <  2:55223/251‑1 (MQ=255)                                                                                                                                                                                    aaaGCGCGCGCATCATGATGGAAGCGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCTCTACCGCTTGGGTGCTGCCAGCCTGCGGGCAAGCCT‑TTGTAAAGCGCAGGGCCACTGTGACGGTGAGAGCGCCAGCAGCTACTAAGTAAGATg                                                         <  2:257429/250‑1 (MQ=255)                                                                                                                                                                                                 catGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTG‑AAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtg                                           <  1:201688/251‑1 (MQ=255)                                                                                                                                                                                                     atgtaaGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTG‑AAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcgct                                       <  1:244635/247‑1 (MQ=255)                                                                                                                                                                                                       ggAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTG‑AAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcgctgc                                     <  1:56678/251‑1 (MQ=255)                                                                                                                                                                                                       ggAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTG‑AAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcgctg                                      <  1:11707/250‑1 (MQ=255)                                                                                                                                                                                                           gTGATCCGGGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTG‑AAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCtt                                 <  2:148676/251‑1 (MQ=255)                                                                                                                                                                                                                        tGGGCGGAGGAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTAATGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTG‑AAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACg                     >  2:24086/1‑250 (MQ=255)                                                                                                                                                                                                                              tAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGAGCTTCCAGCCTGCTGGCAAGCCTGTTG‑ACAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGt               >  1:24665/1‑250 (MQ=255)                                                                                                                                                                                                                                         gTTGGTTTCAACCCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTG‑AAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCGTATCAGGCCTACGAGACGTATCTACCCGTAg   <  1:96969/251‑1 (MQ=255)                                                                                                                                                                                                                                           tgggttCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGAGGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGTCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGTCAAGCCTGTTG‑AAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGTGCGTATCAGGCCTACGAGACGTAGCTACCCGTAgg  <  2:28562/246‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               CCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTG‑AAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGG  >  NC_000913/4617682‑4618165

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