Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,814,080	C→T	R72C (CGC→TGC)	katE →	catalase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,814,080	0	C	T	97.9%	158.6 / ‑4.2	47	R72C (CGC→TGC)	katE	catalase II
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (24/22);  total (25/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 = 9.85e-01

ACTCACTGGCACCTGAGGACGGCTCTCATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTACGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCA  >  NC_000913/1813946‑1814213                                                                                                                                       |                                                                                                                                      aCTCACTGGCACCTGAGGACGGCTCTCATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCaaa                                                                                                                                  >  1:833131/1‑140 (MQ=255)  cTCACTGGCACCTGAGGACGGCTCTCATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAg                                                                                                                                 <  1:1255724/140‑1 (MQ=255)    cACTGGCACCTGAGGACGGCTCTCATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGc                                                                                                                               <  2:176421/140‑1 (MQ=255)    cACTGGCACCTGAGGACGGCTCTCATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGc                                                                                                                               <  2:121092/140‑1 (MQ=255)          cACCTGAGGACGGCTCTCATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGaa                                                                                                                         >  2:127358/1‑140 (MQ=255)               gAGGACGGCTCTCATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTa                                                                                                                    >  2:272702/1‑140 (MQ=255)                 ggACGGCTCTCATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCAGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCa                                                                                                                              >  2:590398/1‑128 (MQ=255)                 ggACGGCTCTCATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCAGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCa                                                                                                                              <  1:590398/128‑1 (MQ=255)                   aCGGCTCTCATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATgcg                                                                                                                >  2:1078973/1‑140 (MQ=255)                   aCGGCTCTCATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATgcg                                                                                                                >  2:1328761/1‑140 (MQ=255)                     ggCTCTCATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCt                                                                                                              <  1:1078973/140‑1 (MQ=255)                     ggCTCTCATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCt                                                                                                              <  2:1050849/140‑1 (MQ=255)                     ggCTCTCATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCt                                                                                                              <  1:272702/140‑1 (MQ=255)                       ctctCATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGa                                                                                                            >  2:1030116/1‑140 (MQ=255)                       ctctCATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGa                                                                                                            >  1:145463/1‑140 (MQ=255)                           cATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCAc                                                                                                        <  1:466243/140‑1 (MQ=255)                           cATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCAc                                                                                                        <  2:741853/140‑1 (MQ=255)                               gTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAAt                                                                                                    <  2:516993/140‑1 (MQ=255)                                tCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATc                                                                                                   <  2:558045/140‑1 (MQ=255)                                   aGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAgg                                                                                                <  1:127358/140‑1 (MQ=255)                                     cGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGc                                                                                              >  2:920929/1‑140 (MQ=255)                                                 cacCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGcc                                                                                  >  2:1250994/1‑140 (MQ=255)                                                      ccAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCcgacga                                                                             <  1:920929/140‑1 (MQ=255)                                                          gTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCaa                                                                         >  1:715581/1‑140 (MQ=255)                                                              acaACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCGAAACt                                                                     <  1:652851/140‑1 (MQ=255)                                                                     aCCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCg                                                              >  2:871138/1‑140 (MQ=255)                                                                     aCCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTAGGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCg                                                              >  2:1184390/1‑140 (MQ=255)                                                                     aCCGCCCCAGAGATCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAGTTATGCGCTGACCACTAATCAGGGCGTGCACATCGCCGACGATCAAAACTCACtaag                                                              >  2:413955/1‑137 (MQ=255)                                                                            cAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGt                                                       >  1:527511/1‑140 (MQ=255)                                                                            cAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGt                                                       <  1:1328761/140‑1 (MQ=255)                                                                            cAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGt                                                       >  2:679059/1‑140 (MQ=255)                                                                                   ccTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTg                                                >  2:551465/1‑140 (MQ=255)                                                                                         aGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCaa                                          >  2:213487/1‑140 (MQ=255)                                                                                            cccTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACgc                                       <  2:715581/140‑1 (MQ=255)                                                                                                   aCGCGTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGa                                <  1:213487/140‑1 (MQ=255)                                                                                                      cgTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGa                             >  2:1267914/1‑140 (MQ=255)                                                                                                      cgTAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGa                             >  1:1133083/1‑140 (MQ=255)                                                                                                      cgTAACGAAAAACTTAATTCTCTGGAAGACGTACGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGa                             >  1:1121103/1‑140 (MQ=255)                                                                                                        tAACGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGAtt                           >  2:1210264/1‑140 (MQ=255)                                                                                                           cGAAAAACTTAATTCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTa                        >  1:666928/1‑140 (MQ=255)                                                                                                                       ttCTCTGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCgaga            >  2:1082113/1‑140 (MQ=255)                                                                                                                           ctGGAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAAt        <  1:1082113/140‑1 (MQ=255)                                                                                                                             ggAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCa      <  1:261982/140‑1 (MQ=255)                                                                                                                             ggAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCa      <  1:551465/140‑1 (MQ=255)                                                                                                                             ggAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGA‑TTTATTCTgcgc                <  1:207928/129‑1 (MQ=255)                                                                                                                             ggAAGACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGA‑TTTATTCTgcgc                >  2:207928/1‑129 (MQ=255)                                                                                                                                 gatgTATGCAAAGGCAGTGAAAATTATGCGATGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCa  <  2:833131/137‑1 (MQ=255)                                                                                                                                 gACGTATGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCa  >  1:504802/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                      ACTCACTGGCACCTGAGGACGGCTCTCATCGTCCAGCGGCTGAACCAACACCGCCAGGTGCACAACCTACCGCCCCAGGGAGCCTGAAAGCCCCTGATACGCGTAACGAAAAACTTAATTCTCTGGAAGACGTACGCAAAGGCAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCA  >  NC_000913/1813946‑1814213

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