Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,162,570	T→A	I572N (ATC→AAC)	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,570	0	T	A	100.0%	176.0 / NA	45	I572N (ATC→AAC)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (21/24);  total (21/24)

CGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGA  >  REL606/4162497‑4162638                                                                          |                                                                     cGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACa                                                                    >  2:133740/1‑76 (MQ=255) cGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACa                                                                    <  1:1492413/76‑1 (MQ=255)      cGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAActct                                                               >  2:1098073/1‑76 (MQ=255)        tACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAActctct                                                             <  2:758803/76‑1 (MQ=255)        tACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAActctct                                                             >  1:657082/1‑76 (MQ=255)         acacCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTg                                                            >  2:599186/1‑76 (MQ=255)              cGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCgt                                                       <  1:1629425/76‑1 (MQ=255)              cGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTAAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCgt                                                       <  2:888466/76‑1 (MQ=255)               gACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCgtg                                                      >  1:914343/1‑76 (MQ=255)                aCTCCCTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCgtgt                                                     <  2:333629/76‑1 (MQ=255)                  tCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTAc                                                   >  1:83120/1‑76 (MQ=255)                   cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACg                                                  <  1:1112723/76‑1 (MQ=255)                     cTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGca                                                >  2:350465/1‑76 (MQ=255)                     cTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGca                                                >  2:1546843/1‑76 (MQ=255)                       aCGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGcaca                                              <  2:237687/76‑1 (MQ=255)                       aCGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGACCGTGTACGcaca                                              >  1:381183/1‑76 (MQ=255)                        cGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAg                                             <  1:1016775/76‑1 (MQ=255)                        cGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAg                                             <  1:644611/76‑1 (MQ=255)                           tCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACt                                          <  2:479243/76‑1 (MQ=255)                            cgcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTa                                         >  1:1383383/1‑76 (MQ=255)                               gTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACg                                      <  1:360899/76‑1 (MQ=255)                                tATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGa                                     >  1:292599/1‑76 (MQ=255)                                 ntGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGaa                                    >  1:231681/2‑76 (MQ=255)                                 aTGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGaa                                    >  1:1403344/1‑76 (MQ=255)                                 aTGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGaa                                    <  2:445803/76‑1 (MQ=255)                                      cAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACgg                               <  2:1372025/76‑1 (MQ=255)                                       aaTCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGc                              >  1:1460121/1‑76 (MQ=255)                                            aaaCCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt                         >  2:1058178/1‑76 (MQ=255)                                            aaaCCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt                         <  1:124130/76‑1 (MQ=255)                                                  cTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGAc                   >  1:1429417/1‑76 (MQ=255)                                                  cTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGAc                   >  1:1080294/1‑76 (MQ=255)                                                    gAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTc                 <  1:1009307/76‑1 (MQ=255)                                                     aaGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGCGACTcc                >  2:1265081/1‑76 (MQ=255)                                                     aaGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTcc                >  2:1265075/1‑76 (MQ=255)                                                          ccGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATc           >  2:1710387/1‑76 (MQ=255)                                                          ccGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATc           >  1:369866/1‑76 (MQ=255)                                                           cGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCg          <  1:1090885/76‑1 (MQ=255)                                                            gAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGt         <  1:1303583/76‑1 (MQ=255)                                                             aaCATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTa        <  2:804183/76‑1 (MQ=255)                                                             aaCATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTa        >  2:161625/1‑76 (MQ=255)                                                             aaCATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTa        <  1:324284/76‑1 (MQ=255)                                                             aaCAGCGGTCTGAACAACTCTCCGTCCGTGTACGCACAGACTAACGAACACGGATTACTTGAGACTCCGTATCGTa        <  2:1307833/76‑1 (MQ=255)                                                              aCATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaa       <  1:892497/76‑1 (MQ=255)                                                              aCATCGGTCTGAACAACGCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaa       <  1:892479/76‑1 (MQ=255)                                                                aTCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAg     <  1:1365785/76‑1 (MQ=255)                                                                 tCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGt    >  2:301860/1‑76 (MQ=255)                                                                 tCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGt    <  1:1189932/76‑1 (MQ=255)                                                                   ggTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGa  <  1:1573722/76‑1 (MQ=255)                                                                          |                                                                     CGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGA  >  REL606/4162497‑4162638

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