Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	535,046	T→A	12.0%	intergenic (+31/‑89)	nusG → / → rplK	transcription antipausing factor NusG/50S ribosomal protein L11

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	535,046	0	T	A	12.0%	141.3 / 5.6	55	intergenic (+31/‑89)	nusG/rplK	transcription antipausing factor NusG/50S ribosomal protein L11
Reads supporting (aligned to +/- strand):  ref base T (23/25);  new base A (4/3);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAATGGTAGTGTCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAGGCCGGCCAG  >  NC_002947/534901‑535185                                                                                                                                                  |                                                                                                                                            cAATGGTAGTGTCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCtgt                                                                                                                                           <  4:233226/148‑1 (MQ=255)         gtgtCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCgg                                                                                                                                  <  5:197965/149‑1 (MQ=255)            tCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCACACAGGGCTGCgggg                                                                                                                                <  8:57235/148‑1 (MQ=255)                 gaagTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTgt                                                                                                                          <  4:331592/149‑1 (MQ=255)                       aaCTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCAc                                                                                                                    >  3:298936/1‑149 (MQ=255)                          tACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTgg                                                                                                                 >  2:215879/1‑149 (MQ=255)                             gAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGa                                                                                                               >  4:47121/1‑148 (MQ=255)                              aaaaGAGTCGCCTGCAGGTTGCGGTGCTAATTTTCGGTCGCTCTACTCCGGTGGAGATCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCAGCAGCCCTGTGTGCTGCGGGGTTTTGTGGTCACTgggtt                                                                                                              >  4:58731/1‑146 (MQ=255)                               aaaGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATa                                                                                                             >  7:302719/1‑148 (MQ=255)                                    gTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGc                                                                                                       <  1:215879/149‑1 (MQ=255)                                     tCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGc                                                                                                       >  1:341999/1‑148 (MQ=255)                                     tCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTCTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCa                                                                                                      <  8:302719/149‑1 (MQ=255)                                       gCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCAGTGAGCTGAGGGGGTATGTCGGCACTGGGATAAAAAGCAAg                                                                                                    >  8:299453/1‑149 (MQ=255)                                       gCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCACACAGGGCTGCggg                                                                                                                                 >  6:247995/1‑120 (MQ=255)                                       gCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCACACAGGGCTGCggg                                                                                                                                 <  5:247995/120‑1 (MQ=255)                                           gCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCAt                                                                                                >  2:315272/1‑149 (MQ=255)                                               gTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCGGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCgg                                                                                            >  8:188022/1‑149 (MQ=255)                                                    ggTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACACCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATACGGGGAGc                                                                                       >  7:38018/1‑149 (MQ=255)                                                           aTTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCACACAGGGCTGCgggg                                                                                                                                >  3:261846/1‑101 (MQ=38)                                                           aTTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCACACAGGGCTGCgggg                                                                                                                                <  4:261846/101‑1 (MQ=38)                                                              ttCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCg                                                                              <  1:144344/148‑1 (MQ=255)                                                                   tCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTa                                                                        <  2:349352/149‑1 (MQ=255)                                                                     gCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTAcc                                                                      >  7:303466/1‑149 (MQ=255)                                                                     gCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTAcc                                                                      <  1:315272/149‑1 (MQ=255)                                                                            tCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCCGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGAtttt                                                               >  3:173610/1‑149 (MQ=255)                                                                              cGGTGGAGCTCGAGTTCGGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGAtttt                                                               <  7:188022/147‑1 (MQ=255)                                                                                     gCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACACCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTgg                                                                                                                 >  4:141243/1‑90 (MQ=255)                                                                                      cTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAg                                                                                 >  2:216139/1‑121 (MQ=255)                                                                                      cTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAg                                                                                 <  1:216139/121‑1 (MQ=255)                                                                                          aGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAGAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATgg                                                 <  1:99268/149‑1 (MQ=255)                                                                                          aGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATgg                                                 <  6:153853/149‑1 (MQ=255)                                                                                            ttCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCt                                               <  1:57023/149‑1 (MQ=255)                                                                                                        gAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATgg                                                 <  7:82982/135‑1 (MQ=255)                                                                                                        gAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATgg                                                 >  8:82982/1‑135 (MQ=255)                                                                                                        gAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAg                                   <  2:341999/149‑1 (MQ=255)                                                                                                        gAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAg                                   <  4:173610/149‑1 (MQ=255)                                                                                                          aaaGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGTAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGc                                 <  8:38018/149‑1 (MQ=255)                                                                                                               tCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACa                            <  3:47121/149‑1 (MQ=255)                                                                                                               tCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACACAAGTCATCCGGGGAGCCATCAGGCGt                                                                             >  8:254305/1‑100 (MQ=255)                                                                                                                 tAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCt                                               <  5:40292/128‑1 (MQ=255)                                                                                                                 tAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCt                                               >  6:40292/1‑128 (MQ=255)                                                                                                                     cGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCTGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCAGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAg                       <  3:58731/148‑1 (MQ=255)                                                                                                                     cGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAg                       >  7:127189/1‑148 (MQ=255)                                                                                                                     cGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAg                       <  8:127189/148‑1 (MQ=255)                                                                                                                      ggACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAg                                                                                                    >  4:74395/1‑70 (MQ=255)                                                                                                                      ggACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAg                                                                                                    <  3:74395/70‑1 (MQ=255)                                                                                                                          gAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGc                                                >  6:181937/1‑118 (MQ=255)                                                                                                                          gAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGc                                                <  5:181937/118‑1 (MQ=255)                                                                                                                                gacaaaaCCCCGCAGCACACAGGGCCGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAg           >  6:231546/7‑149 (MQ=255)                                                                                                                                 cccATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTTCCCGAGTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAgg          >  2:266771/1‑149 (MQ=255)                                                                                                                                   cATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGTTAAAACGCATGTCATCCGGGGAGCCATCAGGCGTTGGTAAAAGATTTTGGAGTAGCTAATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCATGTTAAgt          >  5:161533/1‑146 (MQ=255)                                                                                                                                   cATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTg                                                              >  8:135342/1‑95 (MQ=255)                                                                                                                                   cATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTg                                                              <  7:135342/95‑1 (MQ=255)                                                                                                                                      tctccGCATCCCTGTGTGCTGCGGGGTTTTGTTGTAACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAggccgg      <  5:263876/145‑1 (MQ=255)                                                                                                                                         ccGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAGGCCggccag  <  7:299453/149‑1 (MQ=255)                                                                                                                                          cGCAGCCCTGTGTGCTGCGGGGTTTTGTTGTCACTGGGATAAAACGCAAGTCATCCGGGGAGACATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAgtacgg      >  6:263876/1‑139 (MQ=255)                                                                                                                                                  |                                                                                                                                            CAATGGTAGTGTCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAGGCCGGCCAG  >  NC_002947/534901‑535185

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