Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	535,050	T→G	10.9%	intergenic (+35/‑85)	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,050	0	T	G	10.9%	147.3 / 11.3	55	intergenic (+35/‑85)	nusG/rplK	transcription antipausing factor NusG/50S ribosomal protein L11
Reads supporting (aligned to +/- strand):  ref base T (23/25);  new base G (3/3);  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 = 4.57e-01

GTGTCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAGGCCGGCCAG  >  NC_002947/534909‑535185                                                                                                                                              |                                                                                                                                        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                                                                      <  1:315272/149‑1 (MQ=255)                                                             gCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTAcc                                                                      >  7:303466/1‑149 (MQ=255)                                                                    tCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCCGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGAtttt                                                               >  3:173610/1‑149 (MQ=255)                                                                      cGGTGGAGCTCGAGTTCGGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGAtttt                                                               <  7:188022/147‑1 (MQ=255)                                                                             gCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACACCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTgg                                                                                                                 >  4:141243/1‑90 (MQ=255)                                                                              cTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAg                                                                                 <  1:216139/121‑1 (MQ=255)                                                                              cTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAg                                                                                 >  2:216139/1‑121 (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)                                                                                                                                              |                                                                                                                                        GTGTCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCATCCCATACCCCGCAGCCCTGTGTGCTGCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAGGCCGGCCAG  >  NC_002947/534909‑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: ‑