Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,796,614	G→C	11.8%	intergenic (‑172/+46)	pvdH ← / ← PP_4224	diaminobutyrate‑2‑oxoglutarate transaminase/sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,796,614	0	G	C	11.8%	138.4 / 10.7	51	intergenic (‑172/+46)	pvdH/PP_4224	diaminobutyrate‑2‑oxoglutarate transaminase/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (17/28);  new base C (2/4);  total (19/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

CTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCA  >  NC_002947/4796468‑4796760                                                                                                                                                   |                                                                                                                                                    cTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTa                                                                                                                                                   <  5:132453/149‑1 (MQ=255)      aGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGc                                                                                                                                              <  7:90426/149‑1 (MQ=255)         ggTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTAcgggcggg                                                                                                                                           <  2:244557/149‑1 (MQ=38)           tGGACATGGGTAAATCCTCGCGCGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGtt                                                                                                                                         <  5:243303/149‑1 (MQ=255)           tGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGtt                                                                                                                                         <  2:9332/149‑1 (MQ=255)            ggACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGATACGGTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCTATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCAGCCCGTAGGTGCGGGtt                                                                                                                                         >  7:327400/1‑148 (MQ=255)              aCATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCAc                                                                                                                                      <  6:118775/149‑1 (MQ=255)                aTGGGTATCTCCTCGAGAGCAAGGGGTTGGGCGGTTTTCGCCTTGGAAACGCTTGAGTGGGGGGAGGATTTATCGGTTGTTACGCGTGTTTATTTTTCTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTAccc                                                                                                                                    <  8:110674/149‑1 (MQ=21)                aTGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCAccc                                                                                                                                    <  5:47929/149‑1 (MQ=255)                 tGGGTAAATCCTCGCGAGCAAGCGGATGGGCGTTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCcg                                                                                                                                   <  5:21023/149‑1 (MQ=37)                  gggTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTAGTTAAGCGATTTTATTTTTGTATGTGATGGCCTCTTAGCTGGTAAACCCGCCCGTAGGAGAGGTTTCaaactc                                                                                                                                  >  3:117144/1‑144 (MQ=255)                   ggTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCcgc                                                                                                                                  <  8:323116/148‑1 (MQ=255)                    gTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCCCCcgcg                                                                                                                                 >  1:140807/1‑148 (MQ=255)                       aaTCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAa‑ga                                                                                                                            <  5:290817/149‑1 (MQ=255)                        aTCCTCTCGAGCAATAGGATGGGAGGTTATCGCTTTGGAAAGGGTTGAGTGGGGGGTGGATTTAGGGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAa‑gag                                                                                                                           <  4:117144/149‑1 (MQ=38)                         tCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTa                                                                                                                                                   >  5:373822/1‑125 (MQ=255)                         tCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTa                                                                                                                                                   <  6:373822/125‑1 (MQ=255)                           cTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGc                                                                                                                         <  1:211560/148‑1 (MQ=255)                           cTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGGGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAA‑GAGGc                                                                                                                         <  7:320713/148‑1 (MQ=25)                            tCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCCCCCGCGAA‑GAGGc                                                                                                                         >  2:230176/1‑147 (MQ=255)                                  gCAAGCGGAGGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAg                                                                                                                 <  5:20473/149‑1 (MQ=255)                                    aaGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGc                                                                                                                         >  1:8749/1‑139 (MQ=255)                                    aaGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGc                                                                                                                         <  2:8749/139‑1 (MQ=255)                                       cGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGAt                                                                                                            <  8:179652/149‑1 (MQ=255)                                              gCGGTTATCGCTTTGGAAACGCTTGAGGGGTGGGAGGATTTAGCGGTTGTTCCGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATc                                                                                                     <  8:327400/149‑1 (MQ=255)                                                ggTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCaa                                                                                                   <  4:16622/149‑1 (MQ=255)                                                  ttATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATg                                                                                                 >  6:201777/1‑149 (MQ=255)                                                     tCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCAGAAGGAGCGGGTTCACACGCGAA‑GAGGCCAGCACATGCGATAAACATCAATgac                                                                                               >  2:253172/1‑146 (MQ=255)                                                     tCGCTTTGGAAACGCTTGAGTGGGGGGAGGATGTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGccc                                                                                              <  1:139726/149‑1 (MQ=255)                                                     tCGCTTTGGAAACGCTTGAGTGAGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAA‑GAGGCCAg                                                                                                                      >  7:233723/1‑125 (MQ=21)                                                       gCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAACGAGGCCAGCACAGGCGATAAACATCAATGCCCa                                                                                             <  1:50209/149‑1 (MQ=255)                                                       gCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCa                                                                                             <  3:63586/148‑1 (MQ=255)                                                        cTTTGGAAACGCTTGAGGGGGTGGAGGATTTAGCGTTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGTGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGTA‑GAGGCCAGCACAGGCGATAAACATGAATGCCCAgg                                                                                           <  1:253172/149‑1 (MQ=39)                                                         tttGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGa                                                                                          >  5:290192/1‑149 (MQ=255)                                                         tttGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCCATAAACATAAATTGCCAGGa                                                                                          >  3:356642/1‑149 (MQ=255)                                                              aaaCGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCAGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACgcgc                                                                                     >  1:56618/1‑149 (MQ=255)                                                               aaCGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGcc                                                                                    >  8:182951/1‑149 (MQ=255)                                                                           gggggATGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGcaa                                                                         <  5:201777/148‑1 (MQ=255)                                                                              ggaggaTTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACACGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGcaacaa                                                                      >  1:296495/1‑148 (MQ=255)                                                                                              gttACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCAc                                                     >  5:298412/1‑149 (MQ=255)                                                                                                 aCGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTa                                                   <  6:290192/148‑1 (MQ=255)                                                                                                  cgcgATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCa                                                                                                                       >  7:252098/1‑79 (MQ=255)                                                                                                  cgcgATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCa                                                                                                                       <  8:252098/79‑1 (MQ=255)                                                                                                  cgcgATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACg                                                 <  7:182951/149‑1 (MQ=255)                                                                                                  cgcgATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACg                                                 <  2:336973/149‑1 (MQ=255)                                                                                                            tttttGTCTGTGCTGGCTTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCGGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCaaa                                        <  6:30682/148‑1 (MQ=255)                                                                                                              tttGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAgc                                      <  2:56618/148‑1 (MQ=255)                                                                                                                 gTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCa                                   <  2:296495/148‑1 (MQ=255)                                                                                                                               ctTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACAcc                    >  7:325817/1‑149 (MQ=255)                                                                                                                                ttCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCg                   >  6:196703/1‑149 (MQ=255)                                                                                                                                       gtgaaCCCGCTCCTACGGGCGGGTTTACCCGCGAA‑GAGGCCAGCACAGACGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCAc             >  3:293910/4‑148 (MQ=37)                                                                                                                                            cccgcccgTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGc       >  3:119627/1‑149 (MQ=255)                                                                                                                                            cccgcccgTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGc       >  4:164284/1‑149 (MQ=255)                                                                                                                                              cgcccgTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCt     <  6:298412/149‑1 (MQ=255)                                                                                                                                                 ccgTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGca  >  4:248112/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAA‑GAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCA  >  NC_002947/4796468‑4796760

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