Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,162,569	A→T	I572F (ATC→TTC)	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	T	100.0%	245.4 / NA	60	I572F (ATC→TTC)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (32/28);  total (32/28)

CGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACC  >  REL606/4162497‑4162640                                                                         |                                                                        cGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTc                                                                       <  2:555004/75‑1 (MQ=255)     aCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAAct                                                                   <  2:1307470/75‑1 (MQ=255)      cGTACCCCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAActc                                                                  <  2:1019745/75‑1 (MQ=255)       gTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAActct                                                                 <  1:915132/75‑1 (MQ=255)        tACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAActctc                                                                <  2:1602187/75‑1 (MQ=255)           acCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGt                                                             >  1:39839/1‑75 (MQ=255)           acCCGACTCACTACGGTCGCGTATGGCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGt                                                             >  1:39841/1‑75 (MQ=255)            cccGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTc                                                            >  1:1088548/1‑75 (MQ=255)              cGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCg                                                          <  2:1259586/75‑1 (MQ=255)                  tCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTa                                                      >  1:835481/1‑75 (MQ=255)                  tCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTa                                                      <  1:1467825/75‑1 (MQ=255)                     cTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGc                                                   >  2:131885/1‑75 (MQ=255)                     cTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGc                                                   >  1:922603/1‑75 (MQ=255)                        cGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGcaca                                                <  2:112979/75‑1 (MQ=255)                          gTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGa                                              >  1:938186/1‑75 (MQ=255)                           tCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGAc                                             <  2:1507922/75‑1 (MQ=255)                            cgcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACt                                            <  2:147271/75‑1 (MQ=255)                              cgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTaa                                          <  1:318320/75‑1 (MQ=255)                              cgTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTaa                                          <  1:1399250/75‑1 (MQ=255)                               gTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAAc                                         >  1:862400/1‑75 (MQ=255)                                 aTGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGa                                       >  2:179017/1‑75 (MQ=255)                                  tGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGaa                                      >  1:195869/1‑75 (MQ=255)                                   gTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAAt                                     <  2:147415/75‑1 (MQ=255)                                    tCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                                    <  1:1253176/75‑1 (MQ=255)                                     ccAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                                   >  2:861687/1‑75 (MQ=255)                                      cAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACg                                  <  2:1313192/75‑1 (MQ=255)                                      cAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACg                                  >  1:1226034/1‑75 (MQ=255)                                       aaTCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACgg                                 >  2:932930/1‑75 (MQ=255)                                        aTCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGc                                <  2:1435375/75‑1 (MQ=255)                                         tCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCt                               <  1:776566/75‑1 (MQ=255)                                           gAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTc                             >  2:1320886/1‑75 (MQ=255)                                            aaaCCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTcc                            <  1:645062/75‑1 (MQ=255)                                             aaCCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt                           >  2:956162/1‑75 (MQ=255)                                             aaCCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt                           >  1:1242251/1‑75 (MQ=255)                                             aaCCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt                           >  2:1198035/1‑75 (MQ=255)                                              aCCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                          <  1:426208/75‑1 (MQ=255)                                              aCCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                          >  1:210372/1‑75 (MQ=255)                                                cccTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTga                        >  2:361771/1‑75 (MQ=255)                                                cccTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTga                        >  1:430524/1‑75 (MQ=255)                                                  cTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTgaga                      >  1:1573367/1‑75 (MQ=255)                                                   tGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGAc                     <  2:1300074/75‑1 (MQ=255)                                                   tGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGAc                     <  2:1483524/75‑1 (MQ=255)                                                    gAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACt                    >  1:1615746/1‑75 (MQ=255)                                                    gAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACt                    <  1:840450/75‑1 (MQ=255)                                                     aaGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACGGACTAACGAATACGGCTTCCTTGAGACTc                   >  1:519647/1‑75 (MQ=255)                                                     aaGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTc                   >  1:923594/1‑75 (MQ=255)                                                     aaGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTc                   >  2:1514451/1‑75 (MQ=255)                                                      aGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTcc                  >  1:460117/1‑75 (MQ=255)                                                      aGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTcc                  >  2:208679/1‑75 (MQ=255)                                                      aGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTcc                  >  1:1477732/1‑75 (MQ=255)                                                       ggTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCg                 <  1:426506/75‑1 (MQ=255)                                                        gTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGt                <  1:384599/75‑1 (MQ=255)                                                         tCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTa               <  2:110532/75‑1 (MQ=255)                                                             aaCATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGt           <  1:861687/75‑1 (MQ=255)                                                              aCATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTa          <  2:1088548/75‑1 (MQ=255)                                                              aCATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTa          <  1:752889/75‑1 (MQ=255)                                                                 tCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAg       <  1:592682/75‑1 (MQ=255)                                                                   ggTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAggg     >  1:152993/1‑73 (MQ=255)                                                                   ggTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTg     >  1:815962/1‑75 (MQ=255)                                                                    gTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGa    >  2:950989/1‑75 (MQ=255)                                                                     tCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGAc   >  1:924050/1‑75 (MQ=255)                                                                      cTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGAcc  >  1:1446373/1‑75 (MQ=255)                                                                         |                                                                        CGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACC  >  REL606/4162497‑4162640

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