Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,363,590	G→C	7.1%	S59C (TCC→TGC)	glmM ←	phosphoglucosamine mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,363,590	0	G	C	7.1%	175.7 / 3.7	56	S59C (TCC→TGC)	glmM	phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base G (24/28);  new base C (2/2);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

CATCGTGCGGGTTGTGCGAGGCACTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGC  >  NC_002947/5363447‑5363728                                                                                                                                                |                                                                                                                                             cATCGTGCGGGTGTTGCGAGGCACTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTc                                                                                                                                          <  4:43471/148‑1 (MQ=255)  aTCGTGCGGGTTGTGCGAGGCACTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCg                                                                                                                                         >  1:241702/1‑148 (MQ=255)   tCGTGCGGGTTGTGCGAGGCACTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGaa                                                                                                                                       <  3:904/149‑1 (MQ=255)    cGTGCGGGTTGTGCGAGGCACTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTc                                                                                                                                          <  8:288449/145‑1 (MQ=255)    cGTGCGGGTTGTGCGAGGCACTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTc                                                                                                                                          >  7:288449/1‑145 (MQ=255)     gTGCGGGTTGTGCGAGGCACTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAAc                                                                                                                                      <  5:274860/148‑1 (MQ=255)     gTGCGGGTTGTGCGAGGCACTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTCCGCACTCGaagc                                                                                                                                     <  1:261853/149‑3 (MQ=255)      tGCGGGTTGTGCGAGGCACTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑t                                                                                                                                   <  2:80128/149‑1 (MQ=255)             tgtgCGAGGCACTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCg                                                                                                                            <  5:288623/149‑1 (MQ=255)               tgCGAGGCACTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGCGTGCGGACTCGAACA‑TGTATCCg                                                                                                                            >  3:182325/1‑147 (MQ=255)               tgCGAGGCACTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCgg                                                                                                                           <  1:289671/148‑1 (MQ=255)               tgCGAGGCACTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCgg                                                                                                                           <  5:111945/148‑1 (MQ=255)                 cGAGGCACTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTCCGCACTCGaagc                                                                                                                                     <  6:49748/137‑3 (MQ=255)                 cGAGGCACTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTCCGCACTCGaagc                                                                                                                                     >  5:49748/1‑135 (MQ=255)                  gAGGCACTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑At                                                                                                                      <  2:220432/149‑1 (MQ=255)                        cTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCgt                                                                                                                <  6:77097/149‑1 (MQ=255)                           atgaCGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCgtg                                                                                                               <  2:241702/147‑1 (MQ=255)                           atgaCGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTc                                                                                                             <  5:269244/149‑1 (MQ=255)                                   gCCGGCTTCGGCGTGGAAGGTATGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCcca                                                                                                     >  8:64713/1‑149 (MQ=255)                                      ggCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCccacca                                                                                                  <  1:12749/149‑1 (MQ=255)                                            ggCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCAcgc                                                                                            >  5:254069/1‑149 (MQ=255)                                                  ggaGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCTTGTCCTTGCCCACCAGCACGCGGCAAt                                                                                      <  2:167368/147‑1 (MQ=255)                                                     ggTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGc                                                                                   >  4:80326/1‑149 (MQ=255)                                                     ggTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGc                                                                                   <  6:254069/149‑1 (MQ=255)                                                      gTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGt                                                                                                                                 >  1:185624/1‑103 (MQ=255)                                                      gTACGCGTGAGGTAAGCGATGGACGGTGTCGGCATCTGCCATAGCAGCATCACATCGGCACCTGCAGCGGACAGGCCGGCTTCGAGTGCTGACTCGAACA‑ttt                                                                                                                                 <  2:185624/103‑3 (MQ=255)                                                        aCGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCctt                                                                                >  4:16872/1‑149 (MQ=255)                                                                      cGATGGCCGGTGTCGGCATCGGCCCGAGTAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGc                                                                   >  8:134913/1‑148 (MQ=255)                                                                      cGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGc                                                                   >  1:71454/1‑148 (MQ=255)                                                                             cGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATgccgg                                                           >  6:47618/1‑149 (MQ=255)                                                                             cGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATgccgg                                                           <  4:101518/149‑1 (MQ=255)                                                                                    ggCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCa                                                    >  8:273708/1‑149 (MQ=255)                                                                                     tcATCGGCCCGAGCAGCATCACATCGGCAACTGCCGCGGACAGGCCGTCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCAGCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAg                                                   <  4:182325/148‑1 (MQ=255)                                                                                           gCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGc                                             >  6:72597/1‑149 (MQ=255)                                                                                            cccGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCt                                            >  6:87685/1‑149 (MQ=255)                                                                                                   agcaTCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCAt                                     >  6:268875/1‑149 (MQ=255)                                                                                                        cacaTCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGt                                >  8:149388/1‑149 (MQ=255)                                                                                                         acaTCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATTGTAGCCGGAG‑ATACGCGTGTCCTTGCCCCCCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGt                                >  1:72330/1‑149 (MQ=255)                                                                                                           aTCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAGAATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGt                                >  1:241978/1‑147 (MQ=255)                                                                                                           aTCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAg                             >  3:175110/1‑149 (MQ=255)                                                                                                              ggCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGc                           <  5:87685/148‑1 (MQ=255)                                                                                                              ggCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGc                           <  3:16872/148‑1 (MQ=255)                                                                                                              ggCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGc                           <  7:273708/148‑1 (MQ=255)                                                                                                               gCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGt                         >  2:292366/1‑149 (MQ=255)                                                                                                                   cTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATc                     <  5:72597/149‑1 (MQ=255)                                                                                                                    tGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCAcg                                                                                             <  8:185103/76‑1 (MQ=255)                                                                                                                    tGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCAcg                                                                                             >  7:185103/1‑76 (MQ=255)                                                                                                                     gCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCgg                   <  4:175110/149‑1 (MQ=255)                                                                                                                     gCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGATTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCgg                   >  4:178266/1‑149 (MQ=255)                                                                                                                      ccGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCgg                   <  5:47618/148‑1 (MQ=255)                                                                                                                           gACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTc             >  7:108588/1‑149 (MQ=255)                                                                                                                                   atgttCGAGTCCGCACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGAcgc     >  3:90167/4‑149 (MQ=255)                                                                                                                                     cTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCgg   <  1:292366/149‑1 (MQ=255)                                                                                                                                        cGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGc  <  7:149388/147‑1 (MQ=255)                                                                                                                                           gTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGccc                                                     <  4:259585/93‑1 (MQ=255)                                                                                                                                           gTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGccc                                                     >  3:259585/1‑93 (MQ=255)                                                                                                                                                |                                                                                                                                             CATCGTGCGGGTTGTGCGAGGCACTGATGACGATGCCGGCTTCGGCGTGGAAGGTACGCGTGAGGTAAGCGATGGCCGGTGTCGGCATCGGCCCGAGCAGCATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACA‑TGTAGCCGGAG‑ATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGC  >  NC_002947/5363447‑5363728

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