Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,617,616	T→C	A98A (GCT→GCC)	deoC →	deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,617,616	0	T	C	100.0%	76.7 / NA	25	A98A (GCT→GCC)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/13);  total (12/13)

ACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCC  >  NC_000913/4617369‑4617858                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     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)                                                                                                                                                                                                                              aaCCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGGC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCCGAGAATGTACTGCTG‑AAAGTGATCATCGGAACCGGAGAACTGAAAGACGAAGCGCTGGTCCGTAAAACGTCTGAAATCTCCCTCAAAACGGGTGCGGACTTCATCAAAACCCCTACCCGTaaa                      >  1:96097/1‑250 (MQ=255)                                                                                                                                                                                                                                    ggCGGCAAGCGCCTGCGGGTCCGATGAAGTTGACGTGTTGTTCCCGGACCGCGCGCAGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGCATGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCTTAAAGCGTCTGAAATCTCCATCAAAGCGTGTGCGGACTTCATCAAAACCTCTTCCGGTAAAGTGGCtg               <  2:175052/250‑1 (MQ=255)                                                                                                                                                                                                                                          aaTCGCCTACGGTGCCGATGAAGTTGACGT‑GTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAgc                                                                                            <  1:48977/167‑1 (MQ=255)                                                                                                                                                                                                                                                 tACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTG‑GAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCCTCCAAGCGGGGGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCgccgcc  >  1:354218/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     ACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCC  >  NC_000913/4617369‑4617858

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