Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,169,870	G→A	8.3%	P8L (CCG→CTG)	katG ←	catalase‑peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,169,870	0	G	A	8.3%	134.6 / 2.6	48	P8L (CCG→CTG)	katG	catalase‑peroxidase
Reads supporting (aligned to +/- strand):  ref base G (24/20);  new base A (2/2);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

AAAGTCGAGGCTCTTGAACGCCTTGGCATAGTCGAAGTCCGGGTCGGGGCTGGACTTGGAGGAATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCTCATCAGGTGTACACTGGGCTTGCGAAGCT  >  NC_002947/4169725‑4170006                                                                                                                                                  |                                                                                                                                         aaaGTCGAGGCTCTTGAACGCCTTGGCATAGTCGAAGTCCGGGTCGGGGCTGGACTTGGAGGAATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGc                                                                                                                                       <  4:157854/149‑1 (MQ=255)         ggCTCTTGAACGCCTTGGCATAGTCGAAGTCCGGGTCGGGGCTGGACTTGGAGGAATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGAt                                                                                                                               >  5:174510/1‑149 (MQ=255)                 aaCGCCTTGGCATAGTCGAAGTCCGGGTCGGGGCTGGACTTGGAGGAATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGa                                                                                                                       >  2:248313/1‑149 (MQ=255)                  aCGCCTTGGCATAGTCGAAGTCCGGGTCGGGGCTGGACTTGGAGGAATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACAGGCATTTCGATTCGTTCGa                                                                                                                       >  6:185751/1‑148 (MQ=255)                  aCGCCTTGGCATAGTCGAAGTCCGGGTCGGGGCTGGACTTGGAGGAATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACAGGCATTTCGATTCGTTCGa                                                                                                                       >  2:353743/1‑148 (MQ=255)                      cTTGGCATAGTAGAAGTCCGGGTCGGGGCTGGACTTGGAGGAATGCTGGTGAAGAATTCTCATGTTGAGATGGTCAGGCCACCAGTAACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTAGACATcg                                                                                                                  >  6:348782/1‑148 (MQ=255)                         ggCATAGTCGAAGTCCGGGTCGGGGCTGGACTTGGAGGAATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGc                                                                                                                <  5:186587/148‑1 (MQ=255)                            aTAGTCGAAGTCCGGGTCGGGGCTGGACTTGGAGGAATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCggg                                                                                                             >  3:322665/1‑148 (MQ=255)                             tAGTCGAAGTCCGGGTCGGGGCTGGACTTGGAGGAATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTa                                                                                                           <  1:182453/149‑1 (MQ=255)                                  gAAGTCCGGGTCGGGGCTGGACTTGGAGGAATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCttt                                                                                                      >  5:208207/1‑149 (MQ=255)                                   aaGTCCGGGTCGGGGCTGGACTTGGAGGAATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCAttt                                                                                                                                   <  1:7136/119‑1 (MQ=255)                                   aaGTCCGGGTCGGGGCTGGACTTGGAGGAATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCAttt                                                                                                                                   >  2:7136/1‑119 (MQ=255)                                        cGGGTCGGGGCTGGACTTGGAGGAATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACAAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCg                                                                                                >  1:17700/1‑149 (MQ=255)                                                  cTGGACTTGGAGGAATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACAGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGa                                                                                      <  1:353743/149‑1 (MQ=255)                                                  cTGGACTTGGAGGAATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACAGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGa                                                                                      <  5:185751/149‑1 (MQ=255)                                                          ggaggaATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCg                                                                                <  1:248313/147‑1 (MQ=255)                                                          ggaggaATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGt                                                                              >  3:260081/1‑149 (MQ=255)                                                           gaggaATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTc                                                                             >  8:231140/1‑149 (MQ=255)                                                                aTGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACAAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGc                                                                        <  2:17700/149‑1 (MQ=255)                                                                 tGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGc                                                                        >  7:343081/1‑148 (MQ=255)                                                                             aTTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCa                                                           <  6:193253/149‑1 (MQ=255)                                                                                 tcAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGAttc                                                                                                                             >  3:6270/1‑79 (MQ=255)                                                                                 tcAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGAttc                                                                                                                             <  4:6270/79‑1 (MQ=255)                                                                                 tcAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCgg                                                                                               <  1:128957/109‑1 (MQ=255)                                                                                 tcAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCgg                                                                                               >  2:128957/1‑109 (MQ=255)                                                                                 tcAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGcc                                                        >  4:346781/1‑148 (MQ=255)                                                                                   aGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAg                                                     >  4:229797/1‑149 (MQ=255)                                                                                                cAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTc                                          <  3:271572/147‑1 (MQ=255)                                                                                                cAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTc                                          >  4:271572/1‑147 (MQ=255)                                                                                                 aGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGt                                        >  7:639/1‑148 (MQ=255)                                                                                                        caGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTAccc                                 >  6:229861/1‑148 (MQ=255)                                                                                                          gTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCt                               >  5:168001/1‑148 (MQ=255)                                                                                                           tCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCtca                             >  6:207062/1‑149 (MQ=255)                                                                                                             aCGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCtcat                            <  5:229861/148‑1 (MQ=255)                                                                                                               ggTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCTCATCAg                         <  8:343081/149‑1 (MQ=255)                                                                                                               ggTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCTCATCAg                         >  5:326084/1‑149 (MQ=255)                                                                                                               ggTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCTCATCAg                         >  5:102607/1‑149 (MQ=255)                                                                                                                      ggtgCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCTCATCAGGTGTaca                  >  3:108585/1‑149 (MQ=255)                                                                                                                       gtgCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCTCATCAGGTGTaca                  <  4:260081/148‑1 (MQ=255)                                                                                                                          ccgccACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCTCATCAGGTGTACACTgg              <  3:346781/149‑1 (MQ=255)                                                                                                                              cACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCCACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCTCATCAGGTGTACACTGGGCtt          >  4:91396/1‑149 (MQ=255)                                                                                                                                ccTGCGGTTTGATGGAACGGGCATTTCGATTCTTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCTCATCAGGTGTACACTGGGCTTGc        <  6:208207/149‑1 (MQ=255)                                                                                                                                 cTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCTCATCAGGTGTACACTGGGCTTGCg       <  6:168001/149‑1 (MQ=255)                                                                                                                                     ggTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCTCATCAGGTGTACACTGGGCTTGCGAAGc   <  6:102607/149‑1 (MQ=255)                                                                                                                                      gtgtGATGGAACGGGCAATTCGATTCGTTCGACATCTGCGGGTACATTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCTCATCAGGTGTACACTGGGCTTGCGAAGCt  <  5:348782/146‑1 (MQ=255)                                                                                                                                      gTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCTGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCTCATCAGGTGTACACTGGGCTTGCGAAGCt  <  7:231140/149‑1 (MQ=255)                                                                                                                                      gTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCTCATCAGGTGTACACTGGGCTTGCGAAGCt  <  6:326084/149‑1 (MQ=255)                                                                                                                                                 cGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCTCATCAGGTGTACACTGGGCtg          >  2:50704/1‑129 (MQ=255)                                                                                                                                                  |                                                                                                                                         AAAGTCGAGGCTCTTGAACGCCTTGGCATAGTCGAAGTCCGGGTCGGGGCTGGACTTGGAGGAATGCTGGTGAAGAATTCTCAGGTTGAGCTGGTCAGGCCACCAGTCACGGTTGGTGGTGCCGCCACCTGCGGTTTGATGGAACGGGCATTTCGATTCGTTCGACATCTGCGGGTACCTTTGGGTCGGTTATCCAGATTTCCGGTCTATGCCAGTGTTCTTTCAGCCGAGCGGGTGGGCTCGTACTACCCCTCATCAGGTGTACACTGGGCTTGCGAAGCT  >  NC_002947/4169725‑4170006

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: ‑