Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,526,650	C→G	7.8%	G298G (GGC→GGG)	PP_2215 →	acetyl‑CoA acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,526,650	0	C	G	7.8%	163.8 / 2.5	51	G298G (GGC→GGG)	PP_2215	acetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (21/26);  new base G (2/2);  total (23/28)
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

GCGCTGCGGCGCTGGTACTGATGCGCCGCTCCGAGGCCGACAAACGTGGCCTCAAGCCATTGGCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCATGAAACACCTCGACCTGCCACACGACAAGGTCAATATCCACGGCGGCGC  >  NC_002947/2526523‑2526796                                                                                                                                 |                                                                                                                                                   gcgctgcggcgctgGTACTGATGCGCCGCTCCGAGGCCGACAAACGTGGCCTCAAGCCATTGGCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGaaa                                                                                                                               <  2:275558/149‑1 (MQ=255)                     aTGCGCCGCTCCGAGGCCGACAAACGTGGCCTCAAGCCATTGGCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGTATCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTgg                                                                                                            <  1:248929/149‑1 (MQ=255)                       gcgcCGCTCCGAGGCCGACAAACGTGGCCTCAAGCCATTGGCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCGGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGcc                                                                                                          >  1:100140/1‑148 (MQ=255)                       gcgcCGCTCCGAGGCCGACAAACGTGGCCTCAAGCCATTGGCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGcc                                                                                                          >  7:11605/1‑148 (MQ=255)                       gcgcCGCTCCGAGGCCGACAAACGTGGCCTCAAGCCATTGGCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGcc                                                                                                          <  2:100140/148‑1 (MQ=255)                         gccgcTCCGAGGCCGACAAACGTGGCCTCAAGCCATTGGCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAAc                                                                                                                <  5:206056/140‑1 (MQ=255)                         gccgcTCCGAGGCCGACAAACGTGGCCTCAAGCCATTGGCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAAc                                                                                                                >  6:206056/1‑140 (MQ=255)                          ccgcTCCGAGGCAGACAAACGTGGCCTCAAGCCATTGGCCGTCATCCACGGCCACGCCGCCATTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACAGATCAAACTCACCGCCAGGAACCTGGCCGAAg                                                                                                      >  6:213649/1‑149 (MQ=255)                                cGAGGCCGACAAACGTGGCCTCAAGCCATTGGCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGAc                                                                                                 <  6:255085/148‑1 (MQ=255)                                       gACAAACGTGGCCTCAAGCCATTGGCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTc                                                                                           >  5:62947/1‑147 (MQ=255)                                         cAAACGTGGCCTCAAGCCATTGGCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTc                                                                                           >  7:253907/1‑145 (MQ=255)                                         cAAACGTGGCCTCAAGCCATTGGCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCgaga                                                                                       <  7:22671/149‑1 (MQ=255)                                                 gCCTCAAGCGATTGGCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGGTGT‑TCCCGAGCGCCCCGATCGGCGCGATCGAGTAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGAGCTGTTCGAGATCAACGAg                                                                               <  5:213649/149‑1 (MQ=255)                                                    tCAAGCCATTGGCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGcc                                                                            >  6:143644/1‑149 (MQ=255)                                                          cATTGGCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGcc                                                                      >  3:166203/1‑149 (MQ=255)                                                            ttGGCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGt                                                                    <  6:278325/149‑1 (MQ=255)                                                               gCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGt                                                                 >  6:117954/1‑149 (MQ=255)                                                                     aTCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCt                                                           <  5:143644/149‑1 (MQ=255)                                                                      tCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTg                                                          >  1:181551/1‑149 (MQ=255)                                                                       ccACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTgg                                                         >  6:85094/1‑149 (MQ=255)                                                                        cACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTggc                                                        >  4:234248/1‑149 (MQ=255)                                                                        cACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTggc                                                        <  8:11605/149‑1 (MQ=255)                                                                         aCGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTggcg                                                       >  5:131364/1‑149 (MQ=255)                                                                         aCGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTggcg                                                       <  5:213655/149‑1 (MQ=255)                                                                              cACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCAt                                                  <  8:253907/149‑1 (MQ=255)                                                                                cgccgcCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAATCTGATGAAACGCACCGGCGTGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCATgt                                                <  3:79714/149‑2 (MQ=255)                                                                                   cgccTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCATGAAac                                             <  4:166203/149‑1 (MQ=255)                                                                                      cTTTGCCTACCCCCCCTCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCATGAAac                                             <  6:221184/146‑1 (MQ=255)                                                                                      cTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGCAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCATGAAac                                             >  5:221184/1‑146 (MQ=255)                                                                                            cGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGcc                                                    >  4:307209/1‑133 (MQ=255)                                                                                            cGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGcc                                                    <  3:307209/133‑1 (MQ=255)                                                                                             gACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGaaa                                                                                                                               <  8:280445/57‑1 (MQ=255)                                                                                             gACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGaaa                                                                                                                               >  7:280445/1‑57 (MQ=255)                                                                                               cacCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCATGAAACACCTCGACCTGc                                 <  2:181551/149‑1 (MQ=255)                                                                                                  cccGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGcc                                                    >  2:198406/1‑127 (MQ=255)                                                                                                  cccGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGcc                                                    <  1:198406/127‑1 (MQ=255)                                                                                                      gcgcTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCATGAAACACCTCGACCTGCCACACGa                          <  5:15608/149‑1 (MQ=255)                                                                                                          tGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCATGAAACACCTCGa                                      <  8:201177/133‑1 (MQ=255)                                                                                                          tGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCATGAAACACCTCGa                                      >  7:201177/1‑133 (MQ=255)                                                                                                            t‑tCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCATGAAACACCTCGACCTGCCACACGACAAg                      <  5:85094/147‑1 (MQ=255)                                                                                                              tgtcGATCGCGCCGATCGGGGCGGTCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCATGAAACACCTCGACCTGCCACACGACAAGGt                    >  2:88494/4‑148 (MQ=255)                                                                                                                ccGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCATGAAACACCTCGACCTGCCACACGACAAGGTCaa                 >  3:85774/1‑149 (MQ=255)                                                                                                                 cGATCGCGCCGATCGGGGCGGTCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGcc                                                    >  2:188275/1‑113 (MQ=255)                                                                                                                 cGATCGCGCCGATCGGGGCGGTCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGcc                                                    <  1:188275/113‑1 (MQ=255)                                                                                                                       ccccGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGt                                                                                             <  1:271567/66‑1 (MQ=255)                                                                                                                       ccccGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGt                                                                                             >  2:271567/1‑66 (MQ=255)                                                                                                                       ccccGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCa                                                   >  7:115636/1‑108 (MQ=255)                                                                                                                       ccccGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCa                                                   <  8:115636/108‑1 (MQ=255)                                                                                                                        gccGATCGGGGCGGTCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCATGAAACACCTCGACCTGCCACACGACAAGGTCAATATCCAcg         <  1:88494/148‑1 (MQ=255)                                                                                                                           gATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCATGAAACACCTCGACCTGCCACACGACAAGGTCAATATCCAcggcg      <  6:62947/149‑1 (MQ=255)                                                                                                                               ggCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCATGAAACACCTCGACCTGCCACACGACAAGGTCAATATCCACGGCGgcgc  <  6:131364/149‑1 (MQ=255)                                                                                                                                 |                                                                                                                                                   GCGCTGCGGCGCTGGTACTGATGCGCCGCTCCGAGGCCGACAAACGTGGCCTCAAGCCATTGGCCGTCATCCACGGCCACGCCGCCTTTGCCGACACCCCGGCGCTGT‑TCCCGACCGCCCCGATCGGCGCGATCGACAAACTGATGAAACGCACCGGCTGGAACCTGGCCGAAGTCGACCTGTTCGAGATCAACGAGGCCTTCGCCGTGGTCACCCTGGCGGCCATGAAACACCTCGACCTGCCACACGACAAGGTCAATATCCACGGCGGCGC  >  NC_002947/2526523‑2526796

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑