Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,653,981	T→C	E387E (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,653,981	0	T	C	100.0%	204.6 / NA	60	E387E (GAA→GAG)	rspA	mandelate racemase/muconate lactonizing enzyme family protein RspA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/26);  total (34/26)

TTTCACGTTCGACAATCGCCAGTTGATTCGGTTTTTCAATTAATATGCTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATG  >  NC_000913/1653861‑1654116                                                                                                                         |                                                                                                                                        tttCACGTTCGACAATCGCCAGTTGATTCGGTTTTTCAATTAATATGCTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTcgccg                                                                                                                      >  2:567111/1‑140 (MQ=255)   tCACGTTCGACAATCGCCAGTTGATTCGGTTTTTCAATTAATATGCTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCTACTAGTAGATAAGCAGACTCAAAGGGATATTTCgccgcc                                                                                                                    >  2:285599/1‑140 (MQ=255)   tCACGTTCGACAATCGCCAGTTGATTCGGTTTTTCAATTAATATGCTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCgccgcc                                                                                                                    >  1:680382/1‑140 (MQ=255)         tCGACAATCGCCAGTTGATTCGGTTTTTCAATTAATATGCTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCttt                                                                                                              <  1:734363/140‑1 (MQ=255)                  gCCAGTTGATTCGGTTTTTCAATTAATATGCTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAAt                                                                                                     <  1:566693/140‑1 (MQ=255)                     aGTTGATTCGGTTTTTCAATTAATATGCTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCg                                                                                                  >  1:701637/1‑140 (MQ=255)                          aTTCGGTTTTTCAATTAATATGCTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATAcc                                                                                             <  2:154666/140‑1 (MQ=255)                          aTTCGGTTTTTCAATTAATATGCTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATAcc                                                                                             >  2:421333/1‑140 (MQ=255)                           ttCGGTTTTTCAATTAATATGCTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCa                                                                                            >  1:801679/1‑140 (MQ=255)                            tCGGTTTTTCAATTAATATGCTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCaa                                                                                           >  1:26211/1‑140 (MQ=255)                              ggTTTTTCAATTAATATGCTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGa                                                                                         <  1:95664/140‑1 (MQ=255)                                        ttAATATGCTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGt                                                                               <  2:701637/140‑1 (MQ=255)                                        ttAATATGCTTTTCAGTATCTTTCTGCTTACCAGTTCGACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGt                                                                               <  1:369465/140‑1 (MQ=255)                                               gCTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCgg                                                                        <  1:95109/140‑1 (MQ=255)                                               gCTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCgg                                                                        <  2:801679/140‑1 (MQ=255)                                               gCTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCgg                                                                        <  1:567111/140‑1 (MQ=255)                                                cTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGTTATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGa                                                                       >  1:88902/1‑140 (MQ=255)                                                  tttCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCAGCGAATTCGATAGCAAGACCCGGTTTGTCTCCCGGATg                                                                     >  1:291908/1‑140 (MQ=255)                                                   ttCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATTGAATTCGATACCAACACCCGGGTTGTCTCCCGGATGc                                                                    >  1:368854/1‑140 (MQ=255)                                                   ttCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCgg                                                                        <  2:754619/136‑1 (MQ=255)                                                   ttCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCgg                                                                        >  1:754619/1‑136 (MQ=255)                                                         aTCTTACTCCTTACCAGTTCCACAGCGTGTCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAg                                                              <  1:285599/140‑1 (MQ=255)                                                          tCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGc                                                             >  1:761391/1‑140 (MQ=255)                                                              aCTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATAGCAAGACCGGGTGTGACTCCCGGATGCATATAGCCGtt                                                         >  1:260756/1‑140 (MQ=255)                                                               cTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTa                                                        >  2:568273/1‑140 (MQ=255)                                                                   ttACCAGTTCCACAGCGTGCCAACTTCCATACGTGCGGCTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGa                                                    <  1:313282/140‑1 (MQ=255)                                                                    tACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTTGATACCAAGACCGGGTTTGTCTCCCTGATGGATATAGCCGTTATCGaa                                                   >  2:241275/1‑140 (MQ=255)                                                                    tACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCTGTTTGTCTCCCGGATGCATATAGCCGTTATCGaa                                                   >  1:43362/1‑140 (MQ=255)                                                                    tACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGaa                                                   >  2:89050/1‑140 (MQ=255)                                                                    tACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTATCTACCGGATGCATATAGCCGTTATCGaa                                                   >  1:400245/1‑140 (MQ=255)                                                                     aCCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTAt                                                       >  1:578824/1‑135 (MQ=255)                                                                     aCCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTAt                                                       <  2:578824/135‑1 (MQ=255)                                                                     aCCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGaaa                                                  >  1:617816/1‑140 (MQ=255)                                                                     aCCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGaaa                                                  >  1:202215/1‑140 (MQ=255)                                                                     aCCAGTTCAACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGaaa                                                  >  1:435980/1‑140 (MQ=255)                                                                          ttCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACGAAGACCCGATTTGTCTACCGGATGCATATATCCGTTATCGAAAGTCCa                                             >  2:288410/1‑140 (MQ=255)                                                                          ttCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCa                                             >  1:145461/1‑140 (MQ=255)                                                                               caGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTtgt                                        >  1:190316/1‑140 (MQ=255)                                                                                  cGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCg                                     >  1:424398/1‑140 (MQ=255)                                                                                    tGCCATCTTCCAGACGTTCGACTGGTAGATACGCATGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCggg                                   <  2:260756/140‑1 (MQ=255)                                                                                    tGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCTGGAAGCAAATAGCCGTAATCGAAAGTCCAGTTGTGCggg                                   >  2:270324/1‑140 (MQ=255)                                                                                    tGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCggg                                   >  2:580132/1‑140 (MQ=255)                                                                                    tGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCggg                                   >  1:607463/1‑140 (MQ=255)                                                                                      ccATCTTCCAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGaa                                 <  2:158511/140‑1 (MQ=255)                                                                                             ccAGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCg                          <  2:202215/140‑1 (MQ=255)                                                                                               aGACGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAg                        <  1:580132/140‑1 (MQ=255)                                                                                                  cGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCAt                     <  2:145461/140‑1 (MQ=255)                                                                                                  cGTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCAt                     <  2:617816/140‑1 (MQ=255)                                                                                                   gTGCGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCAtt                    >  2:500818/1‑140 (MQ=255)                                                                                                      cGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGt                 >  1:813995/1‑140 (MQ=255)                                                                                                      cGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGt                 <  2:424398/140‑1 (MQ=255)                                                                                                      cGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGt                 >  2:536708/1‑140 (MQ=255)                                                                                                      cGACTGGTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGt                 >  1:197396/1‑140 (MQ=255)                                                                                                            gTAGATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGaa           <  1:421333/140‑1 (MQ=255)                                                                                                               gATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATaa        <  2:607463/140‑1 (MQ=255)                                                                                                               gATAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCAGATGCATATAGCCGTTATCGAAAGTCCTGTTGTGCGGGACGAATTCGAGAATTTGTTCGGAATaa        >  1:319049/1‑140 (MQ=255)                                                                                                                aTAAGCAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAAc       <  2:146231/140‑1 (MQ=255)                                                                                                                aTAAGCAGGCTCATAGGAATATCTCGCCGCCAGCTTTTCCTCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAAc       <  2:291908/140‑1 (MQ=255)                                                                                                                     cAGGCTCATAGGGATTTTTCGCTGCCAGCATTTCATCAAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATg  <  1:288410/140‑1 (MQ=255)                                                                                                                     cAGGCTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGGAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATg  <  2:88902/140‑1 (MQ=255)                                                                                                                         |                                                                                                                                        TTTCACGTTCGACAATCGCCAGTTGATTCGGTTTTTCAATTAATATGCTTTTCATTATCTTACTCCTTACCAGTTCCACAGCGTGCCATCTTCCAGACGTGCGACTGGTAGATAAGCAGGTTCATAGGGATATTTCGCCGCCAGCTTTTCATCGAATTCGATACCAAGACCCGGTTTGTCTCCCGGATGCATATAGCCGTTATCGAAAGTCCAGTTGTGCGGGAAGACTTCGAGCATTTGTTCGGAATAACCCATG  >  NC_000913/1653861‑1654116

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