Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,162,569	A→C	I572L (ATC→CTC)	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,569	0	A	C	100.0%	150.8 / NA	37	I572L (ATC→CTC)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (18/19);  total (18/19)

GAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTA  >  REL606/4162498‑4162632                                                                        |                                                                gagaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCaa                                                             >  2:973968/1‑76 (MQ=255)   gaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAAct                                                           >  2:819698/1‑76 (MQ=255)     cGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAActct                                                         >  1:462548/1‑76 (MQ=255)       tACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAActctct                                                       <  1:533235/76‑1 (MQ=255)           cccGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTcc                                                   >  2:486498/1‑76 (MQ=255)               acgcACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCgtgt                                               <  1:719493/73‑1 (MQ=255)               aCTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCgtgt                                               >  1:888833/1‑76 (MQ=255)                cTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTa                                              <  2:1402482/76‑1 (MQ=255)                 tCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTAc                                             <  2:918090/76‑1 (MQ=255)                  cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACg                                            >  1:757113/1‑76 (MQ=255)                    cTACGGTCGCGTATGTCCAATCGAACCCCCTGAAGGTCCGAACATCGGTCTGCTCACCTCTCTGTCCGCGTACGca                                          >  1:535460/1‑76 (MQ=255)                    cTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGca                                          >  2:1132505/1‑76 (MQ=255)                    cTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGca                                          >  1:535447/1‑76 (MQ=255)                          tCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACt                                    <  1:1392390/76‑1 (MQ=255)                          tCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACt                                    <  1:522286/76‑1 (MQ=255)                          tCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACt                                    <  1:544398/76‑1 (MQ=255)                          tCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACt                                    <  1:60220/76‑1 (MQ=255)                           cgcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTa                                   <  1:684739/76‑1 (MQ=255)                           cgcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTa                                   <  1:35225/76‑1 (MQ=255)                             cgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAAc                                 <  1:25717/76‑1 (MQ=255)                                aTGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGaa                              >  1:221993/1‑76 (MQ=255)                                aTGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGaa                              >  1:1386634/1‑76 (MQ=255)                                  gTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                            <  1:65774/76‑1 (MQ=255)                                  gTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                            >  1:345207/1‑76 (MQ=255)                                    ccAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACg                          <  1:192855/76‑1 (MQ=255)                                        tCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCtt                      <  2:923470/76‑1 (MQ=255)                                           aaaCCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTctt                   >  2:587620/1‑74 (MQ=255)                                           aaaCCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt                   <  2:1516204/76‑1 (MQ=255)                                            aaCCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                  >  2:1346425/1‑76 (MQ=255)                                               cccTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTgag               <  1:1066761/76‑1 (MQ=255)                                               cccTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTgag               <  1:292062/76‑1 (MQ=255)                                                 cTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGAc             >  2:370102/1‑76 (MQ=255)                                                  tGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACt            <  2:739381/76‑1 (MQ=255)                                                       gTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTa       >  1:1175147/1‑76 (MQ=255)                                                           gAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGt   <  1:879085/76‑1 (MQ=255)                                                            aaCATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTa  >  2:853285/1‑76 (MQ=255)                                                            aaCATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTa  >  2:347084/1‑76 (MQ=255)                                                                        |                                                                GAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTA  >  REL606/4162498‑4162632

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