Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,167,556	G→A	E866K (GAA→AAA)	rpoC →	DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,167,556	0	G	A	100.0%	180.4 / NA	47	E866K (GAA→AAA)	rpoC	DNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/27);  total (20/27)

TCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGTTGTATCTTGTGACA  >  REL606/4167483‑4167628                                                                          |                                                                         tCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACaaa                                                                        <  1:950552/76‑1 (MQ=255)   gtgtAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACa                                                                      <  2:613453/76‑1 (MQ=255)   gtgtAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACa                                                                      <  1:1306895/76‑1 (MQ=255)   gtgtAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACa                                                                      >  1:1327337/1‑76 (MQ=255)    tgtAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAg                                                                     >  1:979603/1‑76 (MQ=255)     gtAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAgt                                                                    <  1:1428389/76‑1 (MQ=255)         ctgctgAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAgtggtg                                                                <  2:380457/76‑1 (MQ=255)           gctgAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGtgtg                                                              <  1:801680/76‑1 (MQ=255)           gctgAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAAAAGTGGtgtg                                                              <  1:92471/76‑1 (MQ=255)             tgAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGAc                                                            >  2:1478485/1‑76 (MQ=255)              gAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAGACAGTGGTGTGAct                                                           >  2:741164/1‑75 (MQ=255)              gAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGAcc                                                           >  2:1725089/1‑76 (MQ=255)               aaGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACct                                                          >  2:284459/1‑76 (MQ=255)               aaGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACct                                                          <  1:974882/76‑1 (MQ=255)                aGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACctg                                                         <  2:1075254/76‑1 (MQ=255)                     ttCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGa                                                    <  1:251773/76‑1 (MQ=255)                       cTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAag                                                  <  1:1417300/76‑1 (MQ=255)                         gAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAagag                                                >  2:1052884/1‑76 (MQ=255)                               gggTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAActct                                          <  1:1079115/76‑1 (MQ=255)                                tgTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTg                                         >  2:616300/2‑76 (MQ=255)                                  tACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTc                                       <  1:5871/76‑1 (MQ=255)                                  tACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTc                                       >  1:1364759/1‑76 (MQ=255)                                    ctgctgATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGa                                     <  1:1850885/76‑1 (MQ=255)                                    ctgctgATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGa                                     <  1:281281/76‑1 (MQ=255)                                      gctgATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGAcg                                   >  2:1453966/1‑76 (MQ=255)                                      gctgATATCCTCGTTCCGCGCAACACACTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGAcg                                   >  2:929974/1‑76 (MQ=255)                                        tgATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGAcgcg                                 <  1:1062379/76‑1 (MQ=255)                                        tgATATCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGAcgcg                                 <  2:700822/76‑1 (MQ=255)                                           tatCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGtt                              <  2:790970/76‑1 (MQ=255)                                           tatCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGtt                              <  2:1685212/76‑1 (MQ=255)                                            atCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTa                             >  1:476672/1‑76 (MQ=255)                                             tCCTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTaa                            <  1:514919/76‑1 (MQ=255)                                              ccTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTaaa                           <  1:940877/76‑1 (MQ=255)                                               cTCGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGCCGACGCGGTTAAAg                          >  2:1383529/1‑76 (MQ=255)                                                 cGTTCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTa                        <  2:1456887/76‑1 (MQ=255)                                                    tCCGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGt                     >  2:1192921/1‑76 (MQ=255)                                                     ccGCGCAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGtt                    >  1:1397210/1‑76 (MQ=255)                                                          cAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCtgtt               >  2:1146537/1‑76 (MQ=255)                                                          cAACACGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCtgtt               >  2:1002922/1‑76 (MQ=255)                                                            acacGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCtgttgt             <  1:1624212/76‑1 (MQ=255)                                                            acacGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCtgttgt             <  2:401441/76‑1 (MQ=255)                                                             cacGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGTTGTa            <  1:1079519/76‑1 (MQ=255)                                                             cacGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACGCTGTCGACGCGGTTAAAGTACGTTCTGTTGTa            <  1:1752474/76‑1 (MQ=255)                                                              acGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGTTGTAt           <  1:320494/76‑1 (MQ=255)                                                               cGCTGCTGCACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGTTGTATc          >  1:234073/1‑76 (MQ=255)                                                                  tgctgcACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGTTGTATCTtg       >  2:162849/1‑76 (MQ=255)                                                                    ctgcACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGTTGTATCTtgtg     >  1:1274733/1‑76 (MQ=255)                                                                    acacacaAACAGGTGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGTTGTATCTtgtg     <  2:1201014/73‑1 (MQ=255)                                                                     tgcACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGTTGTATCTTGTGa    >  2:262806/1‑76 (MQ=255)                                                                      gcACAAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGTTGTATCTTGTGac   <  1:77131/76‑1 (MQ=255)                                                                       cacaAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGTTGTATCTTGTGaca  <  1:342288/76‑1 (MQ=255)                                                                       cacaAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGTTGTATCTTGTGaca  <  1:693521/76‑1 (MQ=255)                                                                          |                                                                         TCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTGTCGACGCGGTTAAAGTACGTTCTGTTGTATCTTGTGACA  >  REL606/4167483‑4167628

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