Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,135,451	A→G	K539K (AAA→AAG)	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,135,451	0	A	G	100.0%	166.8 / NA	46	K539K (AAA→AAG)	katG	catalase‑peroxidase HPI, heme b‑containing
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/18);  total (28/18)

CGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGG  >  NC_000913/4135303‑4135593                                                                                                                                                     |                                                                                                                                               cgcgTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGc                                                                                                                                              <  2:640037/151‑1 (MQ=255)        ggCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCa                                                                                                                                          <  2:422192/148‑1 (MQ=255)        ggCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAgc                                                                                                                                       >  1:617981/1‑151 (MQ=255)         gCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAgc                                                                                                                                       >  1:93003/1‑150 (MQ=255)          cATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAgcgc                                                                                                                                     >  1:315893/1‑151 (MQ=255)              aaTGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCa                                                                                                                                 <  2:66159/151‑1 (MQ=255)                 gCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAgg                                                                                                                               >  2:542595/1‑150 (MQ=255)                       gcgcGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGc                                                                                                                        <  1:349211/151‑1 (MQ=255)                         gcgACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCAt                                                                                                                      >  2:644069/1‑151 (MQ=255)                                  aTGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTAcc                                                                                                             >  2:44170/1‑151 (MQ=255)                                          gccgcAGCCGGGCGTGCGCTGCCTGTTCTGGAGAAAATCCAGAAAGAGGCTGGTAAAGCCTCGCTGGCGGATATCATAGGGCTGGCTGGTGGGGGTGGGTGTGAGAAGGCCGACAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTgcg                                                                                                      >  1:359880/1‑150 (MQ=255)                                                gCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTc                                                                                               <  2:46546/151‑1 (MQ=255)                                                gCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTc                                                                                               <  1:63715/151‑1 (MQ=255)                                                  cGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTcgc                                                                                             <  2:194363/151‑1 (MQ=255)                                                         gCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATg                                                                                      >  2:586748/1‑151 (MQ=255)                                                              gCCTGTTCTGGAGAAAATCCAGGAAGAGTCTGGTAACGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGcgcg                                                                                  >  1:656188/1‑150 (MQ=255)                                                                  gTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCT‑GTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGa                                                                            <  2:520869/151‑1 (MQ=255)                                                                       ggAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCa                                                                         <  1:4850/150‑1 (MQ=255)                                                                            aaaTCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTg                                                                    >  2:453818/1‑150 (MQ=255)                                                                             aaTCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGAc                                                                  >  2:123330/1‑151 (MQ=255)                                                                              aTCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGAc                                                                  >  1:516913/1‑150 (MQ=255)                                                                                ccAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACAtt                                                               >  1:382561/1‑151 (MQ=255)                                                                                 cAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTg                                                              >  1:651765/1‑151 (MQ=255)                                                                                      agggTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATGGAGATg                                                         <  1:211071/148‑1 (MQ=255)                                                                                      agagTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATg                                                         <  2:452529/151‑1 (MQ=255)                                                                                         gTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGttt                                                      >  2:142167/1‑151 (MQ=255)                                                                                             ggTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTg                                                     <  1:15464/148‑1 (MQ=255)                                                                                             ggTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGa                                                    <  2:540624/149‑1 (MQ=255)                                                                                             ggTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgc                                                  >  2:404915/1‑151 (MQ=255)                                                                                              gTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgct                                                 >  1:168258/1‑151 (MQ=255)                                                                                                        gCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGGGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGcc                                        >  1:298415/1‑150 (MQ=255)                                                                                                        ccTGGCGGATACCATAGTGCTTGCTGGTGTTGTTGGTGTTGAGAAGGCCGCATGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGTCGGGTCGCGCTGATGCGCGTCAGCATCAGACTGACATTGAGATGTTTGGGCTGCTGGGGCCa                                       <  1:224653/150‑1 (MQ=255)                                                                                                          tGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAAt                                     >  2:169208/1‑151 (MQ=255)                                                                                                          tGGCGGATATCATAGTGCTGGCTGGCGTGGGTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAAt                                     >  1:431266/1‑151 (MQ=255)                                                                                                              ggATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGc                                  >  1:1812/1‑150 (MQ=255)                                                                                                                atatCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTg                                >  2:448803/1‑150 (MQ=255)                                                                                                                    cATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGt                           >  1:680390/1‑151 (MQ=255)                                                                                                                        gTGCTGGCTGGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTcc                       <  2:319816/151‑1 (MQ=255)                                                                                                                          gctggctgGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGt                     <  1:499593/151‑1 (MQ=255)                                                                                                                          gctggctgGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGt                     >  1:396351/1‑151 (MQ=255)                                                                                                                          gctggctgGTGTGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGt                     <  1:355505/151‑1 (MQ=255)                                                                                                                                  gtgtGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATcg             >  1:273873/1‑151 (MQ=255)                                                                                                                                  gtgtGGTTGGTGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATcg             <  2:622904/151‑1 (MQ=255)                                                                                                                                      ggttggtGTTGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCt         >  2:211500/1‑151 (MQ=255)                                                                                                                                              ttGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTgg  >  1:185369/1‑150 (MQ=255)                                                                                                                                              ttGAGAAGGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTgg  >  2:139533/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                               CGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGG  >  NC_000913/4135303‑4135593

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