Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,163,752	T→G	I966S (ATC→AGC)	rpoB →	DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,163,752	0	T	G	100.0%	92.9 / NA	23	I966S (ATC→AGC)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (12/11);  total (12/11)

ACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCT  >  REL606/4163680‑4163821                                                                         |                                                                      aCAAACGTGCGCTGGAAATCGAAGAAATGTAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCt                                                                    <  1:958300/76‑1 (MQ=255)  cAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTc                                                                   <  2:667987/76‑1 (MQ=255)  cAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTc                                                                   <  1:217100/76‑1 (MQ=255)       gTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGc                                                              <  1:413220/76‑1 (MQ=255)           gcTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGt                                                          <  2:684084/76‑1 (MQ=255)             tGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCt                                                        >  1:498237/1‑76 (MQ=255)              ggAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTg                                                       <  2:339456/76‑1 (MQ=255)               gAAATNGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGt                                                      <  2:778335/76‑1 (MQ=255)                    cGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGc                                                 >  1:223078/1‑76 (MQ=255)                                gCTCAAACAGGCGAAGAAAGACCTGTCTGAAGANCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgct                                     <  1:780655/76‑1 (MQ=255)                                    aaaCAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCNGCCGTATCCgtgctgtgc                                 >  1:751452/1‑76 (MQ=255)                                    aaaCAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtgc                                 >  2:825772/1‑76 (MQ=255)                                           cGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGc                          >  2:569546/1‑76 (MQ=255)                                            gaagaaAGACCTGTCTGAAGAACTGCAGAGCCTCGACGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGcc                         >  1:360722/1‑76 (MQ=255)                                             aagaaagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCg                        >  1:461739/1‑76 (MQ=255)                                               gaaagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGt                      <  2:452077/76‑1 (MQ=255)                                                aaagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTNTCCGTGCTGTGCTGGTAGCCGGTg                     <  2:172052/76‑1 (MQ=255)                                                aaagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTg                     <  2:767412/76‑1 (MQ=255)                                                aaagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTg                     >  2:924287/1‑76 (MQ=255)                                                 aagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCGGTGCTGGTAGCCGGTgg                    <  2:313925/76‑1 (MQ=255)                                                     ccTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGtt                <  2:591625/76‑1 (MQ=255)                                                        gTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGCAGCCGGTGGCGTTGaa             >  2:624942/1‑76 (MQ=255)                                                         tCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAg            <  2:619036/76‑1 (MQ=255)                                                            gaagaaCTGCAGTGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTg         >  2:147833/1‑76 (MQ=255)                                                            gaagaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTg         >  2:523170/1‑76 (MQ=255)                                                                 aCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAg    <  2:287154/76‑1 (MQ=255)                                                                   tGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCt  >  1:421165/1‑76 (MQ=255)                                                                         |                                                                      ACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCT  >  REL606/4163680‑4163821

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