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%	227.4 / 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 (31/22);  total (31/22)

TGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAA  >  REL606/4162399‑4162542                                                                        |                                                                         tGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCg                                                                      <  1:1231590/76‑1 (MQ=255) tGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCg                                                                      >  1:305941/1‑76 (MQ=255)  ggACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGt                                                                     >  2:188028/1‑76 (MQ=255)   gACGAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTg                                                                    <  1:92347/76‑1 (MQ=255)   gACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTg                                                                    >  1:95190/1‑76 (MQ=255)   gACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTg                                                                    <  1:1776148/76‑1 (MQ=255)      cAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAAc                                                                 >  2:482078/1‑76 (MQ=255)       aGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACg                                                                >  1:277705/1‑76 (MQ=255)       aGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCCTGAACg                                                                <  2:693200/76‑1 (MQ=255)         aacaacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTg                                                              <  2:581938/76‑1 (MQ=255)              cccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGc                                                         <  1:1464763/76‑1 (MQ=255)               ccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCt                                                        >  2:547319/1‑76 (MQ=255)                cGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCtt                                                       <  2:1223086/76‑1 (MQ=255)                  cTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCg                                                     >  2:1120717/1‑76 (MQ=255)                   tGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGa                                                    >  1:1122867/1‑76 (MQ=255)                    gTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGaa                                                   >  1:1457497/1‑76 (MQ=255)                    gTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGaa                                                   >  2:763786/1‑76 (MQ=255)                    gTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGaa                                                   >  2:782258/1‑76 (MQ=255)                        gagaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTc                                               >  2:1633100/1‑76 (MQ=255)                             tACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCgaccc                                          >  2:1162164/1‑73 (MQ=255)                             tACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGAc                                          >  2:197706/1‑76 (MQ=255)                                gCACAAGCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTa                                       <  1:1663879/76‑1 (MQ=255)                                gCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTa                                       <  1:1738754/76‑1 (MQ=255)                                gCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTa                                       <  1:874676/76‑1 (MQ=255)                                 cacaAACGTCGTATCTCCGCACTCGGCCCAGTCGGACTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTaa                                      >  1:348592/1‑75 (MQ=255)                                 cacaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTac                                      >  2:1674906/1‑76 (MQ=255)                                   caAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTacac                                    <  2:1822905/76‑1 (MQ=255)                                       cgtcgtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGa                                <  1:1664452/76‑1 (MQ=255)                                        gtcgtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGAc                               >  2:1486525/1‑76 (MQ=255)                                          cgtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTc                             >  1:5181/1‑76 (MQ=255)                                             aTCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACt                          >  1:402351/1‑76 (MQ=255)                                             aTCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACt                          >  1:149223/1‑76 (MQ=255)                                              tctcCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTa                         >  2:1422791/1‑76 (MQ=255)                                              tctcCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTa                         <  1:122161/76‑1 (MQ=255)                                              tctcCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTa                         >  1:1308772/1‑76 (MQ=255)                                                 ccGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACgg                      >  1:62230/1‑76 (MQ=255)                                                 ccGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACgg                      >  2:577105/1‑76 (MQ=255)                                                 ccGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACgg                      >  1:1820270/1‑76 (MQ=255)                                                      cgcGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcgcg                 <  1:757562/74‑1 (MQ=255)                                                        cGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTa               <  2:326518/76‑1 (MQ=255)                                                          gCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATg             <  2:1113404/76‑1 (MQ=255)                                                           cgcAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGt            <  1:932418/74‑1 (MQ=255)                                                           cccAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGt            >  1:1379029/1‑76 (MQ=255)                                                           cccAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGt            >  1:731784/1‑76 (MQ=255)                                                            ccAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTc           >  1:485952/1‑76 (MQ=255)                                                             cAGGTGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTcc          >  2:973048/1‑76 (MQ=255)                                                             cAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTcc          <  2:518364/76‑1 (MQ=255)                                                             cAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTcc          >  2:481302/1‑76 (MQ=255)                                                             cAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTcc          >  2:125284/1‑76 (MQ=255)                                                               ggCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCaa        <  1:1381219/76‑1 (MQ=255)                                                                    tCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaa   <  1:1658241/76‑1 (MQ=255)                                                                    tCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaa   <  1:782258/76‑1 (MQ=255)                                                                    tCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTACGTCCAATCGaa   >  1:1412487/1‑76 (MQ=255)                                                                     cTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaaa  <  2:1602009/76‑1 (MQ=255)                                                                        |                                                                         TGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAA  >  REL606/4162399‑4162542

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