Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,796,627	C→T	16.6%	intergenic (‑185/+33)	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,627	0	C	T	16.6%	112.6 / 13.0	51	intergenic (‑185/+33)	pvdH/PP_4224	diaminobutyrate‑2‑oxoglutarate transaminase/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (17/24);  new base T (5/3);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-01

GGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGCCAGCCCC  >  NC_002947/4796479‑4796774                                                                                                                                                     |                                                                                                                                                        ggACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTc                                                                                                                                                         <  4:266032/149‑1 (MQ=255)   aattGGGTTATTCCTCGCGAGCAAGTGGATGTGCGGTTATGGCTTTGGTTACGCTTGAGTGGGGGGAGGA‑TTAGCGGTTGTTTGGCGATTTTTTTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGTGCGTGTTCAc                                                                                                                                                       <  6:26491/145‑1 (MQ=37)    cATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCAcc                                                                                                                                                      <  3:102097/149‑1 (MQ=255)     aTGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTAccc                                                                                                                                                     >  4:237659/1‑149 (MQ=38)     aTGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTAccc                                                                                                                                                     <  3:237659/149‑1 (MQ=38)     aTGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCAccc                                                                                                                                                     <  3:64820/149‑1 (MQ=255)            aaTCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGaa                                                                                                                                                >  4:9254/1‑147 (MQ=255)            aaTCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                                              <  8:209238/149‑1 (MQ=255)                cTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGc                                                                                                                                           <  8:25034/148‑1 (MQ=255)                  cgcgAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGGTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAg                                                                                                                                        <  3:63646/149‑1 (MQ=255)                       gCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTAACCCGCGAAGAGGCCAGCACAg                                                                                                                                   >  3:250334/1‑149 (MQ=255)                        cAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAgg                                                                                                                                  >  7:237142/1‑149 (MQ=255)                           gCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑a                                                                                                                              <  2:193440/149‑1 (MQ=255)                                 gggCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAAc                                                                                                                         <  8:237142/148‑1 (MQ=255)                                          tCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGccc                                                                                                               <  7:20718/149‑1 (MQ=255)                                             cTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTCCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATa                                                                                                                            <  1:124341/133‑1 (MQ=255)                                              tttGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGGTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCAGCCCGTAGGAGCGGGTTCACCCGCCAAGAGGCCAGCACCGGCG‑ATCAACATCAATGCCCAGGa                                                                                                           >  2:83755/1‑149 (MQ=39)                                                  gataCGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑cgc                                                                                                      <  8:92369/145‑1 (MQ=255)                                                     aCGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCGCAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCg                                                                                                  <  3:203624/149‑1 (MQ=255)                                                        cTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGCATAAACATCAATGCCCAGGA‑‑CGCGCCGGt                                                                                                <  1:22781/149‑1 (MQ=255)                                                                gggggAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGcaa                                                                                        <  3:9254/148‑1 (MQ=255)                                                                                  cgatctgcgATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCCAGCCGT‑Ca                                                                    <  5:238136/143‑1 (MQ=255)                                                                                         cgATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAAAGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑Ca                                                                    >  6:238136/1‑142 (MQ=255)                                                                                          gATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCt                                                            <  7:21660/149‑1 (MQ=255)                                                                                                 tttttGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTGACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGC‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAg                                                     >  5:238140/1‑149 (MQ=255)                                                                                                 ttatcgccTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGACG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCaaa                                                      <  5:122522/141‑1 (MQ=18)                                                                                                 ttatcgccTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGACG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCa                                                        >  6:122522/8‑146 (MQ=18)                                                                                                    ttGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAgcg                                                   <  3:268080/148‑1 (MQ=255)                                                                                                         tgtgCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGc                                              >  1:247181/1‑148 (MQ=255)                                                                                                           tgCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGACCCGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAg                                                <  1:83755/146‑1 (MQ=255)                                                                                                              tGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCgg                                        >  3:110106/1‑149 (MQ=255)                                                                                                                     ttCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCg                                 >  5:64110/1‑149 (MQ=37)                                                                                                                          gggTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCa                            >  8:267097/1‑149 (MQ=255)                                                                                                                           ggTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTc                                                         >  8:10259/1‑119 (MQ=25)                                                                                                                           ggTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTc                                                         <  7:10259/119‑1 (MQ=25)                                                                                                                            gTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCAc                           >  3:256212/1‑148 (MQ=255)                                                                                                                             tgaaCCCGCTCCTACGGGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGa                         >  3:258864/3‑149 (MQ=255)                                                                                                                             tAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGcaa                                                                                        >  7:95748/1‑87 (MQ=255)                                                                                                                             tAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGcaa                                                                                        <  8:95748/87‑1 (MQ=255)                                                                                                                             tAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGcaa                                                                                        <  4:253706/87‑1 (MQ=255)                                                                                                                             tAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGcaa                                                                                        >  3:253706/1‑87 (MQ=255)                                                                                                                              gaaCCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGAt                        >  3:219976/2‑149 (MQ=38)                                                                                                                              gaaCCCGCTCCTACGGGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGc                             <  7:267585/143‑1 (MQ=255)                                                                                                                              gaaCCCGCTCCTACGGGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGc                             >  8:267585/2‑144 (MQ=255)                                                                                                                                aCCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCg                      >  6:99531/1‑149 (MQ=255)                                                                                                                                 cccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATc                       <  4:169292/147‑1 (MQ=255)                                                                                                                                 cccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATc                       >  3:169292/1‑147 (MQ=255)                                                                                                                                  ccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCCACTACGGCTGTCAAAGCGCAGGCTACCGGCacac                                   <  1:196979/135‑1 (MQ=255)                                                                                                                                  ccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCGGCacac                                   >  2:196979/1‑134 (MQ=255)                                                                                                                                  ccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGcc                    >  5:267430/1‑149 (MQ=255)                                                                                                                                     cccgTAGGAGCGGGTTCACCCGCTAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGc                 <  7:244505/149‑1 (MQ=255)                                                                                                                                                   ttCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGCCAGccc   >  8:248766/1‑149 (MQ=255)                                                                                                                                                     cACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGCCAGcccc  <  6:267430/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                        GGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCG‑ATAAACATCAATGCCCAGGA‑‑CGCGCCGGTACCTGCAACAAAACCCGCAAGCCGT‑CACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGCCAGCCCC  >  NC_002947/4796479‑4796774

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