Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,873,999	A→G	8.3%	intergenic (+47/+84)	PP_5147 → / ← PP_5148	haloacid dehalogenase‑like family hydrolase/membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,873,999	0	A	G	8.3%	132.6 / 2.8	48	intergenic (+47/+84)	PP_5147/PP_5148	haloacid dehalogenase‑like family hydrolase/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (23/21);  new base G (2/2);  total (25/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 = 2.55e-01

TTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGCCCCCATCAGCCCGGCAATGGACAGGAAGCAAACCCCACTGAATACC  >  NC_002947/5873854‑5874143                                                                                                                                                  |                                                                                                                                                  ttgCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATc                                                                                                                                                <  3:53968/149‑1 (MQ=255)      gCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATcacca                                                                                                                                            >  3:9956/1‑148 (MQ=255)            gACGCATGTAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCt                                                                                                                                     >  4:249611/1‑149 (MQ=255)               gCATGTAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGa                                                                                                                                  <  3:249611/149‑1 (MQ=255)                 aTGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGGATCACCACAGGCCTTGAgg                                                                                                                                >  1:201951/1‑149 (MQ=255)                 aTGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGGATCACCACAGGCCTTGAgg                                                                                                                                >  1:205320/1‑149 (MQ=255)                   gcagtggcagtgGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGATTTCCCGTGCAGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACCGGCCTTGAGG‑Gc                                                                                                                             >  1:81968/1‑149 (MQ=255)                                 cccGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGt                                                                                                               <  8:177158/149‑1 (MQ=255)                                     gCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAg                                                                                                             >  5:46877/1‑147 (MQ=255)                                              cTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCa                                                                                                  >  5:242052/1‑149 (MQ=255)                                                gCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCAcc                                                                                                >  1:247690/1‑149 (MQ=255)                                                                 aGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGca                                                                                <  2:247690/148‑1 (MQ=255)                                                                   cGTGGGTGGCCGGTGATGTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACaa                                                                             >  7:225034/1‑149 (MQ=255)                                                                     tgggtggCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATc                                                                           <  2:148607/149‑1 (MQ=255)                                                                              ggTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAgg                                                                  <  5:105430/149‑1 (MQ=255)                                                                                tGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGGATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCt                                                                <  2:205320/149‑1 (MQ=255)                                                                                tGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGGATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCt                                                                <  2:201951/149‑1 (MQ=255)                                                                                 gATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTg                                                               >  7:99632/1‑149 (MQ=255)                                                                                       gCTGCGGGATTGAGGTTTTCCCGTGCAGGTCTTATCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCa                                                            <  4:214819/146‑1 (MQ=255)                                                                                       gCTGCGGGATTGAGGGTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGa                                                          <  8:82406/148‑1 (MQ=255)                                                                                           cGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGc                                                     >  7:156770/1‑149 (MQ=255)                                                                                            gggATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCt                                                    >  8:164189/1‑149 (MQ=255)                                                                                                  gAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGcc                                               <  8:225034/148‑1 (MQ=255)                                                                                                   aGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGcccc                                             >  4:63859/1‑149 (MQ=255)                                                                                                    ggTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGccccc                                            <  8:156770/149‑1 (MQ=255)                                                                                                     gTTTTCCGGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACATGCCTTGATG‑GCTGCTCTGTACCTGTAGGAGCGGGCTCACCCGCGAATATGCCGGCACAATCAGCACAAGGCTGTCAGACCTGCTTGgccgcc                                            <  2:81968/148‑1 (MQ=38)                                                                                                     gTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGGCGCTGCGCTGTACCCGTAGGAGCGGGCTc                                                                                                   <  1:152731/94‑1 (MQ=255)                                                                                                     gTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCCGTAGGAGCGGGCTc                                                                                                   >  2:152731/1‑93 (MQ=255)                                                                                                      ttttCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGCCCCCAt                                          >  3:17493/1‑149 (MQ=255)                                                                                                        ttCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGCCCCCATc                                         <  3:63859/148‑1 (MQ=255)                                                                                                         tCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGCCCCCATCa                                        <  7:164189/148‑1 (MQ=255)                                                                                                              tGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGCCCCCATCAGCCCgg                                  <  6:46877/149‑1 (MQ=255)                                                                                                                 cGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGCCCCCATCAGCCCGGCaa                               >  5:93826/1‑149 (MQ=255)                                                                                                                  ggTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTc                                                             >  7:247384/1‑118 (MQ=255)                                                                                                                  ggTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTc                                                             <  8:247384/118‑1 (MQ=255)                                                                                                                   gTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCGCCCGCGAATGGGCCGGCACACGCAGCAGACGGGGGGCAGCCCGGCTTGGCGCGCGGAAGCCCGGCATTg                             >  5:55477/1‑149 (MQ=21)                                                                                                                      ttttCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCATCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGCCCCCATCAGCCCGGCAATGGa                           >  6:173778/1‑148 (MQ=255)                                                                                                                      ttttCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCATCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGCCCCCATCAGCCCGGCAATGGa                           >  8:186108/1‑148 (MQ=255)                                                                                                                      ttttCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGCCCCCATCAGCCCGGCAATGGa                           >  8:46849/1‑148 (MQ=255)                                                                                                                         tCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGGGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGCCCCCATCAGCCCGGCAATGGACAg                        >  7:17871/1‑148 (MQ=255)                                                                                                                           gcgGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGTCCCCCATCAGCCCGGCAATGGACAGGaa                     <  7:46849/149‑1 (MQ=255)                                                                                                                            cgGGCACGCCCGCTCCCCCAGGAATCACCCCAGGCCTGGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACACGGCTGTCAGACCGGCTTGGCCCCCATCACCCCGGCAATGGACAGGAAg                    <  6:249641/149‑1 (MQ=255)                                                                                                                                     ccGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGCCCCCATCCGGCCGGCAATGGACAGGAAGCAAACCCCa           >  4:119152/1‑149 (MQ=255)                                                                                                                                        cTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGccc                                              <  8:82470/111‑1 (MQ=255)                                                                                                                                        cTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGccc                                              >  7:82470/1‑111 (MQ=255)                                                                                                                                           ccACATGAATCACCACAGGCCTTGAGG‑GCTGCGGTGTACCTGTAGGTTGGGGCTCACCCGCGAATATGCCGGCAAAATCATCACAAGGCTGTCAGACCGGCTTGGCCCCCATCAGCCCGGCAATGGACAGGAAGCAAACCCCACTGaa      <  4:9956/148‑1 (MQ=255)                                                                                                                                           ccACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGCCCCCATCAGCCCGGCAATGGACAGGAAGCAAACCCCACTGaa      >  3:179837/1‑148 (MQ=255)                                                                                                                                              caGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGCCCCCATCAGCCCGGCAATGGACAGGAAGCAAACCCCACTGAATAcc  <  6:242052/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  TTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGG‑GCTGCGCTGTACCTGTAGGAGCGGGCTCACCCGCGAATAGGCCGGCACAATCAGCACAAGGCTGTCAGACCGGCTTGGCCCCCATCAGCCCGGCAATGGACAGGAAGCAAACCCCACTGAATACC  >  NC_002947/5873854‑5874143

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