Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,162,521	G→A	G556S (GGT→AGT)	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,162,521	0	G	A	100.0%	157.0 / NA	39	G556S (GGT→AGT)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/12);  total (27/12)

CGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGC  >  REL606/4162448‑4162591                                                                          |                                                                       cGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGt                                                                      <  2:1291551/76‑1 (MQ=255)   cACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTcg                                                                    <  1:1401853/76‑1 (MQ=255)   cACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTcg                                                                    <  1:782746/76‑1 (MQ=255)          cccAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGt                                                             <  1:816288/76‑1 (MQ=255)            cAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTcc                                                           >  1:1573833/1‑76 (MQ=255)            cAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTcc                                                           >  2:1059950/1‑76 (MQ=255)             aGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCa                                                          >  2:1026462/1‑76 (MQ=255)                cGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCAATc                                                       <  1:89055/76‑1 (MQ=255)                cGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCAATc                                                       >  1:1679399/1‑76 (MQ=255)                 ggTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCAATCg                                                      >  1:1366975/1‑76 (MQ=255)                 ggTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCAATCg                                                      >  1:737687/1‑76 (MQ=255)                     tGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAAc                                                  >  2:998575/1‑76 (MQ=255)                        cccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAAcccc                                               <  1:12991/76‑1 (MQ=255)                              aaCGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAgg                                         >  2:672169/1‑76 (MQ=255)                                    cAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGaa                                   >  1:549154/1‑76 (MQ=255)                                     aGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAAc                                  >  2:1246677/1‑76 (MQ=255)                                     aGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAAc                                  >  2:456253/1‑76 (MQ=255)                                     aGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAAc                                  >  1:870586/1‑76 (MQ=255)                                       gCTTCGAAGTTCGAGACGTACCCCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACAt                                <  1:1566537/76‑1 (MQ=255)                                        cTTCGAAGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCCGACGGCCCGAACATc                               >  1:870096/1‑76 (MQ=255)                                         ttCGAAGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCg                              >  1:851217/1‑76 (MQ=255)                                              aGTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTg                         <  1:1292902/76‑1 (MQ=255)                                               gTTCGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGa                        >  2:815350/1‑76 (MQ=255)                                                  ngagaCGTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCa                     >  1:1178250/2‑76 (MQ=255)                                                  cGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCa                     >  2:822600/1‑76 (MQ=255)                                                  cGAGACGTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCa                     >  1:1638542/1‑76 (MQ=255)                                                   gagaCGTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCaa                    >  2:991815/1‑76 (MQ=255)                                                   gagaCGTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCaa                    <  1:207859/76‑1 (MQ=255)                                                       cGTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActct                <  2:277794/76‑1 (MQ=255)                                                        gTACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActctc               <  2:732986/76‑1 (MQ=255)                                                         tACACCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActctct              >  2:391268/1‑76 (MQ=255)                                                          acacCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTg             >  2:611544/1‑76 (MQ=255)                                                          acacCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTg             >  2:1237747/1‑76 (MQ=255)                                                            acCCGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTc           >  1:821930/1‑76 (MQ=255)                                                              ccGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCg         <  2:1234033/76‑1 (MQ=255)                                                               cGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgt        >  1:666870/1‑76 (MQ=255)                                                               cGACTCACTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgt        >  1:750264/1‑76 (MQ=255)                                                                     aCTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGc  >  1:1723610/1‑76 (MQ=255)                                                                     aCTACAGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGc  >  1:375908/1‑76 (MQ=255)                                                                          |                                                                       CGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGC  >  REL606/4162448‑4162591

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: ‑