Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,154,403	G→A	8.1%	S123S (AGC→AGT)	PP_1909 ←	septum formation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,154,403	0	G	A	8.1%	137.1 / 2.7	49	S123S (AGC→AGT)	PP_1909	septum formation protein
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base A (2/2);  total (24/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 = 6.96e-01

ACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGAGCCTATGATCAA  >  NC_002947/2154255‑2154546                                                                                                                                                     |                                                                                                                                                  aCAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGgcg                                                                                                                                                   <  5:104254/149‑1 (MQ=255) aCAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGgcg                                                                                                                                                   <  5:21307/149‑1 (MQ=255)   aaGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCAAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCt                                                                                                                                                 >  7:277540/1‑149 (MQ=255)    ctcTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACTGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTg                                                                                                                                                <  7:108499/147‑1 (MQ=255)       tCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCCGTTCACGCAGTGTCACGGTAAAAGGCAAGCAATCCAGCAGGCATTGCCCGGTGGCGCTg                                                                                                                                                >  8:108499/1‑146 (MQ=255)       tCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGt‑t                                                                                                                                             >  1:106426/1‑148 (MQ=255)       tCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCACTgct                                                                                                                                              >  3:270177/1‑146 (MQ=255)         aCCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑c                                                                                                                                           <  6:299168/147‑1 (MQ=255)           ccccAGCCCCTCCGCCTTGAAACTCCCTGCTCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑cagca                                                                                                                                       <  5:106507/149‑1 (MQ=255)           ccccAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGCTGGCACTgctgaa                                                                                                                                           >  4:144115/1‑142 (MQ=255)              cAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTg                                                                                                                                    <  6:246169/149‑1 (MQ=255)                 ccccTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCACTGC‑T‑GAACAGCGCCa                                                                                                                                 <  4:270177/149‑1 (MQ=255)                   ccTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGCTGGCACTGC‑T‑GAACAGCGCCacc                                                                                                                               <  3:144115/149‑3 (MQ=255)                      ccGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGTGCTGT‑T‑CAGCAGTGCCAGCCCg                                                                                                                            <  7:55723/149‑1 (MQ=255)                      ccGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCg                                                                                                                            <  7:41765/149‑1 (MQ=255)                           ttGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGt                                                                                                                          >  3:126517/1‑146 (MQ=255)                           ttGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGt                                                                                                                          >  6:21962/1‑146 (MQ=255)                            tGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGa                                                                                                                      <  1:134980/149‑1 (MQ=255)                              aaaCTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGa                                                                                                                      >  7:75865/1‑147 (MQ=255)                                   cccTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCa                                                                                                               >  3:196953/1‑149 (MQ=255)                                        cgcAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACctg                                                                                                           <  4:38260/148‑1 (MQ=255)                                         gcAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACctgct                                                                                                         <  8:75865/149‑1 (MQ=255)                                                aGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCt                                                                                                  <  5:161904/149‑1 (MQ=255)                                                        ttCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTc                                                                                            <  8:277540/147‑1 (MQ=255)                                                              tgcCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTc                                                                                            >  8:129810/1‑141 (MQ=255)                                                              tgcCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTc                                                                                            <  7:129810/141‑1 (MQ=255)                                                                ccACGTAACGCTCGACCTGCTCCCTGCTCAGATCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTg                                                                                   <  4:100218/148‑1 (MQ=255)                                                                  aCGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTc                                                                                <  1:110840/149‑1 (MQ=255)                                                                        cGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAgc                                                                          <  5:64592/149‑1 (MQ=255)                                                                                  gctccctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAAc                                                                >  1:146139/1‑149 (MQ=255)                                                                                  gcgccctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAAc                                                                >  4:60386/4‑149 (MQ=255)                                                                                      cctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTAt                                                            >  4:116563/1‑149 (MQ=255)                                                                                       ctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATg                                                           >  2:30643/1‑149 (MQ=255)                                                                                          ctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCgg                                                        <  5:21962/149‑1 (MQ=255)                                                                                           tcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTATGCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCg                                                         >  1:303024/1‑149 (MQ=255)                                                                                              gTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGaa                                                                                                                     <  4:91438/84‑1 (MQ=255)                                                                                              gTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGaa                                                                                                                     >  3:91438/1‑84 (MQ=255)                                                                                              gTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTg                                                                                   >  8:166334/1‑118 (MQ=255)                                                                                              gTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTg                                                                                   <  7:166334/118‑1 (MQ=255)                                                                                                  aCGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGccc                                                 <  5:88956/148‑1 (MQ=255)                                                                                                  aCGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGcc                                                  >  5:292841/1‑147 (MQ=255)                                                                                                               gTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGcc                                   <  2:106426/149‑1 (MQ=255)                                                                                                                tAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGccc                                  <  2:303024/149‑1 (MQ=255)                                                                                                                                cACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATCCCGCAACCTGGt                  >  2:239656/1‑149 (MQ=255)                                                                                                                                cACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGt                  >  2:234440/1‑149 (MQ=255)                                                                                                                                 aCCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTc                 >  5:280045/1‑149 (MQ=255)                                                                                                                                        aTTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTg                                                    >  7:283185/1‑107 (MQ=255)                                                                                                                                        aTTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTg                                                    <  8:283185/107‑1 (MQ=255)                                                                                                                                         ttGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGAGCCTa         >  8:151950/1‑149 (MQ=255)                                                                                                                                            cccGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGAGCCTATGa      >  8:277488/1‑149 (MQ=255)                                                                                                                                                gTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGAGCCTATGATCaa  <  7:151950/149‑1 (MQ=255)                                                                                                                                                    cgcTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGa              >  5:116696/1‑133 (MQ=255)                                                                                                                                                    cgcTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGa              <  6:116696/133‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  ACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGT‑T‑CAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGAGCCTATGATCAA  >  NC_002947/2154255‑2154546

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