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%	212.3 / NA	57	D516G (GAC→GGC)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/25);  total (32/25)

TGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAAT  >  NC_000913/4182663‑4182926                                                                                                                                 |                                                                                                                                        tGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGaacaac                                                                                                                               >  1:1091622/1‑139 (MQ=255) tGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGaacaac                                                                                                                               >  2:178452/1‑139 (MQ=255) tGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGaacaac                                                                                                                               >  1:453046/1‑139 (MQ=255)    ggTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCg                                                                                                                            <  2:1091622/139‑1 (MQ=255)                    tctcTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGc                                                                                                            >  2:82531/1‑139 (MQ=255)                       ctGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGcaca                                                                                                         <  1:1130311/139‑1 (MQ=255)                       ctGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGcaca                                                                                                         <  1:1130959/139‑1 (MQ=255)                        tGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACaa                                                                                                        >  2:1245273/1‑139 (MQ=255)                              aTCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAAcgtcg                                                                                                  >  2:63569/1‑139 (MQ=255)                              aTCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAAcgtcg                                                                                                  >  1:579994/1‑139 (MQ=255)                                      aCCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCg                                                                                          <  1:82531/139‑1 (MQ=255)                                          tGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACt                                                                                      >  2:898964/1‑139 (MQ=255)                                          tGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACt                                                                                      <  1:63569/139‑1 (MQ=255)                                           gATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTc                                                                                     >  1:1400372/1‑139 (MQ=255)                                                 acaGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAACCAGCTGTCTCAGTTTATGGGCCAGAACAACCCg                                                                                                                            >  2:242323/1‑94 (MQ=255)                                                 acaGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAACCAGCTGTCTCAGTTTATGGGCCAGAACAACCCg                                                                                                                            <  1:242323/94‑1 (MQ=255)                                                          gATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTg                                                                      <  2:646535/139‑1 (MQ=255)                                                               aCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCg                                                                                                                            <  1:121271/80‑1 (MQ=255)                                                               aCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCg                                                                                                                            >  2:121271/1‑80 (MQ=255)                                                               aCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACt                                                                                      >  2:747751/1‑118 (MQ=255)                                                               aCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACt                                                                                      <  1:747751/118‑1 (MQ=255)                                                               aCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCg                                                                 >  1:599682/1‑139 (MQ=255)                                                               aCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCg                                                                 >  2:671289/1‑139 (MQ=255)                                                               aCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCg                                                                 >  1:78080/1‑139 (MQ=255)                                                               aCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCg                                                                 >  2:614918/1‑139 (MQ=255)                                                               aCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACACACAAACGTCGTATCTCCGCACt                                                                                      <  1:1132411/118‑1 (MQ=255)                                                               aCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACACACAAACGTCGTATCTCCGCACt                                                                                      >  2:1132411/1‑118 (MQ=255)                                                                     aGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                                                           >  1:260098/1‑139 (MQ=255)                                                                        cGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGc                                                        <  1:921205/139‑1 (MQ=255)                                                                                 cagcagTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa                                               >  2:331889/1‑139 (MQ=255)                                                                                       tGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCTAGAACAACCCg                                                                                                                            <  2:852238/56‑1 (MQ=255)                                                                                       tGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCg                                                                                                                            >  1:852238/1‑56 (MQ=255)                                                                                       tGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCg                                         <  1:671289/139‑1 (MQ=255)                                                                                       tGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCg                                         <  2:599682/139‑1 (MQ=255)                                                                                       tGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCg                                         >  2:387427/1‑139 (MQ=255)                                                                                         aaaGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCgag                                       >  1:572086/1‑139 (MQ=255)                                                                                           agagTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGAc                                     <  2:579994/139‑1 (MQ=255)                                                                                              gTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAAcc                                                                                                                              >  1:595568/1‑47 (MQ=255)                                                                                              gTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAAcc                                                                                                                              <  2:595568/47‑1 (MQ=255)                                                                                              gTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAg                                                                              >  1:1233454/1‑95 (MQ=255)                                                                                              gTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAg                                                                              <  2:1233454/95‑1 (MQ=255)                                                                                              gTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTa                                  >  2:1248799/1‑139 (MQ=255)                                                                                                 cttcGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTTGAGACGTacac                               >  2:1320853/1‑139 (MQ=255)                                                                                                 cttcGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTacac                               >  1:671970/1‑139 (MQ=255)                                                                                                  ttcGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACAcc                              <  2:1400372/139‑1 (MQ=255)                                                                                                   tcGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACAccc                             >  1:1186800/1‑139 (MQ=255)                                                                                                    cGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCg                                         >  2:343987/1‑126 (MQ=255)                                                                                                    cGGTTCCAGCCAGCTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCg                                         <  1:343987/126‑1 (MQ=255)                                                                                                           agccagcTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACt                     <  2:1186800/139‑1 (MQ=255)                                                                                                            gccagcTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTa                    <  2:383527/139‑1 (MQ=255)                                                                                                                 cTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGca                                                                                                           >  1:847324/1‑47 (MQ=255)                                                                                                                 cTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGca                                                                                                           <  2:847324/47‑1 (MQ=255)                                                                                                                 cTGTCTCAGTTTATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTc               >  1:1342250/1‑139 (MQ=255)                                                                                                                          tttATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTc      <  2:260098/139‑1 (MQ=255)                                                                                                                            tATGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCa    <  1:1248799/139‑1 (MQ=255)                                                                                                                              tGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCg                                         <  2:1025394/100‑1 (MQ=255)                                                                                                                              tGGGCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAAt  >  2:951462/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                        TGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAAT  >  NC_000913/4182663‑4182926

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