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

GACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCG  >  REL606/4163679‑4163823                                                                          |                                                                        gACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGcc                                                                       >  1:286230/1‑76 (MQ=255)  aCAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCt                                                                      <  2:1609442/76‑1 (MQ=255)  aCAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCt                                                                      <  2:1065256/76‑1 (MQ=255)   cAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTc                                                                     <  2:21445/76‑1 (MQ=255)   cAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTc                                                                     >  1:1542381/1‑76 (MQ=255)   cAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTc                                                                     >  2:1774177/1‑76 (MQ=255)    aaaCGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCg                                                                    <  2:1789343/76‑1 (MQ=255)        gTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGc                                                                <  1:163797/76‑1 (MQ=255)         tGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCg                                                               <  1:1565754/76‑1 (MQ=255)          gcgcTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCgg                                                              >  1:1489859/1‑76 (MQ=255)           cgcTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCggg                                                             <  2:716205/76‑1 (MQ=255)           cgcTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCggg                                                             >  1:140513/1‑76 (MQ=255)            gcTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGt                                                            >  1:1496541/1‑76 (MQ=255)            gcTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGt                                                            >  1:655747/1‑76 (MQ=255)             cTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTGTGAAGAACTGCAGAGCCTCGAAGCGGGTc                                                           >  2:310312/1‑76 (MQ=255)             cTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTc                                                           >  1:273919/1‑76 (MQ=255)             cTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTc                                                           >  2:1434764/1‑76 (MQ=255)              tGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCt                                                          <  2:1523026/76‑1 (MQ=255)               ggAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTg                                                         >  2:1700645/1‑76 (MQ=255)               ggAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTg                                                         >  2:1409390/1‑76 (MQ=255)                gAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGCAGAACTGCAGAGCCTCGAAGCGGGTCTGt                                                        >  2:323366/1‑76 (MQ=255)                  aaTCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTc                                                      <  2:460004/76‑1 (MQ=255)                         gaaATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTa                                               <  2:738766/76‑1 (MQ=255)                          aaaTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTAt                                              >  2:1501782/1‑76 (MQ=255)                           aaTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGACGCGGGTCTGTTCAGCCGTATc                                             >  2:311783/1‑76 (MQ=255)                             tGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCg                                           <  2:1806803/76‑1 (MQ=255)                               cAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtg                                         >  1:1433372/1‑76 (MQ=255)                                 gCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgct                                       <  2:1085062/76‑1 (MQ=255)                                  cTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctg                                      <  2:1038419/76‑1 (MQ=255)                                   tCAAACAGGCGAAGAAAGACCTGGCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgt                                     <  1:465866/76‑1 (MQ=255)                                    cAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtg                                    >  1:319171/1‑76 (MQ=255)                                    cAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtg                                    >  1:1818039/1‑76 (MQ=255)                                    cAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtg                                    >  2:1497366/1‑76 (MQ=255)                                    cAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtg                                    >  2:1356479/1‑76 (MQ=255)                                       aCAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGtgctgtgctg                                 <  2:216466/76‑1 (MQ=255)                                         aGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGt                               <  2:422197/76‑1 (MQ=255)                                         aGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGt                               >  2:1518341/1‑76 (MQ=255)                                         aGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGt                               >  2:587869/1‑76 (MQ=255)                                          ggCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTa                              <  1:689349/76‑1 (MQ=255)                                             gaagaaAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGcc                           >  2:605254/1‑76 (MQ=255)                                              aagaaagaCCTGTCTGAAGAGCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCg                          <  2:932970/76‑1 (MQ=255)                                               agaaagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCgg                         >  1:1809018/1‑76 (MQ=255)                                               agaaagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGGAGCCgg                         >  1:709632/1‑76 (MQ=255)                                               agaaagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCCGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCgg                         >  1:832061/1‑76 (MQ=255)                                                 aaagcCCTGTCTGAAGACCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTg                       <  2:169475/76‑1 (MQ=255)                                                 aaagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGCAGCCGGTg                       >  2:577640/1‑76 (MQ=255)                                                  aagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTgg                      >  2:359308/1‑76 (MQ=255)                                                  aagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTgg                      >  2:415225/1‑76 (MQ=255)                                                   agaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGc                     >  1:13861/1‑76 (MQ=255)                                                    gaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCg                    >  2:1559211/1‑76 (MQ=255)                                                      ccTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGtt                  >  2:1049382/1‑76 (MQ=255)                                                      ccTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGtt                  >  2:1514688/1‑76 (MQ=255)                                                      ccTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGtt                  >  2:513713/1‑76 (MQ=255)                                                             gaagaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTg           >  2:947535/1‑76 (MQ=255)                                                             gaagaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTg           >  1:1527262/1‑76 (MQ=255)                                                             gaagaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTg           >  2:1651069/1‑76 (MQ=255)                                                              aagacCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTga          <  2:420724/76‑1 (MQ=255)                                                               agaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTgag         <  2:185075/76‑1 (MQ=255)                                                                gaaCTGCGGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTgaga        <  2:139323/76‑1 (MQ=255)                                                                gaaCTGCAGGGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTgaga        <  1:387539/76‑1 (MQ=255)                                                                gaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTgaga        <  2:1035533/76‑1 (MQ=255)                                                                gaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTgaga        >  2:1830943/1‑76 (MQ=255)                                                                 aaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGaa       >  1:197456/1‑76 (MQ=255)                                                                 aaCTGCAGAGCCGCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGaa       <  2:2847/76‑1 (MQ=255)                                                                   cTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGc     >  2:340214/1‑76 (MQ=255)                                                                      cAGAGCCTCGAAGCGGGTCTGTTCAGCTGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCg  >  2:146290/1‑76 (MQ=255)                                                                      cAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCg  >  1:1258250/1‑76 (MQ=255)                                                                          |                                                                        GACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCG  >  REL606/4163679‑4163823

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