Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,104,903	C→A	7.3%	intergenic (‑20/+21)	PP_1879 ← / ← PP_1880	hypothetical protein/outer membrane autotransporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,104,903	0	C	A	7.3%	168.3 / 3.9	55	intergenic (‑20/+21)	PP_1879/PP_1880	hypothetical protein/outer membrane autotransporter
Reads supporting (aligned to +/- strand):  ref base C (28/23);  new base A (2/2);  total (30/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 = 7.08e-01

CAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGACCGCCACCCCTGCACCGA  >  NC_002947/2104757‑2105048                                                                                                                                                   |                                                                                                                                                   cAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCa                                                                                                                                                   <  1:86704/148‑1 (MQ=255) cAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCAt                                                                                                                                                  >  2:155054/1‑149 (MQ=255)   aCATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTg                                                                                                                                                >  1:234765/1‑149 (MQ=255)   aCATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTg                                                                                                                                                <  2:234765/149‑1 (MQ=255)      tCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGcc                                                                                                                                             <  2:81054/149‑1 (MQ=255)            accacGCTTGGTACGGTTGCGGGGACAAGAGAGCCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCt                                                                                                                                       >  6:262362/1‑149 (MQ=255)               acGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGt                                                                                                                                    >  4:153250/1‑149 (MQ=255)                  cTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTc                                                                                                                                 <  8:188844/149‑1 (MQ=255)                   ttGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTc                                                                                                                                          <  3:47507/139‑1 (MQ=255)                   ttGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTc                                                                                                                                          >  4:47507/1‑139 (MQ=255)                     ggTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCt                                                                                                                                <  5:343660/147‑1 (MQ=255)                     ggTACGGTTGGGGGGACAAGAGAGGCCGAGACCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCtt                                                                                                                               >  1:30726/1‑148 (MQ=255)                      gTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTa                                                                                                                             <  3:152885/149‑1 (MQ=255)                                        gagagGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGcc                                                                                                           >  8:300240/1‑149 (MQ=255)                                        gagagGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGGTTCTTTACCAGTTGTAACGGGCGcc                                                                                                           >  8:271801/1‑149 (MQ=255)                                         agagGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCAAGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGcc                                                                                                           >  4:97241/1‑148 (MQ=255)                                             gCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACg                                                                                                                 <  6:210371/138‑1 (MQ=255)                                             gCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACg                                                                                                                 >  5:210371/1‑138 (MQ=255)                                               cGAGTCCAGTTCTGGGGTTGGAGTTCTGCTGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTc                                                                                                   <  3:97241/149‑1 (MQ=255)                                                 aGTCCAGTTCAGGGGTTGGAGTTGTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGTAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCAc                                                                                                 <  7:198911/149‑1 (MQ=255)                                                 aGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCAc                                                                                                 <  3:153250/149‑1 (MQ=255)                                                     cAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGccc                                                                                             <  7:185893/149‑1 (MQ=255)                                                               gTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGtt                                                                                   >  7:216911/1‑149 (MQ=255)                                                                  ggAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAACGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCa                                                                                 <  2:80454/148‑1 (MQ=255)                                                                            ggTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGcc                                                                      <  1:155054/149‑1 (MQ=255)                                                                                 ccGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGc                                                                 <  8:216911/149‑1 (MQ=255)                                                                                 ccGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTAACAGTTGAAACGGGCGCCGAC‑GCTCACGCCCCAAGGCTGATCAAGATTGCTGCCATgg                                                                  >  4:9190/1‑148 (MQ=255)                                                                                    gAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTg                                                               >  2:164552/1‑148 (MQ=255)                                                                                            gTTCTGTGTCATCGTCCATCTCGATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTTACGGGCGCCGAC‑GCTGACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGaa                                                       <  5:262362/148‑1 (MQ=255)                                                                                            gTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGaa                                                       <  2:30726/148‑1 (MQ=255)                                                                                              tCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCTATGGCTGTAATCGAAGTc                                                    >  3:170609/1‑149 (MQ=255)                                                                                                  tgtCATCGTCCTTCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGTCTGTAATCGAAGTCGGc                                                 <  3:9190/148‑1 (MQ=255)                                                                                                       tCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGc                                           >  1:228904/1‑149 (MQ=255)                                                                                                         gTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCt                                         >  6:34985/1‑149 (MQ=255)                                                                                                               ctcCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGcc                                                                      >  6:365769/1‑114 (MQ=255)                                                                                                               ctcCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGcc                                                                      <  5:365769/114‑1 (MQ=255)                                                                                                                 ccATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGTGCGCCGAC‑GCTCACGCCCAATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCAtttt                                 >  5:142184/1‑149 (MQ=255)                                                                                                                 ccATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCAtttt                                 >  6:69190/1‑149 (MQ=255)                                                                                                                 ccATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTAAAGCTTGCTGCCATGGCTGTAATCGAAGTCGCCGTGGGCCTGCCAtttt                                 >  1:288203/1‑149 (MQ=255)                                                                                                                      tgtCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCTATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGcc                            <  4:170609/149‑1 (MQ=255)                                                                                                                      tgtCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGc                             <  7:56107/148‑1 (MQ=255)                                                                                                                           cGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATgg                                                                  >  7:210885/1‑106 (MQ=255)                                                                                                                           cGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATgg                                                                  <  8:210885/106‑1 (MQ=255)                                                                                                                              aTGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGa                    >  3:101273/1‑149 (MQ=255)                                                                                                                                  aTGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTtg                                                                           >  4:23620/1‑90 (MQ=255)                                                                                                                                  aTGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTtg                                                                           <  3:23620/90‑1 (MQ=255)                                                                                                                                   tGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGACCGc                >  6:334760/1‑148 (MQ=255)                                                                                                                                     cAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGACAGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGACCGCCa              >  1:103179/1‑149 (MQ=255)                                                                                                                                     cAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGACCGCCAc             >  4:42242/1‑149 (MQ=255)                                                                                                                                     cAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACCGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTAGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGACCGCCAc             >  7:281540/1‑149 (MQ=255)                                                                                                                                         agggaaggCCAATGCGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGACCGcc               <  3:52198/143‑1 (MQ=255)                                                                                                                                         agggaaggCCAATGCGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGACCGcc               >  4:52198/1‑143 (MQ=255)                                                                                                                                          gggaaggCCAATGCGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGAc                   >  2:256463/1‑138 (MQ=255)                                                                                                                                          gggaaggCCAATGCGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGAc                   <  1:256463/138‑1 (MQ=255)                                                                                                                                           ggaaggCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGACCGCCACCCCTg        >  5:182826/1‑148 (MQ=255)                                                                                                                                                gcgcATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGACCGCCACCCCTGCACCGa  >  4:318855/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGAC‑GCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGACCGCCACCCCTGCACCGA  >  NC_002947/2104757‑2105048

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