Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,134,394	A→G	D187G (GAC→GGC)	katG →	catalase‑peroxidase HPI, heme b‑containing

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,134,394	0	A	G	100.0%	166.9 / NA	45	D187G (GAC→GGC)	katG	catalase‑peroxidase HPI, heme b‑containing
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/22);  total (23/22)

GCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCG  >  NC_000913/4134247‑4134540                                                                                                                                                    |                                                                                                                                                   gCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGt                                                                                                                                                 >  2:39686/1‑151 (MQ=255)         aaaGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAAc                                                                                                                                         >  2:611933/1‑151 (MQ=255)           aGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCg                                                                                                                                       <  1:283927/151‑1 (MQ=255)              gcgTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGa                                                                                                                                     >  2:217987/1‑150 (MQ=255)                           tgGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACt                                                                                                                       >  2:29658/1‑151 (MQ=255)                            ggCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTg                                                                                                                      <  2:242631/151‑1 (MQ=255)                               cAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTgggg                                                                                                                   >  1:175023/1‑151 (MQ=255)                                        aGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGaaaa                                                                                                          <  1:459314/151‑1 (MQ=255)                                         gAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGaaaa                                                                                                          <  2:203498/150‑1 (MQ=255)                                              atGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCt                                                                                                     <  1:562787/150‑1 (MQ=255)                                              atGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTg                                                                                                    <  1:551302/151‑1 (MQ=255)                                                ggTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGc                                                                                                  >  1:532359/1‑151 (MQ=255)                                                        aaTCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCAc                                                                                          <  1:203779/151‑1 (MQ=255)                                                          tctcCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCg                                                                                        >  2:110465/1‑151 (MQ=255)                                                             ccTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGt                                                                                       <  2:194568/149‑1 (MQ=255)                                                             ccTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCa                                                                                     <  1:601706/151‑1 (MQ=255)                                                                    cGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGCTCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGa                                                                               >  2:353757/1‑150 (MQ=255)                                                                        cTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAgcgc                                                                          >  2:170529/1‑151 (MQ=255)                                                                             tATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGc                                                                      <  2:439088/150‑1 (MQ=255)                                                                                ccTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGaa                                                                   >  1:592716/1‑150 (MQ=255)                                                                                       ggTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGc                                                           >  1:5942/1‑151 (MQ=255)                                                                                          aaCGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTgg                                                        <  2:244673/151‑1 (MQ=255)                                                                                            cGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGt                                                      >  2:96497/1‑151 (MQ=255)                                                                                             gTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTg                                                     >  1:142266/1‑151 (MQ=255)                                                                                             gTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTg                                                     >  1:507464/1‑151 (MQ=255)                                                                                                 cgcTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAAc                                                 <  1:24714/151‑1 (MQ=255)                                                                                                     aGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCgag                                             <  2:244162/151‑1 (MQ=255)                                                                                                                gCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAt                                  >  2:645922/1‑151 (MQ=255)                                                                                                                gCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAt                                  <  2:333093/151‑1 (MQ=255)                                                                                                                    ccGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTa                               <  1:105581/150‑1 (MQ=255)                                                                                                                      gTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGt                            >  2:417133/1‑151 (MQ=255)                                                                                                                        aCCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGtt                           >  2:288582/1‑150 (MQ=255)                                                                                                                            tCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAAccc                      >  2:626880/1‑151 (MQ=255)                                                                                                                            tCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAAcc                       <  1:446392/150‑1 (MQ=255)                                                                                                                              ggTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCgg                    >  1:331393/1‑151 (MQ=255)                                                                                                                              ggTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCgg                    <  2:428009/151‑1 (MQ=255)                                                                                                                               gTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGa                   <  2:156676/151‑1 (MQ=255)                                                                                                                                ttttGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACCGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGa                   >  2:177120/1‑150 (MQ=255)                                                                                                                                 tttGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGa                   >  2:667262/1‑149 (MQ=255)                                                                                                                                    ggTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGcc              >  1:251620/1‑151 (MQ=255)                                                                                                                                    ggTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGc               <  2:493154/150‑1 (MQ=255)                                                                                                                                    ggTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGc               <  1:205141/150‑1 (MQ=255)                                                                                                                                         cGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGAt         >  1:431631/1‑151 (MQ=255)                                                                                                                                              gTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAg    <  2:276433/151‑1 (MQ=255)                                                                                                                                                gAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgcg  <  2:430227/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCG  >  NC_000913/4134247‑4134540

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