Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,648,209	G→C	7.3%	intergenic (+16/+8)	elaA → / ← PP_4957	acetyltransferase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,648,209	0	G	C	7.3%	174.5 / 3.9	55	intergenic (+16/+8)	elaA/PP_4957	acetyltransferase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/26);  new base C (2/2);  total (27/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.75e-01

GGCCGGGGGTGCCGGTTTATCTGTCGGCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTCGAACTCGCAGACGATGTC  >  NC_002947/5648063‑5648357                                                                                                                                                   |                                                                                                                                                      ggCCGGGGGTGCCGGTTTATCTGTCGGCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGc                                                                                                                                                      <  1:59182/148‑1 (MQ=255) ggCCGGGGGTGCCGGTGTATCTGTCGGCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGc                                                                                                                                                      <  2:53930/148‑1 (MQ=255)    cGGGGGTGCCGGTTTATCTGTCGGCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGca                                                                                                                                                  <  8:321780/149‑1 (MQ=255)        ggTGCCGGTTTATCTGTCGGCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGCCCGcacagg                                                                                                                                              <  7:204938/149‑3 (MQ=255)          tGCCGGTTTATCTGTCGGCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAg                                                                                                                                            <  5:87598/149‑1 (MQ=255)            ccGGTTTATCTGTCGGCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAggg                                                                                                                                          <  7:237624/149‑1 (MQ=255)                        tCGGCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCa                                                                                                                                  <  8:118850/145‑1 (MQ=255)                             tCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCACACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTcgc                                                                                                                         >  2:194026/1‑149 (MQ=255)                               cGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTcgcg                                                                                                                        <  8:7260/148‑1 (MQ=255)                                 gcgcATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGAt                                                                                                                     >  2:320177/1‑149 (MQ=255)                                          cAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCaa                                                                                                           >  3:193860/1‑149 (MQ=255)                                           aGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAg                                                                                                          >  6:94294/1‑149 (MQ=255)                                           aGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAg                                                                                                          >  2:96279/1‑149 (MQ=255)                                           aGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGAGGATCTGCC‑GCAAg                                                                                                          >  3:101044/1‑149 (MQ=255)                                           aGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAAGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAg                                                                                                          >  7:8303/1‑149 (MQ=255)                                                     tGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGgatga                                                                                                >  5:376233/1‑149 (MQ=255)                                                     tGCCCGCCATGGGCTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGgatga                                                                                                <  7:192790/149‑1 (MQ=255)                                                       cccGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGgatgat                                                                                               <  4:193860/148‑1 (MQ=255)                                                       cccGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATc                                                                                              <  1:320177/149‑1 (MQ=255)                                                        ccGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATc                                                                                              >  1:84412/1‑148 (MQ=255)                                                         cGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCa                                                                                            >  7:29533/1‑149 (MQ=255)                                                          gCCATGGGTTTGCGGTGGTGGGCGAGGAGTCTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTcgc                                                                                                                         <  2:90961/120‑1 (MQ=255)                                                          gCCATGGGTTTGCGGTGGTGGGCGAGGAGTCTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTcgc                                                                                                                         >  1:90961/1‑120 (MQ=255)                                                             aTGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGc                                                                                        <  6:376233/149‑1 (MQ=255)                                                               gggTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCtt                                                                                      >  8:96391/1‑149 (MQ=255)                                                                 gTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGTGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGt                                                                                    <  5:94294/149‑1 (MQ=255)                                                                         tggtggGCGAGGAGTTTTTGGAAGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGcc                                                                            <  8:8303/149‑1 (MQ=255)                                                                                 gaggagTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGt                                                                     >  2:160829/1‑148 (MQ=255)                                                                                      gTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATc                                                                                              >  4:105111/1‑118 (MQ=255)                                                                                      gTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATc                                                                                              <  3:105111/118‑1 (MQ=255)                                                                                       tttttGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCGGTCATGAAt                                                              >  8:52947/1‑149 (MQ=255)                                                                                             gaggaCGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCTGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGa                                                          <  1:96279/148‑1 (MQ=255)                                                                                              aggaCGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATaa                                                       <  8:29533/149‑1 (MQ=255)                                                                                              aggaCGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATa                                                        >  2:157809/1‑148 (MQ=255)                                                                                                   ggCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCAc                                                                                           <  5:275779/108‑1 (MQ=255)                                                                                                   ggCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCAc                                                                                           >  6:275779/1‑108 (MQ=255)                                                                                                     cATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCgg                                                >  4:242915/1‑149 (MQ=255)                                                                                                                   gAATGCGCAAGGGATGAGCCCACCCTGTGCGCGACGCACATCAGGGGTAACCGAGAACTTAGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCATGTCGATCGCCCACCAGTCATGAATGCGATAATGCCCGGCGTGGTtgcgtgcg                                  >  8:196324/1‑149 (MQ=255)                                                                                                                    aaTGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTgcgc                                 <  7:96391/149‑1 (MQ=255)                                                                                                                        cgcAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCGCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCtt                              <  8:318253/148‑1 (MQ=255)                                                                                                                         gcAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTcc                            >  3:365702/1‑149 (MQ=255)                                                                                                                           aaGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCctg                          <  2:84412/149‑1 (MQ=255)                                                                                                                               gATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGAGGTTGCTTGCGCCTTCCTGCTGt                      >  6:121716/1‑149 (MQ=255)                                                                                                                                 tGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTc                    >  7:82738/1‑149 (MQ=255)                                                                                                                                 tGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTc                    <  4:101044/149‑1 (MQ=255)                                                                                                                                  gAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCGGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTc                    <  7:52947/148‑1 (MQ=255)                                                                                                                                    gCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGccccc                                                                         <  4:134463/93‑1 (MQ=255)                                                                                                                                    gCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGccccc                                                                         >  3:134463/1‑93 (MQ=255)                                                                                                                                    gCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTCGaa                 <  8:82738/149‑1 (MQ=255)                                                                                                                                       cACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTAGAACTc              >  8:294286/1‑149 (MQ=255)                                                                                                                                        ctgatgtgCGGCCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTCGAACTCg             >  5:277884/5‑149 (MQ=255)                                                                                                                                        ctgatgtgCGGCCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTCGAACTCg             >  7:161271/5‑149 (MQ=255)                                                                                                                                        aCCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTCGAACTCg             <  1:160829/149‑1 (MQ=255)                                                                                                                                          ccTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTCGAACTCg             >  4:61157/1‑147 (MQ=255)                                                                                                                                                cGGCCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTCGAACTCGCAGACGAt     <  8:161271/149‑1 (MQ=255)                                                                                                                                                 gggCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTCGAACTCGCAGACGATg    <  1:157809/149‑1 (MQ=255)                                                                                                                                                   gCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCCCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTCGAACTAGCAGACGATGTc  >  3:212196/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GGCCGGGGGTGCCGGTTTATCTGTCGGCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCC‑GCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTCGAACTCGCAGACGATGTC  >  NC_002947/5648063‑5648357

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