Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,603,016	C→T	8.9%	A141V (GCC→GTC)	garK →	glycerate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,603,016	0	C	T	8.9%	121.0 / 2.3	45	A141V (GCC→GTC)	garK	glycerate kinase
Reads supporting (aligned to +/- strand):  ref base C (21/19);  new base T (2/2);  total (24/21)
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.92e-01

AGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGCCCATATCGACGCCAGCGGCCTCGACCCGCGCCTGGCCGGGGTGCAGTTTGAAGTC  >  NC_002947/3602869‑3603164                                                                                                                                                     |                                                                                                                                                     aGGCCAGTGGCCTGCAAGAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGcc                                                                                                                                                     <  2:102940/149‑1 (MQ=255)        tGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATCCCGGCAGTg                                                                                                                                              >  1:285031/1‑149 (MQ=255)            cTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCAt                                                                                                                                          >  3:103645/1‑149 (MQ=255)                        ccAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACt                                                                                                                              <  7:79194/149‑1 (MQ=255)                          aaCAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCCGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTg                                                                                                                            >  8:136040/1‑149 (MQ=255)                            cAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGgc                                                                                                                          >  8:145945/1‑149 (MQ=255)                            cAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGgc                                                                                                                          <  2:285031/149‑1 (MQ=255)                                 cAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGgcctgcg                                                                                                                     >  2:149971/1‑149 (MQ=255)                                       gATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCACCGAAGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCcggct                                                                                                               >  8:14025/1‑149 (MQ=255)                                       gATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCctgct                                                                                                               <  1:149971/149‑1 (MQ=255)                                           cTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGa                                                                                                            <  6:144025/148‑1 (MQ=255)                                           cTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGAt                                                                                                           <  6:240443/149‑1 (MQ=255)                                             tGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGc                                                                                                         >  6:267950/1‑149 (MQ=255)                                                cGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCTACGACGCCGCCCCCCGGCGTGTGGTGATGGCCAACGGCGGCAGCGCCCCCAACGATGCCGGCAGGGGCTTGTTGCGAGCCCTTGCCGTGCGCCTGCTGGATGGTg                                                                                                       >  1:62284/1‑149 (MQ=255)                                                   agcagcaCCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGa                                                                                                                                   <  2:276427/117‑1 (MQ=255)                                                   agcagcaCCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGa                                                                                                                                   >  1:276427/1‑117 (MQ=255)                                                     cagcaCCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGc                                                                                                 <  4:43943/149‑1 (MQ=255)                                                         aCCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGCCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAgc                                                                                             >  4:44716/1‑149 (MQ=255)                                                               ggTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGa                                                                                       <  5:267950/149‑1 (MQ=255)                                                                 tACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTgg                                                                                        >  3:98632/1‑146 (MQ=255)                                                                    cGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGTCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAgg                                                                                  >  4:118248/1‑149 (MQ=255)                                                                    cGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGTCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAgg                                                                                  >  2:24771/1‑149 (MQ=255)                                                                      gTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAgcgc                                                                                           >  2:270215/1‑138 (MQ=255)                                                                      gTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAgcgc                                                                                           <  1:270215/138‑1 (MQ=255)                                                                      gTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGggc                                                                                >  5:295884/1‑149 (MQ=255)                                                                                  ccGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTg                                                                    >  7:217207/1‑149 (MQ=255)                                                                                   cGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTgg                                                                   <  4:103645/149‑1 (MQ=255)                                                                                         tct‑cGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGc                                                             >  5:198656/1‑149 (MQ=255)                                                                                           t‑cGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCCGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTg                                                           >  5:37447/1‑149 (MQ=255)                                                                                               aCGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGcc                                                        <  6:295884/149‑1 (MQ=255)                                                                                                 gCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGGCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGccc                                                       >  1:234860/1‑148 (MQ=255)                                                                                                 gCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGccc                                                       <  2:234860/148‑1 (MQ=255)                                                                                                    ggCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGCCCatat                                                   >  8:219958/1‑149 (MQ=255)                                                                                                     gcgcCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGCCCATATc                                                  >  1:303098/1‑149 (MQ=255)                                                                                                         ccGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGCCCATATCGACg                                              <  5:47017/149‑1 (MQ=255)                                                                                                          cGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGCCCATATCGACGc                                             <  8:217207/149‑1 (MQ=255)                                                                                                                           cATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGCCCATATCGACGCCAGCGGCCTCGACCcgc                            >  8:292542/1‑149 (MQ=255)                                                                                                                                gcggcAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGCCCATATCGACGCCAGCGGCCTCGACCCGCGCCTg                       <  2:303098/149‑1 (MQ=255)                                                                                                                                   gcAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGCCCATATCGACGCCAGCGGCCTCGACCCGCGCCTGGcc                    >  6:263873/1‑149 (MQ=255)                                                                                                                                    cAGCGCCACCAACGATGTCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGCCCATATCGACGCCAGCGGCCTCGACCCGCGCCTGGCCg                   <  3:118248/149‑1 (MQ=255)                                                                                                                                    cAGCGCCACCAACGATGTCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGCCCATATCCACGCCAGCGGCCTCGACCCGCGCCTGGCCg                   <  1:24771/149‑1 (MQ=255)                                                                                                                                      gcgcCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGCCCATATCGACGCCAGCGGCCTCGACCCGCGCCTGGCCggg                 <  3:44716/149‑1 (MQ=255)                                                                                                                                      gcgcCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGCCCATATCGACGCCAGCGGCCTCGACCCGCGCCTGGCCggg                 <  7:145945/149‑1 (MQ=255)                                                                                                                                         ccaccaACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGCCCATATCGACGCCAGCGGTCTCGACCCGCGCCTGGCCGGGGTg              <  7:136040/149‑1 (MQ=255)                                                                                                                                           accaACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGCCCATATCCACGCCAGCGGCCTCGACCCGCGCCTGGCCGGGGTGCa            >  3:65842/1‑149 (MQ=255)                                                                                                                                                  aTGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGCCCATATCGACGCCAGCGGCCTCGACCCGCGCCTGGCCGGGGTGCAGTTTGaa     <  2:62284/149‑1 (MQ=255)                                                                                                                                                     ccGTCAGTGGAATGGTGCGCGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCGTGCATTGGCCAAGCTGGCCCATATCGACGCCAGCGGCCTCGACCCGCGCATGTCCGGGGTGCAGTTTGAAGTc  <  7:14025/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AGGCCAGTGGCCTGCAACAGGTGCCAACAGGCCAACGCGATGCTTGCCGCAGCAGCACCTGGGGTACCGGTGAACTGATTGCCGCAGCTCT‑CGACGCCGGCGCCCGGCGTGTGGTGCTGGCCATCGGCGGCAGCGCCACCAACGATGCCGGCAGTGGCATGTTGCGAGCACTTGGCCTGCGCCTGCTGGATGCTGACGGGCAAGCGCTGGAAGAGGGCGGCCTGGCCTTGGCCAAGCTGGCCCATATCGACGCCAGCGGCCTCGACCCGCGCCTGGCCGGGGTGCAGTTTGAAGTC  >  NC_002947/3602869‑3603164

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