Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	474,195	Δ1 bp	11.8%	coding (793/1983 nt)	dnaG ←	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	474,195	0	C	.	11.8%	97.2 / 8.9	34	coding (793/1983 nt)	dnaG	DNA primase
Reads supporting (aligned to +/- strand):  ref base C (16/14);  new base . (2/2);  total (18/16)
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.82e-01

AAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGG  >  NC_002947/474058‑474332                                                                                                                                          |                                                                                                                                           aaaaCAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑tcgtc                                                                                                                               >  1:7420/1‑150 (MQ=255)             cACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGA‑GATCATCGTCGTCGAg                                                                                                                            >  2:350044/1‑141 (MQ=255)             cACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGA‑GATCATCGTCGTCGAg                                                                                                                            <  1:350044/141‑1 (MQ=255)              aCGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTtg                                                                                                                  <  2:147884/150‑1 (MQ=255)                         cacGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTTCCCGGCCAGGGCAATGACGTACATGTAGTGCCTGGCGACGATTATC‑TCGTCGAGGTTGCGGTTGTGCTTgcgcgc                                                                                                       >  1:233533/1‑150 (MQ=255)                                                    ttCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTgc                                                                                                           <  1:102759/119‑1 (MQ=255)                                                    ttCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTgc                                                                                                           >  2:102759/1‑119 (MQ=255)                                                      cGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCt                                                                          >  1:183380/1‑150 (MQ=255)                                                        cTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTg                                                                                                                        >  2:483587/1‑102 (MQ=255)                                                        cTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTg                                                                                                                        <  1:483587/102‑1 (MQ=255)                                                                   gCCAAGGGTGTCCACGGCATTGCGCAGGCCTTGCTGGGTCAGGGCAATTACGTCCTTTTAGCCCTCGACGACCATGATC‑TCGTCGAGGTTGCAGTTGTGCTTCCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCgg                                                             <  2:96724/150‑1 (MQ=255)                                                                   gCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCgg                                                             <  1:294558/150‑1 (MQ=255)                                                                    ccAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGCCGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCggg                                                            <  2:38585/150‑1 (MQ=255)                                                                      aaGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGt                                                          >  2:474851/1‑150 (MQ=255)                                                                      aaGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGt                                                          >  1:19453/1‑150 (MQ=255)                                                                               cACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGa                                                 <  2:377482/150‑1 (MQ=255)                                                                                      ttGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAgg                                          >  2:369697/1‑150 (MQ=255)                                                                                            aGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGa                                                 <  1:205135/137‑1 (MQ=255)                                                                                            aGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGa                                                 >  2:205135/1‑137 (MQ=255)                                                                                              gCCTTGCTGGGCCAGGGCAATGACGTCCATATAGACCTCGACGACGAGGATC‑GCGTCGGGGTTGCGGTTAGTTTTGCGGGCATAAAGCAGCTCGGGCAGCTCCTCCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTGCAAGTACTTggg                                  >  1:403979/1‑150 (MQ=255)                                                                                                 ttGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCTCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCtt                               >  1:29806/1‑150 (MQ=255)                                                                                                 ttGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTACTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCtt                               >  1:29807/1‑150 (MQ=255)                                                                                                         ccAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGc                       >  2:464230/1‑150 (MQ=255)                                                                                                          cAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCgg                                                    >  2:227706/1‑120 (MQ=255)                                                                                                          cAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCgg                                                    <  1:227706/120‑1 (MQ=255)                                                                                                          cAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGcc                      <  2:403979/150‑1 (MQ=255)                                                                                                           aGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCg                     >  1:243333/1‑150 (MQ=255)                                                                                                            gggCAATGACGTCCATGTAGCCCTCGACGACGATTATC‑TCGTCGAGGTTGCGGTTGTCGATGCCCGCATCGTACAGCCCGTACAGTCCCTGCCCCCTGTAGAAAACCGGGGGTTGCGGGGAGTTCAGTTATTTGCGCTTGTTGTCGctta                    <  2:302716/150‑4 (MQ=255)                                                                                                                    aCGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCgg            <  1:272479/150‑1 (MQ=255)                                                                                                                      gTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGcc          >  1:95113/1‑150 (MQ=255)                                                                                                                              aGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAgg  <  2:7420/150‑1 (MQ=255)                                                                                                                              aGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAgg  <  2:95113/150‑1 (MQ=255)                                                                                                                                 ccTCGACGA‑GATCATCGTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAgg  <  2:230768/147‑1 (MQ=255)                                                                                                                                 ccTCGACGA‑GATCATCGTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAgg  >  1:230768/1‑147 (MQ=255)                                                                                                                                          |                                                                                                                                           AAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGG  >  NC_002947/474058‑474332

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