Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,252,213	T→G	8.9%	intergenic (‑181/+59)	PP_3726 ← / ← rocE	enoyl‑CoA hydratase/isomerase family protein/amino acid permease RocE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,252,213	0	T	G	8.9%	150.3 / 5.2	54	intergenic (‑181/+59)	PP_3726/rocE	enoyl‑CoA hydratase/isomerase family protein/amino acid permease RocE
Reads supporting (aligned to +/- strand):  ref base T (25/24);  new base G (3/2);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GGACTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGCACAGGGCAATGAAGCCA  >  NC_002947/4252066‑4252359                                                                                                                                                    |                                                                                                                                                   ggaCTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTc                                                                                                                                                   <  3:76732/149‑1 (MQ=255)                   aCTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGGTCGAGCGCCGCCTACGCGGTGCGCGATCTca                                                                                                                                            <  4:153041/138‑1 (MQ=37)                   aCTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGGTCGAGCGCCGCCTACGCGGTGCGCGATCTca                                                                                                                                            >  3:153041/1‑138 (MQ=37)                      aaaCGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAAt                                                                                                                              >  7:23876/1‑149 (MQ=255)                       aaCGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGaaaa                                                                                                                               <  7:263917/147‑1 (MQ=255)                        aCGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGc                                                                                                                            >  6:305276/1‑149 (MQ=255)                            gtgtTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTGTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGc                                                                                                                         >  4:344234/1‑148 (MQ=255)                                  aGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCaa                                                                                                                   >  7:244706/1‑148 (MQ=255)                                   gCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAaca                                                                                                                 >  7:147837/1‑149 (MQ=255)                                    cAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAaca                                                                                                                 >  3:256154/1‑148 (MQ=255)                                     aaGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACAcc                                                                                                               >  6:51463/1‑149 (MQ=255)                                     aaGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACCGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACAcc                                                                                                               >  7:16232/1‑149 (MQ=255)                                            gATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCcga                                                                                                         <  2:299321/148‑1 (MQ=255)                                            gATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCcga                                                                                                         >  1:299321/1‑148 (MQ=255)                                                 gCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCcgacgac                                                                                                     >  5:130368/1‑147 (MQ=255)                                                  cTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGCCGACGAg                                                                                                  >  4:182521/1‑149 (MQ=255)                                                  cTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAg                                                                                                  >  4:187169/1‑149 (MQ=255)                                                    gTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCcgacga                                                                                                      >  7:88257/1‑143 (MQ=255)                                                    gTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCcgacga                                                                                                      <  8:88257/143‑1 (MQ=255)                                                     ttACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCt                                                                                                <  5:316050/148‑1 (MQ=255)                                                       aCAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTgc                                                                                              <  6:130368/148‑1 (MQ=255)                                                          aCTTGGCTCAACTGGGTGTCTGGCGTTACTCAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGTTGGCGCACCGCGTAGGCGGCGCTGGACCTCACAGGCGCTGAAAATGAAGCGACCAACACCTGCCGACGAGGAGCTgcggc                                                                                           <  8:16232/148‑1 (MQ=255)                                                           cTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGt                                                                                         <  8:23876/149‑1 (MQ=255)                                                                cacaACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAAc                                                                                    <  8:147837/149‑1 (MQ=255)                                                                   aaCTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAAcccc                                                                                 >  7:206568/1‑149 (MQ=255)                                                                        ggTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGgt                                                                             <  5:305276/148‑1 (MQ=255)                                                                          tgtCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGtgtg                                                                          >  2:198310/1‑149 (MQ=255)                                                                               ggCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTggg                                                                     >  5:333199/1‑149 (MQ=255)                                                                                   ttACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCg                                                                 <  5:51463/149‑1 (MQ=255)                                                                                       acaAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAgg                                                             >  1:113158/1‑149 (MQ=255)                                                                                       acaAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGGCAAACCCCTGGTGTGCTGGGTGCGAAgg                                                             <  6:300500/149‑1 (MQ=255)                                                                                                 tCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGAc                                                                                                                      >  8:219539/1‑82 (MQ=255)                                                                                                 tCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGAc                                                                                                                      <  7:219539/82‑1 (MQ=255)                                                                                                   gTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGa                                                  <  2:113158/148‑1 (MQ=255)                                                                                                       tGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGAcggcg                                             <  8:244706/149‑1 (MQ=255)                                                                                                            gCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGa                                        <  1:236650/149‑1 (MQ=255)                                                                                                                     ccTGTGAGGTCGAGCGCCGCCTACGCGGTGCGCGATCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGaagaag                                >  1:60068/1‑148 (MQ=37)                                                                                                                      cTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTa                              <  6:333199/149‑1 (MQ=255)                                                                                                                      cTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTa                              >  2:348972/1‑149 (MQ=255)                                                                                                                       tgtgAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTAc                             >  1:234177/1‑149 (MQ=255)                                                                                                                         tgAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTAcgc                           <  1:144305/149‑1 (MQ=255)                                                                                                                             gtCGAGCGCCGCCTACGCGGTGCGCGATCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCa                       >  1:166445/2‑149 (MQ=38)                                                                                                                              tCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCa                       >  5:356763/1‑148 (MQ=255)                                                                                                                                gcgcACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATa                     <  3:344234/148‑1 (MQ=255)                                                                                                                                 cgcACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGc                   <  6:356763/149‑1 (MQ=255)                                                                                                                                    aCCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGca                  >  7:212820/1‑147 (MQ=255)                                                                                                                                    aCCGCGTAGGCGGCGCGCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGcaca                <  7:297013/149‑1 (MQ=255)                                                                                                                                      cgcgTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGa                                        <  7:136895/123‑1 (MQ=255)                                                                                                                                      cgcgTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGa                                        >  8:136895/1‑123 (MQ=255)                                                                                                                                       gcgTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGCACAggg             <  4:256154/149‑1 (MQ=255)                                                                                                                                            ggcggcGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGCACAGGGCAATg        <  1:348972/149‑1 (MQ=255)                                                                                                                                             gcggcgCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGCACAGGGCAATg        >  5:308490/1‑148 (MQ=255)                                                                                                                                               ggcgCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGa                                        <  5:244358/114‑1 (MQ=255)                                                                                                                                               ggcgCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGa                                        >  6:244358/1‑114 (MQ=255)                                                                                                                                                    tCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGCACAGGGCAATGAAGCCa  >  3:315084/1‑147 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GGACTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGCACAGGGCAATGAAGCCA  >  NC_002947/4252066‑4252359

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