Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,783,197:1	+T	8.5%	coding (226/927 nt)	PP_2436 ←	LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,783,197	1	.	T	8.5%	173.9 / 12.2	47	coding (226/927 nt)	PP_2436	LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (22/21);  new base T (2/2);  total (24/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 = 8.57e-01

GGCGCCAGGCGGCGCAGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTGGGTAATGCCGACGATCTG  >  NC_002947/2783052‑2783341                                                                                                                                                   |                                                                                                                                                 ggCGCCAGGCGGCGCAGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑Gtt                                                                                                                                               <  3:293787/149‑1 (MQ=255)    gcCAGGCGGCGCAGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGACAGCGC‑GTTgcg                                                                                                                                            >  8:133543/1‑149 (MQ=255)    gcCAGGCGGCGCAGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAAC‑GCGCTGGCGCTGTTgcg                                                                                                                                            <  4:45612/149‑1 (MQ=255)          cggcgCAGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTGCGCAAC‑GCGCTGGCGCTGTTGCGgacatc                                                                                                                                      <  4:33605/149‑3 (MQ=255)                aGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCgat                                                                                                                                <  2:217486/149‑1 (MQ=255)                    gcTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGatgatg                                                                                                                            >  1:159325/1‑149 (MQ=255)                                      aGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGtgtg                                                                                                          >  1:39814/1‑149 (MQ=255)                                       gATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTgg                                                                                                         <  2:39814/149‑1 (MQ=255)                                        aTGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTggg                                                                                                        <  4:297666/149‑1 (MQ=255)                                           aCCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCAc                                                                                                     <  2:159325/149‑1 (MQ=255)                                                  cGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCt                                                                                              >  1:285370/1‑149 (MQ=255)                                                    accgaggcggtcATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCtt                                                                                             <  1:328104/141‑1 (MQ=255)                                                      cAGGTCGGTCATACTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTgcg                                                                                          >  5:77178/1‑149 (MQ=255)                                                          tcggtcATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGt                                                                                      <  5:64558/149‑1 (MQ=255)                                                           cggtcATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTa                                                                                     >  1:13651/1‑149 (MQ=255)                                                            ggtcATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTa                                                                                     <  6:216252/148‑1 (MQ=255)                                                              tcATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACg                                                                                   >  2:168696/1‑148 (MQ=255)                                                                     gATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGaa                                                                             >  8:102677/1‑147 (MQ=255)                                                                         cGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCgg                                                                       <  7:173581/149‑1 (MQ=255)                                                                         cGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCgg                                                                       <  4:288648/149‑1 (MQ=255)                                                                          ggAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCggg                                                                      >  4:164387/1‑149 (MQ=255)                                                                            aaGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTc                                                                    <  6:77178/149‑1 (MQ=255)                                                                              ggTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCAt                                                                  >  3:111436/1‑149 (MQ=255)                                                                                  ttgttgGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTg                                                                >  4:133329/1‑147 (MQ=255)                                                                                     ttgGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGt                                                           >  1:131876/1‑149 (MQ=255)                                                                                             ctgcGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCAc                                                                                <  5:204291/120‑1 (MQ=255)                                                                                             ctgcGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCAc                                                                                >  6:204291/1‑120 (MQ=255)                                                                                                       aaCGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAg                                         >  5:251881/1‑149 (MQ=255)                                                                                                        aCGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAgg                                        <  4:82203/149‑1 (MQ=255)                                                                                                        aCGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAgg                                        <  7:102677/149‑1 (MQ=255)                                                                                                          gTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGgcgc                                              <  7:176773/141‑1 (MQ=255)                                                                                                          gTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGgcgc                                              >  8:176773/1‑141 (MQ=255)                                                                                                           tGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCa                                     >  3:266175/1‑149 (MQ=255)                                                                                                              ggCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTg                                  <  2:285370/149‑1 (MQ=255)                                                                                                                cTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGa                                <  2:131876/149‑1 (MQ=255)                                                                                                                   tCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGac                               <  2:13651/147‑1 (MQ=255)                                                                                                                     cTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTg                           >  4:177740/1‑149 (MQ=255)                                                                                                                      tGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGc                          >  6:222081/1‑149 (MQ=255)                                                                                                                           cGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCt                     >  6:143419/1‑149 (MQ=255)                                                                                                                            gATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTg                    <  4:266175/149‑1 (MQ=255)                                                                                                                             ctgtgCGCAAC‑GCGCTGGCGCTGTTGCGGACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTgg                   >  5:301037/2‑149 (MQ=255)                                                                                                                             aTGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTgg                   >  8:181732/1‑149 (MQ=255)                                                                                                                              tgtgCGCAAC‑GCGCTGGCGCTGTTGCGGACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTggg                  >  4:248116/1‑149 (MQ=255)                                                                                                                                tCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTGGGt                 >  6:200063/1‑148 (MQ=255)                                                                                                                                  cGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTGGGTAAt              >  1:174441/1‑149 (MQ=255)                                                                                                                                  cGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTGGGTAAt              <  5:267087/149‑1 (MQ=255)                                                                                                                                    cAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTGGGTAATGc            <  1:168696/149‑1 (MQ=255)                                                                                                                                              agccc‑GTTTCCCACAGGCCCGAGCATGTTCTCGGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGAGGGTCATAGAAGGATTCGCGCAGGCGCGACAGTGCATTCGACACTGCCGGCTGGGTAATGCCGACGATCTg  <  6:70417/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GGCGCCAGGCGGCGCAGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGC‑GTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTGGGTAATGCCGACGATCTG  >  NC_002947/2783052‑2783341

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