Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,358,589	A→G	V86A (GTC→GCC)	rhmA ←	2‑keto‑3‑deoxy‑L‑rhamnonate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,358,589	0	A	G	100.0%	101.9 / NA	30	V86A (GTC→GCC)	rhmA	2‑keto‑3‑deoxy‑L‑rhamnonate aldolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/12);  total (18/12)

TCCCCAGCGCGCAGCCCGTGCCACACTGGCCCCGACACCACGCTCACCGTAGGGAGGATAGCGCGTGGCAGACACCACCTGACGTGCCTGTTCGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGACTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCAGCAACCAGTCATAACCAGAA  >  NC_000913/2358447‑2358735                                                                                                                                               |                                                                                                                                                   tCCCCAGCGCGCAGCCCGTGCCACACTGGCCCCGACACCACGCTCACCGTAGGGAGGATAGCGCGTGGCAGACACCACCTGACGTGCCTGTTCGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGttt                                                                                                                                             >  2:624229/1‑150 (MQ=255)    ccAGCGCGCAGCCCGTGCCACACTGGCCCCGACACCACGCTCACCGTAGGGAGGATAGCGCGTGGCAGACACCACCTGACGTGCCTGTTCGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATc                                                                                                                                         <  1:437323/151‑1 (MQ=255)             aGCCCGTGCCACACTGGCCCCGACACCACGCTCACCGTAGGGAGGATAGCGCGTGGCAGACACCACCTGACGTGCCTGTTCGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGttt                                                                                                                                 >  1:371796/1‑150 (MQ=255)               cccGTGCCACACTGGCCCCGACACCACGCTCACCGTAGGGAGGATAGCGCGTGGCAGACACCACCTGACGTGCCTGTTCGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACt                                                                                                                              >  2:550556/1‑151 (MQ=255)                      cacaCTGGCCCCGACACCACGCTCACCGTAGGGAGGATAGCGCGTGGCAGACACCACCTGACGTGCCTGTTCGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTcc                                                                                                                       >  2:363205/1‑151 (MQ=255)                          cTGGCCCCGACACCACGCTCACCGTAGGGAGGATAGCGCGTGGCAGACACCACCTGACGTGCCTGTTCGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCAcc                                                                                                                    <  1:437929/150‑1 (MQ=255)                               cccGACACCACGCTCACCGTAGGGAGGATAGCGCGTGGCAGACACCACCTGACGTGCCTGTTCGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACcggacgg                                                                                                              <  1:388240/151‑1 (MQ=255)                               cccGACACCACGCTCACCGTAGGGAGGATAGCGCGTGGCAGACACCACCTGACGTGCCTGTTCGACAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACcggacgg                                                                                                              <  1:24036/151‑1 (MQ=255)                                  gACACCACGCTCACCGTAGGGAGGATAGCGCGTGGCAGACACCACCTGACGTGCCTGTTCGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATc                                                                                                           >  1:46521/1‑151 (MQ=255)                                     accacGCTCACCGTAGGGAGGATAGCGCGTGGCAGACACCACCTGACGTGCCTGTTCGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACg                                                                                                        >  2:560522/1‑151 (MQ=255)                                       cacGCTCACCGTAGGGAGGATAGCGCGTGGCAGACACCACCTGACGTGCCTGTTCGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATGTCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACggg                                                                                                      <  2:142294/151‑1 (MQ=255)                                         cGCTCACCGTAGGGAGGATAGCGCGTGGCAGACACCACCTGACGTGCCTGTTCGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGtt                                                                                                    >  1:271319/1‑151 (MQ=255)                                                 gTAGGGAGGATAGCGCGTGGCAGACACCACCTGACGTGCCTGTTCGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCa                                                                                            <  1:253964/151‑1 (MQ=255)                                                                       gACACCCCCTGACGTGCCTGTTCGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAg                                                                        <  2:210481/149‑1 (MQ=255)                                                                          accaccTGACGTGCCTGTTCGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAGCtgat                                                                   >  1:47273/1‑151 (MQ=255)                                                                                             cGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAAt                                                   <  2:268749/148‑1 (MQ=255)                                                                                             cGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAAt                                                   >  1:268756/1‑148 (MQ=255)                                                                                              ggCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAATGgtg                                               >  2:257325/1‑151 (MQ=255)                                                                                              ggCAGTATCGACCATCGCGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCCCCTGTAGCTGATGATAAAGATCCTGAATGgt                                                >  2:309798/1‑150 (MQ=255)                                                                                                  gTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTg                                           >  1:11648/1‑151 (MQ=255)                                                                                                     tCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGgcg                                        >  1:372690/1‑151 (MQ=255)                                                                                                         ccATCGGGATCAGTAGAGTTTGCGCGCCAATCTCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCCCCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTg                                    <  2:329931/151‑1 (MQ=255)                                                                                                             cGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTcc                                >  2:295416/1‑151 (MQ=255)                                                                                                              gggATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTccc                               >  1:199962/1‑151 (MQ=255)                                                                                                                gATCAGTAGAGTTTGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCg                             <  1:463460/151‑1 (MQ=255)                                                                                                                             tGCGCGCCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCAGCAAcc                >  2:360295/1‑151 (MQ=255)                                                                                                                               cgcgcCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCAGCAAcc                >  1:72973/1‑149 (MQ=255)                                                                                                                                gcgcCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCAGCAACCAGt             >  1:312253/1‑151 (MQ=255)                                                                                                                                  gcCAATATCCAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCAGCAACCAGTc            <  1:93571/150‑1 (MQ=255)                                                                                                                                           cAGGGCTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCAGCAACCAGTCATAACCAGaa  <  2:3231/151‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                   TCCCCAGCGCGCAGCCCGTGCCACACTGGCCCCGACACCACGCTCACCGTAGGGAGGATAGCGCGTGGCAGACACCACCTGACGTGCCTGTTCGGCAGTATCGACCATCGGGATCAGTAGAGTTTGCGCGCCAATATCCAGGACTTGTTTAATCAGCGGTTTACTGCCTTCCACCGGACGGATCACGGGTTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCAGCAACCAGTCATAACCAGAA  >  NC_000913/2358447‑2358735

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