Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,235,966	Δ1 bp	coding (61/243 nt)	yjbE →	uncharacterized protein YjbE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,235,964	0	G	.	100.0%	147.4 / NA	33	coding (59/243 nt)	yjbE	uncharacterized protein YjbE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (19/14);  total (19/14)

CATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGA  >  NC_000913/4235884‑4236038                                                                                 |                                                                           cattcaATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTg                                                                            >  1:793865/1‑81 (MQ=255)  attcaATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTgg                                                                           <  2:793814/81‑1 (MQ=255)     caATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTggctg                                                                        >  1:793838/1‑78 (MQ=255)       aTGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCa                                                                     >  2:793844/1‑81 (MQ=255)           aGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTg                                                                 >  2:793840/1‑81 (MQ=255)                       aTGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGc                                                      >  1:793871/1‑80 (MQ=255)                            aaaaGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCg                                                >  2:793856/1‑81 (MQ=255)                             aaaGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCgg                                               >  2:793842/1‑81 (MQ=255)                              aaGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGc                                              <  1:793815/81‑1 (MQ=255)                                gTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAg                                            >  1:793868/1‑81 (MQ=255)                                 ttCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAgg                                           <  1:793816/81‑1 (MQ=255)                                  tCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAggg                                          >  2:793846/1‑81 (MQ=255)                                     gTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCg                                       >  2:793851/1‑81 (MQ=255)                                      tATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTc                                        <  2:793817/79‑1 (MQ=255)                                         ggCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGgcag                                   >  1:793850/1‑81 (MQ=255)                                          gCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGgcag                                   >  1:793864/1‑80 (MQ=255)                                           cATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGgcagca                                 <  1:793818/81‑1 (MQ=255)                                                ttGCCATATCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGt                             <  1:793819/80‑1 (MQ=255)                                                 tGCCATATCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTc                            <  2:793820/80‑1 (MQ=255)                                                    cATATCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGtt                        <  1:793821/81‑1 (MQ=255)                                                     atatCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGttt                       >  2:793847/1‑81 (MQ=255)                                                      tatCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTc                      <  1:793823/81‑1 (MQ=255)                                                      tatCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTc                      <  2:793822/81‑1 (MQ=255)                                                       atCTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTcg                     >  1:793857/1‑81 (MQ=255)                                                         cTGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTcggc                   >  2:793858/1‑81 (MQ=255)                                                          tGCGCTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTcggcg                  >  2:793872/1‑81 (MQ=255)                                                            cgcTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTcggc                   <  1:793824/78‑1 (MQ=255)                                                             gcTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTcggcgtg                >  2:793870/1‑78 (MQ=255)                                                             gcTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCggg                >  2:793853/1‑80 (MQ=255)                                                              cTTGCGGCGACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCggg                >  1:793876/1‑79 (MQ=255)                                                                 gcggcgACTCCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGt            <  2:793825/80‑1 (MQ=255)                                                                  cggcgACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTtc          <  2:793826/81‑1 (MQ=255)                                                                  cggcgACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTtc          <  1:793827/81‑1 (MQ=255)                                                                      gACTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTtcctcc      >  1:793855/1‑81 (MQ=255)                                                                        cTTCTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCcgc    >  1:793859/1‑81 (MQ=255)                                                                           cTGCGT‑GGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGa  <  1:793828/80‑1 (MQ=255)                                                                                 |                                                                           CATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGA  >  NC_000913/4235884‑4236038

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