Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,336,499	A→G	8.1%	intergenic (‑43/+37)	PP_2052 ← / ← PP_5500	bifunctional sugar‑phosphatase/mannitol‑1‑phosphate 5‑dehydrogenase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,336,499	0	A	G	8.1%	135.1 / 2.8	49	intergenic (‑43/+37)	PP_2052/PP_5500	bifunctional sugar‑phosphatase/mannitol‑1‑phosphate 5‑dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/18);  new base G (2/2);  total (29/20)
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.86e-01

GCTTGCTTGAGTGTCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAGGCGGGGCGGAGCGGGGCGTGCAGAGGAGGTTG  >  NC_002947/2336356‑2336623                                                                                                                                                |                                                                                                                               gcttgcttGAGTGTCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGc                                                                                                                           >  8:246336/1‑149 (MQ=255)      cttgAGTGTCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGTTGCAGGC‑gagat                                                                                                                     <  5:120881/149‑2 (MQ=255)      cttgAGTGTCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑gagat                                                                                                                     <  4:173484/149‑2 (MQ=255)      cttgAGTGTCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGaa                                                                                                                     >  8:59777/1‑149 (MQ=255)            tgtCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGgtgt                                                                                                               <  7:246336/149‑1 (MQ=255)              tCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGgtgtgt                                                                                                             >  6:348590/1‑149 (MQ=255)                cTGATGCTCAGGCGTTCGTGGTGTAATTGGGGGCAATCAATATCGAACTGGGGGGTTAAATTTTGTTTCCCGTTATTTTAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGa                                                                                                           <  1:332301/149‑1 (MQ=255)                 tGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGgt                                                                                                                 >  4:38367/1‑142 (MQ=255)                 tGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGgt                                                                                                                 <  3:38367/142‑1 (MQ=255)                 tGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAg                                                                                                          >  6:266084/1‑149 (MQ=255)                         aGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCt                                                                                                        <  4:79014/143‑1 (MQ=255)                         aGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCt                                                                                                        >  3:79014/1‑143 (MQ=255)                           gcgcTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTggg                                                                                                <  3:35428/149‑1 (MQ=255)                               tCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCa                                                                                            <  7:59777/149‑1 (MQ=255)                                  gtgtCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTggg                                                                                                >  8:123928/1‑142 (MQ=255)                                  gtgtCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTggg                                                                                                <  7:123928/142‑1 (MQ=255)                                     tCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGCAGAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAg                                                                                      <  1:60409/149‑1 (MQ=255)                                       gAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGTTGCAGGC‑gaga                                                                                                                      <  4:70009/115‑1 (MQ=255)                                       gAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGTTGCAGGC‑gaga                                                                                                                      >  3:70009/1‑115 (MQ=255)                                       gAAGAGCGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTGGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGt                                                                                    <  5:292876/148‑1 (MQ=255)                                           agGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCgggg                                                                               <  5:348590/149‑1 (MQ=255)                                           agGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCgggg                                                                               <  5:266084/149‑1 (MQ=255)                                               tGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTTAGATGCAGTGGGGGCA‑GTCAGGTCGGGGGca                                                                            >  7:277792/1‑148 (MQ=255)                                                gCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAgg                                                                                                                            >  1:70119/1‑101 (MQ=255)                                                gCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAgg                                                                                                                            <  2:70119/101‑1 (MQ=255)                                                   atccatATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGt                                                                       >  3:134356/1‑149 (MQ=255)                                                          tCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAg                                                                >  3:339204/1‑149 (MQ=255)                                                           cGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGa                                                                   >  8:171412/1‑145 (MQ=255)                                                           cGAAGATGGCGATGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGa                                                                   <  7:171412/145‑1 (MQ=255)                                                            gactaTGGCGCGGAAAATCTGTTTACCGTGAACGTACAGCCTTGGCGGTTGCCCCTTGTCGTTGCGGCAGGCATCTCCCCTGCAACAGGC‑GAGAACGGTGTGTGCGCTGCAGTGGGGGCA‑GTCAGGGCGGGGGCACACGTTCGACAggg                                                              <  8:286626/145‑1 (MQ=255)                                                             aaGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAgggg                                                             >  1:10841/1‑149 (MQ=255)                                                                         aaaaTCTGTTTACCGTGATAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTAGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTAGACAGGGGGACGCTCAGGAg                                                 >  4:308545/1‑149 (MQ=255)                                                                            aTCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAGGGGGAc                                                          >  2:310759/1‑137 (MQ=255)                                                                            aTCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAGGGGGAc                                                          <  1:310759/137‑1 (MQ=255)                                                                                    aCCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCa                                          >  1:60265/1‑145 (MQ=255)                                                                                    aCCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTAGCCAGGCGGAGGATCAGGAGTGCCGCAGAgg                                      >  7:317650/1‑149 (MQ=255)                                                                                        tGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAGGGGGACGc                                                        <  4:68017/127‑1 (MQ=255)                                                                                        tGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAGGGGGACGc                                                        >  3:68017/1‑127 (MQ=255)                                                                                        tGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAg                                  >  4:27148/1‑149 (MQ=255)                                                                                                aaGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAggcggggcg                          >  5:228326/1‑149 (MQ=255)                                                                                                         cTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAGGCGGGGCGGAGCGGGGc                 >  6:149035/1‑149 (MQ=255)                                                                                                         cTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGGGAGCTGCAGTGGGGGGA‑GTCAGCTCCGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAGGCGGGGCGGAGCGGGGc                 >  1:110911/1‑149 (MQ=255)                                                                                                               ccccTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCagag                                       >  4:45027/1‑121 (MQ=255)                                                                                                               ccccTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCagag                                       <  3:45027/121‑1 (MQ=255)                                                                                                                 ccTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAGGCGGGGCGGAGCGGGGCGTGCagag         >  7:48921/1‑149 (MQ=255)                                                                                                                   ttGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCATGTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAggcgggg                            >  1:122448/1‑129 (MQ=255)                                                                                                                         ttGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGa                                                  >  1:148554/1‑100 (MQ=255)                                                                                                                                accgttctcGCCTGTTGCAGGC‑GAGATCGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCag                                         >  3:316204/6‑102 (MQ=38)                                                                                                                                              cAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAGGCGGGGCGGAGCGGGGCGTGCAGAGGAGGTTg  >  3:144109/1‑127 (MQ=255)                                                                                                                                                |                                                                                                                               GCTTGCTTGAGTGTCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCAGGCATCTCGCCTGCAACAGGC‑GAGAACGGTGTGTGAGCTGCAGTGGGGGCA‑GTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAGGCGGGGCGGAGCGGGGCGTGCAGAGGAGGTTG  >  NC_002947/2336356‑2336623

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