Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,903,363	T→C	S490P (TCG→CCG)	ygcE →	putative sugar kinase YgcE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,903,363	0	T	C	100.0%	102.5 / NA	30	S490P (TCG→CCG)	ygcE	putative sugar kinase YgcE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/13);  total (17/13)

ATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCC  >  NC_000913/2903232‑2903490                                                                                                                                    |                                                                                                                                aTGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATAGGGGAATGTGGTTTTGATCCAGAGCTGGCCGAATGCg                                                                                                                         >  1:282232/1‑140 (MQ=255)            tACGACGCGATGAATAAAGGCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATGCAGAGCGGGCCGAATGCGTAAGTGAGGaa                                                                                                              >  1:318867/1‑140 (MQ=255)               gacgCGATGAATAAAGGCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGc                                                                                                           <  1:167456/140‑1 (MQ=255)                      tGAATAAAGGCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCcgggcgg                                                                                                    >  1:634545/1‑140 (MQ=255)                         aTAAAGGCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGaa                                                                                                 <  2:634545/140‑1 (MQ=255)                             aGGCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTg                                                                                             >  2:482671/1‑140 (MQ=255)                                  ttttCAACGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCgg                                                                                        >  2:500900/1‑140 (MQ=255)                                   tttCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGc                                                                                       <  1:482671/140‑1 (MQ=255)                                   tttCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGc                                                                                       <  1:500900/140‑1 (MQ=255)                                               aaCCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCt                                                                           >  1:60215/1‑140 (MQ=255)                                               aaCCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGAAAGGCCGCGCGGGATACTGCCCGGCCTGAACATACCt                                                                           >  1:352588/1‑140 (MQ=255)                                                          cTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGAATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATc                                                                >  1:547030/1‑140 (MQ=255)                                                           tGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGGAGGCCGGGAGGGGAACGGCCCGGCCTGAACATACCTGAATGGTTATcc                                                               >  2:321751/1‑140 (MQ=255)                                                                 gATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTg                                                         >  1:128106/1‑140 (MQ=255)                                                                  aTCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGa                                                        >  1:76881/1‑140 (MQ=255)                                                                      tGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGAcgcg                                                    <  1:95670/140‑1 (MQ=255)                                                                      tGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGAATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGAcgcg                                                    <  2:547030/140‑1 (MQ=255)                                                                       gAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCgg                                                   >  2:22214/1‑140 (MQ=255)                                                                                      aGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGGCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGa                                    <  1:501553/140‑1 (MQ=255)                                                                                               tggggaattggggaatGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTaaa                           >  2:428839/1‑140 (MQ=255)                                                                                               tggggaattggggaatGCGGATGCGATCCAGCGCTAGCCCAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATGCCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTaaa                           <  2:352588/140‑1 (MQ=255)                                                                                                ggggaattggggaatGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTaaaa                          >  2:831166/1‑140 (MQ=255)                                                                                                     attggggaatGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGa                     >  2:577708/1‑140 (MQ=255)                                                                                                        ggggaatGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATaa                  <  1:683616/140‑1 (MQ=255)                                                                                                             atGTGGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGt             <  2:60215/140‑1 (MQ=255)                                                                                                              tgtgGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTc            <  1:577708/140‑1 (MQ=255)                                                                                                              tgtgGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTc            <  2:76881/140‑1 (MQ=255)                                                                                                               gtgGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGt             <  1:683661/138‑1 (MQ=255)                                                                                                               gtgGATTCGATCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGt             >  2:683661/1‑138 (MQ=255)                                                                                                                        aTCCAGAGCTGGCCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTACGACAACTTAAGAAATGATAAAAAGTCCGGTTTATcc  >  1:413801/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                ATGCCAAACGCTACGACGCGATGAATAAAGGCATTTTCAAAGACCTAACCAAACATACTGATGTGATCCTGAAAAAATCGTATGAAGTGATGCATGGGGAATTGGGGAATGTGGATTCGATCCAGAGCTGGTCGAATGCGTAAGTGAGGAAGGCCGGGCGGGAAACTGCCCGGCCTGAACATACCTGAATGGTTATCCCCGCTGACGCGGGGAACATAAGTTCGACAACTTAAAAAATGATAACAAGTCCGGTTTATCC  >  NC_000913/2903232‑2903490

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