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%	127.8 / NA	33	I966S (ATC→AGC)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (17/16);  total (17/16)

ACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCG  >  REL606/4163680‑4163823                                                                         |                                                                        aCAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCt                                                                      <  2:130564/76‑1 (MQ=255)   aaaCGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCg                                                                    <  2:280053/76‑1 (MQ=255)     aCGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGaa                                                                  >  1:98760/1‑76 (MQ=255)       gTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGc                                                                <  1:119714/76‑1 (MQ=255)           gcTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGt                                                            <  1:619680/76‑1 (MQ=255)              ggAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTg                                                         <  2:550725/76‑1 (MQ=255)               gAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGt                                                        <  2:1262549/76‑1 (MQ=255)               gAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGt                                                        <  1:490414/76‑1 (MQ=255)                 aaTCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTc                                                      <  1:294709/76‑1 (MQ=255)                  aTCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCa                                                     <  2:65120/76‑1 (MQ=255)                  aTCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCa                                                     >  1:140764/1‑76 (MQ=255)                  aTCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCa                                                     >  2:215810/1‑76 (MQ=255)                    cGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACATGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGc                                                   >  1:246516/1‑76 (MQ=255)                       agaaATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGt                                                <  2:335220/76‑1 (MQ=255)                               aGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgc                                        >  1:1309813/1‑76 (MQ=255)                                 cTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctg                                      <  2:942120/76‑1 (MQ=255)                                  acAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgt                                     >  1:859075/2‑76 (MQ=255)                                   cAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtg                                    <  2:148325/76‑1 (MQ=255)                                    caaCAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCctgctgtgc                                   >  1:817770/2‑76 (MQ=255)                                       cAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTgg                                >  2:105897/1‑76 (MQ=255)                                           cGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGc                            <  1:1468115/76‑1 (MQ=255)                                            gaagaaAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGcc                           <  2:814019/76‑1 (MQ=255)                                            gaagaaAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGcc                           <  2:602905/76‑1 (MQ=255)                                                aaagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTg                       >  2:1457774/1‑76 (MQ=255)                                                aaagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTg                       >  2:378620/1‑76 (MQ=255)                                                 aagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCcggcgg                      >  2:677012/1‑76 (MQ=255)                                                 aagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTgg                      >  1:1531098/1‑76 (MQ=255)                                                  agaCGTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGc                     <  2:611968/76‑1 (MQ=255)                                                    aCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCGGGTGGCGt                   >  2:81325/1‑76 (MQ=255)                                                     ccTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGtt                  >  2:995552/1‑76 (MQ=255)                                                             aagaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTga          >  2:755822/1‑76 (MQ=255)                                                               gaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTgaga        >  2:551956/1‑76 (MQ=255)                                                                     cAGAGCCTCGAAGCGGGTCTGTTCAGCCTTAGCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCg  <  1:187743/76‑1 (MQ=255)                                                                     cAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCg  >  2:118184/1‑76 (MQ=255)                                                                         |                                                                        ACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCG  >  REL606/4163680‑4163823

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