Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,146,012	A→G	D50G (GAC→GGC)	rpoD →	RNA polymerase sigma factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,146,012	0	A	G	100.0%	168.2 / NA	41	D50G (GAC→GGC)	rpoD	RNA polymerase sigma factor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/21);  total (20/21)

GAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACG  >  REL606/3145939‑3146083                                                                          |                                                                        gAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGGCa                                                                       >  1:223520/1‑76 (MQ=255) gAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGGCa                                                                       >  1:431041/1‑76 (MQ=255)    gTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGGCATgg                                                                    <  2:296381/76‑1 (MQ=255)          gACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGGCATGGGCATTc                                                              <  2:828034/76‑1 (MQ=255)           aCCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGGCATGGGCATTCa                                                             >  1:434898/1‑76 (MQ=255)             cATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGGCATGGGCATTCAgg                                                           <  1:218800/76‑1 (MQ=255)              aTCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGGCATGGGCACTTAGGt                                                          <  1:69913/76‑1 (MQ=255)                 tGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCTAGGGCATGGGCATTCAGGTGAt                                                       <  1:6957/76‑1 (MQ=255)                   ccGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGGCATGGGCATTCAGGTGATgg                                                     <  2:918602/76‑1 (MQ=255)                       aaGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGGCATGGGCATTCAGGTGATGgaaga                                                 >  2:763705/1‑76 (MQ=255)                       aaGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGGCATGGGCATTCAGGTGATGgaaga                                                 <  1:257602/76‑1 (MQ=255)                         gATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGGCATGGGCATTCAGGTGATGGaagaag                                               >  1:216583/1‑76 (MQ=255)                           tatCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGGCATGAGCATTCAGGTGATGGAAGAAGCa                                             >  1:312419/1‑76 (MQ=255)                               gtcgATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCgg                                         <  1:359983/76‑1 (MQ=255)                                tcgATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGa                                        <  1:544595/76‑1 (MQ=255)                                tcgATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGa                                        >  1:1066692/1‑76 (MQ=255)                                 cgATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGAt                                       <  1:699092/76‑1 (MQ=255)                                    ttCAGATCAGATCGAAGACATCATCCAAATGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGcc                                    >  2:78289/1‑76 (MQ=255)                                     tcagatcagaTCGAAGACATCATCCAAATGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCg                                   <  1:543617/76‑1 (MQ=255)                                     tcagatcagaTCGAAGACATCATCCAAATGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCg                                   <  1:862018/76‑1 (MQ=255)                                      cagatcagatCGAAGACATCATCCAAATGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCga                                  >  2:439273/1‑76 (MQ=255)                                       agatcagatcGAAGACATCATCCAAATGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCgat                                 >  1:828593/1‑76 (MQ=255)                                            agatcGAAGACATCATCCAAATGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCt                            >  2:485685/1‑76 (MQ=255)                                             gatcGAAGACATCATCCAAATGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTg                           <  1:229827/76‑1 (MQ=255)                                                cGAAGACATCATCCAAATGATCAACGGCATGGNCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTNATg                        >  2:1121975/1‑76 (MQ=255)                                                 gAAGACATCATCCAAATGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATgc                       >  1:570930/1‑76 (MQ=255)                                                  aaGACATCATCCAAATGATCAACGGCATGGGCTTTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATgct                      >  2:935296/1‑76 (MQ=255)                                                   aGACATCATCCAAATGATCAACGGCATGGACATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATgctg                     <  1:66166/76‑1 (MQ=255)                                                     aCATCATCCAAATGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATgctggc                   >  1:650486/1‑76 (MQ=255)                                                        tcatcCAAATGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGa                <  1:573766/76‑1 (MQ=255)                                                        tcatcCAAATGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGa                >  1:1180580/1‑76 (MQ=255)                                                            ccAAATGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAac            <  1:799451/76‑1 (MQ=255)                                                                aTGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACcg        >  1:815317/1‑76 (MQ=255)                                                                 tGATCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACcgc       >  1:1183202/1‑76 (MQ=255)                                                                 tGATCAACGGCATAGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACcgc       >  2:8372/1‑76 (MQ=255)                                                                    tCAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGa    <  1:358593/76‑1 (MQ=255)                                                                     cAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGAc   >  1:887143/1‑76 (MQ=255)                                                                     cAACGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGAc   <  2:884466/76‑1 (MQ=255)                                                                      aaCNGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGNTGATCTGATGCTGGCTGAAAACACCGCGGACg  <  2:1130085/76‑1 (MQ=255)                                                                      aaCGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGTTCTGATGCTGGCTGACAACACCGCGGACg  <  1:78289/76‑1 (MQ=255)                                                                      aaCGGCATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACg  <  1:666804/76‑1 (MQ=255)                                                                          |                                                                        GAGGTCAATGACCATCTGCCGGAAGATATCGTCGATTCAGATCAGATCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACG  >  REL606/3145939‑3146083

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