Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,342,683	T→C	G137G (GGT→GGC)	kdsC →	3‑deoxy‑D‑manno‑octulosonate 8‑phosphate phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,342,683	0	T	C	100.0%	186.1 / NA	49	G137G (GGT→GGC)	kdsC	3‑deoxy‑D‑manno‑octulosonate 8‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/23);  total (26/23)

CACTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGTTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTGGCGCAGGGCAAACTGGATG  >  NC_000913/3342548‑3342816                                                                                                                                        |                                                                                                                                      cactcactTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTa                                                                                                                                    >  2:214278/1‑139 (MQ=255) cactcactTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTa                                                                                                                                    >  1:126585/1‑139 (MQ=255)  actcactTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGACCTGCGGGGAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAATTGGGCTTaa                                                                                                                                   >  1:507099/1‑139 (MQ=255)     cactTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCg                                                                                                                                <  2:535002/139‑1 (MQ=255)     cactTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCg                                                                                                                                <  1:498064/139‑1 (MQ=255)      actTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGt                                                                                                                               >  1:677185/1‑139 (MQ=255)        ttGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCACCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCg                                                                                                                             >  1:764436/1‑139 (MQ=255)           tATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCg                                                                                                                          >  1:1031236/1‑139 (MQ=255)                ggggCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGcc                                                                                                                     <  2:840900/139‑1 (MQ=255)                ggggCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGcc                                                                                                                     >  2:378616/1‑139 (MQ=255)                ggggCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGcc                                                                                                                     <  2:677185/139‑1 (MQ=255)                  ggCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATTGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGa                                                                                                                   >  1:310029/1‑139 (MQ=255)                  ggCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGa                                                                                                                   >  2:600283/1‑139 (MQ=255)                  ggCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGa                                                                                                                   >  2:214592/1‑139 (MQ=255)                       tCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATgcgc                                                                                                              <  2:795025/139‑1 (MQ=255)                       tCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATgcgc                                                                                                              >  1:84804/1‑139 (MQ=255)                       tCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATgcgc                                                                                                              <  1:151453/139‑1 (MQ=255)                           acaaacTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATcc                                                                                                          <  2:691505/139‑1 (MQ=255)                                     cGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATc                                                                                                >  1:98104/1‑139 (MQ=255)                                         tttAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCcgc                                                                                            >  2:1111221/1‑139 (MQ=255)                                                aTCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGa                                                                                     >  2:55094/1‑139 (MQ=255)                                                  ctgctgGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATcc                                                                                                          >  2:1079017/1‑116 (MQ=255)                                                  ctgctgGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATcc                                                                                                          <  1:1079017/116‑1 (MQ=255)                                                  ctgctgGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGAtt                                                                                   >  2:254790/1‑139 (MQ=255)                                                       ggAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTg                                                                              <  1:562866/139‑1 (MQ=255)                                                               tGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGcc                                                                                                                     <  2:774666/92‑1 (MQ=255)                                                               tGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGcc                                                                                                                     >  1:774666/1‑92 (MQ=255)                                                                ggCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCAtt                                                                     >  2:428149/1‑139 (MQ=255)                                                                       gCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCg                                                              >  1:612194/1‑139 (MQ=255)                                                                                 aTGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGgcgc                                                    <  2:891726/139‑1 (MQ=255)                                                                                 aTGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGgcgc                                                    <  2:84804/139‑1 (MQ=255)                                                                                  tgtgGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCa                                                   <  2:612194/139‑1 (MQ=255)                                                                                  tgtgGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTAGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCa                                                   >  2:122931/1‑139 (MQ=255)                                                                                                gCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATgc                                                                                                                >  2:298003/1‑64 (MQ=255)                                                                                                gCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATgc                                                                                                                <  1:298003/64‑1 (MQ=255)                                                                                                gCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTg                                     <  1:428149/139‑1 (MQ=255)                                                                                                gCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTg                                     >  1:869128/1‑139 (MQ=255)                                                                                                        ctcATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACttat                             <  1:1066181/139‑1 (MQ=255)                                                                                                        ctcATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACttat                             <  1:378616/139‑1 (MQ=255)                                                                                                                   gCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTg                                     <  2:161411/120‑1 (MQ=255)                                                                                                                   gCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTg                                     >  1:161411/1‑120 (MQ=255)                                                                                                                   gCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTGgcg                  >  1:809136/1‑139 (MQ=255)                                                                                                                       gTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTGGCGCAgg              <  2:507099/139‑1 (MQ=255)                                                                                                                                  aGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCa                                                                       >  2:464961/1‑71 (MQ=255)                                                                                                                                  aGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCa                                                                       <  1:464961/71‑1 (MQ=255)                                                                                                                                  aGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGcgcg                                            >  1:902722/1‑98 (MQ=255)                                                                                                                                  aGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGcgcg                                            <  2:902722/98‑1 (MQ=255)                                                                                                                                  aGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTGGCGCAGGGCAAACTGGAt   >  2:971268/1‑139 (MQ=255)                                                                                                                                   gTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTGGCGCAGGGCAAACTGGATg  <  2:485865/139‑1 (MQ=255)                                                                                                                                   gTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTGGCGCAGGGCAAACTGGATg  <  1:600283/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      CACTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGTTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTGGCGCAGGGCAAACTGGATG  >  NC_000913/3342548‑3342816

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