Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	557,336:1	+A	11.4%	intergenic (+79/‑130)	rpsN → / → rpsH	30S ribosomal protein S14/30S ribosomal protein S8

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	557,336	1	.	A	11.4%	148.7 / 16.9	44	intergenic (+79/‑130)	rpsN/rpsH	30S ribosomal protein S14/30S ribosomal protein S8
Reads supporting (aligned to +/- strand):  ref base . (19/20);  new base A (2/3);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01

GGCCGTAACATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAG  >  NC_002947/557189‑557477                                                                                                                                                     |                                                                                                                                              ggCCGTAACATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑g                                                                                                                                              <  3:117596/149‑1 (MQ=255)  gCCGTAACATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑gg                                                                                                                                             <  5:277601/149‑1 (MQ=255)         cATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCTGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTg                                                                                                                                       >  1:188296/1‑148 (MQ=255)         cATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCTGGTCTGGTCAAGGCATCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTg                                                                                                                                       <  2:188296/148‑1 (MQ=255)         cATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCC‑ATATAGGGGCTTgt                                                                                                                                      <  6:248616/149‑1 (MQ=255)          aTGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTgtt                                                                                                                                     >  6:50469/1‑149 (MQ=255)           tgctgcGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTgttt                                                                                                                                    <  7:311/149‑1 (MQ=255)             ctgcGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGCCTTgtttcg                                                                                                                                  >  4:184417/1‑149 (MQ=255)             ctgcGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCC‑ATATAGGGGCTTgtttg                                                                                                                                   <  2:315867/148‑2 (MQ=255)              tgcGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCC‑ATATAGGGGCTTgttt                                                                                                                                    >  4:213336/1‑146 (MQ=255)              tgcGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCC‑ATATAGGGGCTTgttt                                                                                                                                    <  3:213336/146‑1 (MQ=255)                    aaGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGc                                                                                                                           >  3:31361/1‑149 (MQ=255)                     aGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGcc                                                                                                                          >  5:265788/1‑149 (MQ=255)                     aGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGcc                                                                                                                          <  6:118662/149‑1 (MQ=255)                                               gTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGTGTGGAAGGGGAGCAATCTTCCGTCCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCATGTTTCGATTCTGCCTTGTGACTAGAAAGACCGCCTCAcat                                                                                                >  5:225668/1‑147 (MQ=255)                                                    gTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTgtttcgttt                                                                                                                               <  1:8410/113‑1 (MQ=255)                                                    gTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTgtttcgttt                                                                                                                               >  2:8410/1‑113 (MQ=255)                                                     tCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCt                                                                                          >  1:319049/1‑149 (MQ=255)                                                        aGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGTAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTggg                                                                                       <  4:31361/149‑1 (MQ=255)                                                                            cttgCAAGGCGCGGGTGGAAGGGGAGGAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCGTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCgg                                                                   <  6:225668/146‑1 (MQ=255)                                                                                  aGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGtt                                                             <  2:319049/149‑1 (MQ=255)                                                                                   ggCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTg                                                            <  4:90043/149‑1 (MQ=255)                                                                                   ggCGCGGGTGCAAGGGGAGCAATCTTCCGACCGCCGGTGAACTTGCCAACAAACAAGCCGCTATAT‑GGGGCTTGTTTCGTTTCAGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTg                                                            <  3:184417/149‑1 (MQ=255)                                                                                    gcgcGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGt                                                           <  2:7327/149‑1 (MQ=255)                                                                                      gcgGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCggg                                                                  >  1:278313/1‑140 (MQ=255)                                                                                      gcgGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCggg                                                                  <  2:278313/140‑1 (MQ=255)                                                                                               aGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGc                                                  >  7:244235/1‑147 (MQ=255)                                                                                                      cAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCtaat                                         <  5:50469/149‑1 (MQ=255)                                                                                                           ttCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGTGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATa                                                       <  1:34748/130‑1 (MQ=255)                                                                                                           ttCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGAGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATa                                                       >  2:34748/1‑130 (MQ=255)                                                                                                             ccGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCt                                   <  6:265788/148‑1 (MQ=255)                                                                                                              cGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCtt                                  <  3:30275/148‑1 (MQ=255)                                                                                                                  cgccgGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGttttt                                                                                  <  1:90563/96‑1 (MQ=255)                                                                                                                  cgccgGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGttttt                                                                                  >  2:90563/1‑96 (MQ=255)                                                                                                                             ttGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATc                    >  5:109344/1‑147 (MQ=255)                                                                                                                             ttGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATc                    >  5:330962/1‑147 (MQ=255)                                                                                                                             ttGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATc                    >  3:256574/1‑147 (MQ=255)                                                                                                                                     aaacaaGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGa          <  6:109344/149‑1 (MQ=255)                                                                                                                                       acaaGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCtt                                  >  2:274192/1‑123 (MQ=255)                                                                                                                                       acaaGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCtt                                  <  1:274192/123‑1 (MQ=255)                                                                                                                                         aaGCCCC‑ATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTAt      >  6:5217/1‑149 (MQ=255)                                                                                                                                          aGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATg     >  2:136498/1‑149 (MQ=255)                                                                                                                                           gCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATg     >  3:45197/1‑148 (MQ=255)                                                                                                                                             cccTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCAt            >  4:161372/1‑139 (MQ=255)                                                                                                                                             cccTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCAt            <  3:161372/139‑1 (MQ=255)                                                                                                                                             cccTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAg  <  4:256574/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              GGCCGTAACATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAG  >  NC_002947/557189‑557477

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