Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,973,487	G→C	8.7%	G129G (GGC→GGG) ‡	PP_2600 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,973,487	0	G	C	8.7%	137.0 / 3.7	46	G129G (GGC→GGG) ‡	PP_2600	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/23);  new base C (2/2);  total (21/25)
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.02e-01

AAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCAAGGAAGGCATCGACGCGCTCGGCATTGTAGTTTGCCCCCAGCTCATGG  >  NC_002947/2973344‑2973634                                                                                                                                                    |                                                                                                                                                    aaGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTgg                                                                                                                                                <  5:136428/149‑1 (MQ=255) aaGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTgg                                                                                                                                                >  4:90194/1‑149 (MQ=255) aaGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTgg                                                                                                                                                <  4:23615/149‑1 (MQ=255)  aGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGc                                                                                                                                               <  4:273613/149‑1 (MQ=255)    gcggcAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGt                                                                                                                                             >  1:260729/1‑149 (MQ=255)     cggcAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAGGCGCAGCCCTTGGCg                                                                                                                                              >  8:195193/1‑147 (MQ=255)     cggcAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAGGCGCAGCCCTTGGCg                                                                                                                                              <  7:195193/147‑1 (MQ=255)           ggTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCca                                                                                                                                      >  6:390305/1‑149 (MQ=255)               ggCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCcatgc                                                                                                                                   <  5:390305/148‑1 (MQ=255)                gCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCcatgcat                                                                                                                                 >  5:182551/1‑149 (MQ=255)                gCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCcatgcat                                                                                                                                 <  2:107254/149‑1 (MQ=255)                 cACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGCTCTAGCGCAAGCGTGGCAGCCGGCCGCAAACCTCGCGGGCATCGTAGCG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTCCGCCTCGCGTATCGCcatgcatg                                                                                                                                >  5:269752/1‑149 (MQ=255)                    cGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGAGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTcc                                                                                                                             >  4:147072/1‑149 (MQ=255)                     ggCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCg                                                                                                                            <  3:90194/149‑1 (MQ=255)                      gCCATGATCGGCATGCCCATGAACGAGGGTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCgg                                                                                                                           <  6:294385/149‑1 (MQ=255)                                      ccATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGAtgccgttg                                                                                                                 >  5:369580/1‑143 (MQ=255)                                      ccATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGAtgccgttg                                                                                                                 <  6:369580/143‑1 (MQ=255)                                           aaCGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCaaa                                                                                                      >  2:55148/1‑149 (MQ=255)                                              gAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTg                                                                                                   >  4:227859/1‑149 (MQ=255)                                                      gcgcATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTcc                                                                                                                             <  5:360978/115‑1 (MQ=255)                                                      gcgcATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTcc                                                                                                                             >  6:360978/1‑115 (MQ=255)                                                      gcgcATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCg                                                                                           >  8:115643/1‑149 (MQ=255)                                                      gcgcATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCg                                                                                           <  2:260729/149‑1 (MQ=255)                                                          aTGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTa                                                                                       <  1:55148/149‑1 (MQ=255)                                                                  gcgcAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGTCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTtctc                                                                                 <  5:255352/147‑1 (MQ=255)                                                                  gcgcAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTtctc                                                                                 <  7:115643/147‑1 (MQ=255)                                                                    gcAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGaa                                                                              >  4:225897/1‑148 (MQ=255)                                                                      aaGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGaaaaa                                                                           >  3:349355/1‑149 (MQ=255)                                                                      aaGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGaaaaa                                                                           <  3:227859/149‑1 (MQ=255)                                                                      aaGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGaaaaa                                                                           <  3:225897/149‑1 (MQ=255)                                                                      aaGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGaaaaa                                                                           <  6:182551/149‑1 (MQ=255)                                                                        gCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGTGCCTCCTCAGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGa                                                                               <  1:126790/147‑1 (MQ=255)                                                                                  agGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATg                                                               >  1:330315/1‑149 (MQ=255)                                                                                              tCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCggg                                                           <  1:141943/141‑1 (MQ=255)                                                                                              tCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCggg                                                           >  2:141943/1‑141 (MQ=255)                                                                                                    gCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCaagg                                               >  8:189564/1‑147 (MQ=255)                                                                                                     cATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCaaggaag                                            <  3:147072/149‑1 (MQ=255)                                                                                                      aTCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCaaggaagg                                           >  5:348749/1‑149 (MQ=255)                                                                                                         gTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCAAGGAAGGCAt                                        <  4:349355/149‑1 (MQ=255)                                                                                                                    gcggcgATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCAAGGAAGGCATCGAcgcg                                 <  2:330315/149‑1 (MQ=255)                                                                                                                          aTCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCAGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGccac                                                                     <  6:184406/107‑1 (MQ=255)                                                                                                                          aTCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCAGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGccac                                                                     >  5:184406/1‑107 (MQ=255)                                                                                                                                   gatacgccAGGCGCAGCCCTTGGGGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCAAGGAAGGCATCGACGCGCTCGGCATTGTAGtt                  >  5:77742/5‑149 (MQ=255)                                                                                                                                    atacgccAGGCGCAGCCCTTGGTGTATCGCCATGCATGTACGGTATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCTAAAAACTCCACCATATGCGGGCACTCCTCAAGGAAGGCATCGACGCTCTCGGCATTGTAGttt                 <  6:77742/146‑1 (MQ=255)                                                                                                                                       cgccAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCAAGGAAGGCATCGACGCGCTCGGCATTGTAGTTTGcc              >  7:168147/1‑149 (MQ=255)                                                                                                                                          cAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCAAGGAAGGCATCGACGCGCTCGGCATTGTAGTTTGccccc           <  2:133190/149‑1 (MQ=255)                                                                                                                                                   cgcCTGGCGTATCGCCTTGCATGTCCGGGAGGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCAAGGAAGGCATCGACGCGCTCGGCATTGTAGTTTGCCCCCAGCTCATg   <  6:269752/148‑1 (MQ=255)                                                                                                                                                   cgcCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCAAGGAAGGCATCGACGCGCTCGGCATTGTAGTTTGCCCCCAGCTCATgg  <  6:348749/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    AAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCAAGGAAGGCATCGACGCGCTCGGCATTGTAGTTTGCCCCCAGCTCATGG  >  NC_002947/2973344‑2973634

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