Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	226,066	G→A	8.1%	D98N (GAT→AAT)	PP_0174 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	226,066	0	G	A	8.1%	141.4 / 2.4	49	D98N (GAT→AAT)	PP_0174	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/24);  new base A (2/2);  total (23/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 = 6.45e-01

CGCTGCTGGTCAACCGCCTGGATGAGGCCAACGACCTGCAGGCGGTGTACGTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCCAGCAGCACTACTTCGGCAAAGGCATCGACCCTTCCTATGTGCAGGCGGCGGAGTGTTACCGGCAGGCCTTCGA  >  NC_002947/225926‑226214                                                                                                                                             |                                                                                                                                                     cGCTGCTGGTCAACCGCCTGGATGAGGCCAACGACCTGCAGGCGGTGTACGTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTg                                                                                                                                              <  4:345605/149‑1 (MQ=255)  gctgctGGTCAACCGCCTGGATGAGGCCAACGACCTGCAGGCGGTGTACGTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTg                                                                                                                                              <  6:267809/148‑1 (MQ=255)         gTCAACCGCCTGGATGAGGCCAACGACCTGCAGGCGGTGTACGTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGTTCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCg                                                                                                                                       <  7:43121/148‑1 (MQ=255)         gTCAACCGCCTGGATGAGGCCAACGACCTGCAGGCGGTGTACGTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCGGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCg                                                                                                                                       <  7:48991/148‑1 (MQ=255)            aaCCGCCTGGATGAGGCCAACGACCTGCAGGCGGTGTACGTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCa                                                                                                                                   >  7:269507/1‑149 (MQ=255)            aaCCGCCTGGATGAGGCCAACGACCTGCAGGCGGTGTACGTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCa                                                                                                                                   >  8:228968/1‑149 (MQ=255)               cGCCTGGATGAGGCCAACGACCTGCAGGCGGTGTACGTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCt                                                                                                                                >  3:37327/1‑149 (MQ=255)                       tGAGGCCAACGACCTGCAGGCGGTGTACGTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTgc                                                                                                                           >  1:219046/1‑146 (MQ=255)                       tGAGGCCAACGACCTGCAGGCGGTGTACGTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTgc                                                                                                                           <  2:219046/146‑1 (MQ=255)                           gCCAACGACCTGCAGGCGGTGTACGTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATggc                                                                                                                    >  7:95998/1‑149 (MQ=255)                             cAACGACCTGCAGGCGGTGTACGTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATggcgg                                                                                                                  >  7:33436/1‑149 (MQ=255)                                  aCCTGCAGGCGGTGTACGTCGCCCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGCCCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTg                                                                                                             >  6:21969/1‑149 (MQ=255)                                  aCCTGCAGGCGGTGTACGTCGACCTGCGGCGTTAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTg                                                                                                             <  2:73014/149‑1 (MQ=255)                                   ccTGCAGGCGGTGTACGTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCGCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGc                                                                                                            <  2:253181/149‑1 (MQ=255)                                       cAGGCGGTGTACGTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAggg                                                                                                        <  7:214386/149‑1 (MQ=255)                                         ggCGGTGTACGTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCa                                                                                                      <  2:170510/149‑1 (MQ=255)                                             gtgtACGTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAAtgc                                                                                                  >  4:290548/1‑149 (MQ=255)                                                   gTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAAtgctgctgc                                                                                            >  5:338544/1‑149 (MQ=255)                                                         cTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGtt                                                                                      >  3:235696/1‑149 (MQ=255)                                                            cggcgTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCaa                                                                                   <  8:269507/149‑1 (MQ=255)                                                                 tCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTg                                                                              <  7:228968/149‑1 (MQ=255)                                                                   aGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGc                                                                            <  4:37327/149‑1 (MQ=255)                                                                     gCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGc                                                                            >  6:199052/1‑147 (MQ=255)                                                                     gCATTACAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCATGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCca                                                                          >  7:237788/1‑149 (MQ=255)                                                                        ttGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCcagca                                                                       <  5:199052/149‑1 (MQ=255)                                                                          gCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCcagcag                                                                      <  8:95998/148‑1 (MQ=255)                                                                            aaGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGc                                                                            <  6:66196/140‑1 (MQ=255)                                                                            aaGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGc                                                                            >  5:66196/1‑140 (MQ=255)                                                                                 aGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCAATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCCAGCAGCACTACtt                                                              >  7:204014/1‑149 (MQ=255)                                                                                  gCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCCAGCAGCACTACTTc                                                             <  3:290548/149‑1 (MQ=255)                                                                                   cGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCAATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCCAGCAGCACTACTTCg                                                            <  8:204014/149‑1 (MQ=255)                                                                                                aaCGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCAcc                                                                                                                                 >  4:160246/1‑67 (MQ=255)                                                                                                aaCGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCAcc                                                                                                                                 <  3:160246/67‑1 (MQ=255)                                                                                                aaCGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGtt                                                                                      <  4:260007/110‑1 (MQ=255)                                                                                                aaCGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGtt                                                                                      >  3:260007/1‑110 (MQ=255)                                                                                                aaCGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCCAGCAGCACTACTTCGGCAAAGGCATCGa                                               >  1:232076/1‑149 (MQ=255)                                                                                                         ggCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCAATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAggg                                                                                                        >  5:357198/1‑83 (MQ=255)                                                                                                         ggCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCAATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAggg                                                                                                        <  6:357198/83‑1 (MQ=255)                                                                                                          gCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCCAGCAGCACTACTTCGGCAAAGGCATCGACCCTTCCTAt                                     <  4:235696/149‑1 (MQ=255)                                                                                                              gATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCCAGCAGCACTACTTCGGCAAAGGCATCGACCCTTCCTAtgtg                                  >  1:9954/1‑148 (MQ=255)                                                                                                              gATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCCAGCAGCACTACTTCGGCAAAGGCATCGACCCTTCCTAtgtg                                  >  4:311038/1‑148 (MQ=255)                                                                                                                      tCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCCAGCAGCACTACTTCGGCAAAGGCATCGACCCTTCCTATGTGCAGgcggcg                         >  3:244610/1‑149 (MQ=255)                                                                                                                      tCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCCAGAAGCACTACTTCGGCAAAGGCATCGACCCTTCCTATGTGCAggcgg                           >  3:336481/1‑147 (MQ=255)                                                                                                                         aCGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCCAGCAGCACTACTTCGGCAAAGGCATCGACCCTTCCTATGTGCAGGCGGCGGAg                      >  2:188061/1‑149 (MQ=255)                                                                                                                            gCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCCAGCAGCACTACTTCGGCAAAGGCATCGACCCTTCCTATGTGCAGGCGGCGGAgtgt                   <  1:188061/149‑1 (MQ=255)                                                                                                                                gTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCCAGCAGCACTACTTCGGCAAAGGCATCGACCCTTCCTATGTGCAGGCGGCGGAGTGTTAcc               <  3:311038/149‑1 (MQ=255)                                                                                                                                      gCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCCAGCAGCACTACTTCGGCAAATGCATCGACCCTTCCTATGTGCAGGCGGGAGAGTGTTACCGGCAg          <  4:336481/148‑1 (MQ=255)                                                                                                                                          gCCGAGACCGTGCTGGCCGGGCACCGGGTGCGCATGGCGGCCCTGCAGGGCAATGCTGCGGCCTGGTTCAACCTTGGCCAGCAGCACTACTTCGGCAAAGGCATCGACCCTTCCTATGTGCAGGCGGCGGAGTGTTACCGGCAGGCCtt     <  8:237788/149‑1 (MQ=255)                                                                                                                                           ccGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCCAGCAGCACTACTTCGGCAAAGGCATCGACCCTTCCTATGTGCAGGCGGCGGAGTGTTACCGGCAGGCCTTc    <  6:338544/149‑1 (MQ=255)                                                                                                                                             gATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCCAGCAGCACTACTTCGGCAAAGGCATCGACCCTTCCTATGTGCAGGCGGCGGAGTGTTACCGGCAGGCCTTCGa  >  4:299353/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                     CGCTGCTGGTCAACCGCCTGGATGAGGCCAACGACCTGCAGGCGGTGTACGTCGACCTGCGGCGTCAGGCATTGCAAGGCAGCGTGTCCGCATTGAACGACCTCGGCTGGATCTGGCTCAACGGCAAGTACTGGCGGGCCGATACCGTGCTGGCCGGGCACCTGTTGCGCATGGCGGCCCTGCAGGGCAATGCTGCTGCCTGGTTCAACCTTGGCCAGCAGCACTACTTCGGCAAAGGCATCGACCCTTCCTATGTGCAGGCGGCGGAGTGTTACCGGCAGGCCTTCGA  >  NC_002947/225926‑226214

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