Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,715,326	Δ1 bp	intergenic (‑6/‑97)	yfiE ← / → eamB	putative LysR‑type DNA‑binding transcriptional regulator YfiE/cysteine/O‑acetylserine exporter EamB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,715,323	0	C	.	100.0%	97.4 / NA	23	intergenic (‑3/‑100)	yfiE/eamB	putative LysR‑type DNA‑binding transcriptional regulator YfiE/cysteine/O‑acetylserine exporter EamB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (6/17);  total (6/17)

GCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGG  >  NC_000913/2715243‑2715391                                                                                 |                                                                     gCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAAc                                                                      >  2:485289/1‑81 (MQ=255)   acaATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAAcc                                                                     >  1:485294/1‑80 (MQ=255)           tGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAA‑CCCAAACCTAt                                                           >  1:485293/1‑81 (MQ=255)              gAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAA‑CCCAAACCTATCGa                                                        <  2:485257/81‑1 (MQ=255)               aaGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAA‑CCCAAACCTATCGaa                                                       <  1:485259/81‑1 (MQ=255)               aaGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAA‑CCCAAACCTATCGaa                                                       <  1:485258/81‑1 (MQ=255)                      cAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAA‑CCCAAACCTATCGAAAATATCg                                                <  2:485260/81‑1 (MQ=255)                              aaCCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAA‑CCCAAACCTATCGAAAATATCGAATCTAg                                         >  1:485297/1‑80 (MQ=255)                              aaCCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAA‑CCCAAACCTATCGAAAATATCGAATCTAGa                                        <  1:485261/81‑1 (MQ=255)                              aaCCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAA‑CCCAAACCTATCGAAAATATCGAATCTAGa                                        <  2:485262/81‑1 (MQ=255)                                       ccaccaCGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAA‑CCCAAACCTATCGAAAATATCGAATCTAGAATATaaaaa                               >  2:485302/1‑81 (MQ=255)                                               gTTTTAAGCGTAATAAAACGGCGCAGATCCATAA‑CCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAAcattcatt                       >  2:485292/1‑81 (MQ=255)                                                   tAAGCGTAATAAAACGGCGCAGATCCATAA‑CCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAttttt                    <  2:485263/80‑1 (MQ=255)                                                   tAAGCGTAATAAAACGGCGCAGATCCATAA‑CCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAttttt                    <  1:485264/80‑1 (MQ=255)                                                     aGCGTAATAAAACGGCGCAGATCCATAA‑CCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAttttttt                  >  2:485286/1‑80 (MQ=255)                                                      gCGTAATAAAACGGCGCAGATCCATAA‑CCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTaa                <  1:485265/81‑1 (MQ=255)                                                      gCGTAATAAAACGGCGCAGATCCATAA‑CCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTaa                <  2:485266/81‑1 (MQ=255)                                                        gTAATAAAACGGCGCAGATCCATAA‑CCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAAt              <  1:485267/81‑1 (MQ=255)                                                           ataaAACGGCGCAGATCCATAA‑CCTAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGtt           <  1:485269/81‑1 (MQ=255)                                                           ataaAACGGCGCAGATCCATAA‑CCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGtt           <  2:485270/81‑1 (MQ=255)                                                           ataaAACGGCGCAGATCCATAA‑CCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGtt           <  1:485268/81‑1 (MQ=255)                                                             aaaaCGGCGCAGATCCATAA‑CCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTcc         >  1:485296/1‑81 (MQ=255)                                                              aaaCGGCGCAGATCCATAA‑CCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCA‑TTTTTTAAATGTTCCg        <  2:485271/80‑1 (MQ=255)                                                               aaCGGCGCAGATCCATAA‑CCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCgt       <  2:485272/81‑1 (MQ=255)                                                                    cgcAGATCCATAA‑CCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg  <  1:485273/81‑1 (MQ=255)                                                                                 |                                                                     GCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGG  >  NC_000913/2715243‑2715391

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