Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,556,558	C→T	7.7%	G307S (GGT→AGT)	gntZ ←	6‑phosphogluconate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,556,558	0	C	T	7.7%	148.6 / 2.3	52	G307S (GGT→AGT)	gntZ	6‑phosphogluconate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (26/22);  new base T (2/2);  total (28/24)
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.97e-01

ACCAGGTCGCCATTGGCGCCGAACAGGAACAGGGTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCTTCGCCACTGTCGGACACCGAGCCACTGAATTGCGACAGCT  >  NC_002947/4556426‑4556704                                                                                                                                      |                                                                                                                                                   aCCAGGTCGCCATTGGCGCCGAACAGGAACAGGGTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACTGTAgg                                                                                                                                               >  7:78956/1‑138 (MQ=255) aCCAGGTCGCCATTGGCGCCGAACAGGAACAGGGTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACTGTAgg                                                                                                                                               >  7:84571/1‑138 (MQ=255) aCCAGGTCGCCATTGGCGCCGAACAGGAACAGGGTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACTGTAgg                                                                                                                                               <  8:84571/138‑1 (MQ=255) aCCAGGTCGCCATTGGCGCCGAACAGGAACAGGGTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACTGTAgg                                                                                                                                               <  8:78956/138‑1 (MQ=255)  ccAGGTCGCCATTGGCGCCGAACAGGAACAGGGTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGc                                                                                                                                    <  2:366315/148‑1 (MQ=255)  ccAGGTCGCCATTGGCGCCGAACAGGAACAGGGTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCt                                                                                                                                   >  5:381768/1‑149 (MQ=255)   cAGGTCGCCATTGGCGCCGAACAGGAACAGGGTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTg                                                                                                                                  >  5:371504/1‑149 (MQ=255)             ttGGCGCCGAACAGGAACAGGGTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACgg                                                                                                                        <  7:26067/149‑1 (MQ=255)                gcgcCGAACAGGAACAGGGTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAg                                                                                                                     >  5:239095/1‑149 (MQ=255)                       acaggaacaggGTACAAGGGGGGCCGGCCGCAAGGGTCAGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGcccgcg                                                                                                               >  8:142453/1‑148 (MQ=255)                        caggaacaggGTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGaa                                                                                                             <  6:381768/149‑1 (MQ=255)                         aggaacaggGTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGgc                                                                                                                                >  2:299422/1‑129 (MQ=255)                         aggaacaggGTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGgc                                                                                                                                <  1:299422/129‑1 (MQ=255)                           gaacaggGTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAg                                                                                                          >  3:137418/1‑149 (MQ=255)                           gaacaggGTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAg                                                                                                          >  5:277164/1‑149 (MQ=255)                                gggTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAg                                                                                                                                                <  7:338175/106‑1 (MQ=255)                                gggTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAg                                                                                                                                                >  8:338175/1‑106 (MQ=255)                                gggTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTgcgc                                                                                                     >  6:345764/1‑149 (MQ=255)                                  gTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTg                                                                                                   <  6:10364/149‑1 (MQ=255)                                       aGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCCCCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGacaa                                                                                              >  3:298264/1‑149 (MQ=255)                                        gggggggCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGacaac                                                                                             <  6:371504/149‑1 (MQ=255)                                                   cGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAAcaccggccccgg                                                                                  >  4:286116/1‑149 (MQ=255)                                                    gCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAAcaccggcaccgg                                                                                  >  2:2540/1‑148 (MQ=255)                                                     cAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACAccggcaccggca                                                                                >  1:134351/1‑149 (MQ=255)                                                         tgTCTGTTTAGTCATTCGCCTTTCTTCTAGACGTGGACACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGc                                                                            <  4:298264/148‑1 (MQ=255)                                                                  aGTCATTCGCCTTTGTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAgc                                                                                                                    >  7:352198/1‑100 (MQ=255)                                                                  aGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAgc                                                                                                                    <  8:352198/100‑1 (MQ=255)                                                                          gCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGt                                                           >  7:382518/1‑149 (MQ=255)                                                                          gCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGt                                                           >  7:138019/1‑149 (MQ=255)                                                                             tttCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGa                                                        >  1:39317/1‑149 (MQ=255)                                                                               tcttctCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGa                                                        <  5:345764/147‑1 (MQ=255)                                                                                 ttctCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGt                                                    >  1:316477/1‑149 (MQ=255)                                                                                   ctcGACGTGGCCACCGAAGCCCAGGCGCATGGCCGAAGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAacac                                                                                           >  1:67941/1‑109 (MQ=255)                                                                                   ctcGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAacac                                                                                           <  2:67941/108‑1 (MQ=255)                                                                                       aCGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCg                                              <  1:2540/149‑1 (MQ=255)                                                                                       aCGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCg                                              >  3:287216/1‑149 (MQ=255)                                                                                        cGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCgg                                             <  6:239095/149‑1 (MQ=255)                                                                                           ggCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCt                                          >  2:1822/1‑149 (MQ=255)                                                                                              cACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCTTCg                                       >  6:402336/1‑149 (MQ=255)                                                                                                     gCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCTTCGCCACTGt                                <  8:382518/149‑1 (MQ=255)                                                                                                       ccAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCTTCGCCACTGTCg                              <  6:2885/149‑1 (MQ=255)                                                                                                             gcATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCTTCGCCACTGTCGGACAcc                        >  7:18229/1‑149 (MQ=255)                                                                                                                   ccG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACTGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCTTCGCCACTGTCGGACACCGAGCCa                  <  8:138019/149‑1 (MQ=255)                                                                                                                        gagGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCTTCGCCACTGTCGGACACCGAGCCACTGa              <  6:277164/149‑1 (MQ=255)                                                                                                                          ggATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGAcgg                                                                 >  7:65319/1‑96 (MQ=255)                                                                                                                          ggATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGAcgg                                                                 <  8:65319/96‑1 (MQ=255)                                                                                                                           gATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCTTCGCCACTGTCGGACACCGAGCCACTGAAtt           >  3:337239/1‑149 (MQ=255)                                                                                                                           gATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCTTCGCCACTGTCGGACACCGAGCCACTGAAtt           <  2:316477/149‑1 (MQ=255)                                                                                                                                 gTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCTTCGCCACTGTCGGACACCGAGCCACTGAATTGCGa       >  8:230396/1‑147 (MQ=255)                                                                                                                                 gTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCTTCGCCACTGTCGGACACCGAGCCACTGAATTGCGa       <  4:137418/147‑1 (MQ=255)                                                                                                                                 gTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCTTCGCCACTGTCGGACACCGAGCCACTGAATTGCGa       <  7:230396/147‑1 (MQ=255)                                                                                                                                  tCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCTTCGCCACTGTCGGACACCGAGCCACTGAATTGCGACAg    <  2:39317/149‑1 (MQ=255)                                                                                                                                     ccGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCTTCGCCa                                    <  8:190017/114‑1 (MQ=255)                                                                                                                                     ccGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCTTCGCCa                                    >  7:190017/1‑114 (MQ=255)                                                                                                                                      cGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCTTCGCCACTGTCGGACACCGAGCCACTGAATTGCGACAGCt  >  8:362495/1‑147 (MQ=255)                                                                                                                                      |                                                                                                                                                   ACCAGGTCGCCATTGGCGCCGAACAGGAACAGGGTACAAGGGGGGGCGGCCGCAAGGGTCTGTTTAGTCATTCGCCTTTCTTCTCGACGTGGCCACCGAAGCCCAGGCGCATGGCCG‑AGAGGATCTTGTCACCGTAGGTGCCCTGTTGCTGGCGCGAACGGAAGCGCGCGAACAGTGCGCTGGACAACACCGGCACCGGCACCGCCTGCTCGACGGCGGCGTCGATGGTCCAGCGGCCTTCGCCACTGTCGGACACCGAGCCACTGAATTGCGACAGCT  >  NC_002947/4556426‑4556704

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