Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,858,432	G→A	7.7%	noncoding (62/97 nt) noncoding (63/97 nt)	PP_mr49 ← PP_mr50 ←	Bacteria_small_SRP ncRNA Bacteria_large_SRP ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,858,432	0	G	A	7.7%	148.8 / 2.5	52	noncoding (62/97 nt) noncoding (63/97 nt)	PP_mr49 PP_mr50	Bacteria_small_SRP ncRNA Bacteria_large_SRP ncRNA
Reads supporting (aligned to +/- strand):  ref base G (22/26);  new base A (2/2);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.12e-01

CGAAATTCCAGCAGCTGATGAAGCTGTCGGCGGCGGGAGCATTGGGCGAAAAGCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAA  >  NC_002947/4858284‑4858577                                                                                                                                                     |                                                                                                                                                  cGAAATTCCAGCAGCTGATGAAGCTGTCGGGGGCGGGAGCATTGGGCGAAAAGCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCg                                                                                                                                                   >  1:31195/1‑149 (MQ=255) cGAAATTCCAGCAGCTGATGAAGCTGTCGGCGGCGGGAGCATTGGGCGAAAAGCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCg                                                                                                                                                   <  4:103582/149‑1 (MQ=255)    aaTTCCAGCAGCTGATGAAGCTGTCGGCGGCGGGAGCATTGGGCGAAAAGCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTa                                                                                                                                                <  2:99171/149‑1 (MQ=255)        ccAGCAGCTGATGAAGCTGTCGGCGGCGGGAGCATTGGGCGAAAAGCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATgg                                                                                                                                             <  2:303868/148‑1 (MQ=255)             agcTGATGAAGCTGTCGGCGGCGGGAGCATTGGGCGAAAAGCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAGATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTc                                                                                                                                       <  3:220941/149‑1 (MQ=255)             agcTGATGAAGCTGTCGGCGGCGGGAGCATTGGGCGAAAAGCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTc                                                                                                                                       >  1:95660/1‑149 (MQ=255)                     aaGCTGTCGGCGGCGTGAGCATTGGGCGAAAAGCAGGCAAGAACGAAGGCGTTCGTTACACATTGAGGTTTTCCGACCCCAGAAATGCAGGCGGCACCCCCAGCCACACCCCGGCACTCGATGTAACCGCTATGGCTGCAccaagaccg                                                                                                                               >  2:287792/1‑142 (MQ=255)                      aGCTGTCGGCGGCGGGAGCATTGGGCGAAAAGCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCgg                                                                                                                               <  7:128066/148‑1 (MQ=255)                                  cggGAGCATTGGGCGAAAAGCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCACTATGGCTGCTCCCTTCCGGGc                                                                                                                             <  6:6728/138‑1 (MQ=255)                                  cggGAGCATTGGGCGAAAAGCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCACTATGGCTGCTCCCTTCCGGGc                                                                                                                             <  6:12043/138‑1 (MQ=255)                                  cggGAGCATTGGGCGAAAAGCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCACTATGGCTGCTCCCTTCCGGGCCAGatcggg                                                                                                                    >  5:6728/1‑142 (MQ=255)                                  cggGAGCATTGGGCGAAAAGCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCACTATGGCTGCTCCCTTCCGGGCCAGatcggg                                                                                                                    >  5:12043/1‑142 (MQ=255)                                     gAGCATTGGGCGAAAAGCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAAc                                                                                                               <  2:95660/149‑1 (MQ=255)                                      aGCATTGGGCGAAAAGCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGCCCTGCCCAGGTTAACg                                                                                                              >  4:72853/1‑149 (MQ=255)                                                gAAAAGCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAAt                                                                                                    >  6:41933/1‑149 (MQ=255)                                                gAAAAGCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAAt                                                                                                    >  2:205064/1‑149 (MQ=255)                                                     gCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCggg                                                                                               <  5:41933/149‑1 (MQ=255)                                                      cAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCgggg                                                                                              <  6:346471/149‑1 (MQ=255)                                                       aGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCggggg                                                                                             <  7:303708/149‑1 (MQ=255)                                                         gCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAAc                                                                                                               >  2:45891/1‑129 (MQ=255)                                                         gCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAAc                                                                                                               <  1:45891/129‑1 (MQ=255)                                                                    ggCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGt                                                                                                                   <  8:198758/114‑1 (MQ=255)                                                                    ggCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATAGCTGCTCCCTTCCGGGCCTGACCAGGt                                                                                                                   >  7:198758/1‑114 (MQ=255)                                                                      cgttcgttACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGccac                                                                              <  2:31195/149‑1 (MQ=255)                                                                         tcgttACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGgct                                                                                                                                           >  3:296157/1‑85 (MQ=255)                                                                         tcgttACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGgct                                                                                                                                           <  4:296157/85‑1 (MQ=255)                                                                           gttACAGATTGATGGTTTCCGACCCCACAAATGGAGGCGGCCCCACCACCCACCCCCCGGCACACGATGTTCCCGCTTTGGCTGGTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACt                                                                          <  3:72853/148‑1 (MQ=255)                                                                             tACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACa                                                                       >  6:18552/1‑149 (MQ=255)                                                                                         ttttCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGag                                                           >  2:229866/1‑149 (MQ=255)                                                                                            tCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGagca                                                         <  2:339285/148‑1 (MQ=255)                                                                                                    cAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGgc                                                >  8:91088/1‑149 (MQ=255)                                                                                                       aaaTGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGgcgc                                              >  8:292054/1‑148 (MQ=255)                                                                                                       aaaTGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGcc                                             <  3:257145/149‑1 (MQ=255)                                                                                                                cggcCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCg                                    <  5:18552/149‑1 (MQ=255)                                                                                                                  gcCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACTTGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGt                                  <  1:44937/149‑1 (MQ=255)                                                                                                                  gcCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCAATACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGt                                  >  1:313068/1‑149 (MQ=255)                                                                                                                   ccccACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTc                                 >  4:235721/1‑149 (MQ=255)                                                                                                                       accaGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTg                             <  5:240631/149‑1 (MQ=255)                                                                                                                       accaGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTg                             >  6:240631/1‑149 (MQ=255)                                                                                                                       accaGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTg                             <  3:235721/149‑1 (MQ=255)                                                                                                                        ccaGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTg                             >  7:223106/1‑148 (MQ=255)                                                                                                                           gCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGc                         <  5:153295/149‑1 (MQ=255)                                                                                                                           gCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGc                         <  1:315161/149‑1 (MQ=255)                                                                                                                            ccACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAg                          >  6:258982/1‑147 (MQ=255)                                                                                                                              agacCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGATCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGa                      <  6:124494/147‑1 (MQ=255)                                                                                                                                  cccGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTAc                                                                        >  6:128399/1‑95 (MQ=255)                                                                                                                                  cccGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTAc                                                                        <  5:128399/95‑1 (MQ=255)                                                                                                                                    cGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAAc                <  7:91088/149‑1 (MQ=255)                                                                                                                                    cGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAAc                <  1:205064/149‑1 (MQ=255)                                                                                                                                          tCGAAGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACAACATGTCTCTCTGAGCAGCGAGAGGCGCAATTGTACCGGTCTTTGGAGCAGATACAACAACTGc          >  1:256222/1‑149 (MQ=255)                                                                                                                                            gATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCgg        >  6:224189/1‑149 (MQ=255)                                                                                                                                            gATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCgg        >  1:42717/1‑149 (MQ=255)                                                                                                                                               gTTCCCGATATGGCTGCTCCCTTCCGGGCATGACCAGGTTTACGGGTAATCAATGCGGGGGGTCCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCa      <  2:256222/148‑1 (MQ=255)                                                                                                                                                   ccGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATAAAACCACTGCGGCAGaaa  >  5:38138/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CGAAATTCCAGCAGCTGATGAAGCTGTCGGCGGCGGGAGCATTGGGCGAAAAGCAGGCAAGAACGAAGGCGTTCGTTACAGATTGAGGTTTTCCGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAA  >  NC_002947/4858284‑4858577

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