Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,182,881	A→C	E546A (GAA→GCA)	rpoB →	RNA polymerase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,182,881	0	A	C	100.0%	114.6 / NA	33	E546A (GAA→GCA)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (16/17);  total (16/17)

CGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTA  >  NC_000913/4182741‑4183021                                                                                                                                             |                                                                                                                                             cGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCgag                                                                                                                                      <  1:49973/149‑1 (MQ=255)                         ttCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTc                                                                                                              <  8:14699/149‑1 (MQ=255)                               ccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTc                                                                                                              <  8:18325/143‑1 (MQ=255)                               ccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTc                                                                                                              >  7:18325/1‑143 (MQ=255)                                      tctcAGTTTATGGACCAGAACAACCCGCTCTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAAt                                                                                                 >  6:88419/1‑149 (MQ=255)                                            tttATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAAc                                                                                           >  4:26572/1‑149 (MQ=255)                                                       gAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTc                                                                                >  2:38169/1‑149 (MQ=255)                                                       gAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTc                                                                                >  2:70623/1‑149 (MQ=255)                                                                      tGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGa                                                                 >  5:69028/1‑149 (MQ=255)                                                                        agaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGa                                                                 <  5:88419/147‑1 (MQ=255)                                                                           ttACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGtt                                                                                                                                          <  6:96510/71‑1 (MQ=255)                                                                           ttACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGtt                                                                                                                                          >  5:96510/1‑71 (MQ=255)                                                                                   aaaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActctct                                                       >  8:74599/1‑146 (MQ=255)                                                                                   aaaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActctct                                                       <  7:74599/146‑1 (MQ=255)                                                                                               tcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACg                                                                                                                                   >  1:19987/1‑58 (MQ=255)                                                                                               tcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACg                                                                                                                                   <  2:19987/58‑1 (MQ=255)                                                                                                     cTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTaa                                  >  7:46395/1‑149 (MQ=255)                                                                                                      tCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACt                                                                                                                        <  3:37806/62‑1 (MQ=255)                                                                                                      tCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACt                                                                                                                        >  4:37806/1‑62 (MQ=255)                                                                                                           ccAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                            >  8:103165/1‑149 (MQ=255)                                                                                                            cAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                           <  8:46395/149‑1 (MQ=255)                                                                                                             aGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACg                          <  6:69028/149‑1 (MQ=255)                                                                                                              ggCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACgg                         <  1:70623/149‑1 (MQ=255)                                                                                                                       aCCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTga                <  3:26572/149‑1 (MQ=255)                                                                                                                       aCCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTga                <  7:103165/149‑1 (MQ=255)                                                                                                                       aCCCGTGAACGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTga                <  1:38169/149‑1 (MQ=255)                                                                                                                              aaCGTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCg         >  6:13281/1‑149 (MQ=255)                                                                                                                                 gTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGt                                                                     >  2:140571/1‑86 (MQ=255)                                                                                                                                 gTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGt                                                                     <  1:140571/86‑1 (MQ=255)                                                                                                                                 gTGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTAt      <  5:13281/149‑1 (MQ=255)                                                                                                                                  tGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTga                <  1:89383/138‑1 (MQ=255)                                                                                                                                  tGCAGGCTTCGCAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTga                >  2:89383/1‑138 (MQ=255)                                                                                                                                       gCTTCGCAGTTCGCGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCgga  >  5:72119/1‑145 (MQ=255)                                                                                                                                             |                                                                                                                                             CGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTA  >  NC_000913/4182741‑4183021

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