Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,164,563	C→A	N1236K (AAC→AAA)	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,164,563	0	C	A	100.0%	113.6 / NA	31	N1236K (AAC→AAA)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (15/16);  total (15/16)

TGTACGATGGTCGCACTGGTGAACAGTTCGAGCGTCCGGTAACCGTTGGTTACATGTACATGCTGAAACTGAACCACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTG  >  REL606/4164490‑4164632                                                                          |                                                                      tGTACGATGGTCGCACTGGTGAACAGTTCGAGCGTCCGGTAACCGTTGGTTACATGTACATGCTGAAACTGAaaca                                                                     <  1:576775/76‑1 (MQ=255)   tACGATGGTCGCACTGGTGAACAGTTCGAGCGTCCGGTAACCGTTGGTTACATGTACATGCTGAAACTGAaacacc                                                                   >  2:565222/1‑76 (MQ=255)   tACGATGGTCGCACTGGTGAACAGTTCGAGCGTCCGGTAACCGTTGGTTACATGTACATGCTGAAACTGAaacacc                                                                   >  2:1099522/1‑76 (MQ=255)    aCGATGGTCGCACTGGTGAACAGTTCGAGCGTCCGGTAACCGTTGGTTACATGTACATGCTGAAACTGAAACACCt                                                                  <  1:669244/76‑1 (MQ=255)     cGATGGTCGCACTGGTGAACAGTTCGAGCGTCCGGTAACCGTTGGTTACATGTACATGCTGAAACTGAAACACCTg                                                                 <  1:295638/76‑1 (MQ=255)     cGATGGTCGCACTGGTGAACAGTTCGAGCGTCCGGTAACCGTTGGTTACATGTACATGCTGAAACTGAAACACCTg                                                                 <  1:650399/76‑1 (MQ=255)       aTGGTCGCACTGGTGAACAGTTCGAGCGTCCGGTAACCGTTGGTTACATGTACATGCTGAAACTGAAACACCTGGt                                                               >  1:712736/1‑76 (MQ=255)         ggTCGCAGTGGTGAACAGTTCGAGCGTCCGGTAACCGTTGGTTACATGTACATGCTGAAACTGAAACACCTGGTcg                                                             >  1:721467/1‑76 (MQ=255)         ggTCGCACTGGTGAACAGTTCGAGCGTCCGGTAACCGTTGGTTACATGTACATGCTGAAACTGAAACACCTGGTcg                                                             <  1:847002/76‑1 (MQ=255)             gCACTGGTGAACAGTTCGAGCGTCCGGTAACCGTTGGTTACATGTACATGCTGAAACTGAAACACCTGGTcgacga                                                         >  2:288621/1‑76 (MQ=255)                cTGGTGAACAGTTCGAGCGTCCGGTAACCCTTGGTTACATGTACATGCTGAAACTGAAACACCTGGTCGACGACaa                                                      <  1:1044676/76‑1 (MQ=255)                  ggTGAACAGTTCGAGCGTCCGGTAACCGTTGGTTACATGTACATGCTGAAACTGAAACACCTGGTCGACGACAAGa                                                    <  2:566017/76‑1 (MQ=255)                     gAACAGTTCGAGCGTCCGGTAACCGTTGGTTACATGTACATGCTGAAACTGAAACACCTGGTCGACGACAAGATGc                                                 >  2:742422/1‑76 (MQ=255)                          gTTCGAGCGTCCGGTAACCGTTGGTTACATGTACATGCTGAAACTGAAACACCTGGTCGACGACAAGATGCAcgcg                                            <  1:502093/76‑1 (MQ=255)                                 cGTCCGGTAACCGTTGGTTACATGTACATGCTGAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCGTTCCa                                     <  2:632790/76‑1 (MQ=255)                                  gttcGGTAACCGTTGGTTACATGTACATGCTGAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCGTTCCAc                                    <  1:537117/73‑1 (MQ=255)                                              ttGGTTACATGTACATGCTGAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAg                        <  1:157211/76‑1 (MQ=255)                                               tGGTTACATGTACATGCTGAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGc                       >  1:666055/1‑76 (MQ=255)                                                 nttACATGTACATGCTGAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCt                     <  1:979994/75‑1 (MQ=255)                                                  ttACATGTACATGCTGAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTg                    >  1:326746/1‑76 (MQ=255)                                                  tgacatgtacatgCTGAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCCTTCCACCGGTTCTTACAGCCTg                    >  1:835620/3‑76 (MQ=255)                                                   tacatgtacatgCTGAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTgg                   >  1:951660/1‑76 (MQ=255)                                                    acatgtacatgCTGAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGt                  >  1:587652/1‑76 (MQ=255)                                                     catgtacatgCTGAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGtt                 >  2:1139624/1‑76 (MQ=255)                                                     catgtacatgCTGAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGtt                 >  1:314298/1‑76 (MQ=255)                                                       tgtacatgCTGAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTAc               >  1:47897/1‑76 (MQ=255)                                                         tacatgCTGAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTc             <  2:431195/76‑1 (MQ=255)                                                           catgCTGAAACTGAAACATCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTcag           <  1:883280/76‑1 (MQ=255)                                                            atgCTGAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTANTcagc          >  2:929166/1‑76 (MQ=255)                                                             tgCTGAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTcagca         <  2:326746/76‑1 (MQ=255)                                                              gCTGAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTATTACAGCCTGGTTACTcagcag        >  1:798315/1‑76 (MQ=255)                                                               ctgaaactgaaaCACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCagcagc       >  1:2578/1‑76 (MQ=255)                                                                 gAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCg     <  1:716763/76‑1 (MQ=255)                                                                 gAAACTGAAACACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTCACTACAGCCTGGTTACTCAGCAGCCg     >  2:847376/1‑76 (MQ=255)                                                                    acggAAACACCTGGTCGACGTNAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTg  <  2:969150/73‑1 (MQ=255)                                                                          |                                                                      TGTACGATGGTCGCACTGGTGAACAGTTCGAGCGTCCGGTAACCGTTGGTTACATGTACATGCTGAAACTGAACCACCTGGTCGACGACAAGATGCACGCGCGTTCCACCGGTTCTTACAGCCTGGTTACTCAGCAGCCGCTG  >  REL606/4164490‑4164632

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