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%	279.8 / NA	67	I572L (ATC→CTC)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (38/29);  total (38/29)

GAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCG  >  REL606/4162498‑4162641                                                                        |                                                                         gagaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCCGAAGGTCCGAACATCCGTCTGCTCaa                                                                      >  2:1384837/1‑76 (MQ=255)  agaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAAc                                                                     >  2:1199253/1‑76 (MQ=255)    aCGTACAACCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAActc                                                                   >  1:1276918/1‑76 (MQ=255)     cGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAActct                                                                  <  2:294330/76‑1 (MQ=255)      gTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAActctc                                                                 >  1:1232863/1‑76 (MQ=255)       tccacCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAActctct                                                                <  1:1926232/74‑1 (MQ=255)        acacCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTg                                                               >  2:1677512/1‑76 (MQ=255)         cacCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGt                                                              >  2:1582808/1‑76 (MQ=255)         cacCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGt                                                              <  1:1783902/76‑1 (MQ=255)           cccGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTcc                                                            <  1:1880539/76‑1 (MQ=255)             cGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCgt                                                          >  2:644616/1‑76 (MQ=255)              gACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCgtg                                                         >  2:1430904/1‑76 (MQ=255)                  cACTACTGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACg                                                     <  2:1363401/76‑1 (MQ=255)                  cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACg                                                     <  1:341809/76‑1 (MQ=255)                  cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACg                                                     <  2:1104419/76‑1 (MQ=255)                   aCTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGc                                                    >  1:1728629/1‑76 (MQ=255)                    cTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGca                                                   <  2:390736/76‑1 (MQ=255)                    cTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGca                                                   <  2:110938/76‑1 (MQ=255)                     tACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGcac                                                  >  2:1848719/1‑76 (MQ=255)                        ggTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGa                                               <  1:533548/76‑1 (MQ=255)                        ggTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGa                                               >  1:1058645/1‑76 (MQ=255)                        ggTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGa                                               >  2:420115/1‑76 (MQ=255)                        ggTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGa                                               >  2:70127/1‑76 (MQ=255)                        ggTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGa                                               >  2:754874/1‑76 (MQ=255)                        ggTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGa                                               <  1:1767857/76‑1 (MQ=255)                        ggTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGa                                               <  1:494347/76‑1 (MQ=255)                         gTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGAc                                              <  2:1871835/76‑1 (MQ=255)                           cgcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTa                                            >  1:945909/1‑76 (MQ=255)                            gcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTaa                                           <  2:118818/76‑1 (MQ=255)                             ngTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAAc                                          >  1:1512721/2‑76 (MQ=255)                             cgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCTGTGTACGCACAGACTAAc                                          <  1:572797/76‑1 (MQ=255)                             cgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAAc                                          >  2:864124/1‑76 (MQ=255)                                aTGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGaa                                       >  2:675863/1‑76 (MQ=255)                                aTGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGaa                                       >  2:186164/1‑76 (MQ=255)                                 tGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCGCAGACTAACGAAt                                      <  1:1306454/76‑1 (MQ=255)                                 tGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAAt                                      <  2:1674685/76‑1 (MQ=255)                                  ggccAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                                     <  2:49819/74‑1 (MQ=255)                                   tCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                                    >  2:1767214/1‑76 (MQ=255)                                    ccAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACg                                   >  2:909899/1‑76 (MQ=255)                                     cAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACgg                                  >  2:366094/1‑76 (MQ=255)                                         cGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTc                              >  2:998098/1‑76 (MQ=255)                                         cGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTc                              >  2:207189/1‑76 (MQ=255)                                          naaaCCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTcc                             >  1:1708676/2‑76 (MQ=255)                                            aaCCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                           <  1:1858373/76‑1 (MQ=255)                                            aaCCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                           >  2:1230983/1‑76 (MQ=255)                                            aaCCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                           >  1:1484334/1‑76 (MQ=255)                                             aCCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                          <  1:1741127/76‑1 (MQ=255)                                             aCCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                          >  2:931426/1‑76 (MQ=255)                                             aCCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                          <  2:1831689/76‑1 (MQ=255)                                             aCCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                          >  1:1884291/1‑76 (MQ=255)                                               cccTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTgag                        <  1:1044650/76‑1 (MQ=255)                                               cccTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTgag                        <  2:1119920/76‑1 (MQ=255)                                                 cTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGAc                      >  2:1325333/1‑76 (MQ=255)                                                   gAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTc                    >  1:577332/1‑76 (MQ=255)                                                   gAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTc                    >  1:1855281/1‑76 (MQ=255)                                                    aaGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTcc                   >  1:1105982/1‑76 (MQ=255)                                                      ggTCCGAACATCGGTCTGCTAAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGt                 <  2:1318873/76‑1 (MQ=255)                                                               aTCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAg        <  2:1963162/76‑1 (MQ=255)                                                               aTCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAg        >  2:1884466/1‑76 (MQ=255)                                                               aTCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAg        <  1:1285172/76‑1 (MQ=255)                                                                tCGGTCTGCTCAACNCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGt       >  1:1079971/1‑76 (MQ=255)                                                                 cGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTg      <  1:763038/76‑1 (MQ=255)                                                                  ggTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGa     <  1:884410/76‑1 (MQ=255)                                                                  ggTCTGCTCAACGCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGa     <  2:248929/76‑1 (MQ=255)                                                                    tCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGAcc   <  1:347060/76‑1 (MQ=255)                                                                    tCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGAcc   >  1:1832836/1‑76 (MQ=255)                                                                     cTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTTAGACTCCGTATCGTAAAGAGACCg  >  2:691842/1‑76 (MQ=255)                                                                     cTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCg  >  2:1429995/1‑76 (MQ=255)                                                                        |                                                                         GAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCG  >  REL606/4162498‑4162641

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