Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_012967	3,339,158	A→C	100%	intergenic (+22/‑4)	dusB → / → fis	tRNA dihydrouridine synthase DusB/DNA‑binding transcriptional regulator Fis

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_012967	3,339,158	0	A	C	100.0%	44.3 / NA	16	intergenic (+22/‑4)	dusB/fis	tRNA dihydrouridine synthase DusB/DNA‑binding transcriptional regulator Fis
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (5/11);  total (5/11)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTTGCGTAAACAGAAATAAAGAGCTGACAGAACTATGTTCGAACAACGCGTAAATTC  >  NC_012967/3339128‑3339184                               |                           tttGCGTAAACAGAAATAAAGAGCTGACAGCACTAt                       <  1:768788/36‑1 (MQ=38) tttGCGTAAACAGAAATAAAGAGCTGACAGCACTAt                       >  1:89192/1‑36 (MQ=38)   tGCGTAAACAGAAATAAAGAGCTGACAGCACTATGt                     <  1:339432/36‑1 (MQ=39)    acGTAAACAGAAATAAAGAGCTGACAGCACTATGtt                    <  1:919511/35‑1 (MQ=38)     cGTAAACAGAAATAAAGAGCTGACAGCACCATGtcc                   >  1:366208/1‑34 (MQ=25)      gTAAACAGAAATAAAGAGCTGACAGCACTATGTTCg                  <  2:133926/36‑1 (MQ=39)       tAAACAGAAATAAAGAGCTGACAGCACTATGTTCGa                 <  2:413743/36‑1 (MQ=39)        aaaCAGAAATAAAGAGCTGACAGCACTATGTTCGaa                <  2:152870/36‑1 (MQ=39)        aaaCAGAAATAAAGAACTCACAGCACTATGTTCGaa                <  1:466994/36‑1 (MQ=16)           cAGAAATAAAAAACTGACAGCACTATGTTCGaacaa             <  1:551724/36‑1 (MQ=19)            aGAAATAAAGAGCTGACAGCACTATGTTCGaacaac            >  1:194947/1‑36 (MQ=39)                aTAAAGAGCTGACAGCACTATGTTCGAACAACGCGt        <  1:256298/36‑1 (MQ=39)                  aaaGTGCTGACAGCACTATGTTCGAACAACGCGTaa      <  1:713281/36‑1 (MQ=37)                  aaaGGGCTGACAGCACTATGTTCGAACAACGCGTaa      <  1:517003/36‑1 (MQ=37)                     gagCTGACAGCACTATGTTCGAACAACGCGTAAAtt   >  2:166040/1‑36 (MQ=39)                      agCTGACAGCACTATGTTCGAACAACGCGTAAATTc  >  1:800911/1‑36 (MQ=39)                               |                           TTTGCGTAAACAGAAATAAAGAGCTGACAGAACTATGTTCGAACAACGCGTAAATTC  >  NC_012967/3339128‑3339184

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