Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,654,101	T→C	E347E (GAA→GAG)	rspA ←	mandelate racemase/muconate lactonizing enzyme family protein RspA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,654,101	0	T	C	100.0%	177.7 / NA	49	E347E (GAA→GAG)	rspA	mandelate racemase/muconate lactonizing enzyme family protein RspA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/28);  total (21/28)

AGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGCCAGTACGTACCTGATACAGCGAAGC  >  NC_000913/1653975‑1654232                                                                                                                               |                                                                                                                                    aGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCa                                                                                                                        <  1:1077307/140‑1 (MQ=255)   cAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATg                                                                                                                      <  1:583632/140‑1 (MQ=255)    aGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAAccc                                                                                                                         >  2:98235/1‑136 (MQ=255)    aGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAAccc                                                                                                                         <  1:98235/136‑1 (MQ=255)     ggTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGaa                                                                                                                               <  1:432264/129‑1 (MQ=255)     ggTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGaa                                                                                                                               >  2:432264/1‑129 (MQ=255)     ggTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTa                                                                                                                    <  2:912558/140‑1 (MQ=255)     ggTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTa                                                                                                                    <  2:341491/140‑1 (MQ=255)       ttCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTAtt                                                                                                                  >  1:390535/1‑140 (MQ=255)       ttCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTAtt                                                                                                                  >  1:478249/1‑140 (MQ=255)             gggATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGa                                                                                                            <  1:435385/140‑1 (MQ=255)                atatTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGacac                                                                                                         >  2:201701/1‑140 (MQ=255)                           ccAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTgggg                                                                                              <  2:37415/140‑1 (MQ=255)                           ccAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTgggg                                                                                              <  1:627634/140‑1 (MQ=255)                             aGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGAc                                                                                            <  1:201701/140‑1 (MQ=255)                             aGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGAc                                                                                            <  2:621175/140‑1 (MQ=255)                                ttttCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCca                                                                                         <  2:567121/140‑1 (MQ=255)                                  ttCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCcaca                                                                                       <  2:140032/140‑1 (MQ=255)                                                  ccAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAgcgc                                                                       <  2:853921/140‑1 (MQ=255)                                                     aGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGcgccgc                                                                    >  2:444272/1‑140 (MQ=255)                                                     aGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGcgccgc                                                                    >  1:58615/1‑140 (MQ=255)                                                     aGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATATATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGcgccgc                                                                    >  2:64491/1‑140 (MQ=255)                                                       aCCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAg                                                                  <  1:444272/140‑1 (MQ=255)                                                        cccGGTTTGTCTCCCTGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGc                                                                 >  2:433629/1‑140 (MQ=255)                                                            gTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATg                                                             >  1:409388/1‑140 (MQ=255)                                                              ttGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCa                                                           <  2:58615/140‑1 (MQ=255)                                                              ttGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCa                                                           <  1:64491/140‑1 (MQ=255)                                                                          tGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACaa                                               <  2:478249/140‑1 (MQ=255)                                                                                  gCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAg                                       >  2:631920/1‑140 (MQ=255)                                                                                   ccGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAg                                       >  2:877713/1‑139 (MQ=255)                                                                                   ccGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAg                                       <  1:877713/139‑1 (MQ=255)                                                                                    cGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGa                                     <  1:631920/140‑1 (MQ=255)                                                                                    cGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGa                                     <  1:343928/140‑1 (MQ=255)                                                                                      ttATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGAcc                                   >  2:756172/1‑140 (MQ=255)                                                                                      ttATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGAcc                                   >  2:804030/1‑140 (MQ=255)                                                                                       tATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCg                                  >  2:510694/1‑140 (MQ=255)                                                                                        aTCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCgt                                 <  1:804030/140‑1 (MQ=255)                                                                                                tCCAGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGCCa                         >  2:703943/1‑140 (MQ=255)                                                                                                   aGTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGCCAgta                      >  2:221918/1‑140 (MQ=255)                                                                                                    gTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGCCAgtac                     <  2:409388/140‑1 (MQ=255)                                                                                                    gTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGCCAgtac                     <  2:953429/140‑1 (MQ=255)                                                                                                    gTTGTGCGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGCCAgtac                     <  1:703943/140‑1 (MQ=255)                                                                                                          cGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGCCAgtac                     <  1:69174/134‑1 (MQ=255)                                                                                                          cGGGAAGACTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGCCAgtac                     >  2:69174/1‑134 (MQ=255)                                                                                                                 aCTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGCCAgt                       >  1:229731/1‑125 (MQ=255)                                                                                                                 aCTTCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGCCAgt                       <  2:229731/125‑1 (MQ=255)                                                                                                                   ttCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGCCAgt                       <  1:980659/123‑1 (MQ=255)                                                                                                                   ttCGAGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGCCAgt                       >  2:980659/1‑123 (MQ=255)                                                                                                                       aGCATTTGCTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGATCCAGTACGTACCTGATACAGCGAAGc  >  1:759824/1‑140 (MQ=255)                                                                                                                               |                                                                                                                                    AGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATGTATTCCTGGACACCGAAATTGGGGACCCACAGATCAAAGTGCAGCGCCGCAGCCATGCAGACTGGTGACAAATCGGAAGGACCGTGTGAGCCAGTACGTACCTGATACAGCGAAGC  >  NC_000913/1653975‑1654232

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: ‑