Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,816,002	C→T	D712D (GAC→GAT)	katE →	catalase HPII

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,816,002	0	C	T	100.0%	99.1 / NA	30	D712D (GAC→GAT)	katE	catalase HPII
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/17);  total (13/17)

CTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTACCTTTGCCGGTGCACCTTCGCTGACGGTCGATGCGGTCATTGTCCCTTGCGGCAATATCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGACCAGGGTGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTG‑ACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCA  >  NC_000913/1815765‑1816235                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           ctACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTACCTTTGCCGGTGCACCTTCGCTGACGGTCGATGCGGTCATTGTCCCTTGCGGCAATATCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAAACAAGATCGCTGATCAGGGGgaaga                                                                                                                                                                                                                                 >  2:319321/1‑249 (MQ=255)      cccGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTACCTTTGCCGGTGCACCTTCGCTGACGGTCGATGCGGTCATTGTCCCTTGCGGCAATATCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGAAGGGAt                                                                                                                                                                                                                           <  2:197441/250‑1 (MQ=255)         gaaggggTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTACCTTTGCCGGTGCACCTTCGCTGACGGTCGATGCGGTCATTGTCCCTTGCGGCAATATCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGAAGGGATtgtg                                                                                                                                                                                                                       <  1:110197/247‑1 (MQ=255)                 gAAGTGACTGCGGATGACTGTACGGTGTTGCCTATAGCCGCTACCTTTGCCGGTGCACCCTCGCTGACGGTCGATGCGGTCATTGTCCCTTGCGGCAATATCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGGAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGAAGGGATTGTGGAAGCTGa                                                                                                                                                                                                               <  2:30717/251‑1 (MQ=255)                      gACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTACCTTTGCCGGTGCACCTTCGCTGACGGTCGATGCGGTCATTGTCCCTTGCGGCAATATCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGAAACAATCAAGATCGCTGATCAGGGGGAAGAAGGGATTGTGGAAGCTGACAgcg                                                                                                                                                                                                          >  2:116467/1‑251 (MQ=255)                             gATGACGGTACGGTGTTGCCTATAGCCGCTACCTTTGCCGGTGCACCTTCGCTGACGGTCGATGCGGTCATTGTCCCTTGCGGCAATATCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGGTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGCAGGGATTGTGGAAGCTGACAgcggcgacgg                                                                                                                                                                                                   >  2:357830/1‑244 (MQ=255)                                              gCCTATAGCCGCTACCTTTGCCGGTGCGCCTTCGCTGACGGTCGATGCGGTCATTGTCCCTTGCGGCAATATCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAAc                                                                                                                                                                                  <  1:266894/251‑1 (MQ=255)                                                           aCCTTTGCCGGTGCACCTTCGCTGACGGTCGATGCGGTCATTGTCCCTTGCGGCAATATCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGAt                                                                                                                                                                     <  1:116467/251‑1 (MQ=255)                                                             cTTTGCCGGTGCACCTTCGCTGACGGTCGATGCGGTCATTGTCCCTTGCGGCAATATCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAAAAATCAAGATAGCTGATCAGGGGGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATgg                                                                                                                                                                   >  1:393427/1‑251 (MQ=255)                                                                                  gACGGTCGATGCGGTCATTGTCCGTTGCGGCAATATCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCGTACAAACACCTTAAGCCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGAAGGGATTGTGGAAGCTGAAAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCAc                                                                                                                                              <  2:3558/251‑1 (MQ=255)                                                                                  gACGGTCGATGCGGTCATTGTCCCTTGCGGCAATATCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTc                                                                                                                                                <  2:78642/249‑1 (MQ=255)                                                                                            gCGGTCATTGTCCCTTGCGGCAATATCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTa                                                                                                                                     <  1:319321/250‑1 (MQ=255)                                                                                               gTCATTGTCCCTTGCGGCAATATCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACTCCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGAt                                                                                                                                 <  2:150413/251‑1 (MQ=255)                                                                                                  aTTGTCCCTTGCGGCAATATCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACgg                                                                                                                                                                                                   <  2:265977/182‑1 (MQ=255)                                                                                                  aTTGTCCCTTGCGGCAATATCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGGGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACgg                                                                                                                                                                                                   >  1:265977/1‑182 (MQ=255)                                                                                                            gCGGCAATATCGCGGATATCGCTGACAACCGCGATGCCAGCTACTACCTGATTGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGTTgaag                                                                                                                                                                                                                                  <  2:401451/141‑1 (MQ=255)                                                                                                                    atCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACCCCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTG‑ACAAAATTCCTGCCTGAt                                                                                                            <  1:349853/250‑1 (MQ=255)                                                                                                                             atCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTG‑ACAAAATTCCTGCCTGATGGGAggcggc                                                                                                  <  1:298723/251‑1 (MQ=255)                                                                                                                              tCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCCCGCCAGTTTCAAGCAACAATCAAAATCGCTGATCAGGGTGAAGAAGGGAATGTGGAAGCTGACAGCGCTGACGGGAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACAACGCGTGTGGTTCCGGCTTCCTAAGATTGGACAAA‑TTCCTGCCTGATGGGGGgcggcg                                                                                                 >  2:261225/1‑251 (MQ=255)                                                                                                                                 cTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTCCAAACACCTTAAACCGATTGCGCTGGCGGGCGACGCGCGCAAGTTTTAAGCAACAATCAAGATTGCTGATCAGGGGGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGGAGGTTTATGGATGAACTGCTAAAGCTGATGGCAGCACACCGCGTGTGGTCCCGCATTCCTAAGCTTG‑GCAAAATTCCTGCCTGATTGGAggcgggcca                                                                                               >  2:160678/1‑246 (MQ=255)                                                                                                                                            gATGCCAACTACTACCTGATAGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGGGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTG‑ACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCt                                                                                   <  2:303483/251‑1 (MQ=255)                                                                                                                                                 cAACTACTACCTGATGGAAGCCTACAAACTCCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTG‑ACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAAt                                                                              <  1:197736/251‑1 (MQ=255)                                                                                                                                                   actactACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGGGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTG‑ACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATg                                                                             >  2:102874/1‑250 (MQ=255)                                                                                                                                                    ctactaCCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTG‑ACAAAATTCCTGCCTGATGGGAGGCGGCGCCATTGCGCCGCCTCAAATGAt                                                                           >  2:33923/1‑251 (MQ=255)                                                                                                                                                      actacCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCCGGGTGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACCCCGCGTGTGGTCACGCATTCCTAAGATTT‑ACAAAATTTCTGCCTGATTGGAGGCGGCGCCATTGTGCCGCCTCCAATTAttt                                                                         >  2:309961/1‑251 (MQ=255)                                                                                                                                                         acCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTG‑ACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTa                                                                        <  1:160678/249‑1 (MQ=255)                                                                                                                                                              aTGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGGGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTG‑ACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTTCGCCGCCTCAAATGATTTAAATAGTg                                                                 >  2:54654/1‑251 (MQ=255)                                                                                                                                                                                          ttGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGGGAAGAAGGGATTGTGGAAGCTTACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACAACGCGTGTGGTCACGCATTCCTAAGATTG‑ACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACCTAGTTGGCTTTGTTTATGCCCGGTGCGGCGTgg                                     >  2:208380/1‑250 (MQ=255)                                                                                                                                                                                                               gcAAGTTTAAAGCAACAATCAAGATCGCTGATCAGGGTGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTG‑ACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaa                 >  1:23746/1‑250 (MQ=255)                                                                                                                                                                                                                             acaaTCAAGATCGCTGATCAGGGTGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTG‑ACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCa  <  1:102874/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           CTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTACCTTTGCCGGTGCACCTTCGCTGACGGTCGATGCGGTCATTGTCCCTTGCGGCAATATCGCGGATATCGCTGACAACGGCGATGCCAACTACTACCTGATGGAAGCCTACAAACACCTTAAACCGATTGCGCTGGCGGGTGACGCGCGCAAGTTTAAAGCAACAATCAAGATCGCTGACCAGGGTGAAGAAGGGATTGTGGAAGCTGACAGCGCTGACGGTAGTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTG‑ACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCA  >  NC_000913/1815765‑1816235

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