Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,118,141	Δ1 bp	intergenic (‑51/‑374)	glpF ← / → zapB	glycerol facilitator/cell division factor ZapB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,118,140	0	C	.	100.0%	125.0 / NA	28	intergenic (‑50/‑375)	glpF/zapB	glycerol facilitator/cell division factor ZapB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (16/12);  total (16/12)

ACTGGCCTTTCAAGGTTGATGTTTGACTCATAATCGGATCCTGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGACCTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAACGAGCAAAAACGAGAAATATCGAACTTAAAATGTG  >  NC_000913/4118060‑4118215                                                                                 |                                                                            aCTGGCCTTTCAAGGTTGATGTTTGACTCATAATCGGATCCTGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGAc                                                                             <  2:764287/81‑1 (MQ=255)  cTGGCCTTTCAAGGTTGATGTTTGACTCATAATCGGATCCTGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGAct                                                                            <  1:764288/81‑2 (MQ=255)     gCCTTTCAAGGTTGATGTTTGACTCATAATCGGATCCTGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGActcc                                                                          <  2:764289/80‑4 (MQ=255)           cAAGGTTGATGTTTGACTCATAATCGGATCCTGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGA‑CTCCACGATGc                                                                  >  2:764336/1‑81 (MQ=255)            aaGGTTGATGTTTGACTCATAATCGGATCCTGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGA‑CTCCACGATGCt                                                                 >  1:764326/1‑81 (MQ=255)                 tGATGTTTGACTCATAATCGGATCCTGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAg                                                            >  1:764342/1‑81 (MQ=255)                      tttGACTCATAGTCGGATCCTGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATg                                                       <  2:764290/81‑1 (MQ=255)                      tttGACTCATAATCGGATCATGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCAt                                                        <  2:764291/80‑1 (MQ=255)                           cTCATAATCGGATCCTGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTa                                                  <  2:764293/81‑1 (MQ=255)                           cTCATAATCGGATCCTGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTa                                                  <  2:764292/81‑1 (MQ=255)                              aTAATCGGATCCTGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAAc                                               >  2:764339/1‑81 (MQ=255)                                aaTCGGATCCTGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACt                                              >  1:764340/1‑80 (MQ=255)                                 aTCGGATCCTGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTa                                            <  1:764294/81‑1 (MQ=255)                                     gATCCTGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGt                                        >  2:764338/1‑81 (MQ=255)                                           gAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAAc                                    >  1:764349/1‑79 (MQ=255)                                            aaGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAACGAg                                 <  1:764295/81‑1 (MQ=255)                                                  tAATGTTTGTTGTATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAACGAGCaaaaa                           >  1:764332/1‑81 (MQ=255)                                                   aaTGTTTGTTGTATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAACGAGCaaaaa                           <  1:764296/80‑1 (MQ=255)                                                   aaTGTTTGTTGTATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAACGAGCAAAAAc                          >  1:764344/1‑81 (MQ=255)                                                           ttgtATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAACGAGCAAAAACGAGAAata                  <  1:764297/81‑1 (MQ=255)                                                            tgtATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAACGAGCAAAAACGAGAAatat                 >  1:764343/1‑81 (MQ=255)                                                            tgtATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAACGAGCAAAAACGAGAAatat                 >  2:764352/1‑81 (MQ=255)                                                             gtATGCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAACGAGCAAAAACGAGAAatat                 <  1:764298/80‑1 (MQ=255)                                                                 gCGTGAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAACGAGCAAAAACGAGAAATATCGAAc            >  1:764335/1‑81 (MQ=255)                                                                     gAAAGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAACGAGCAAAAACGAGAAATATCGAACTTa         >  1:764348/1‑80 (MQ=255)                                                                      aaaGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAACGAGCAAAAACGAGAAATATCGAACTTaaa       <  2:764299/81‑1 (MQ=255)                                                                      aaaGTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAACGAGCAAAAACGAGAAATATCGAACTTaa        >  2:764350/1‑80 (MQ=255)                                                                         gTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTATCTATCAAAAACGAGAAATATCGAACTTATAAtg    >  1:885511/1‑81 (MQ=255)                                                                         gTCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAACGAGCAAAAACGAGAAATATCGAACTTAAAAtg    <  2:764300/81‑1 (MQ=255)                                                                          tCACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAACGAGCAAAAACGAGAAATATCGAACTTAAAAtgt   <  2:764301/81‑1 (MQ=255)                                                                           cACGGA‑CTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAACGAGCAAAAACGAGAAATATCGAACTTAAAAtgtg  >  2:764369/1‑81 (MQ=255)                                                                                 |                                                                            ACTGGCCTTTCAAGGTTGATGTTTGACTCATAATCGGATCCTGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGACCTCCACGATGCTTGTAGGCATGCTGTAAACTTATCGTTAACGAGCAAAAACGAGAAATATCGAACTTAAAATGTG  >  NC_000913/4118060‑4118215

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