Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,622,934	G→T	8.7%	A268A (GCC→GCA)	PP_4941 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,622,934	0	G	T	8.7%	133.3 / 4.1	46	A268A (GCC→GCA)	PP_4941	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/20);  new base T (2/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01

GCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTGGTA  >  NC_002947/5622788‑5623080                                                                                                                                                   |                                                                                                                                                   gCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggc                                                                                                                                                  <  4:61265/149‑1 (MQ=255)     gCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGTGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcaggc                                                                                                                                              >  2:378700/1‑149 (MQ=255)          tCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCCATGCAGGCAGcac                                                                                                                                         <  8:92152/149‑3 (MQ=255)          tCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCCATGCAGGCAGcac                                                                                                                                         <  7:141316/149‑3 (MQ=255)               cagcagGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACCCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTCCAgg                                                                                                                                     >  7:108590/1‑148 (MQ=255)                     gTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAg                                                                                                                               >  5:233322/1‑148 (MQ=255)                        gtgCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAg                                                                                                                            <  5:86637/148‑1 (MQ=255)                        gtgCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCCTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAg                                                                                                                            <  3:72238/148‑1 (MQ=255)                             aGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGc                                                                                                                      >  6:191412/1‑149 (MQ=255)                                  ggtgCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTAt                                                                                                                 <  5:327490/149‑1 (MQ=255)                                    tgCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATg                                                                                                                >  8:335639/1‑148 (MQ=255)                                     gcgcCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGc                                                                                                               >  1:140282/1‑148 (MQ=255)                                               ggTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTcc                                                                                                    <  7:335639/149‑1 (MQ=255)                                                     ccGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACg                                                                                              >  8:164535/1‑149 (MQ=255)                                                     ccGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACg                                                                                              >  8:310714/1‑149 (MQ=255)                                                     ccGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACg                                                                                              >  5:149738/1‑149 (MQ=255)                                                                 acgggTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAAcgc                                                                                  <  4:276229/149‑1 (MQ=255)                                                                  cgggTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAAcgcg                                                                                 <  6:233322/149‑1 (MQ=255)                                                                        gAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCacga                                                                           >  8:252505/1‑149 (MQ=255)                                                                              gCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCGAACGCGGCACGACGATGc                                                                     >  2:293769/1‑149 (MQ=255)                                                                                  aaCACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTc                                                                 >  5:354832/1‑149 (MQ=255)                                                                                     acGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCa                                                              <  1:21288/149‑1 (MQ=255)                                                                                             atCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACg                                                        >  7:223905/1‑147 (MQ=255)                                                                                                  ccAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGgcgc                                                 <  7:252505/149‑1 (MQ=255)                                                                                                    aGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGccc                                               <  1:378700/149‑1 (MQ=255)                                                                                                         ggCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCa                                          <  7:232083/149‑1 (MQ=255)                                                                                                          gCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAg                                         <  2:140282/149‑1 (MQ=255)                                                                                                          gCTTCACGTGAGGCGGCCGGCTGGTCGGGCGATGCCAGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGGCCCAGGCACGCAGGCGCCCCTTCAg                                         <  5:112372/149‑1 (MQ=255)                                                                                                            ttCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGgc                                       >  2:362561/1‑149 (MQ=255)                                                                                                            ttCACGGGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGgc                                       <  8:108590/149‑1 (MQ=255)                                                                                                                cGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGgcgcgc                                   >  3:251831/1‑149 (MQ=255)                                                                                                                 gTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggca                                                                                                                                                 >  3:278657/1‑38 (MQ=255)                                                                                                                 gTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggca                                                                                                                                                 <  4:278657/38‑1 (MQ=255)                                                                                                                  tGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGgcgc                                                 <  8:306519/133‑1 (MQ=255)                                                                                                                  tGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGgcgc                                                 >  7:306519/1‑133 (MQ=255)                                                                                                                    aGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGTTGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGcggcc                               >  5:318738/1‑149 (MQ=255)                                                                                                                      gcggcCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGcggcctg                             >  2:119588/1‑148 (MQ=255)                                                                                                                      gcggcCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGT‑CCAGGCACGCAGGCGCCCCTTCAGGCGCGcggcctcg                            <  6:59855/149‑1 (MQ=255)                                                                                                                          ccggcTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAg                                         <  1:334191/133‑1 (MQ=255)                                                                                                                          ccggcTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAg                                         >  2:334191/1‑133 (MQ=255)                                                                                                                           cggcTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGcggcctcggcct                        >  6:201474/1‑149 (MQ=255)                                                                                                                                 gTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCa                  <  6:149738/149‑1 (MQ=255)                                                                                                                                 gTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCAGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGc                    <  7:164535/147‑1 (MQ=255)                                                                                                                                   ctgcagctgcCTGCCATGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAgg                >  3:184081/6‑149 (MQ=255)                                                                                                                                         ctgcctgcCATGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATc          >  2:19105/1‑149 (MQ=255)                                                                                                                                         ctgcctgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATc          <  4:251831/149‑1 (MQ=255)                                                                                                                                           gcctgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATctg        >  8:228790/1‑149 (MQ=255)                                                                                                                                            cctgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATctgc       >  2:4651/1‑149 (MQ=255)                                                                                                                                              tgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATctgctg     >  7:144677/1‑149 (MQ=255)                                                                                                                                                 aTGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTGGTa  <  5:191412/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTGGTA  >  NC_002947/5622788‑5623080

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