Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,162,470	A→C	T539P (ACC→CCC)	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,470	0	A	C	100.0%	126.4 / NA	31	T539P (ACC→CCC)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (19/12);  total (19/12)

ATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAA  >  REL606/4162398‑4162542                                                                         |                                                                         aTGGACCAGAACAACCCGCTGTCTGAGATTATGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGcccc                                                                       >  2:1154466/1‑76 (MQ=255)        aGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACg                                                                >  2:726890/1‑76 (MQ=255)            caacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCa                                                            <  1:97591/76‑1 (MQ=255)                  ncTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAATGTGCAGGCTTc                                                      >  1:1157050/2‑76 (MQ=255)                  gCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTc                                                      >  2:1061993/1‑76 (MQ=255)                      tCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAg                                                  <  2:980866/76‑1 (MQ=255)                         gagaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTc                                               >  1:1087120/1‑76 (MQ=255)                          agaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCg                                              <  1:1073248/76‑1 (MQ=255)                          agaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCg                                              >  2:799470/1‑76 (MQ=255)                              tACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGAc                                          >  1:1044791/1‑76 (MQ=255)                               aCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACg                                         <  2:339076/76‑1 (MQ=255)                               aCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACg                                         >  2:233700/1‑76 (MQ=255)                                   acaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTaca                                     >  2:926837/1‑76 (MQ=255)                                      aaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACAccc                                  >  2:526215/1‑76 (MQ=255)                                       aCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCg                                 <  1:634055/76‑1 (MQ=255)                                       aCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCg                                 <  2:1138401/76‑1 (MQ=255)                                        cgtcgcATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGa                                >  1:106755/1‑76 (MQ=255)                                          tcgtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACt                              <  1:1097951/76‑1 (MQ=255)                                           cgtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTc                             >  1:187002/1‑76 (MQ=255)                                              aTCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACt                          >  2:579395/1‑76 (MQ=255)                                                    gCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTc                    >  2:584844/1‑76 (MQ=255)                                                    gCACCCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTc                    >  2:886218/1‑76 (MQ=255)                                                       cTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcgcg                 >  2:1010390/1‑76 (MQ=255)                                                         cGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTa               >  1:1030149/1‑76 (MQ=255)                                                               aGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCa         <  2:1081175/76‑1 (MQ=255)                                                                ggCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCaa        <  1:83016/76‑1 (MQ=255)                                                                   ggTGTGCCCCGGGAACAGGGAGGCCGCGAAGGGCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAGTCg     <  1:1061911/76‑1 (MQ=255)                                                                   ggTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCg     >  2:538303/1‑76 (MQ=255)                                                                   ggTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCg     >  2:376009/1‑76 (MQ=255)                                                                   ggTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCg     <  1:1061918/76‑1 (MQ=255)                                                                   ggTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCg     <  1:440728/76‑1 (MQ=255)                                                                     tCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaa   >  2:516080/1‑76 (MQ=255)                                                                      cTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaaa  <  1:1061004/76‑1 (MQ=255)                                                                         |                                                                         ATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAA  >  REL606/4162398‑4162542

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