Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,530,025	T→A	14.9%	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	14.9%	85.4 / 3.1	37	intergenic (+31/‑137)	PP_4018/topB	acetyltransferase/DNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base T (19/12);  new base A (4/2);  total (23/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGT  >  NC_002947/4529879‑4530171                                                                                                                                                   |                                                                                                                                                   gcgcCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTa                                                                                                                                                  >  2:206219/1‑149 (MQ=255)       aGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCATCGCGGGGTTACCagcg                                                                                                                                            >  3:15501/1‑149 (MQ=255)         gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa                                                                                                                                          <  1:206219/149‑1 (MQ=255)         gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa                                                                                                                                          <  5:212209/149‑1 (MQ=255)         gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa                                                                                                                                          <  4:419931/149‑1 (MQ=255)         gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa                                                                                                                                          <  3:411342/149‑1 (MQ=255)          tGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCATCGCGGGGTTAACCGCGAag                                                                                                                                         >  3:229031/1‑149 (MQ=255)            cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                        <  8:248874/148‑1 (MQ=255)            cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                        <  3:367581/148‑1 (MQ=255)                    gccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCt                                                                                                                                >  6:225235/1‑148 (MQ=255)                    gccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCt                                                                                                                                >  5:127301/1‑148 (MQ=255)                                            tCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCCCCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGAc                                                                                                        <  5:432349/148‑1 (MQ=255)                                              gcgTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACaa                                                                                                      >  2:157899/1‑148 (MQ=255)                                                  aGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACaaaaaa                                                                                                  <  6:127301/148‑1 (MQ=255)                                                        cGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAACAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAAACTGTa                                                                                            >  4:268472/1‑148 (MQ=255)                                                         gAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAca                                                                                          <  5:225235/149‑1 (MQ=255)                                                               ttCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCATAGCGGGTTTAACCGCGAAGAGGCCGGCTCAGGCGATATAAATCCCCCCAGAAAAAAAACCTGTAaatccata                                                                                    >  6:19365/1‑148 (MQ=255)                                                                          tGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCAGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCACATAAAACt                                                                         >  5:136234/1‑149 (MQ=255)                                                                           ggCATCCCGCATGACCACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAACACACCTGAACATCCATCCAGATAAAAcat                                                                        >  6:103249/1‑147 (MQ=255)                                                                             cATCCCGCAGGTCCACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCGCTTCGCGGGTTTCCCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGc                                                                      <  4:112065/149‑1 (MQ=255)                                                                                cccGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCctc                                                                   >  2:105015/1‑149 (MQ=255)                                                                                 ccGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGCTATACATCCCCGCTGACAAAAACCCTGTCCATCCATCCAGATAAAACTTGCctct                                                                  >  4:41589/1‑149 (MQ=255)                                                                                  cGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCCCCTGCCCTGGCCTCTTCGCGGGTATCCCAGCGAAGAGGCCGGCTCCGGCGATATACATCCCCGCTGACAAAAACCCTGTACAGCCCGCCAGATAAACCTTGCCACTg                                                                 >  5:76577/1‑149 (MQ=255)                                                                                         cGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGc                                                                      >  3:371364/1‑137 (MQ=255)                                                                                                gTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTgc                                                    <  2:82740/148‑1 (MQ=255)                                                                                                    gggagtCGGGCGCCTGATCCAGCCACTGCAGTGGCTGGTTGGCTGGTTTCGCCGGGATGAGGGCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAATCTTGCCTCTGCGCCCCCGCTTGCGCATc                                               <  6:76577/144‑1 (MQ=255)                                                                                                    gcgACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATc                                               >  6:162322/1‑149 (MQ=255)                                                                                                           gcgcgcCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTgc                                                                >  1:98282/1‑125 (MQ=255)                                                                                                                      ccATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCAcccc                              >  3:206468/1‑148 (MQ=255)                                                                                                                                  ccggCCTCTTCGCGGGTAAACCCGCGAAGAGGCCAGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTccc                 >  7:133631/3‑149 (MQ=255)                                                                                                                                     gCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAGATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGGATCCTATGGCCCTCCCCCCCGCTGAAGCGTCCCCgg              >  2:112780/1‑149 (MQ=255)                                                                                                                                     gCCTCTTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCCCCCCGATGAAGCGTCCCCgg              >  6:312455/1‑149 (MQ=255)                                                                                                                                     gCCTCTTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCGGACAAAAAACCTGTACATACATCCAGATAAAACTTGGCTATGCGCCCCCGCTTGCGCATGCTATGGCCCACCCCCCCCATGAAGCGTcccccg              >  4:6034/1‑147 (MQ=255)                                                                                                                                         ctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCTCCCCCTCTTGCGCATACTATTGCCAAACACCCCGATTAAGCTTCACCGGCCa           >  4:198087/1‑148 (MQ=255)                                                                                                                                            cgcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc        <  5:162322/148‑1 (MQ=255)                                                                                                                                              cgtgTTTAACCGCGTAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGcc     <  6:136234/146‑1 (MQ=255)                                                                                                                                                 gTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGt  >  5:247501/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGT  >  NC_002947/4529879‑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: ‑