Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,167,458	A→C	E833A (GAG→GCG)	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,458	0	A	C	100.0%	201.9 / NA	48	E833A (GAG→GCG)	rpoC	DNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (31/17);  total (31/17)

GTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCC  >  REL606/4167385‑4167529                                                                          |                                                                        gtTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAgcgc                                                                       <  2:724622/76‑1 (MQ=255)  ttACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGcc                                                                      >  1:629458/1‑76 (MQ=255)   tACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCg                                                                     >  1:921902/1‑76 (MQ=255)       gAAGACGATTGTGGTACNCATGAAGGTATCATGATGACTCCGGTTATCGAGGCTGGTGACGTTCAAACGCCGCTgc                                                                 >  1:226418/1‑76 (MQ=255)        aaGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTgcg                                                                >  2:1297143/1‑76 (MQ=255)          gACGATTGTGCTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAGAGCGCCGCTGcgcg                                                              >  1:152249/1‑76 (MQ=255)            cGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGAt                                                            >  1:702753/1‑76 (MQ=255)              aTTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATcg                                                          >  1:946918/1‑76 (MQ=255)              aTTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATcg                                                          <  2:1175961/76‑1 (MQ=255)               ttGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATcgc                                                         <  2:1149298/76‑1 (MQ=255)                 gtgGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGt                                                       <  1:1467801/76‑1 (MQ=255)                     tACCCATGAGGGTATCATGATGACTCCGGTTNACGAAGGTGGTGACGTTAAAGCGCCGCACCGANATCGCGTACTg                                                   >  1:218112/1‑76 (MQ=255)                     tACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTg                                                   >  1:943762/1‑76 (MQ=255)                     tACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTg                                                   >  1:354428/1‑76 (MQ=255)                      aCCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTgg                                                  <  1:425724/76‑1 (MQ=255)                      aCCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTgg                                                  >  1:591934/1‑76 (MQ=255)                         cATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTc                                               <  1:974373/76‑1 (MQ=255)                          aTGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCg                                              <  1:814331/76‑1 (MQ=255)                            gAAGGTAACATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCgtg                                            <  1:32787/76‑1 (MQ=255)                            gAAGGCATCATGATGACTCCGGTTATCGAGGGTGCTGACGTNAAAGCGCCGCTGCGCGATNGCNTACTGGGTCgtg                                            >  1:1495325/1‑76 (MQ=255)                              aGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTa                                          >  1:1633019/1‑76 (MQ=255)                                gTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAGAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAAc                                        <  1:28594/76‑1 (MQ=255)                                 tATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAAct                                       <  1:1022227/76‑1 (MQ=255)                                  aTCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAActg                                      >  2:1365298/1‑76 (MQ=255)                                    cATGATGACGCCGGTTAGCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCAATCGCGTACTGNGTNGTGTAActgct                                    >  2:279935/1‑76 (MQ=255)                                     atgatgACTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAActgctg                                   >  1:956513/1‑76 (MQ=255)                                      tgatgaCTCCGNTTNGCGAGGGTGGTTACGTTAAAGCGCCGCTGCGCGATCGCGTACACGATCGTGTAAGTGCTGa                                  <  2:281952/76‑1 (MQ=255)                                       gatgaCTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGaa                                 <  2:1037469/76‑1 (MQ=255)                                        atgaCTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAg                                <  1:497394/76‑1 (MQ=255)                                         tgaCTCCNGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAgc                               >  2:1468355/1‑75 (MQ=255)                                         tgaCTCCGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGa                               >  2:848089/1‑76 (MQ=255)                                              ccGGTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTc                          >  2:527538/1‑76 (MQ=255)                                                 gTTATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGa                       <  1:938186/76‑1 (MQ=255)                                                  ttATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGaa                      >  2:176296/1‑76 (MQ=255)                                                   tATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAg                     >  1:499377/1‑76 (MQ=255)                                                   tATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAg                     <  2:730955/76‑1 (MQ=255)                                                   tATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAg                     >  2:153463/1‑76 (MQ=255)                                                   tATCGAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACCNTCTGAAg                     <  2:1479958/76‑1 (MQ=255)                                                       gAGGGTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCgg                 >  2:907954/1‑76 (MQ=255)                                                         gggTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCgggn               <  1:397742/76‑2 (MQ=255)                                                         gggTGGTGACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGt               <  1:1358919/76‑1 (MQ=255)                                                            tggtgACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTAct            >  1:1230491/1‑76 (MQ=255)                                                             ggtgACGTTAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTActg           >  2:876077/1‑76 (MQ=255)                                                                 actttAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGa       >  2:219782/4‑76 (MQ=255)                                                                    ttAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGatat    >  2:1395989/1‑76 (MQ=255)                                                                     tAAAGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATc   >  2:516409/1‑76 (MQ=255)                                                                      aaaGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATcc  >  2:249973/1‑76 (MQ=255)                                                                      aaaGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATcc  >  2:1566714/1‑76 (MQ=255)                                                                      aaaGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATcc  >  2:1228946/1‑76 (MQ=255)                                                                      aaaGCGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATcc  >  2:1213775/1‑76 (MQ=255)                                                                          |                                                                        GTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCC  >  REL606/4167385‑4167529

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