Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,149,461	C→T	7.8%	S98N (AGC→AAC)	PP_1008 ←	RNA polymerase subunit sigma‑24

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,149,461	0	C	T	7.8%	143.6 / 2.5	51	S98N (AGC→AAC)	PP_1008	RNA polymerase subunit sigma‑24
Reads supporting (aligned to +/- strand):  ref base C (25/22);  new base T (2/2);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

AAACACTCAGTTGTTCGGCGATTTCCAGATAGGTGGCGCCATGCAGGCGGCT‑GAGGATGAAAATGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTGCGCGCGGTCTGAAACACGAACGGCTCGATGTGCTCGATGG  >  NC_002947/1149315‑1149605                                                                                                                                                    |                                                                                                                                                 aaaCACTCAGTTGTTCGGCGATTTCCAGATAGGTGGCGCCATGCAGGGGGCT‑GAGGATGAAAATGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCGGAGGTGCTTGAGCAACTGCTCGGCGGGGGCGGCGTCCTCGgtgctg                                                                                                                                                >  4:281979/1‑149 (MQ=255) aaaCACTCAGTTGTTCGGCGATTTCCAGATAGGTGGCGCCATGCAGGCGGCT‑GAGGATGAAAATGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGctgctg                                                                                                                                                <  6:218526/149‑1 (MQ=255)     acTCAGTTGTTCGGCGATTTCCAGATAGGTGGCGCCATGCAGGCGGCTGGAGGATGAAAATGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCtggtg                                                                                                                                             >  1:82888/1‑149 (MQ=255)       tCAGTTGTTCGGCGATTTCCAGATAGGTGGCGCCATGCAGGCGGCT‑GAGGATGAAAATGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGtggct                                                                                                                                          >  1:127550/1‑149 (MQ=255)         aGTTGTTCGGCGATTTCCAGATAGGTGGCGCCATGCAGGCGGCT‑GAGGATGAAAATGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGtggctgg                                                                                                                                        >  8:119756/1‑149 (MQ=255)              ttCGGCGATTTCCAGATAGGTGGCGCCATGCAGGCGGCT‑GAGGATGAAAATGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCgcgg                                                                                                                                   <  8:155820/149‑1 (MQ=255)                cGGCGATTTCCAGATAGGTGGCGCCATGCAGGCGGCT‑GAGGATGAAAATGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCgcggt                                                                                                                                  <  2:127550/148‑2 (MQ=255)                          cAGATAGGTGGGGCCATGCAGGCGGCT‑GAGGATGAAAATGCGTTGTTGGCGTTCGCTTAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTgc                                                                                                                        <  7:247151/148‑1 (MQ=255)                                 gTGGCGCCATGCAGGCGGCT‑GAGGATGAAAATGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTtc                                                                                                                <  2:82888/149‑1 (MQ=255)                                   ggCGCCATGCAGGCGGCT‑GAGGATGAAAATGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTtcgt                                                                                                              >  8:69626/1‑149 (MQ=255)                                     cgcCATGCAGGCGGCT‑GAGGATGAAAATGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTtcgtcg                                                                                                            <  8:314491/149‑1 (MQ=255)                                        cATGCAGGCGGCT‑GAGGATGAAAATGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACt                                                                                                                            <  8:287222/130‑1 (MQ=255)                                        cATGCAGGCGGCT‑GAGGATGAAAATGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACt                                                                                                                            >  7:287222/1‑130 (MQ=255)                                           gCAGGCGGCT‑GAGGATGAAAATGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCAcg                                                                                                      >  5:192796/1‑149 (MQ=255)                                                            gAAAATGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCgcgg                                                                                                                                   >  4:149051/1‑104 (MQ=255)                                                            gAAAATGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCgcgg                                                                                                                                   <  3:149051/104‑1 (MQ=255)                                                              aaaTGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGTTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGt                                                                                    >  6:152357/1‑149 (MQ=255)                                                              aaaTGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGTTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGt                                                                                    >  6:157601/1‑149 (MQ=255)                                                                 tGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCgg                                                                                                           <  3:147028/123‑1 (MQ=255)                                                                 tGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGAAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCgg                                                                                                           >  4:147028/1‑123 (MQ=255)                                                                     ttgttgGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCa                                                                             >  2:315825/1‑149 (MQ=255)                                                                             gTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCATGCACCCGGCGt                                                                     >  2:222013/1‑149 (MQ=255)                                                                                 gCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCAc                                                                 <  3:281979/149‑1 (MQ=255)                                                                                   tCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGc                                                                       >  6:176436/1‑141 (MQ=255)                                                                                   tCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGc                                                                       <  5:176436/141‑1 (MQ=255)                                                                                   tCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACg                                                                >  3:331009/1‑148 (MQ=255)                                                                                   tCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACg                                                                >  6:58524/1‑148 (MQ=255)                                                                                      gTTGATTGAGGCTGACGCTGAAGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCAGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAgg                                                            >  8:52279/1‑149 (MQ=255)                                                                                       ttgattgaGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAggt                                                           >  5:233779/1‑149 (MQ=255)                                                                                         gattgaGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACtgtg                                                                                                                         <  1:145324/85‑1 (MQ=255)                                                                                         gattgaGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACtgtg                                                                                                                         >  2:145324/1‑85 (MQ=255)                                                                                         gattgaGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAggtg                                                          <  5:58524/148‑1 (MQ=255)                                                                                               ggCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAgc                                                   <  6:192796/149‑1 (MQ=255)                                                                                                    aCGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCa                                               <  1:222013/148‑1 (MQ=255)                                                                                                      gCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGt                                            >  6:213702/1‑149 (MQ=255)                                                                                                      gCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGt                                            >  6:218866/1‑149 (MQ=255)                                                                                                          aGGTGCTTGAGCAAGTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTgc                                         <  7:69626/148‑1 (MQ=255)                                                                                                              gCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTgcgc                                       >  1:268916/1‑146 (MQ=255)                                                                                                              gCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTgcgc                                       <  2:268916/146‑1 (MQ=255)                                                                                                                    gCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTGCGCGCGGTCTGa                              <  1:315825/149‑1 (MQ=255)                                                                                                                      aaCTGCTCGGCGTGGGCGGCATCCTCGCTGTTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTGCGCGCGGTCTGaaa                            <  5:152357/149‑1 (MQ=255)                                                                                                                      aaCTGCTCGGCGTGGGCGGCATCCTCGCTGTTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTGCGCGCGGTCTGaaa                            <  5:157601/149‑1 (MQ=255)                                                                                                                      aaCTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTGCGCGCGGTCTGaaa                            <  2:125598/149‑1 (MQ=255)                                                                                                                           cTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAggctg                                                         <  2:70393/115‑4 (MQ=255)                                                                                                                           cTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAggctg                                                         >  1:70393/1‑112 (MQ=255)                                                                                                                            tCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCAt                                                                                        >  8:73022/1‑83 (MQ=255)                                                                                                                            tCGGCGTGGGCGGCATCCTCACTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTGCGCGCGGTCTGAAACACGaa                      >  3:245580/1‑149 (MQ=255)                                                                                                                                       ggcATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTGCGCGCGGTCTGAAACACGAACGGCTCGa              >  2:18386/1‑146 (MQ=255)                                                                                                                                       ggcATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTGCGCGCGGTCTGAAACACGAACGGCTCGa              <  1:18386/146‑1 (MQ=255)                                                                                                                                         cATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTGCGCGCGGTCTGAAaca                          >  7:331602/1‑132 (MQ=255)                                                                                                                                         cATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTGCGCGCGGTCTGAAaca                          <  8:331602/132‑1 (MQ=255)                                                                                                                                              tCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTGCGCGCGGTCTGAAACACGAACGGCt                 <  3:215438/136‑1 (MQ=255)                                                                                                                                              tCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTGCGCGCGGTCTGAAACACGAACGGCt                 >  4:215438/1‑136 (MQ=255)                                                                                                                                                gctgctGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTGCGCGCGGTCTGAAACACGAACGGCTCGATGTGCTCGATgg  >  1:240913/1‑149 (MQ=255)                                                                                                                                                gctgctGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTGCGCGCGGTCTGAAACACGAACGGCTCGATGTGCTCGATg   <  4:331009/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 AAACACTCAGTTGTTCGGCGATTTCCAGATAGGTGGCGCCATGCAGGCGGCT‑GAGGATGAAAATGCGTTGTTGGCGTTCGCTCAGTTGATTGAGGCTGACGCTGAGGTGCTTGAGCAACTGCTCGGCGTGGGCGGCATCCTCGCTGCTGGTGGCTGGCGCGGCGACACTGTGCAGCACTTCGTCGGGCACGTCGTCGACCAGCATGCGTGCCTGCACCCGGCGTGCACGCAGGTGGTCCAGCGCCAGGTTGCGCGCGGTCTGAAACACGAACGGCTCGATGTGCTCGATGG  >  NC_002947/1149315‑1149605

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