Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	388,654	Δ1 bp	intergenic (+8/+99)	tauD → / ← hemB	taurine dioxygenase/porphobilinogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	388,654	0	G	.	100.0%	63.5 / NA	14	intergenic (+8/+99)	tauD/hemB	taurine dioxygenase/porphobilinogen synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (9/5);  total (9/5)

GATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA  >  NC_000913/388610‑388695                                             |                                          aatttccaggtcagacgggctgcaagttgcagaccgttaCGAGA‑TGGAc                                      <  1:1730751‑M2/10‑1 (MQ=255)  atttccaggtcagacgggctgcaagttgcagaccgttaCGAGA‑TGGACg                                     >  1:109416‑M2/39‑49 (MQ=255)    ttccaggtcagacgggctgcaagttgcagaccgttaCGAGA‑TGGACGGt                                   >  1:2074275‑M2/37‑49 (MQ=255)        aggtcagacgggctgcaagttgcagaccgttaCGAGA‑TGGACGGTcccc                               <  1:207698‑M2/17‑1 (MQ=255)           tcagacgggctgcaagttgcagaccgttaCGAGA‑TGGACGGTCCCCTCg                            >  1:721071‑M2/30‑49 (MQ=255)           tcagacgggctgcaagttgcagaccgttaCGAGA‑TGGACGGTCCCCTCg                            >  1:820941‑M2/30‑49 (MQ=255)                        aagttgcagaccgttaCGAGA‑TGGACGGTCCCCTCGCCCCCTTggg                  >  1:656560‑M2/17‑46 (MQ=255)                           ttgcagaccgttaCGAGA‑TGGACGGTCCCCTCGCCCCCTTGGGGAGAgg            >  1:1353254‑M2/14‑49 (MQ=255)                              cagaccgttaCGAGA‑TGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt          <  1:1134194‑M2/38‑1 (MQ=255)                              cagaccgttaCGAGA‑TGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt          >  1:215392‑M2/11‑48 (MQ=255)                               agaccgttaCGAGA‑TGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt          >  1:1751908‑M2/10‑47 (MQ=255)                                    gttaCGAGA‑TGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTg   <  1:288851‑M2/45‑1 (MQ=255)                                     ttaCGAGA‑TGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGt    >  1:1412954‑M2/4‑47 (MQ=255)                                     ttaCGAGA‑TGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGa  <  1:884083‑M2/46‑1 (MQ=255)                                             |                                          GATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA  >  NC_000913/388610‑388695

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