Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,997,189:1	+G	8.7%	intergenic (‑158/+42)	PP_5252 ← / ← PP_5253	hypothetical protein/arylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,997,189	1	.	G	8.7%	168.1 / 11.4	46	intergenic (‑158/+42)	PP_5252/PP_5253	hypothetical protein/arylesterase
Reads supporting (aligned to +/- strand):  ref base . (17/25);  new base G (2/2);  total (19/27)
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.36e-01

CTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATC  >  NC_002947/5997044‑5997336                                                                                                                                                   |                                                                                                                                                    cTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACA‑ACg                                                                                                                                                  >  7:197867/1‑149 (MQ=255)      ggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGt                                                                                                                                             <  7:215779/149‑1 (MQ=255)       gCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACagacgttg                                                                                                                                               >  7:76859/1‑140 (MQ=255)       gCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑Ac                                                                                                                                                   >  8:240205/1‑142 (MQ=255)       gCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑Ac                                                                                                                                                   <  7:240205/142‑1 (MQ=255)        ccGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAg                                                                                                                                           <  3:27173/149‑1 (MQ=255)              gggggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGT‑TGTAGGAGCgg                                                                                                                                     <  8:76859/149‑1 (MQ=255)                  gACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCAACAAAA‑ACGGCTGTAGGAGCGGGttt                                                                                                                                 >  4:58230/1‑149 (MQ=255)                     tGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTAcc                                                                                                                              >  2:54152/1‑149 (MQ=255)                      gCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTAccc                                                                                                                             >  2:3175/1‑149 (MQ=255)                       cGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTCAATCGATGTGCGGTGGGTTGAAGAAGGTTAGACCTGGACAAGCTTTAATCTGCCTCGGACGGCCTCTTCGCGGCTAAACCCCCTCCTACA‑ACGTCTGTAGTAGCCGGTTTaaacg                                                                                                                            >  5:295930/1‑145 (MQ=255)                                 cGGTTGATGCAAGGCCTGGGCTTGAATGGATGGGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTc                                                                                                                  <  3:111242/149‑1 (MQ=255)                                        tGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAAt                                                                                                           <  1:3175/149‑1 (MQ=255)                                                cTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTg                                                                                                   <  3:58230/149‑1 (MQ=255)                                                         aaTGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGc                                                                                          >  1:193898/1‑149 (MQ=255)                                                               tgtgCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGc                                                                                      >  4:81652/1‑147 (MQ=255)                                                                ttgCGGTGGGTGTAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGTAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCaa                                                                                    <  8:197867/147‑1 (MQ=255)                                                                  gctgTGGTTTGTAGTAGGTTAGCCCTGGTCAATCTTTTTGGTGGCTGTGCCGGCGTGTTGGGGTGTAAACCCGCTCCTAGA‑ACGTCTGTTGGTGAGGGTTTACCCGCGAATGTTTCCGACAATTACCTCTGTATGTACGCCAGCAAGt                                                                                  <  6:295930/145‑1 (MQ=255)                                                                     gTGGGTTTAAGAAGTTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTc                                                                              <  1:274841/149‑1 (MQ=255)                                                                      tGGGTTGAAGAAGGTTAGCCCTGGGCAAGCTTTTTGCGGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCa                                                                             <  4:93649/149‑1 (MQ=255)                                                                          ttGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCa                                                                         <  2:193898/149‑1 (MQ=255)                                                                             aagaagGTTAGCCTTGGTCAAGCTTTTTGCTGCGTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCACACAATTACCGCTGAAGGAACGCCAGCaa                                                                                    <  5:38062/135‑1 (MQ=38)                                                                                aagGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAgctgct                                                                   <  3:81652/149‑1 (MQ=255)                                                                                 agGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGctgctg                                                                  <  1:54152/149‑1 (MQ=255)                                                                                    atAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTAAGCTGCTgcg                                                                >  1:109532/2‑148 (MQ=255)                                                                                      aGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTc                                                                                 >  5:243628/1‑129 (MQ=255)                                                                                      aGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTc                                                                                 <  6:243628/129‑1 (MQ=255)                                                                                        cccTGTCCAAGGTTTTTGCGCCCTGGGCCGGCCTCTGCGCGGGAAACCCCGCCCCAACA‑CCGTCTTTCGTCGGGGTTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGa                                                           <  5:39336/149‑1 (MQ=255)                                                                                             gTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCAc                                                      >  2:53208/1‑149 (MQ=255)                                                                                             gTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCAc                                                      >  6:80910/1‑149 (MQ=255)                                                                                                   cTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCAc                                                      <  4:114574/143‑1 (MQ=255)                                                                                                   cTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCAc                                                      >  3:114574/1‑143 (MQ=255)                                                                                                    tttttGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCAt                                                             <  3:293849/135‑1 (MQ=255)                                                                                                    tttttGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCCCTCCTACA‑ACGTGTGTAGGAGCGGGTTTACCCGCGAATGTTTACCACAATTACCGCGGAAGGAACGCCCGCAAGTCCTCATTAAGCTTCTGCGCAt                                                             >  4:293849/1‑135 (MQ=37)                                                                                                             ccTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGTAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgcg                                        <  3:170120/147‑1 (MQ=255)                                                                                                               tGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTTCCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCaa                                                 >  5:325320/1‑136 (MQ=255)                                                                                                               tGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCaa                                                 <  6:325320/136‑1 (MQ=255)                                                                                                                     ggccTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGc                              <  2:318544/149‑1 (MQ=255)                                                                                                                      gactctTCGCGCGTATACCCGCTCCTACT‑ACGTCTGTAGGTGCGGGTTTACCCGCGAATGTTTCCCACCATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGGCACTGCTAAACCGTGCGGCGCCCCGGCa                             <  6:321759/147‑1 (MQ=255)                                                                                                                                gggTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgc                                         <  3:178128/127‑1 (MQ=255)                                                                                                                                gggTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgc                                         >  4:182389/1‑127 (MQ=255)                                                                                                                                gggTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgc                                         >  4:178128/1‑127 (MQ=255)                                                                                                                                gggTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgc                                         <  3:182389/127‑1 (MQ=255)                                                                                                                                        ccGCTCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCGGCGCATGAGTCACTGCAAACCCGTGCGGCGCCACGGCCTACACCTTCAGCTCggag           >  5:182491/1‑147 (MQ=38)                                                                                                                                          gCTCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgc            >  7:205013/1‑146 (MQ=255)                                                                                                                                          gCTCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgc            <  8:205013/146‑1 (MQ=255)                                                                                                                                                 cA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATc  <  2:109532/149‑1 (MQ=255)                                                                                                                                                  a‑aCGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAAt   >  2:8146/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATC  >  NC_002947/5997044‑5997336

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