Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,812,482	T→C	F198L (TTC→CTC)	katE →	catalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,812,482	0	T	C	100.0%	228.1 / NA	86	F198L (TTC→CTC)	katE	catalase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (47/39);  total (47/39)

GGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACTTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTA  >  NC_000913/1812410‑1812553                                                                         |                                                                        ggCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCg                                                                      <  2:490709/76‑1 (MQ=255)   ctctGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGtt                                                                    >  2:1645086/1‑76 (MQ=255)   ctctGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGACGGTATTTTTGACCTCGtt                                                                    <  1:3265253/76‑1 (MQ=255)    tctGCTGATACCGTGCGTGATATCCGTGTCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTg                                                                   >  1:4376919/1‑76 (MQ=255)    tctGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTg                                                                   <  2:3577943/76‑1 (MQ=255)    tctGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTg                                                                   >  2:3793387/1‑76 (MQ=255)      tgctgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGc                                                                 <  1:450210/76‑1 (MQ=255)       gctgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCa                                                                >  1:1653160/1‑76 (MQ=255)       gctgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCa                                                                >  2:363640/1‑76 (MQ=255)        ctgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCaa                                                               >  2:401477/1‑76 (MQ=255)         tgATGCCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAAt                                                              <  2:3109843/76‑1 (MQ=255)         tgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAAt                                                              >  1:2540460/1‑76 (MQ=255)         tgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAAt                                                              <  2:4097721/76‑1 (MQ=255)          gATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATa                                                             >  2:1754978/1‑76 (MQ=255)          gATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATa                                                             >  1:3444920/1‑76 (MQ=255)           atcccGTGCGTGATATCCGTGTCTTTGCGACCACGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATaa                                                            <  2:3824449/73‑1 (MQ=255)           aTACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATaa                                                            <  1:3228415/76‑1 (MQ=255)            tACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAac                                                           <  1:4432432/76‑1 (MQ=255)            tACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAac                                                           >  1:4162601/1‑76 (MQ=255)            tACCGTGCGGGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAac                                                           >  1:957072/1‑76 (MQ=255)              ccGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCACTAacac                                                         >  2:245050/1‑76 (MQ=255)                 tgcgtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGcc                                                      >  2:3306318/1‑76 (MQ=255)                  gcgtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCa                                                     <  1:3970757/76‑1 (MQ=255)                  gcgtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCa                                                     <  1:1395592/76‑1 (MQ=255)                    gtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAt                                                   >  2:3343340/1‑76 (MQ=255)                    gtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAt                                                   >  1:1662457/1‑76 (MQ=255)                    gtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAt                                                   >  1:1412029/1‑76 (MQ=255)                     tgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGTCCTCGTTGGCAATAACCCGCCAAtc                                                  <  1:1942707/76‑1 (MQ=255)                      gATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAtct                                                 <  1:472110/76‑1 (MQ=255)                      gATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAtct                                                 <  1:1367097/76‑1 (MQ=255)                         atCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAtcttct                                              >  2:4701597/1‑76 (MQ=255)                          tCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATcttctt                                             >  2:4101236/1‑76 (MQ=255)                            cGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTa                                           >  1:1265844/1‑76 (MQ=255)                             gTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTAt                                          <  2:2836400/76‑1 (MQ=255)                             gTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTAt                                          >  1:3158689/1‑76 (MQ=255)                              tGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATc                                         >  2:1208305/1‑76 (MQ=255)                               ggCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATcc                                        >  1:1946524/1‑76 (MQ=255)                                 cTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAg                                      >  1:1111398/1‑76 (MQ=255)                                 cTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAg                                      >  2:1486837/1‑76 (MQ=255)                                  tttGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAgg                                     <  1:3626917/76‑1 (MQ=255)                                  tttGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAgg                                     <  1:3638619/76‑1 (MQ=255)                                  tttGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAgg                                     >  2:2735900/1‑76 (MQ=255)                                  tttGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAgg                                     >  2:1281094/1‑76 (MQ=255)                                    tGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGAt                                   >  2:1755968/1‑76 (MQ=255)                                    tGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGAt                                   >  1:4016952/1‑76 (MQ=255)                                    tGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGAt                                   <  1:10868/76‑1 (MQ=255)                                       caccaAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATgcg                                <  1:1979572/76‑1 (MQ=255)                                       caccaAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATgcg                                >  2:3429499/1‑76 (MQ=255)                                        cccaAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATgcgc                               <  2:1069596/75‑1 (MQ=255)                                         ccaAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCa                              >  1:2549474/1‑76 (MQ=255)                                         ccaAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCa                              <  2:3326302/76‑1 (MQ=255)                                         ccaAGTTCTATACCGAAGAGGGTACTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCa                              >  2:787510/1‑76 (MQ=255)                                          caAGTTGTATACCGAAGAGGGGATTTTTGACCTCGTTGGCAGTAACACGCCAATCTTCTTTATCCAGGAGGCGCAt                             <  2:1191569/76‑1 (MQ=255)                                           aaGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATa                            <  1:611488/76‑1 (MQ=255)                                           aaGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATa                            <  1:125805/76‑1 (MQ=255)                                           aaGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATa                            >  2:3137446/1‑76 (MQ=255)                                              ttCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAAt                         >  2:2793984/1‑76 (MQ=255)                                               tCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAAtt                        >  1:2207368/1‑76 (MQ=255)                                                cTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTc                       >  1:3577822/1‑76 (MQ=255)                                                 tataCCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTcc                      >  2:3786979/1‑76 (MQ=255)                                                 tataCCGAAGAGGGTAGTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTcc                      <  2:517867/76‑1 (MQ=255)                                                  ataCCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTccc                     <  1:1478711/76‑1 (MQ=255)                                                   taCCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTcc                      <  1:2320080/74‑1 (MQ=255)                                                   taCCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTcc                      >  2:2320079/1‑74 (MQ=255)                                                    aCCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCg                   <  1:3044294/76‑1 (MQ=255)                                                     ccGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGa                  >  1:2851320/1‑76 (MQ=255)                                                        aaGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGAttt               >  1:1172125/1‑76 (MQ=255)                                                          gagGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTg             >  2:2449823/1‑76 (MQ=255)                                                          gagGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTg             <  1:1930107/76‑1 (MQ=255)                                                          gagGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTg             >  2:2158739/1‑76 (MQ=255)                                                           agGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGt            <  2:3269804/76‑1 (MQ=255)                                                           agGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGt            <  1:1615531/76‑1 (MQ=255)                                                             ggTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTc          <  1:4780335/76‑1 (MQ=255)                                                             ggTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTc          >  1:1793207/1‑76 (MQ=255)                                                              gTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCa         <  1:1101031/76‑1 (MQ=255)                                                              gTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCa         <  1:1631291/76‑1 (MQ=255)                                                               tATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCAt        <  2:3181878/76‑1 (MQ=255)                                                               tATTTTTGACCGCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCAt        <  1:1196462/76‑1 (MQ=255)                                                                aTTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATg       <  1:2793984/76‑1 (MQ=255)                                                                  ttttGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCg     >  1:1134243/1‑76 (MQ=255)                                                                  ttttGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCg     >  2:2182115/1‑76 (MQ=255)                                                                   tttGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCgg    >  1:1322688/1‑76 (MQ=255)                                                                   tttGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCgg    >  2:395591/1‑76 (MQ=255)                                                                   tttGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCgg    <  2:1172124/76‑1 (MQ=255)                                                                    ttGACCTCGTTGGCAATAACACGCCAATCTTCGTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGt   <  1:1280062/76‑1 (MQ=255)                                                                     tGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTa  >  1:1531750/1‑76 (MQ=255)                                                                         |                                                                        GGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACTTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTA  >  NC_000913/1812410‑1812553

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