Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,655,044	T→C	D33G (GAT→GGT)	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,655,044	0	T	C	100.0%	108.9 / NA	32	D33G (GAT→GGT)	rspA	mandelate racemase/muconate lactonizing enzyme family protein RspA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/15);  total (17/15)

CGACGCCAGTAAGCACCTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAA  >  NC_000913/1654909‑1655178                                                                                                                                        |                                                                                                                                       cgACGCCAGTAAGCACCTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACcccc                                                                                                                                    >  2:429743/1‑140 (MQ=255)  gACGCCAGTAAGCACCTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCa                                                                                                                                   <  2:587663/140‑1 (MQ=255)    cGCCAGTAAGCACCTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAg                                                                                                                                 <  1:642470/140‑1 (MQ=255)    cGCCAGTAAGCACCTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAg                                                                                                                                 <  1:429743/140‑1 (MQ=255)    cGCCAGTAAGCACCTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAg                                                                                                                                 <  1:494216/140‑1 (MQ=255)              cACCTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATa                                                                                                                       <  1:613076/140‑1 (MQ=255)                           aaaCTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATAcc                                                                                                                     >  1:707769/1‑129 (MQ=255)                           aaaCTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATAcc                                                                                                                     <  2:707769/129‑1 (MQ=255)                              cTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGa                                                                                                       >  1:494758/1‑140 (MQ=255)                                               tACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGaaa                                                                                      >  2:594334/1‑140 (MQ=255)                                                      cgcATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGcc                                                                               >  1:575198/1‑140 (MQ=255)                                                       gcATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGccc                                                                              <  1:123205/140‑1 (MQ=255)                                                          tCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCgg                                                                           <  1:637055/140‑1 (MQ=255)                                                           cgcgACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGa                                                                          <  2:494758/140‑1 (MQ=255)                                                            gcgACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGAc                                                                         >  1:94693/1‑140 (MQ=255)                                                                           cGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACtt                                                          >  2:631340/1‑140 (MQ=255)                                                                           cGGACAAAGGTGATCCTGCAAAGAAAAGGCCACGGAATGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACtt                                                          <  1:275268/140‑1 (MQ=255)                                                                            ggACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTc                                                         >  1:556335/1‑140 (MQ=255)                                                                             gACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCa                                                        <  2:449060/140‑1 (MQ=255)                                                                                     tGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTAc                                                >  2:481956/1‑140 (MQ=255)                                                                                      gATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACg                                               <  1:393568/140‑1 (MQ=255)                                                                                       aTCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTc                                                         <  2:412812/129‑1 (MQ=255)                                                                                       aTCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTc                                                         >  1:412812/1‑129 (MQ=255)                                                                                                       gCCACGGAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCACGGACAGGAAACAAAAACTTCAGCCTTTACGTTCTTCATGTTCGATTc                              >  2:340258/1‑140 (MQ=255)                                                                                                            ggAAAGCTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGc                         >  2:720151/1‑140 (MQ=255)                                                                                                               aaGCTCACGTCCATTTAGGGTGGCACCCCGAAGGCCCATAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATc                      <  1:353405/140‑1 (MQ=255)                                                                                                                 gCTCACGTCCATTGACGGTGGCGCCGCCATGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGc                    <  1:340258/140‑1 (MQ=255)                                                                                                                  cTCACGTCCATTGAGGGTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTGCAGGCGTTTCGATCTTGATGTGAGATTCCTTGCATCGCt                   >  2:353405/1‑140 (MQ=255)                                                                                                                    cACGTCCATTGAGGGTGGCACCCACAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTg                 >  2:508590/1‑140 (MQ=255)                                                                                                                      cGTCCATTGAGGGTGGTACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTc               >  2:402587/1‑140 (MQ=255)                                                                                                                                 ggTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGCTTCGTTCCATCGCTTGTCGTGATGCATGa    >  1:281721/1‑140 (MQ=255)                                                                                                                                  gTGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGaa   <  2:584720/140‑1 (MQ=255)                                                                                                                                   tGGCACCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGaaa  >  2:640444/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       CGACGCCAGTAAGCACCTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAA  >  NC_000913/1654909‑1655178

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