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%	216.6 / NA	53	T539P (ACC→CCC)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (26/27);  total (26/27)

ATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAA  >  REL606/4162398‑4162542                                                                         |                                                                         aTGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGcccc                                                                       <  1:194944/76‑1 (MQ=255) aTGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGcccc                                                                       <  1:945000/76‑1 (MQ=255)    gACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTg                                                                    <  2:1560452/76‑1 (MQ=255)     aCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGa                                                                   >  2:693731/1‑76 (MQ=255)       cAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAAc                                                                 >  1:226294/1‑76 (MQ=255)         gAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGt                                                               <  1:210325/76‑1 (MQ=255)          aacaacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTg                                                              <  2:538989/76‑1 (MQ=255)           acaacCCGCTGTCTTAGATTACGCACAAACGTCGTATCTCCGCACGCGGCCCAGGCGGTCTGCCCCGTGAACGTGc                                                             <  2:268524/76‑1 (MQ=255)             aacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAg                                                           >  2:1534316/1‑76 (MQ=255)              acCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAgg                                                          >  2:1865317/1‑76 (MQ=255)                ccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCt                                                        >  2:1842445/1‑76 (MQ=255)                  gCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTc                                                      <  2:1714436/76‑1 (MQ=255)                     gTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGaa                                                   >  2:1038467/1‑76 (MQ=255)                        tGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGtt                                                >  2:645368/1‑76 (MQ=255)                         gagaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTc                                               <  1:314259/76‑1 (MQ=255)                          agaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCg                                              <  1:1070085/76‑1 (MQ=255)                          agaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCg                                              <  2:2017718/76‑1 (MQ=255)                          agaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCg                                              >  2:814400/1‑76 (MQ=255)                          agaTTACCCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCg                                              >  2:933054/1‑76 (MQ=255)                             ttACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCgaga                                           >  2:1091134/1‑76 (MQ=255)                              tACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGAc                                          >  2:2063370/1‑76 (MQ=255)                                 gCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTa                                       <  2:712180/76‑1 (MQ=255)                                  cacaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTac                                      <  1:203817/76‑1 (MQ=255)                                   acaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTaca                                     >  2:1978294/1‑76 (MQ=255)                                   acaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTaca                                     <  2:667340/76‑1 (MQ=255)                                    caAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTacac                                    <  2:384917/76‑1 (MQ=255)                                     aaaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACAcc                                   >  1:773046/1‑76 (MQ=255)                                     aaaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACAcc                                   <  1:2006609/76‑1 (MQ=255)                                     aaaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACAcc                                   >  2:1713874/1‑76 (MQ=255)                                      aaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACAccc                                  >  2:95568/1‑76 (MQ=255)                                        cgtcgtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGCGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGa                                >  2:1022865/1‑76 (MQ=255)                                            gtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCa                            <  2:1159415/76‑1 (MQ=255)                                             tATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCAc                           >  2:610081/1‑76 (MQ=255)                                              aTCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGCACGAGACGTACACCCGACTCACt                          >  2:1462331/1‑76 (MQ=255)                                              aTCCCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACt                          <  1:1297717/76‑1 (MQ=255)                                                ctcCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTAc                        <  1:961186/76‑1 (MQ=255)                                                ctcCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTAc                        >  1:957692/1‑76 (MQ=255)                                                 tcCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACg                       <  2:1651274/76‑1 (MQ=255)                                                 tcCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACg                       >  1:246092/1‑76 (MQ=255)                                                  ccGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACgg                      >  2:679193/1‑76 (MQ=255)                                                   cGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACggc                     >  1:646800/1‑75 (MQ=255)                                                    gcgcTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCGGGCTTCGAAGTTCGAGACGGACACCCGACTCACTACGGTc                    <  1:1529891/73‑1 (MQ=255)                                                     cACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcg                   <  1:775395/76‑1 (MQ=255)                                                      aCTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcgc                  >  2:541313/1‑76 (MQ=255)                                                        tCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGt                >  2:995769/1‑76 (MQ=255)                                                        tCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGt                >  1:1108577/1‑76 (MQ=255)                                                            cccAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGt            <  1:2182/76‑1 (MQ=255)                                                            cccAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGt            <  1:1460229/76‑1 (MQ=255)                                                               aGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGCCTCACTACGGTCGCGGATGTCCa         >  2:1599783/1‑76 (MQ=255)                                                               aGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCa         >  2:16262/1‑76 (MQ=255)                                                               aGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCa         <  2:1463057/76‑1 (MQ=255)                                                                 gCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAAt       <  2:964620/76‑1 (MQ=255)                                                                  cGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATc      <  1:1780454/76‑1 (MQ=255)                                                                   ggTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCg     <  2:169498/76‑1 (MQ=255)                                                                     tCTGCCCCGTGAACGTGCAGGCTTCGAAGTCCGAGACGTACTCCCGACTCACTACGGTCGCGTATGTCCAATCGaa   <  1:1501049/76‑1 (MQ=255)                                                                      cTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaaa  <  1:1953107/76‑1 (MQ=255)                                                                         |                                                                         ATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAA  >  REL606/4162398‑4162542

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