Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,233,543	T→C	11.5%	V136A (GTC→GCC)	PP_4611 →	ECF family RNA polymerase sigma‑70 factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,233,543	0	T	C	11.5%	128.3 / 8.4	52	V136A (GTC→GCC)	PP_4611	ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base T (29/17);  new base C (4/2);  total (33/19)
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.06e-01

CTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCG  >  NC_002947/5233401‑5233685                                                                                                                                                 |                                                                                                                                               cTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACgg                                                                                                                                          >  6:348010/1‑149 (MQ=255)    cAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCt                                                                                                                                       >  8:185284/1‑149 (MQ=255)       gcgcTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTggg                                                                                                                                    <  3:173842/149‑1 (MQ=255)       gcgcTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTggg                                                                                                                                    <  7:185284/149‑1 (MQ=255)        cgcTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGc                                                                                                                                   <  3:221854/149‑1 (MQ=255)         cctcgcacggcgcgacctcttccgggacctgtgcgagGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGc                                                                                                                                   >  1:364120‑M2/38‑148 (MQ=255)           tcgcacggcgcgacctcttccgggacctgtgcgagGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGcc                                                                                                                                  >  5:244031‑M2/36‑147 (MQ=255)             GCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGgccagca                                                                                                                              >  2:12039‑M1/1‑146 (MQ=255)                  ggTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGCCagca                                                                                                                              >  3:151845/1‑141 (MQ=255)                  ggTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGCCagca                                                                                                                              <  4:151845/144‑4 (MQ=255)                    tCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtc                                                                                                                       >  3:115201/1‑149 (MQ=255)                    gcgacctcttccgggacctgtgcgagGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtc                                                                                                                       <  2:364120‑M2/123‑1 (MQ=255)                        ggAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcgctc                                                                                                                   >  6:119135/1‑149 (MQ=255)                        ggAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcgctc                                                                                                                   >  4:267536/1‑149 (MQ=255)                        cctcttccgggacctgtgcgagGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcgctc                                                                                                                   <  6:244031‑M2/127‑1 (MQ=255)                                 cgcgCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTgg                                                                                                          <  5:119135/149‑1 (MQ=255)                                    gcCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTGGGAGCCCTTGGTAAAAACCGCCATAAATCTGTCCCACCTCGACCGCCAGGGGCAAAGCGAAATAGCTAGCGAACTgtg                                                                                                         >  1:294070/1‑147 (MQ=255)                                           gagGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCt                                                                                                >  3:52603‑M2/4‑149 (MQ=255)                                                  aGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACg                                                                                         >  4:129183/1‑149 (MQ=255)                                                       gCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCaa                                                                                    <  3:267536/149‑1 (MQ=255)                                                       gCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCaa                                                                                    <  4:52603/149‑1 (MQ=255)                                                             atcTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGt                                                                                       >  4:321032/1‑140 (MQ=255)                                                             atcTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGt                                                                                       <  3:321032/140‑1 (MQ=255)                                                                    agaCATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGTTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAAc                                                                       <  4:115201/149‑1 (MQ=255)                                                                                           ttGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACtt                                                >  4:350440/1‑149 (MQ=255)                                                                                               tGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTcgcg                                            >  3:348176/1‑149 (MQ=255)                                                                                                  aCGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCtg                                         >  1:362084/1‑149 (MQ=255)                                                                                                   cGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCtgt                                        <  3:129183/149‑1 (MQ=255)                                                                                                      ccTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATgcgc                                                        <  1:283006/130‑1 (MQ=255)                                                                                                      ccTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATgcgc                                                        >  2:283006/1‑130 (MQ=255)                                                                                                       cTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAAc                                    >  2:242089/1‑149 (MQ=255)                                                                                                             a‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTc                               >  1:14713/1‑149 (MQ=255)                                                                                                             a‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTc                               >  5:52794/1‑149 (MQ=255)                                                                                                                     gTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTcc                        >  2:298919/1‑149 (MQ=255)                                                                                                                      tCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGGAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTccc                       >  2:97994/1‑149 (MQ=255)                                                                                                                      tCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTccc                       >  1:345164/1‑149 (MQ=255)                                                                                                                        aaaCGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTcccc                      <  5:348010/148‑1 (MQ=255)                                                                                                                              gcCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCtggtg               >  7:199845/1‑149 (MQ=255)                                                                                                                                cTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGc             >  7:239425/1‑149 (MQ=255)                                                                                                                                cTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGc             >  7:336729/1‑149 (MQ=255)                                                                                                                                 ttCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGTGCATTTCGCTGGCCACGGTCAAGCGCt                                                                               <  2:61354/82‑1 (MQ=255)                                                                                                                                 ttCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCt                                                                               >  1:61354/1‑82 (MQ=255)                                                                                                                                  tCCTTCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCCGCCCTGCGCTGCTACTTGGCGCTGTGAACGCCTCGTTTGCCCCGCTGGTGGCg            >  2:233072/1‑148 (MQ=255)                                                                                                                                   ccTGCTGGCCCGGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACtt                                                >  6:246570/1‑111 (MQ=255)                                                                                                                                   ccTGCTGGCCCGGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACtt                                                <  5:246570/111‑1 (MQ=255)                                                                                                                                   ccTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCg            <  3:350440/147‑1 (MQ=255)                                                                                                                                      cttgGCCCAGGCCGTCGACCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGaa         >  5:26685/3‑147 (MQ=255)                                                                                                                                        tgGCCCAGGCCGTCGACCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAgg     >  3:341295/1‑149 (MQ=255)                                                                                                                                          gCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGc    <  2:294070/148‑1 (MQ=255)                                                                                                                                           cccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCg  >  6:92283/1‑149 (MQ=255)                                                                                                                                                gTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCgg           >  7:300725/1‑135 (MQ=255)                                                                                                                                                gTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCgg           <  8:300725/135‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                               CTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTT‑CCA‑CCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCG  >  NC_002947/5233401‑5233685

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