Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,169,556	C→T	17.4%	A113T (GCC→ACC)	katG ←	catalase‑peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,169,556	0	C	T	17.4%	46.1 / 3.7	23	A113T (GCC→ACC)	katG	catalase‑peroxidase
Reads supporting (aligned to +/- strand):  ref base C (11/8);  new base T (2/2);  total (13/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

AGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTGTTGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGCCACCAATCCTGCGAGTCGGTCATCA  >  NC_002947/4169412‑4169696                                                                                                                                                 |                                                                                                                                             aGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTGTTGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAg                                                                                                                                          <  7:64676/149‑1 (MQ=255)                                                                          cTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGc                                                                 >  1:146897/1‑149 (MQ=255)                                                                             cTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGTGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGa                                                              >  4:155873/1‑149 (MQ=255)                                                                             cTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGTGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGa                                                              >  4:150894/1‑149 (MQ=255)                                                                             cTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGa                                                              >  4:224589/1‑149 (MQ=255)                                                                                  tCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGaa                                                          >  8:245876/1‑148 (MQ=255)                                                                                     gCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGc                                                      >  6:221171/1‑149 (MQ=255)                                                                                               gcgcCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGt                                            >  1:96630/1‑149 (MQ=255)                                                                                                   cGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGcc                                        >  2:54963/1‑149 (MQ=255)                                                                                                       ccTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGCTCCCGGCGCTATGCCAGGCCATGCGGATGAACAGCCGGCCGTAGTGGCCGAAg                                    >  8:55003/1‑149 (MQ=255)                                                                                                           gTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTcg                                 <  5:221171/148‑1 (MQ=255)                                                                                                              cAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGTGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTcggccg                             <  3:155873/149‑1 (MQ=255)                                                                                                              cAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGTGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTcggccg                             <  3:150894/149‑1 (MQ=255)                                                                                                               aGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTcggccgg                            >  7:200959/1‑149 (MQ=255)                                                                                                                 gCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGACTTcggccacta                          <  7:245876/149‑5 (MQ=255)                                                                                                                     aCGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGccac                        <  3:224589/147‑1 (MQ=255)                                                                                                                           tCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCTGGCCACCAATCCt                 <  1:54963/148‑1 (MQ=255)                                                                                                                              ggCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGG‑CGAAGTCGGCCGGCCACCAATCCTGCGAg            >  4:285156/1‑149 (MQ=255)                                                                                                                                ccAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGCCACCAATCCTGCGAGt           <  2:96630/149‑1 (MQ=255)                                                                                                                                    gtGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGCCACCAATCCTGCGAGTCg         <  8:200959/146‑1 (MQ=255)                                                                                                                                     tGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGCCACCAATCCTGCGAGTCGGt       >  7:126648/1‑148 (MQ=255)                                                                                                                                     tGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGCCACCAATCCTGCGAGTCGGt       >  7:296194/1‑148 (MQ=255)                                                                                                                                          aGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGCCACCAATCCTGCGAGTCGGtcatca  <  2:248023/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                             AGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTGTTGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGCCACCAATCCTGCGAGTCGGTCATCA  >  NC_002947/4169412‑4169696

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