Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,617,976	T→C	I218I (ATT→ATC)	deoC →	deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,617,976	0	T	C	100.0%	76.6 / NA	29	I218I (ATT→ATC)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/14);  total (15/14)

GTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTG‑AAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGG  >  NC_000913/4617734‑4618220                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      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)                                                                                                                                                                                                            cgtgcgtACTGCGGAAGAGGCGCAGAAAGATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTTCTTCCAGCCTGCTGGCAAGCCTGTTG‑AAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATc                                   <  2:64804/251‑1 (MQ=255)                                                                                                                                                                                                                 gtACTGCGGAAGATGCGCAGAAATATCTCGCCATCGGAAATGAACTGTTCCGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTG‑AAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACTCCTGATTCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCTTAGGCCGGATAAGGCGTAGACGCATCCGGCa                              <  1:24663/251‑1 (MQ=255)                                                                                                                                                                                                                 gtACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATCAAATGTACGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTAGGTGCTTCCAGCCTGATGGCAAGCCTGTTG‑CAAGCGCTGGGCCACGGAGACGGTAAGAGCGCCAGCCGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCCGGCCTACGAGACGTATCTACCCGTCGGCCGGATAAGGCGTAGACGCATCCGGc                               >  1:160517/1‑250 (MQ=255)                                                                                                                                                                                                                         gAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTG‑AAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGCCGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAgccgc                      >  1:118536/1‑251 (MQ=255)                                                                                                                                                                                                                                       atatCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTG‑AAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGACGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTc        >  1:227305/1‑251 (MQ=255)                                                                                                                                                                                                                                        tatCTCGCCATCGCAGATGACCTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTG‑ACAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCAGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCCAAAGCCGCCTCATACTCTTTTcc       >  2:192610/1‑251 (MQ=255)                                                                                                                                                                                                                                             cGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTG‑AAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCggg  >  1:5179/1‑251 (MQ=255)                                                                                                                                                                                                                                               ccATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTG‑AAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCggg  >  1:172686/1‑249 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      GTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTG‑AAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGG  >  NC_000913/4617734‑4618220

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