Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,182,791	A→G	D516G (GAC→GGC)	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,791	0	A	G	100.0%	151.0 / NA	41	D516G (GAC→GGC)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/17);  total (24/17)

TGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTC  >  NC_000913/4182663‑4182922                                                                                                                                 |                                                                                                                                    tGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGaacaac                                                                                                                           >  2:391170/1‑139 (MQ=255) tGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGaacaac                                                                                                                           >  1:128745/1‑139 (MQ=255)    ggTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCg                                                                                                                        <  2:339994/139‑1 (MQ=255)    ggTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCg                                                                                                                        <  1:1215847/139‑1 (MQ=255)    ggTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCg                                                                                                                        <  1:959115/139‑1 (MQ=255)                  tgtctCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTAc                                                                                                          >  2:1011349/1‑139 (MQ=255)                       ctGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGcaca                                                                                                     <  1:1011349/139‑1 (MQ=255)                                    aTACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTAtctc                                                                                        >  1:174788/1‑139 (MQ=255)                                        ccTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCa                                                                                    <  2:174788/139‑1 (MQ=255)                                        ccTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCa                                                                                    >  1:222119/1‑139 (MQ=255)                                          tGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGAtt                                                                                                            <  2:699147/113‑1 (MQ=255)                                          tGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGAtt                                                                                                            >  1:699147/1‑113 (MQ=255)                                          tGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACt                                                                                  >  2:938707/1‑139 (MQ=255)                                           gATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTc                                                                                 >  1:519706/1‑139 (MQ=255)                                             tGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCgg                                                                               >  1:5781/1‑139 (MQ=255)                                             tGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCgg                                                                               >  1:679054/1‑139 (MQ=255)                                              gCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGc                                                                              <  2:679054/139‑1 (MQ=255)                                                  caGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAg                                                                          <  2:222119/139‑1 (MQ=255)                                                  caGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAg                                                                          >  2:119498/1‑139 (MQ=255)                                                        atGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTc                                                                    >  1:540268/1‑139 (MQ=255)                                                               aCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCg                                                             >  1:529761/1‑139 (MQ=255)                                                               aCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCg                                                             >  2:728555/1‑139 (MQ=255)                                                               aCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCg                                                             >  1:103268/1‑139 (MQ=255)                                                               aCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCg                                                             >  2:1289405/1‑139 (MQ=255)                                                                     aGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                                                       >  2:731272/1‑139 (MQ=255)                                                                      gCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAAcgt                                                                                                >  2:1288737/1‑97 (MQ=255)                                                                      gCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAAcgt                                                                                                <  1:1288737/97‑1 (MQ=255)                                                                         gATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCa                                                   <  1:119498/139‑1 (MQ=255)                                                                                      gTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGcaca                                                                                                     >  1:284965/1‑76 (MQ=255)                                                                                      gTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGcaca                                                                                                     <  2:284965/76‑1 (MQ=255)                                                                                       tGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCg                                     >  2:26206/1‑139 (MQ=255)                                                                                       tGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCg                                     >  2:676794/1‑139 (MQ=255)                                                                                                 cttcGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGa                                                          <  2:120306/108‑1 (MQ=255)                                                                                                 cttcGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGa                                                          >  1:120306/1‑108 (MQ=255)                                                                                                 cttcGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCg                                            <  2:824723/122‑1 (MQ=255)                                                                                                 cttcGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGCCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCg                                            >  1:824723/1‑122 (MQ=255)                                                                                                     ggTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGa                       <  2:128745/139‑1 (MQ=255)                                                                                                            gccagcTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTa                >  1:1002518/1‑139 (MQ=255)                                                                                                              cagcTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACg              <  1:391170/139‑1 (MQ=255)                                                                                                                          tttATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTc  <  1:664116/139‑1 (MQ=255)                                                                                                                          tttATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTc  <  1:1289405/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                    TGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTC  >  NC_000913/4182663‑4182922

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑