Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,003,238	G→A	intergenic (‑45/+50)	yigG ← / ← rarD	inner membrane protein/putative transporter RarD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,003,238	0	G	A	100.0%	41.9 / NA	13	intergenic (‑45/+50)	yigG/rarD	inner membrane protein/putative transporter RarD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (4/9);  total (4/9)

CCGTTAGATGTTGGGATAAATATACGTAACATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATGCCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTACGCTGCGTA  >  NC_000913/4003162‑4003312                                                                             |                                                                           ccGTTAGATGTTGGGATAAATATACGTAACATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATAccc                                                                         <  2:735352/80‑1 (MQ=255) ccGTTAGATGTTGGGATAAATATACGTAACATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATAccc                                                                         <  2:735353/80‑1 (MQ=255)      aGATGTTGGGATAAATATACGTAACATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATACCCActct                                                                    >  1:735373/1‑80 (MQ=255)       gATGTTGGGATAAATATACGTAACATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATACCCACtctc                                                                   <  2:735354/80‑1 (MQ=255)        aTGTTGGGATAAATATACGTAACATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATACCCACTCTCAt                                                                 <  2:735355/81‑1 (MQ=255)                 tAAATATACGTAACATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATACCCACTCTCATGCAAGAGTc                                                        <  1:735357/81‑1 (MQ=255)                 tAAATATACGTAACATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATACCCACTCTCATGCAAGAGTc                                                        <  2:735357/81‑1 (MQ=255)                      ataCGTAACATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATACCCACTCTCATGCAAGAGTCGGCa                                                    >  2:735373/1‑80 (MQ=255)                                   aTTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATACCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAat                                      <  1:735359/81‑1 (MQ=255)                                                             gTCAGCATGGCAATATACCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTc             >  1:735380/1‑80 (MQ=255)                                                              tCAGCATGGCAATATACCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTa           <  2:735359/81‑1 (MQ=255)                                                                aGCATGGCAATATACCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTAc          <  1:735361/80‑1 (MQ=255)                                                                       cAATATACCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTACGCTGCGta  >  2:735385/1‑81 (MQ=255)                                                                             |                                                                           CCGTTAGATGTTGGGATAAATATACGTAACATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATGCCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTACGCTGCGTA  >  NC_000913/4003162‑4003312

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