Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,530,025	T→A	9.6%	intergenic (+31/‑137)	PP_4018 → / → topB	acetyltransferase/DNA topoisomerase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,530,025	0	T	A	9.6%	107.7 / 3.8	42	intergenic (+31/‑137)	PP_4018/topB	acetyltransferase/DNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base T (19/19);  new base A (2/2);  total (21/21)
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.66e-01

CCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGT  >  NC_002947/4529883‑4530171                                                                                                                                                |                                                                                                                                                   ccAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcg                                                                                                                                              <  7:174167/149‑1 (MQ=255)  cAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcgc                                                                                                                                             <  2:221612/149‑1 (MQ=255)  cAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGCTTACCcgc                                                                                                                                             <  4:103764/149‑1 (MQ=255)         acgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAGTTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAgg                                                                                                                                      <  6:72724/149‑1 (MQ=255)          cgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGc                                                                                                                                     <  7:201630/149‑1 (MQ=255)          cgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGc                                                                                                                                     <  6:18561/149‑1 (MQ=255)            ccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCg                                                                                                                                   >  8:202471/1‑149 (MQ=255)                           aGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcgc                                                                                                                                             <  8:216641/124‑1 (MQ=255)                           aGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcgc                                                                                                                                             >  7:216641/1‑124 (MQ=255)                           aGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCTAAGAGGCCGGCTCAGGCGATATAc                                                                                                                    <  7:93601/149‑1 (MQ=255)                           aGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGAtata                                                                                                                     >  1:122522/1‑148 (MQ=255)                                 tGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGCGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATcccc                                                                                                              >  6:237/1‑149 (MQ=255)                                              aGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAAc                                                                                                 >  7:177641/1‑149 (MQ=255)                                                    cGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTa                                                                                            >  6:309056/1‑148 (MQ=255)                                                       ggAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGGTCTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcat                                                                                         >  1:224825/1‑149 (MQ=255)                                                              cTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc                                                                                    >  1:54464/1‑146 (MQ=255)                                                              cTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACAtccatcca                                                                                  <  5:309056/148‑1 (MQ=255)                                                              cTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATAAATCCCCGCAGACAAAAAACCTGt                                                                                             >  5:33136/1‑137 (MQ=255)                                                               tCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAagag                                                                                                                                       <  8:182652/94‑1 (MQ=255)                                                               tCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAagag                                                                                                                                       >  7:182652/1‑94 (MQ=255)                                                               tCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTAGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACAAACCCGCTGACAAAAAACCTGTAaatccaac                                                                                    >  7:256558/1‑143 (MQ=255)                                                                    gCTGCCACCCCGCAGGGCTACTTGGTGCGC‑GCCGCGCGCGCCTTATCCATCCCGTCCAGTGGCCTCTTCGGGGTTTTACCCGCTATGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATaaaa                                                                           <  8:177641/149‑1 (MQ=255)                                                                             ccGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATACATCCAGATAAAACTTGCctct                                                                  >  4:34896/1‑149 (MQ=255)                                                                                 ttGTCGACTGGGTGCGC‑GACTGGCGCTCCTGATCCATCACCTGCACTGGCCTCTTCGCGGTTTTACCCGCTAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTgcgc                                                              <  5:237/148‑1 (MQ=255)                                                                                      gACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAAc                                                                                                 <  5:278570/109‑1 (MQ=255)                                                                                      gACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAAc                                                                                                 >  6:278570/1‑109 (MQ=255)                                                                                         aTGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCtt                                                      <  7:202471/149‑1 (MQ=255)                                                                                            gTGCGC‑GACGCGCGCGCCTTATCCATCACCTCCAGTGGCCTCTTCGCGGTTTAACCCGCTAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTgc                                                    <  8:49295/148‑1 (MQ=255)                                                                                             tGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTgcg                                                   >  5:150421/1‑148 (MQ=255)                                                                                                 c‑gACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATcc                                              <  6:150421/149‑1 (MQ=255)                                                                                                                gATCCATCACCTGCACTGGCCTCTTAGCGGGTTTTCCCGCGAAGAGGCCGGCTCAGGCGCTATACATCCCCGCTGACAAAAAACCTGTTCATCCCTCCAGATTAAACTTGCCTCTGCGCCCCCGCTTGCGCATGCTATGGGCCaccacc                                >  3:155309/1‑149 (MQ=255)                                                                                                                gATCCATCACCTGCACTGGCCTCATAGCGGGTTTACCCGCGAAGAGGGCGGCTCAGGCGATATACATACCCCCGGGCAAAAAACCTGTACATCCAGGCAGATAAAAATTgg                                                                      >  8:199251/1‑110 (MQ=37)                                                                                                                     aTCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCAt                                                >  7:281153/1‑128 (MQ=255)                                                                                                                     aTCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCGTGCGCCGCCTATGGCCCCCCACCCCGAt                           >  2:242027/1‑149 (MQ=255)                                                                                                                     aTCACCTGCACTGCCCTCTGCGCGGGTTTCCCCGCGAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCAt                                                <  8:281153/128‑1 (MQ=255)                                                                                                                        aCCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGa                         <  2:54464/148‑1 (MQ=255)                                                                                                                          cTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCCTCCAGATAAAACTTGCCTCTGCGCCCCCCCTTGCGCAACCTATGGCCCACCACCCCGATGAAGc                      >  3:77265/1‑149 (MQ=255)                                                                                                                          cTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGc                      >  7:75034/1‑149 (MQ=255)                                                                                                                                 ggCCTCTTCGCGGGTTTCCCGGCGACGCGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCg               <  1:242027/149‑1 (MQ=255)                                                                                                                                       ttCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGc                      >  3:178056/1‑136 (MQ=255)                                                                                                                                       ttCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGc                      <  4:178056/136‑1 (MQ=255)                                                                                                                                              gTTTACCCGCTAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATACATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGt  <  3:34896/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   CCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGC‑GACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGT  >  NC_002947/4529883‑4530171

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