Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,434,404	C→T	Q98* (CAG→TAG)	ytfH →	putative transcriptional regulator YtfH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,434,404	0	C	T	100.0%	134.6 / NA	38	Q98* (CAG→TAG)	ytfH	putative transcriptional regulator YtfH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/21);  total (17/21)

CAGGATGGTTTTCTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGG  >  NC_000913/4434326‑4434475                                                                               |                                                                        cAGGATGGTTTTCTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAg                                                                       <  1:816278/81‑1 (MQ=255)  aGGATGGTTTTCTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAgg                                                                      >  1:816328/1‑81 (MQ=255)    gATGGTTTTCTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTg                                                                    <  1:816279/81‑1 (MQ=255)     aTGGTTTTCTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGa                                                                   <  2:816280/81‑1 (MQ=255)      tGGTTTTCTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAg                                                                  <  2:816281/81‑1 (MQ=255)        gTTTTCTTAACCGTATCGCGTATCCGGTGGAGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGc                                                                 <  2:816282/80‑1 (MQ=255)            tCTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGaaa                                                             <  2:816283/80‑1 (MQ=255)             cTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGaaaa                                                            <  1:816284/80‑1 (MQ=255)             cTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAg                                                           >  2:816332/1‑81 (MQ=255)             cTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAg                                                           >  2:816318/1‑81 (MQ=255)               tAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGt                                                         >  2:816336/1‑81 (MQ=255)               tAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGt                                                         >  2:816330/1‑81 (MQ=255)                aaCCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGtt                                                        >  1:816320/1‑81 (MQ=255)                aaCCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGt                                                         <  2:816285/80‑1 (MQ=255)                 aCCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGtt                                                        <  2:816286/80‑1 (MQ=255)                       tCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCa                                                  <  1:816287/80‑1 (MQ=255)                        cgcgTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACt                                                >  1:816337/1‑81 (MQ=255)                         gcgTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACt                                                <  1:816288/80‑1 (MQ=255)                         gcgTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTg                                               >  1:816329/1‑81 (MQ=255)                         gcgTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTg                                               >  1:816323/1‑81 (MQ=255)                           gTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTGGc                                             >  1:816341/1‑81 (MQ=255)                            tATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTGGCg                                            <  2:816289/81‑1 (MQ=255)                                cGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACt                                        <  1:816290/81‑1 (MQ=255)                                    ggtgCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGAt                                    >  2:816344/1‑81 (MQ=255)                                       gccgccGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGAtt                                   >  2:816350/1‑79 (MQ=255)                                        ccgccgCATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAg                                >  2:816343/1‑81 (MQ=255)                                          gccgcATGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGtt                              <  2:816291/81‑1 (MQ=255)                                               aTGTGGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAAtt                         >  2:816345/1‑81 (MQ=255)                                                  tgGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGGGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGc                      <  1:816292/81‑1 (MQ=255)                                                  tgGAATATAGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGc                      <  1:816293/81‑1 (MQ=255)                                                       tataGCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAg                 <  1:816294/81‑1 (MQ=255)                                                           gCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGgtgt             >  2:816342/1‑81 (MQ=255)                                                           gCCTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGgtg              <  1:816295/80‑1 (MQ=255)                                                            ccTCACGCCGCTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGtt            <  1:816296/81‑1 (MQ=255)                                                                    cgcTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGc    >  1:816356/1‑81 (MQ=255)                                                                    cgcTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGc    <  2:816297/81‑1 (MQ=255)                                                                     gcTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCg   <  1:816298/81‑1 (MQ=255)                                                                      cTGGGCGAGTAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCgg  >  2:816347/1‑81 (MQ=255)                                                                               |                                                                        CAGGATGGTTTTCTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGG  >  NC_000913/4434326‑4434475

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