Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,592,798	G→C	10.2%	K331N (AAG→AAC) ‡	PP_2270 →	DNA primase/helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,592,798	0	G	C	10.2%	139.9 / 6.6	49	K331N (AAG→AAC) ‡	PP_2270	DNA primase/helicase
Reads supporting (aligned to +/- strand):  ref base G (23/21);  new base C (2/3);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

CCGGGAAACGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTGGCGAACGAGCTGGGAATCATC  >  NC_002947/2592659‑2592942                                                                                                                                            |                                                                                                                                                  ccGGGAAACGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTa                                                                                                                                           >  2:96252/1‑148 (MQ=255)     gAAACGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCACCTTGCGCTTCacaa                                                                                                                                      <  1:155097/149‑4 (MQ=255)      aaaCGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGc                                                                                                                                     <  8:123700/149‑1 (MQ=255)        aCGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTg                                                                                                                                   >  6:122673/1‑149 (MQ=255)         cGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCACCTTGCGCTTCa                                                                                                                                         >  8:211721/1‑142 (MQ=255)         cGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCACCTTGCGCTTCa                                                                                                                                         >  4:245377/1‑142 (MQ=255)         cGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCACCTTGCGCTTCa                                                                                                                                         <  7:211721/142‑1 (MQ=255)         cGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCACCTTGCGCTTCa                                                                                                                                         <  3:245377/142‑1 (MQ=255)          gCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGtc                                                                                                                                 >  8:320596/1‑149 (MQ=255)              ccACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTgcgc                                                                                                                                             >  7:127423/1‑133 (MQ=255)              ccACGTTCCCAAGGACACCGCTGGCTGGACCGATAAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTgcgc                                                                                                                                             <  8:127423/133‑1 (MQ=255)                aCGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATgg                                                                                                                                      >  5:986/1‑138 (MQ=255)                     cccAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCAt                                                                                                                     >  8:170805/1‑149 (MQ=255)                      ccAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCAt                                                                                                                     >  8:327952/1‑148 (MQ=255)                       cAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCa                                                                                                                   <  8:290646/149‑1 (MQ=255)                        aaGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATc                                                                                                                    >  4:263898/1‑147 (MQ=255)                          ggACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGc                                                                                                                <  3:263898/149‑1 (MQ=255)                           gACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATgg                                                                                                                                      <  7:311799/127‑1 (MQ=255)                           gACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATgg                                                                                                                                      >  8:311799/1‑127 (MQ=255)                            acacCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCtt                                                                                                              <  5:122673/149‑1 (MQ=255)                              acCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTAt                                                                                                            <  7:274294/149‑1 (MQ=255)                                         gACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGac                                                                                                  <  7:170805/148‑1 (MQ=255)                                                gaggagCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACa                                                                                                       >  1:13094/1‑136 (MQ=255)                                                gaggagCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACa                                                                                                       <  2:13094/136‑1 (MQ=255)                                                             cAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGac                                                                             <  7:327952/149‑1 (MQ=255)                                                                       tGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCTCCCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGa                                                                     <  1:96252/148‑1 (MQ=255)                                                                          gcTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATgag                                                                 >  7:245792/1‑148 (MQ=255)                                                                                cGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCt                                                          >  8:171656/1‑149 (MQ=255)                                                                                 gAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTc                                                         >  5:198425/1‑149 (MQ=255)                                                                                 gAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTc                                                         >  1:155289/1‑149 (MQ=255)                                                                                          gTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGc                                                 <  2:155289/148‑1 (MQ=255)                                                                                             aTGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTa                                                                                                             <  7:19530/85‑1 (MQ=255)                                                                                             aTGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTa                                                                                                             >  8:19530/1‑85 (MQ=255)                                                                                                    acgCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCTATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGaaa                                      <  7:320596/149‑1 (MQ=255)                                                                                                     cgCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGa                                        >  2:26013/1‑146 (MQ=255)                                                                                                       cATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAagag                                   <  7:171656/149‑1 (MQ=255)                                                                                                        aTTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGCTGAGAAGGGCTCGCTTGAGCAGATCATGaaa                                      >  2:87908/1‑145 (MQ=255)                                                                                                        aTTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGaaa                                      >  1:62866/1‑145 (MQ=255)                                                                                                                    cGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTgg                      <  6:198425/149‑1 (MQ=255)                                                                                                                               ttGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATga                                                                  <  7:257950/94‑1 (MQ=255)                                                                                                                               ttGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATga                                                                  >  8:257950/1‑94 (MQ=255)                                                                                                                                tgtgAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTGGCGAACGAGCTgg          <  8:245792/149‑1 (MQ=255)                                                                                                                                tgtgAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTGGCGAACGAGCTgg          >  7:101730/1‑149 (MQ=255)                                                                                                                                    aaGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGCTGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTGGCGAACGAGCTGGGaa       <  1:87908/148‑1 (MQ=255)                                                                                                                                    aaGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAg                                                  <  3:266216/105‑1 (MQ=255)                                                                                                                                    aaGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAg                                                  >  4:266216/1‑105 (MQ=255)                                                                                                                                    aaGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTGGCGAACGAGCTGGGAat      >  2:115262/1‑149 (MQ=255)                                                                                                                                       cgcAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTGGCGCACGAGCTGGGAatcat   >  4:253819/1‑149 (MQ=255)                                                                                                                                        gcAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTGGCGAACGAGCTGGGAatcatc  <  1:115262/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                  CCGGGAAACGCTTCCACGTTCCCAAGGACACCGCTGGCTGGACCGATGAGGAGCTGGATGCAGCCGTAGATGCGCTGGGCGAAAACCTGGTCATGTACGACGCATTTGGTGAAACCGAGTGGGACATTGTGAAGCGCAAGGTGCGCTACATGGCTGTCT‑CCGAGGGCATCAAGCTTATCTACATTGACCACCTCACTGCGATGGCTGACACCGCAGATGAGAAGGGCTCGCTTGAGCAGATCATGAAAGAGATGGCGGGCCTGGCGAACGAGCTGGGAATCATC  >  NC_002947/2592659‑2592942

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