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

ACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGA  >  REL606/4163680‑4163824                                                                         |                                                                         aCAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCt                                                                       <  2:103480/76‑1 (MQ=255) aCAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCt                                                                       <  2:1735017/76‑1 (MQ=255)  cAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTc                                                                      >  1:1221743/1‑76 (MQ=255)  cAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTc                                                                      <  2:207081/76‑1 (MQ=255)     aCGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAAACCTGTCTGAAGAACTGCAGAGCCTCGaa                                                                   <  1:1605909/76‑1 (MQ=255)           gcTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGt                                                             >  1:1716868/1‑76 (MQ=255)            cTGGAAATCGAAGAAATGCAGCTCAAACAGTCGAAGAAATACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTc                                                            <  2:563274/76‑1 (MQ=255)             tGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCt                                                           <  2:1089048/76‑1 (MQ=255)             tGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCt                                                           <  2:437855/76‑1 (MQ=255)              ggAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTg                                                          >  1:1899329/1‑76 (MQ=255)               gAAATCGAAGAAATGCAGCTCAAGTGGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGt                                                         <  2:257894/76‑1 (MQ=255)               gAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGt                                                         >  1:1436763/1‑76 (MQ=255)               gAAATCGAAGAAATGCAGCTAAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGt                                                         <  2:1549655/76‑1 (MQ=255)                aaaTCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGtt                                                        <  2:397635/76‑1 (MQ=255)                 aaTCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGCTCTGTTc                                                       <  2:563812/76‑1 (MQ=255)                  aTCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCa                                                      <  1:1222991/76‑1 (MQ=255)                    cGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGc                                                    >  1:1660934/1‑76 (MQ=255)                      aagaaATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCg                                                  <  1:1246546/76‑1 (MQ=255)                        gaaATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTa                                                <  2:1284486/76‑1 (MQ=255)                        gaaATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTa                                                >  2:879312/1‑76 (MQ=255)                        gaaATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTa                                                >  1:713103/1‑76 (MQ=255)                        gaaATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTa                                                >  1:1149131/1‑76 (MQ=255)                        gaaATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTa                                                >  1:1344628/1‑76 (MQ=255)                           aTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATcc                                             >  2:543607/1‑76 (MQ=255)                           aTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATcc                                             >  2:50266/1‑76 (MQ=255)                           aTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATcc                                             <  2:1045009/76‑1 (MQ=255)                           aTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATcc                                             >  2:982129/1‑76 (MQ=255)                                   cAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtg                                     >  1:812859/1‑76 (MQ=255)                                   cAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtg                                     >  2:1706737/1‑76 (MQ=255)                                    aaaCAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAg                                                     <  1:1676597/59‑1 (MQ=255)                                    aaaCAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAg                                                     >  2:1676597/1‑59 (MQ=255)                                      aCAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGtgctgtgctg                                  <  2:587373/76‑1 (MQ=255)                                      aCAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGtgctgtgctg                                  <  2:860801/76‑1 (MQ=255)                                      aCAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGtgctgtgctg                                  <  1:1947998/76‑1 (MQ=255)                                       cAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTgg                                 <  1:1719511/76‑1 (MQ=255)                                       cAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTgg                                 >  2:1481228/1‑76 (MQ=255)                                        aGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGt                                <  1:1785079/76‑1 (MQ=255)                                        aGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGt                                >  1:852962/1‑76 (MQ=255)                                         ggCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTa                               <  2:536175/76‑1 (MQ=255)                                          gCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAg                              >  1:1530444/1‑76 (MQ=255)                                             aagaaagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCg                           >  1:1469107/1‑76 (MQ=255)                                              agaaagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCgg                          >  1:1426932/1‑76 (MQ=255)                                               gaaagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGt                         >  2:967181/1‑76 (MQ=255)                                               gaaagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGt                         <  2:1071957/76‑1 (MQ=255)                                               gaaagaCCTGTCGGAAGAACTGCGGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGt                         <  2:1978318/76‑1 (MQ=255)                                                 aagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTgg                       >  2:1788169/1‑76 (MQ=255)                                                 aagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCAGGCTGTGCTGGTAGCCGGTgg                       >  2:1864983/1‑76 (MQ=255)                                                  agcccTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGc                      <  2:1488566/73‑1 (MQ=255)                                                   gccctgcctgAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCg                     <  2:1506214/74‑1 (MQ=255)                                                   gaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCg                     <  2:1113335/76‑1 (MQ=255)                                                    aCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGt                    >  2:878072/1‑76 (MQ=255)                                                      cTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTg                  >  2:42234/1‑76 (MQ=255)                                                        gTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGaa                >  1:634903/1‑76 (MQ=255)                                                           tGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCTGTGCTGTGCTGGTGGCCGGTGGCGATGAAGCt             >  1:1664407/1‑76 (MQ=255)                                                             aagaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTga           >  2:922255/1‑76 (MQ=255)                                                             aagaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTga           >  1:1633711/1‑76 (MQ=255)                                                              agaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTCGTGTGCTGGTAGCCGGTGGCGTTGAAGCTgag          >  2:851633/1‑76 (MQ=255)                                                               gaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTgaga         >  1:1301230/1‑76 (MQ=255)                                                                aaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGaa        >  2:601081/1‑76 (MQ=255)                                                                aaCGGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGaa        <  1:1094745/76‑1 (MQ=255)                                                                  cTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGc      >  1:1306534/1‑76 (MQ=255)                                                                   tGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCt     >  2:904568/1‑76 (MQ=255)                                                                    gCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTc    >  1:151628/1‑76 (MQ=255)                                                                      agagCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGa  >  2:827534/1‑76 (MQ=255)                                                                      agagCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGCAGCCGGTGGCGTTGCAGCCGAGAAGCTGGa  >  2:339160/1‑76 (MQ=255)                                                                         |                                                                         ACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGA  >  REL606/4163680‑4163824

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