Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,350,937	T→C	7.0%	intergenic (+30/‑128)	PP_2945 → / → arfB	two‑component system sensor histidine kinase/response regulator/peptidyl‑tRNA hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,350,937	0	T	C	7.0%	163.7 / 2.5	57	intergenic (+30/‑128)	PP_2945/arfB	two‑component system sensor histidine kinase/response regulator/peptidyl‑tRNA hydrolase
Reads supporting (aligned to +/- strand):  ref base T (27/26);  new base C (2/2);  total (29/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 = 7.26e-01

TCCCTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCT  >  NC_002947/3350791‑3351069                                                                                                                                                   |                                                                                                                                     tCCCTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCAt                                                                                                                                      <  7:119733/147‑1 (MQ=255)  cccTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATg                                                                                                                                     >  1:167203/1‑147 (MQ=255)   ccTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAg                                                                                                                                  <  4:28581/149‑1 (MQ=255)   ccTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAg                                                                                                                                  >  3:28581/1‑149 (MQ=255)           aTGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGc                                                                                                                          >  5:280506/1‑149 (MQ=255)            tGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTCCCGCTGGCTGGCATGAAGTGCACGGcc                                                                                                                         >  8:290717/1‑149 (MQ=255)             ggCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGccc                                                                                                                        >  8:188329/1‑149 (MQ=255)              gCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCt                                                                                                                       <  2:8841/149‑1 (MQ=255)               cATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGaa                                                                                                                                   >  2:32977/1‑136 (MQ=255)               cATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGaa                                                                                                                                   <  1:32977/136‑1 (MQ=255)                    tAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAg                                                                                                                 >  3:308495/1‑149 (MQ=255)                      aaGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTg                                                                                                               >  3:255443/1‑149 (MQ=255)                           acacGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAggg                                                                                                          >  1:185550/1‑149 (MQ=255)                                         tAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCgg                                                                                            <  2:167203/149‑1 (MQ=255)                                            ccGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCtt                                                                                         >  6:44083/1‑149 (MQ=255)                                                  cAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCa                                                                                    <  6:280506/148‑1 (MQ=255)                                                   aaGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAg                                                                                                            >  2:206048/1‑123 (MQ=255)                                                   aaGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAg                                                                                                            <  1:206048/123‑1 (MQ=255)                                                   aaGCCAGGTCGGAGTGCGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCAc                                                                                   <  7:290717/148‑1 (MQ=255)                                                    aGCCAGGTCGGAGTGCGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCAAGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCAGGCTTGGGTCGGCATGCGCAcgg                                                                                 >  4:48691/1‑147 (MQ=255)                                                           tCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCaaa                                                                          <  7:188329/149‑1 (MQ=255)                                                                   ggAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGATCCCTg                                                                  >  4:212336/1‑149 (MQ=255)                                                                   ggAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTg                                                                  <  5:44083/149‑1 (MQ=255)                                                                   ggAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTg                                                                  >  6:229603/1‑149 (MQ=255)                                                                    gAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCACACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTg                                                                  >  7:189174/1‑148 (MQ=255)                                                                     aacaacGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGt                                                                                                                >  6:14081/1‑101 (MQ=255)                                                                     aacaacGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGt                                                                                                                <  5:14081/101‑1 (MQ=255)                                                                        aacGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATtg                                                             >  7:189471/1‑149 (MQ=255)                                                                        aacGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATtg                                                             >  7:75842/1‑149 (MQ=255)                                                                         acGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATtgt                                                            <  2:185550/149‑1 (MQ=255)                                                                               cGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTAt                                                       <  6:12589/148‑1 (MQ=255)                                                                                    aCAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGa                                                                                                             >  3:299147/1‑89 (MQ=255)                                                                                    aCAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGa                                                                                                             <  4:299147/89‑1 (MQ=255)                                                                                    aCAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATgggg                                                   >  6:209818/1‑147 (MQ=255)                                                                                        cTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCaaa                                                                          >  2:100271/1‑120 (MQ=255)                                                                                        cTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTACTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCaaa                                                                          <  1:100271/120‑1 (MQ=255)                                                                                              gTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGt                                       >  6:211381/1‑149 (MQ=255)                                                                                               tAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTc                                      <  5:229603/149‑1 (MQ=255)                                                                                                  aGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTCCACGGCCCTGCCGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCtgt                                   >  3:220834/1‑149 (MQ=255)                                                                                                   gCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCACGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATgggg                                                   <  5:172611/132‑1 (MQ=255)                                                                                                   gCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCACGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATgggg                                                   >  6:172611/1‑132 (MQ=255)                                                                                                   gCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCACGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATgggg                                                   <  5:176951/132‑1 (MQ=255)                                                                                                   gCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCACGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATgggg                                                   >  6:176951/1‑132 (MQ=255)                                                                                                    cAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCCCGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGCACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGa                                 >  5:285759/1‑149 (MQ=255)                                                                                                    cAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGa                                 <  4:304391/149‑1 (MQ=255)                                                                                                         cgcAACGCGGTGATGGCGCACGCGTTCTACCTCTGGCTGGCATGCAGTGCAAGGCCCTGACGAGGGGAGGGCTCGCTTGGGTGGGCTTGCGGACTGATTCAAAGAACCCTGCTTTGTGCTATGGGGGAGCGGCACAGTCTGTGATAAt                             <  4:220834/148‑1 (MQ=255)                                                                                                            aaCGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCtg                                    <  5:45134/138‑1 (MQ=255)                                                                                                            aaCGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCtg                                    >  6:45134/1‑138 (MQ=255)                                                                                                                   tgATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGAc                   >  1:8463/1‑148 (MQ=255)                                                                                                                        tcgcCCGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGA‑AATGGGCAATGACAAGGGcc            <  8:274955/148‑1 (MQ=255)                                                                                                                        gcgcACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGc             <  4:255443/149‑1 (MQ=255)                                                                                                                              gcgTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATc       <  2:8463/149‑1 (MQ=255)                                                                                                                              gcgTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATc       <  8:189174/149‑1 (MQ=255)                                                                                                                                 ttCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTc                                      >  3:277678/1‑115 (MQ=255)                                                                                                                                 ttCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTc                                      <  4:277678/115‑1 (MQ=255)                                                                                                                                   cTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGc   <  8:75842/148‑1 (MQ=255)                                                                                                                                      ccGCTGGCTGGCATGAAGTGCACGGCCCTGACGATGGGAGGTCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCt  <  3:212336/146‑1 (MQ=255)                                                                                                                                      ccGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCt  <  5:211381/146‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                     TCCCTTTCCTATGGCATTGTAAAGAAACACGGTGGGGAGATAACCGTATCAAGCCAGGTCGGAGTGGGAACAACGTTTCGTGTACAACTGCCGGTAAAGCAGAGCGCAACGCTGTGATGGCGCACGCGTTCTACCGCTGGCTGGCATGAAGTGCACGGCCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCT  >  NC_002947/3350791‑3351069

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