Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,161,978	C→T	P375S (CCG→TCG)	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,978	0	C	T	100.0%	129.5 / NA	37	P375S (CCG→TCG)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/23);  total (14/23)

ACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACCTGTCTG  >  REL606/4161908‑4162050                                                                       |                                                                         aCGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGtcgccg                                                                     <  1:885075/76‑1 (MQ=255)    tgtCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGCCt                                                                  >  2:29122/1‑76 (MQ=255)    tgtCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACt                                                                  <  2:733541/76‑1 (MQ=255)     gtCGACCCAACTAACGACCGTCTTAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTc                                                                 >  1:690759/1‑76 (MQ=255)     gtCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTc                                                                 <  2:478062/76‑1 (MQ=255)        gACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTg                                                              <  2:825092/76‑1 (MQ=255)        gACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTg                                                              >  1:825093/1‑76 (MQ=255)         aCCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGa                                                             <  1:228791/76‑1 (MQ=255)              aCTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAagcag                                                        <  2:990919/76‑1 (MQ=255)                   cGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGaa                                                   <  1:1246214/76‑1 (MQ=255)                        gTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCt                                              >  2:600253/1‑76 (MQ=255)                        gTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCt                                              <  2:343806/76‑1 (MQ=255)                        gTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCt                                              >  1:584949/1‑76 (MQ=255)                         tCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTg                                             <  1:968470/76‑1 (MQ=255)                          cTGAGCGCACTGGTAGAAATCTCCCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGt                                            <  1:957608/76‑1 (MQ=255)                          cTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGt                                            >  2:369029/1‑76 (MQ=255)                            gAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCTGAAAGCCTGTTc                                          <  2:345762/76‑1 (MQ=255)                            gAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTc                                          <  2:268674/76‑1 (MQ=255)                             aGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCTTGAAGCAGCGGAAAGCCTGTTCg                                         <  2:1004537/76‑1 (MQ=255)                               ggcACTGGAAGAAATCTACCGCATGAGGCGCCGTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCgag                                       <  1:1032935/75‑1 (MQ=255)                               cgcACTGGTAGAAATCTCCCGCATGATGCGCCCGGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCgag                                       <  1:868222/76‑1 (MQ=255)                                    tGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCt                                  >  2:308557/1‑76 (MQ=255)                                     ggTAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTg                                 >  2:484698/1‑76 (MQ=255)                                       tAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGtt                               >  1:1062206/1‑76 (MQ=255)                                       tAGAAATCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGtt                               >  2:167892/1‑76 (MQ=255)                                          aaaTCTACCGCATGATGCTCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGttctt                            <  1:207912/76‑1 (MQ=255)                                          aaaTCTACCGAATGGTGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGttctt                            <  1:1490414/76‑1 (MQ=255)                                            aTCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTtcttct                          <  2:1234545/76‑1 (MQ=255)                                            aTCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTtcttct                          >  1:739652/1‑76 (MQ=255)                                            aTCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTtcttct                          <  2:949909/76‑1 (MQ=255)                                             tCTACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTtctc                         <  1:397302/76‑1 (MQ=255)                                               tACCGCATGATGCGCCCTGGCGAGTCGCCGACTCGTGTAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCg                       <  2:346357/76‑1 (MQ=255)                                                   gCATGATGCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGa                   >  1:1358684/1‑76 (MQ=255)                                                      tgatgCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCg                       >  2:270685/1‑69 (MQ=255)                                                      tgatgCGCCCTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCg                       <  1:270686/69‑1 (MQ=255)                                                             cccTGGCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGAc         <  2:1212257/76‑1 (MQ=255)                                                                 ggCGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACTTGt     >  2:1406735/1‑76 (MQ=255)                                                                   cGAGTCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACCTGTCt   <  1:372462/76‑1 (MQ=255)                                                                    gAGCCCCCGACTCCTGAAGCAGCAGAAAGCCTGTTCGAGAACATGTTCTTCTCCGAAGACCGTTATCACCTGTCTg  <  2:1318942/76‑1 (MQ=255)                                                                       |                                                                         ACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGTTCTTCTCCGAAGACCGTTATGACCTGTCTG  >  REL606/4161908‑4162050

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