Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,161,969	C→T	P372S (CCT→TCT)	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,161,969	0	C	T	100.0%	209.5 / NA	56	P372S (CCT→TCT)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (39/17);  total (39/17)

TCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATG  >  REL606/4161897‑4162041                                                                         |                                                                         tctGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTg                                                                       >  1:585124/1‑76 (MQ=255) tctGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTg                                                                       >  1:986145/1‑76 (MQ=255) tctGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTg                                                                       >  1:1333789/1‑76 (MQ=255) tctGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTg                                                                       <  2:570470/76‑1 (MQ=255) tctGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTg                                                                       >  2:979362/1‑76 (MQ=255)  ctGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTgg                                                                      <  1:92344/76‑1 (MQ=255)  ctGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTgg                                                                      >  1:1702698/1‑76 (MQ=255)    gAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCg                                                                    >  2:1443350/1‑76 (MQ=255)     aaaCCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGa                                                                   >  1:1865762/1‑76 (MQ=255)        ccTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAgcc                                                                >  1:115098/1‑76 (MQ=255)         cTTACGTGCCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAgccg                                                               <  2:898829/76‑1 (MQ=255)          ttACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAgccgc                                                              >  2:344307/1‑76 (MQ=255)          ttACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAgccgc                                                              >  2:1630926/1‑76 (MQ=255)           tACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAgccgcc                                                             >  1:1290207/1‑76 (MQ=255)               tgtCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACt                                                         >  1:1076845/1‑76 (MQ=255)                gtCGACCCAACTAACGACCGGCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTc                                                        >  2:1536208/1‑76 (MQ=255)                  cGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGt                                                      >  1:1045477/1‑76 (MQ=255)                   gACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTg                                                     >  2:127724/1‑76 (MQ=255)                   gACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTg                                                     >  2:267156/1‑76 (MQ=255)                   cccccAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTg                                                     <  2:1199373/74‑1 (MQ=255)                       cAACTAACGACCGTCTTAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAagc                                                 <  1:508490/76‑1 (MQ=255)                          cTAACGCCCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAagcagc                                              <  1:1169900/76‑1 (MQ=255)                          cTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAagcagc                                              >  1:1243121/1‑76 (MQ=255)                           tAACGACCGTCTGAGCGCACTGGTAGAAATCTCCCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCg                                             <  1:42108/76‑1 (MQ=255)                           tAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCg                                             <  2:1712159/76‑1 (MQ=255)                            aaCGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCgg                                            >  1:2024155/1‑76 (MQ=255)                            aaCGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCgg                                            >  2:1233096/1‑76 (MQ=255)                            aaCGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCgg                                            >  2:799552/1‑76 (MQ=255)                               gACCGTCTGAGCGCACTGGTAGAAATCTCCCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGaaa                                         <  1:1024480/76‑1 (MQ=255)                               gACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGCGCAGCCGCGGaaa                                         >  2:732368/1‑76 (MQ=255)                               caCCGTGTGAGCGCACTGGTAGAAATCCCCCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGaaa                                         <  2:727371/75‑1 (MQ=255)                                  cGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGcc                                      >  2:1656965/1‑76 (MQ=255)                                   gTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCt                                     <  1:2048975/76‑1 (MQ=255)                                      tGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGtt                                  >  2:720596/1‑76 (MQ=255)                                      tGAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGtt                                  >  2:1027053/1‑76 (MQ=255)                                       gAGCGCACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTc                                 <  2:893891/76‑1 (MQ=255)                                          cgcAGTGGTAGAAATCTCCCACATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCgag                              <  2:354800/76‑1 (MQ=255)                                            cACTGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGaa                            >  1:58882/1‑76 (MQ=255)                                               tGGTAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCt                         <  2:1699460/76‑1 (MQ=255)                                                  tAGAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGtt                      >  2:1491017/1‑76 (MQ=255)                                                    gAAATCTCCCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGttct                    <  2:600639/76‑1 (MQ=255)                                                    gAAATCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGttct                    >  2:1735056/1‑76 (MQ=255)                                                      naTCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGttcttc                  >  1:1249609/2‑76 (MQ=255)                                                       aTCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTtcttct                 >  2:1664132/1‑76 (MQ=255)                                                        tCTACCGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTtctc                >  1:1106591/1‑76 (MQ=255)                                                             cGCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAg           >  1:599803/1‑76 (MQ=255)                                                              gCATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGa          <  1:704687/76‑1 (MQ=255)                                                               cATGATGCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGAc         <  2:165890/76‑1 (MQ=255)                                                                atgatgCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGAcc        >  2:1860539/1‑76 (MQ=255)                                                                atgatgCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGAcc        >  1:537980/1‑76 (MQ=255)                                                                atgatgCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGAcc        >  2:1319173/1‑76 (MQ=255)                                                                atgatgCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGAcc        >  2:1151272/1‑76 (MQ=255)                                                                atgatgCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGAcc        >  1:794358/1‑76 (MQ=255)                                                                   atgCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGtt     >  1:1999617/1‑76 (MQ=255)                                                                    tgCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTGTTCTCCGAAGACCGTTa    <  2:991768/76‑1 (MQ=255)                                                                    tgCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTa    <  2:1339449/76‑1 (MQ=255)                                                                    tgCGCTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTa    >  2:960760/1‑76 (MQ=255)                                                                      cgcTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATg  <  1:1568274/76‑1 (MQ=255)                                                                      cgcTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATg  <  1:1771795/76‑1 (MQ=255)                                                                      cgcTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATg  <  2:841115/76‑1 (MQ=255)                                                                      cgcTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATg  >  1:416629/1‑76 (MQ=255)                                                                      cgcTCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATg  <  1:762098/76‑1 (MQ=255)                                                                         |                                                                         TCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATG  >  REL606/4161897‑4162041

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