Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,161,306	C→T	R151C (CGT→TGT)	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,161,306	0	C	T	100.0%	113.5 / NA	31	R151C (CGT→TGT)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/16);  total (15/16)

AATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGC  >  REL606/4161233‑4161376                                                                          |                                                                       aaTTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACTGt                                                                      <  2:731088/76‑1 (MQ=255)  aTTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACTGTa                                                                     <  1:403538/76‑1 (MQ=255)     cggCTCATGACAGACAACGGTGCCTTTGTGATAAACGGTATTGAGCGTATCATCGTTTCCCAGCTGCACTGTAGTc                                                                  <  1:131363/74‑1 (MQ=255)       gCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACTGTAGTCCg                                                                <  2:529250/76‑1 (MQ=255)        cTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACTGTAGTCCgg                                                               >  2:913119/1‑76 (MQ=255)        cTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACTGTAGTCCgg                                                               >  1:627884/1‑76 (MQ=255)                agacaACGGTACCTTTGTTATCAGCGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACTGTAGTCCGGGCGtcttc                                                       <  1:1147467/76‑1 (MQ=255)                     aCGGTACCTTTGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACTGTAGTCCGGGCGTCTTCTTTGa                                                  >  1:705275/1‑76 (MQ=255)                       ggTACCTTTGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACTGTAGTCCGGGCGTCTTCTTTGACt                                                >  1:537853/1‑76 (MQ=255)                       ggTACCTTTGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACTGTAGTCCGGGCGTCTTCTTTGACt                                                <  2:72508/76‑1 (MQ=255)                         tACCTTTGTTATCAACGGTACTGAGCGTGTTATCGTATCCCAGCTGCACTGTAGTCCGGGCGTCTTCTTTGACTcc                                              >  2:739926/1‑76 (MQ=255)                            cTTTGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACTGTAGTCCGGGCGTCTTCTTAGACTCCGAc                                           >  2:21233/1‑76 (MQ=255)                             tttGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACTGTAGTCCGGGCGTCTTCTTTGACTCCGACa                                          <  1:350454/76‑1 (MQ=255)                              ttGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCATCTGCACTGTAGTCCGGGCGTCTTCTTTGACTCCGACaa                                         <  1:876036/76‑1 (MQ=255)                              ttGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACTGTAGTCCGGGCGTCTTCTTTGACTCCGACaa                                         <  2:65739/76‑1 (MQ=255)                               tGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACTGTAGTCCGGGCGTCTTCTTTGACTCCGACaaa                                        >  1:1137287/1‑76 (MQ=255)                                     cAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACTGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTaaa                                  <  2:101628/76‑1 (MQ=255)                                             cTGAGCGTGTTATCGTTTCCCAGCTGCACTGTAGTCCGGGCGTCCTCTTTGACTCCGACAAAGGTAAAACCCActc                          >  1:843430/1‑76 (MQ=255)                                              tGAGCGTGTTATCGTTTCCCAGCTGCACTGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCActct                         >  2:777074/1‑76 (MQ=255)                                                aGCGTGTTATCGTTTCCCAGCTGCACTGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTc                       >  1:1150251/1‑76 (MQ=255)                                                aGCGTGTTATCGTTTCCCAGCTGCACTGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTc                       <  2:584798/76‑1 (MQ=255)                                                    tgtTATCGTTTTCCATCTGCACTGTAGTCCTGGCGGCTTCTTNGNCTCCGACNAAGGTAAAACCCACTCTTCGGGt                   <  1:1097437/76‑1 (MQ=255)                                                      ttATCGTTTACCAGCTGCACTGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTaa                 <  1:226310/76‑1 (MQ=255)                                                       tATCGTTTCCCAGCTGCACTGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGCAAAACCCACTCTTCGGGTaaa                >  1:444369/1‑76 (MQ=255)                                                         tCGTTTCCCAGCTGCACTGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGt              >  2:574112/1‑76 (MQ=255)                                                          cGTTTCCCAGCTGCACTGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTg             <  1:733653/76‑1 (MQ=255)                                                             ttCCCAGCTGCACTGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTg          >  2:527149/1‑76 (MQ=255)                                                                  aGCTGCACTGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATaa     <  1:753323/76‑1 (MQ=255)                                                                   gCTGCACTGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAAc    <  2:705275/76‑1 (MQ=255)                                                                    cTGCACTGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAAcg   <  1:769334/76‑1 (MQ=255)                                                                    cTGCACTGTAGTCCGGGCGTCTACTTTGATTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGAATAGcg   >  1:1031445/1‑76 (MQ=255)                                                                     tGCACTGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAAcgc  >  1:504051/1‑76 (MQ=255)                                                                          |                                                                       AATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGC  >  REL606/4161233‑4161376

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