Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	128,152	C→T	8.0%	S210S (AGC→AGT)	polA →	DNA polymerase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	128,152	0	C	T	8.0%	143.1 / 2.7	50	S210S (AGC→AGT)	polA	DNA polymerase I
Reads supporting (aligned to +/- strand):  ref base C (24/22);  new base T (2/2);  total (26/24)
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.43e-01

GCGTGCACGAGAAATTCGGCGTCGGCCCGGAACACATTATCGACTTCCTGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTGGCCACCATCAAGGTTGACGTGCCGCTGGATGTGGAGG  >  NC_002947/128004‑128300                                                                                                                                                     |                                                                                                                                                      gcGTGCACGAGAAATTCGGCGTCGGCCCGGAACACATTATCGACTTCCTGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGc                                                                                                                                                       <  8:130745/149‑1 (MQ=255)       aCGAGAAATTCGGCGTCGGCCCGGAACACATTATCGACTTCCTGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCtt                                                                                                                                                 >  3:173220/1‑149 (MQ=255)         gagaAATTCGGCGTCGGCCCGGAACACATTATCGACTTCCTGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTa                                                                                                                                               >  6:230221/1‑149 (MQ=255)          agaAATTCGGCGTCGGCCCGGAACACATTATCGACTTCCTGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTAc                                                                                                                                              <  1:130644/149‑1 (MQ=255)              aTTCGGCGTCGGCCCGGAACACATTATCGACTTCCTGGCCCTGATGGGTGACAAGGTCGGCAACATCCCCGGCGTGCCGGGCGTGGTCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACgac                                                                                                                                           >  5:270187/1‑146 (MQ=255)              aTTCGGCGTCGGCCCGGAACACATTATCGACTACCTGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGAGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTCTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGcc                                                                                                                                           >  7:23747/1‑148 (MQ=255)               ttCGGCGTCGGCCCGGAACACATTATCGACTTCCTGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCaa                                                                                                                                         >  7:209254/1‑149 (MQ=255)                     gTCGGCCCGGAACACATTATCGACTTCCTGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGa                                                                                                                                   <  5:4027/149‑1 (MQ=255)                     gTCGGCCCGGAACACATTATCGACTTCCTGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGa                                                                                                                                   <  8:252214/149‑1 (MQ=255)                        ggCCCGGAACACATTATCGACTTCCTGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACa                                                                                                                                 >  1:76374/1‑148 (MQ=255)                             ggAACACATTATCGACTTCCTGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTAACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGc                                                                                                                           >  6:150270/1‑149 (MQ=255)                               aaCACATTATCGACTTCCTGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCa                                                                                                                         <  6:250636/149‑1 (MQ=255)                                  acaTTATCGACTTCCTGGCCCTGATGGGTGACAAGGTCGACAACGTCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAgcg                                                                                                                      >  1:56004/1‑149 (MQ=255)                                     ttATCGACTTCCTGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTg                                                                                                                   <  8:209254/149‑1 (MQ=255)                                     ttATCGACTTCCTGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTg                                                                                                                   <  5:122070/149‑1 (MQ=255)                                     ttATCGACTTCCTGGCCCTGATGGGTGACAAGGGCGACAACATCCGCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTg                                                                                                                   <  8:23747/149‑1 (MQ=255)                                             ttCCTGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCa                                                                                                                         >  8:102506/1‑135 (MQ=255)                                             ttCCTGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCa                                                                                                                         <  7:102506/135‑1 (MQ=255)                                                 tGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCgcggc                                                                                                       >  1:132533/1‑149 (MQ=255)                                                       tGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGcc                                                                                                    >  4:29930/1‑146 (MQ=255)                                                              tGACAAGGTCGACAACATCCCCGGCGGGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCGGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGc                                                                                         <  4:173220/149‑1 (MQ=255)                                                                  aaGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGc                                                                                     >  4:97268/1‑149 (MQ=255)                                                                             caTCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTggagg                                                                          >  7:17483/1‑149 (MQ=255)                                                                                  ccGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCAc                                                                     >  4:85321/1‑149 (MQ=255)                                                                                  ccGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCAc                                                                     >  6:275285/1‑149 (MQ=255)                                                                                         gCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGAcg                                                              >  8:104108/1‑149 (MQ=255)                                                                                          ccGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGAcgc                                                             >  3:80418/1‑149 (MQ=255)                                                                                           cGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGTCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGGCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGAcgcg                                                            <  2:56004/149‑1 (MQ=255)                                                                                           cGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTAACCGGCATCGGCGGCGGGCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGACCAAGCTGGAGGAGCACCGCGAcgcg                                                            <  5:150270/149‑1 (MQ=255)                                                                                           cGGGCGTGGGCGAAAAAACGGCGGGGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGAcgcg                                                            <  7:175689/149‑1 (MQ=255)                                                                                                 tGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTc                                                      <  8:17483/149‑1 (MQ=255)                                                                                                 ggggCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCATGACCt                                                                                             <  8:103249/110‑1 (MQ=255)                                                                                                      gAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTc                                                 <  6:10179/149‑1 (MQ=255)                                                                                                             cggcggTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCCAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGaa                                          >  7:223266/1‑149 (MQ=255)                                                                                                               gcggTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACt                                        >  7:185823/1‑149 (MQ=255)                                                                                                                cggTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTg                                       <  3:29930/149‑1 (MQ=255)                                                                                                                    ggggCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTGGcca                                   <  5:230221/149‑1 (MQ=255)                                                                                                                        cTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTGGccacca                                >  1:48417/1‑148 (MQ=255)                                                                                                                             gACCGGCATCGGCGGCGGTCTCAGTGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTGGccac                                  >  1:64197/1‑141 (MQ=255)                                                                                                                             gACCGGCATCGGCGGCGGTCTCAGTGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTGGccac                                  >  1:67584/1‑141 (MQ=255)                                                                                                                             gACCGGCATCGGCGGCGGTCTCAGTGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTGGccac                                  <  2:64197/141‑1 (MQ=255)                                                                                                                             gACCGGCATCGGCGGCGGTCTCAGTGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTGGccac                                  <  2:67584/141‑1 (MQ=255)                                                                                                                             gACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTGGCCACCATCAAg                           <  7:104108/148‑1 (MQ=255)                                                                                                                               ccGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTGGCCACCATCAAGGtt                        >  6:10349/1‑149 (MQ=255)                                                                                                                                cGGTATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTGGCCACCATCAAGGtt                        >  4:246336/1‑148 (MQ=255)                                                                                                                                cGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTGGccacc                                 >  5:20021/1‑139 (MQ=255)                                                                                                                                cGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTGGccacc                                 <  6:20021/139‑1 (MQ=255)                                                                                                                                   cATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTGGCCACCATCAAGGTTGACg                    <  5:10349/149‑1 (MQ=255)                                                                                                                                              gTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTGGCCACCATCAAGGTTGACGTGCCGCTGGAt         <  2:48417/149‑1 (MQ=255)                                                                                                                                                 tcAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCTAAGCCGG‑CAAGCTGGAGGAGCACCCCGACGCGGCATTACAGTCCGACGAAAAGGGCCCCAAAAAGGTTGACCGGCCGCTTGAtgtt      >  1:215430/1‑148 (MQ=255)                                                                                                                                                   aGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTGGCCACCATCAAGGTTGACGTGCCGCTGGATGTGGa    <  5:275285/149‑1 (MQ=255)                                                                                                                                                     cGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTGGCCACCATCAAGGTTGACGTGCCGCTGGATGTGGAgg  <  5:68065/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GCGTGCACGAGAAATTCGGCGTCGGCCCGGAACACATTATCGACTTCCTGGCCCTGATGGGTGACAAGGTCGACAACATCCCCGGCGTGCCGGGCGTGGGCGAAAAAACGGCGGTGGGGCTGTTGACCGGCATCGGCGGCGGTCTCAGCGACCTTTACGCCAACCTGGACAAAGTGCCAGCGCTGGCCATTCGCGGCGCCAAGACCT‑TGCCGGCCAAGCTGGAGGAGCACCGCGACGCGGCATTCCTGTCCTACGAACTGGCCACCATCAAGGTTGACGTGCCGCTGGATGTGGAGG  >  NC_002947/128004‑128300

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