Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,742,829	T→C	F33L (TTC→CTC)	yiaK →	2,3‑diketo‑L‑gulonate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,742,829	0	T	C	100.0%	154.0 / NA	43	F33L (TTC→CTC)	yiaK	2,3‑diketo‑L‑gulonate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/23);  total (20/23)

AAAAAGCGATCAAAAAAACAAGGAAGCCTGGGATGAAAGTGACATTTGAGCAGTTAAAAGCAGCCTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATT  >  NC_000913/3742701‑3742960                                                                                                                                 |                                                                                                                                    aaaaaGCGATCAAAAAAACAAGGAAGCCTGGGATGAAAGTGACATTTGAGCAGTTAAAAGCAGCCTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGcacc                                                                                                                          >  2:62876/1‑140 (MQ=255)  aaaaGCGATCAAAAAAACAAGGAAGCCTGGGATGAAAGTGACATTTGAGCAGTTAAAAGCAGCCTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCC‑GTGCAGAGATGCTCGCCCGcaccac                                                                                                                        >  2:571408/1‑140 (MQ=255)         aTCAAAAAAACAAGGAAGCCTGGGATGAAAGTGACATTTGAGCAGTTAAAAGCAGCCTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATc                                                                                                                  >  2:255805/1‑140 (MQ=255)               aaaaCAAGGAAGCCTGGGATGAAAGTGACATTTGAGCAGTTAAAAGCAGCCTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGt                                                                                                            <  1:255805/140‑1 (MQ=255)                        aaGCCTGGGATGAAAGTGACATTTGAGCAGTTAAAAGCAGCCTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCa                                                                                                   <  1:487693/140‑1 (MQ=255)                          gCCTGGGATGAAAGTGACATTTGAGCAGTTAAAAGCAGCCTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACg                                                                                                 <  1:62876/140‑1 (MQ=255)                              gggATGAAAGTGACATTTGAGCAGTTAAAAGCAGCCTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGt                                                                                             <  1:164505/140‑1 (MQ=255)                              gggATGAAAGTGACATTTGAGCAGTTAAAAGCAGCCTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGt                                                                                             >  2:366802/1‑140 (MQ=255)                              gggATGAAAGTGACATTTGAGCAGTTAAAAGCAGCCTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGt                                                                                             <  2:169708/140‑1 (MQ=255)                                 aTGAAAGTGACATTTGAGCAGTTAAAAGCAGCCTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTaa                                                                                          <  2:823614/140‑1 (MQ=255)                                      aGTGACATTTGAGCAGTTAAAAGCAGCCTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGtt                                                                                     >  1:659824/1‑140 (MQ=255)                                        tGACATTTGAGCAGTTAAAAGCAGCCTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCC‑GTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTcc                                                                                  <  1:571408/140‑1 (MQ=255)                                          aCATTTGAGCAGTTAAAAGCAGCCTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTccc                                                                                 <  1:823671/140‑1 (MQ=255)                                          aCATTTGAGCAGTTAAAAGCAGCCTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTccc                                                                                 <  1:802948/140‑1 (MQ=255)                                                      ttAAAAGCAGCCTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTttcattca                                                                     >  1:100877/1‑140 (MQ=255)                                                                ccTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGaa                                                           >  2:178256/1‑140 (MQ=255)                                                                ccTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGaa                                                           >  1:79023/1‑140 (MQ=255)                                                                ccTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGaa                                                           >  1:782279/1‑140 (MQ=255)                                                                     aaTCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACgg                                                      <  1:178256/140‑1 (MQ=255)                                                                        cGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGa                                                   <  2:100877/140‑1 (MQ=255)                                                                          ggTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGc                                                     >  2:282539/1‑136 (MQ=255)                                                                          ggTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGc                                                     <  1:282539/136‑1 (MQ=255)                                                                          ggTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGata                                                 >  2:489778/1‑140 (MQ=255)                                                                            tCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATc                                               >  1:655146/1‑140 (MQ=255)                                                                             cTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCa                                              <  2:659824/140‑1 (MQ=255)                                                                                  tttCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCt                                         <  2:498571/140‑1 (MQ=255)                                                                                     cACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGAt                                      >  1:219486/1‑140 (MQ=255)                                                                                          ggcgTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGccca                                 >  1:743584/1‑140 (MQ=255)                                                                                          ggcgTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGccca                                 >  1:816258/1‑140 (MQ=255)                                                                                          ggcgTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGccca                                 >  2:304876/1‑140 (MQ=255)                                                                                           gcgTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaa                                <  1:304876/140‑1 (MQ=255)                                                                                           gcgTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaa                                >  2:309577/1‑140 (MQ=255)                                                                                               tGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGTGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaaccca                            <  2:782279/140‑1 (MQ=255)                                                                                               tGACAGCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaaccca                            <  1:59913/140‑1 (MQ=255)                                                                                                    gCGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGt                       <  2:655146/140‑1 (MQ=255)                                                                                                     cGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaaccca                            >  2:464169/1‑134 (MQ=255)                                                                                                     cGAAACGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaaccca                            <  1:464169/134‑1 (MQ=255)                                                                                                        aaCGGGTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATaa                   >  2:339331/1‑140 (MQ=255)                                                                                                         aCGGCTGACGCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAAc                  <  1:569444/140‑1 (MQ=255)                                                                                                                  gCCTGTGCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCgg         <  2:743584/140‑1 (MQ=255)                                                                                                                     tgtgCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGgcgc      >  2:513781/1‑140 (MQ=255)                                                                                                                       tgCAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGTGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCaa    <  1:299694/140‑1 (MQ=255)                                                                                                                         cAGAGATGCTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAAtt  <  1:438296/140‑1 (MQ=255)                                                                                                                                 |                                                                                                                                    AAAAAGCGATCAAAAAAACAAGGAAGCCTGGGATGAAAGTGACATTTGAGCAGTTAAAAGCAGCCTTTAATCGGGTCTTAATTTCACGCGGCGTTGACAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATT  >  NC_000913/3742701‑3742960

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