Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	914,644:1	+C	9.8%	coding (1355/1743 nt)	fruA →	PTS fructose transporter subunit IIBC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	914,644	1	.	C	9.8%	190.2 / 14.6	51	coding (1355/1743 nt)	fruA	PTS fructose transporter subunit IIBC
Reads supporting (aligned to +/- strand):  ref base . (25/21);  new base C (2/3);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-01

ATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCAA  >  NC_002947/914497‑914779                                                                                                                                                        |                                                                                                                                        aTGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑g                                                                                                                                        <  5:143105/149‑1 (MQ=255)       accTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTAGGCGTTCTCGGTGGGGCTGCTGGCCGCGTCGAGCTACGCGCCAAAC‑GGCGGC‑ACCATGCCGGCGGC‑CATGGTCgccgcc                                                                                                                                   <  1:234079/148‑1 (MQ=37)          tgctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGAtt                                                                                                                               >  1:241373/1‑149 (MQ=255)            ctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTgg                                                                                                                             >  2:163969/1‑149 (MQ=255)            ctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGC‑ACCATGCCGGCGGC‑CATGGTCgccgccattg                                                                                                                               <  5:211059/147‑5 (MQ=38)               gggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCtt                                                                                                                          >  8:92499/1‑149 (MQ=255)               gggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCtt                                                                                                                          <  6:195889/149‑1 (MQ=255)               gggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCtt                                                                                                                          >  5:193922/1‑149 (MQ=255)               gggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCtt                                                                                                                          <  1:163969/149‑1 (MQ=255)                 gCTGTTGCTGGGCGGCATGATGTGCGTAGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTAGGTTGGGCTGCAGGGCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCCG‑CAAGGG‑GCCGGCGAttcggatgg                                                                                                                        >  8:333047/1‑142 (MQ=255)                     ttGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTAt                                                                                                                    >  5:305315/1‑149 (MQ=255)                                   gatgTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTgg                                                                                                      >  4:191459/1‑149 (MQ=255)                                         cGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGcc                                                                                                                >  6:227139/1‑133 (MQ=255)                                         cGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGcc                                                                                                                <  5:227139/133‑1 (MQ=255)                                              aCCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGTGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAg                                                                                           >  2:226178/1‑149 (MQ=255)                                                cTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGtt                                                                                         >  7:71543/1‑149 (MQ=255)                                                          cGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagc                                                                               >  4:151684/1‑149 (MQ=255)                                                          cGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagc                                                                               >  8:296304/1‑149 (MQ=255)                                                                acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAgcgc                                                                         <  7:92499/149‑1 (MQ=255)                                                                acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAgcgc                                                                         <  3:191459/149‑1 (MQ=255)                                                                acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCCAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAgcgc                                                                         >  1:42111/1‑149 (MQ=255)                                                                       cgccTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCgg                                                                  <  4:214909/149‑1 (MQ=255)                                                                          cTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGc                                                                >  1:282831/1‑148 (MQ=255)                                                                                     gggTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCgg                                                                  <  6:313021/134‑1 (MQ=255)                                                                                      cgTGGGGCGGCTGGCCGCGGCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATCGT‑GCCGCCGATTGGCGTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCggcgctgg                                                   <  5:345637/148‑1 (MQ=255)                                                                                            gctgctGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGcgctggcgctgg                                             >  5:103108/1‑149 (MQ=255)                                                                                             ctgctgGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTggg                                            <  3:169858/149‑1 (MQ=255)                                                                                             ctgctgGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTggg                                            >  7:193104/1‑149 (MQ=255)                                                                                                ctgGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTgggg                                           <  6:305315/147‑1 (MQ=255)                                                                                                ctgGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCt                                         >  4:139962/1‑149 (MQ=255)                                                                                                ctgGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCt                                         >  5:278485/1‑149 (MQ=255)                                                                                                    ccTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGc                                     >  1:46850/1‑149 (MQ=255)                                                                                                    ccTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCAGGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGc                                     <  2:282831/149‑1 (MQ=255)                                                                                                        gtcgAGCTACGCGCC‑AAT‑GGCTGAGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGTGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCa                                 >  8:162960/1‑149 (MQ=255)                                                                                                        gtcgAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGc                                                                <  2:24594/118‑1 (MQ=255)                                                                                                        gtcgAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGc                                                                >  1:24594/1‑118 (MQ=255)                                                                                                        gtcgAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTc                                  >  2:241025/1‑148 (MQ=255)                                                                                                        gtcgAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTc                                  <  7:296304/148‑1 (MQ=255)                                                                                                        gtcgAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTc                                  <  1:241025/148‑1 (MQ=255)                                                                                                              cttcgcgcg‑ATT‑GGCTGGGACCATGGCGGCGTG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCTCCTTCCTGGCGCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCAtctc                             <  7:283308/144‑1 (MQ=255)                                                                                                                ccaagcC‑AATCGGCGGC‑ACCATGCCGGCGGC‑CATGGTCGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGCGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTcc                            >  3:257347/5‑147 (MQ=25)                                                                                                                    gcC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGTGTATTGCCACCTTCCTGGCCCGTCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAggggg                     <  2:241373/149‑1 (MQ=255)                                                                                                                      c‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGc                    >  2:319086/1‑148 (MQ=255)                                                                                                                               ggcgACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGt            <  2:42111/149‑1 (MQ=255)                                                                                                                                ggcaCCATGCCGGCGGC‑CATGGTCGCCGCC‑ATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGtt           >  5:132805/4‑149 (MQ=38)                                                                                                                                gcgACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGtt           <  3:173678/149‑1 (MQ=255)                                                                                                                                      catgccggcGGCCATGGTCGCCGCCAATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCggg                                                                 <  6:191596/88‑1 (MQ=37)                                                                                                                                      aTGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCGCCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTgccgc     >  1:248257/1‑149 (MQ=255)                                                                                                                                         gccgccGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCaa  <  3:139962/149‑1 (MQ=255)                                                                                                                                         gccgccGG‑CATGGT‑GCCGCCGATTGGCTGGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCaa  <  6:193922/149‑1 (MQ=255)                                                                                                                                                  cATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGa                  >  2:179368/1‑125 (MQ=255)                                                                                                                                                  cATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGa                  <  1:179368/125‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                        ATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCC‑AAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCAA  >  NC_002947/914497‑914779

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