Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,617,574	C→A	D84E (GAC→GAA)	deoC →	deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,617,574	0	C	A	100.0%	66.2 / NA	23	D84E (GAC→GAA)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (10/13);  total (10/13)

CTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCAT  >  NC_000913/4617332‑4617819                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        cTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCgcgc                                                                                                                                                                                                                                                  >  1:366178/1‑250 (MQ=255)                   cGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGc                                                                                                                                                                                                                               >  2:112977/1‑251 (MQ=255)                       cACTGAAATTGAGGGCCCTGCCCCCCCTGAAGGACGACGACACCGACGAGAATGTAATTGCTTTGTTTCTCCGGGCCAAAACCCCGG‑TCGGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc                                                                                                                                                                                                                           <  1:221741/251‑1 (MQ=255)                                      ccatgcccccccTGAATGACGACGACCCCGACGAAAAAGTAATTGCTCTGTGTCATCGGGCCAAACCCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCcgat                                                                                                                                                                                                            <  1:112977/242‑1 (MQ=255)                                                                cccGACGATAAAGTAATTGCGCTGTGACATCAGCCCAAACCCCCGG‑TCGCAAATCCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACGTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCTTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACcg                                                                                                                                                                                  <  1:152271/250‑1 (MQ=255)                                                                  cgacgaGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACcgcg                                                                                                                                                                                >  1:271164/1‑251 (MQ=255)                                                                     cgaGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCgcgcgc                                                                                                                                                                             >  1:344025/1‑251 (MQ=255)                                                                                            caccagtcAAAACCCCGGGTCGGCAATTCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAggt                                                                                                                                                       <  1:262823/244‑1 (MQ=255)                                                                                                                 ggCAATACCGCCGCTATCTGTATCTATCCTAGATTTCTCCCGTTTGCTCGCAAAACACTGCACGAGCAGGCCCCCCAGGAAATCCGTATTCCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGAGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGCCGTTGACGTGGTGTTCCCGTACCGCGCGCTGCTGGCGGTTTACGCGCAggttcgcttttacctctgcgactc                                                                                                                                  >  2:182439/1‑231 (MQ=255)                                                                                                                                atcggTATCTAACCTCGCTTTATCCCGCTTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAAACCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAAGCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGAGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTg                                                                                                                  <  2:262054/247‑1 (MQ=255)                                                                                                                                      atctatcCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGATGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTGAAGCGCAGGTTGGGTTT‑ACACTGGTGAAAGCCttgtaagaggcgaggcgagg                                                                                                            >  2:64970/1‑232 (MQ=255)                                                                                                                                               cGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCggaagt                                                                                                   >  1:382374/1‑246 (MQ=255)                                                                                                                                                   ttATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAATGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTAct                                                                                                <  1:338816/250‑1 (MQ=255)                                                                                                                                                           aatgctcgataaccactGAAAGAGCAGGGCACCCCAGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGGAACGATGACATCGAAAGCGCGCTGGCAGAAACCCGTGCGGCAAGCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCATCGAATGTACTGCTG‑AAAga                                                                                      <  2:120141/236‑2 (MQ=255)                                                                                                                                                              gCTCGCAAAACACTGAAAGAGCCGGGCACCCCGGAAATCCGTATTGCGACGGTAACCAACTGCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTG‑AAAGTGa                                                                                    <  2:245686/250‑1 (MQ=255)                                                                                                                                                                     aaaCACTGAAAGAGCAGGGCACCCCGGAAATCCGTAGTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTG‑AAAGTGATCATCGaa                                                                            <  1:182439/251‑1 (MQ=255)                                                                                                                                                                      aacactGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTG‑AAAGTGATCATCGaaa                                                                           <  1:231837/247‑1 (MQ=255)                                                                                                                                                                                             cGGAAATCCGTTTTTCTACGGTAACCAACCTCCCACACGGTTACGAGGACCTCGAAATCGCTCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAgcg                                                    <  1:263416/251‑1 (MQ=255)                                                                                                                                                                                                 aaTCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAAAGAGCAGGCTGGTTTTGAC‑CTGGGGAAAGCCTGTAAAGAGGCTTGCGCGGCCGCGCATGTACTGATG‑GAAGGGAACATCGAAACCGGCGCACTggaacacgcacggctg                                                 >  1:367953/1‑235 (MQ=255)                                                                                                                                                                                                  aTCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGAt                                               <  1:211004/251‑1 (MQ=255)                                                                                                                                                                                                    ccgtactgcgtcGGTAACCAACTTCCCACACGGTAATGATGACATCGAAATCGCGCTGGCAGAAACACTTGCGGCAATCTCCTGCTGTTCCGATGAAGTTGACGTGGTGTTGCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGTGTTGGC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGGTG‑AAAGTGATCATCGAAACCGGCGAACTGAAAGGCGAATCGCTGATcc                                             <  2:293604/240‑1 (MQ=255)                                                                                                                                                                                                                         cTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTG‑AAAGTGATCATCGAAAACGGCGCACTGAAAGACGAAACGCTGATCCGTAAAGCGACTGAAATCTcc                         >  1:53209/1‑250 (MQ=255)                                                                                                                                                                                                                             ccACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAA‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGGCGCGAATGTACTGCTG‑AAAGAGGTCATCGGAAACGGCGAAATGAAAAACGAAGCGCTGATCCGGAAAGCGGCTGAAATCTCCCTCa                     >  1:239719/1‑250 (MQ=255)                                                                                                                                                                                                                                        cgatgacatcgaaatcGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGGTTTGAA‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGGACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTATGAAATCTCCATCAAAGCGGGGGCgg         >  1:72879/5‑251 (MQ=255)                                                                                                                                                                                                                                               atcgaaatcgCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGCAAAG‑GGTCATCGAAACCGGCGCACTGAAAAACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCACTCCACGCGGGTGAGGACTtcat  >  1:257366/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        CTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCAT  >  NC_000913/4617332‑4617819

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