Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,184,140	A→G	I966V (ATC→GTC)	rpoB →	RNA polymerase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,184,140	0	A	G	100.0%	158.5 / NA	43	I966V (ATC→GTC)	rpoB	RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/21);  total (22/21)

GAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACA  >  NC_000913/4184062‑4184215                                                                               |                                                                            gAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTc                                                                           <  1:829973/81‑1 (MQ=255)  aaaaaGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTcc                                                                          <  1:829974/81‑1 (MQ=255)   aaaaGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCt                                                                         >  2:830012/1‑81 (MQ=255)    aaaGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTc                                                                        <  1:829975/81‑1 (MQ=255)    aaaGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTc                                                                        <  1:829976/81‑1 (MQ=255)    aaaGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTc                                                                        >  1:830035/1‑81 (MQ=255)     aaGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTc                                                                        >  2:830016/1‑80 (MQ=255)     aaGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTc                                                                        >  1:830011/1‑80 (MQ=255)       gACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCg                                                                       >  1:830030/1‑79 (MQ=255)           aaCGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCgg                                                                 <  1:829978/81‑1 (MQ=255)           aaCGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCgg                                                                 <  2:829977/81‑1 (MQ=255)                   cTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGtt                                                          >  2:830021/1‑80 (MQ=255)                   cTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGtt                                                          >  1:830026/1‑80 (MQ=255)                   cTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTACTCGAAGCGGGTCTGTTc                                                         <  2:829979/81‑1 (MQ=255)                    tGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGtt                                                          >  1:830037/1‑79 (MQ=255)                       aaaTCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGcc                                                     <  2:829980/81‑1 (MQ=255)                            gaagaaATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATc                                                <  2:829981/81‑1 (MQ=255)                                aaaTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtg                                            <  1:829982/81‑1 (MQ=255)                                aaaTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtg                                            >  2:830015/1‑81 (MQ=255)                                 aaTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgc                                           >  2:830031/1‑81 (MQ=255)                                  aTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgct                                          >  2:830048/1‑81 (MQ=255)                                    gCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgt                                        >  1:830029/1‑81 (MQ=255)                                     cAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtg                                       >  1:830039/1‑81 (MQ=255)                                     cAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtg                                       <  2:829983/81‑1 (MQ=255)                                      aGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtgc                                      >  2:830038/1‑81 (MQ=255)                                       gCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtgct                                     <  1:829984/81‑1 (MQ=255)                                           aaaCAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTa                                 <  2:829985/81‑1 (MQ=255)                                            aaCAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAg                                >  1:830032/1‑81 (MQ=255)                                             aCAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGc                               >  1:830023/1‑81 (MQ=255)                                               aGGCGAAGAAAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCg                             <  1:829986/81‑1 (MQ=255)                                                   gaagaaAGACCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTgg                         <  1:829987/81‑1 (MQ=255)                                                    aagaaagaCCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGc                        >  2:830042/1‑81 (MQ=255)                                                       aaagaCCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGtt                     >  1:830028/1‑81 (MQ=255)                                                         agaCCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTg                    >  1:830033/1‑80 (MQ=255)                                                           aCCTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAg                 >  1:830041/1‑81 (MQ=255)                                                            ccTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGc                <  1:829988/81‑1 (MQ=255)                                                             cTGTCTGAAGAACTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCt               >  1:830046/1‑81 (MQ=255)                                                                   gaagaaCTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAg         >  2:830044/1‑81 (MQ=255)                                                                    aagaaCTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTgaga           <  2:829989/78‑1 (MQ=255)                                                                      gaaCTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTc      <  2:829991/81‑1 (MQ=255)                                                                      gaaCTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTc      <  1:829990/81‑1 (MQ=255)                                                                        aCTGCAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGa    <  1:829992/81‑1 (MQ=255)                                                                            cAGGTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACa  <  2:829993/79‑1 (MQ=255)                                                                               |                                                                            GAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACA  >  NC_000913/4184062‑4184215

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