Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_012967	3,893,550:1	+G	100%	intergenic (+5/‑51)	kup → / → ECB_RS25215	low affinity potassium transporter Kup/IS3 family transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_012967	3,893,550	1	.	G	90.9%	19.4 / ‑2.7	11	intergenic (+5/‑51)	kup/ECB_RS25215	low affinity potassium transporter Kup/IS3 family transposase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (7/3);  total (7/4)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA‑G‑ACTGTACTGCACCCATTTTGTTGGACGAT  >  NC_012967/3893515‑3893580                                     |                                tATCGAACTGGGTACTCAGGTCGAAATCTAACGCCa                                  >  1:149672/1‑36 (MQ=255)  aTCGAACTGGGTACTCAAGTCGAAATCTAACGCCA‑g                                <  1:188005/36‑1 (MQ=255)     gAACTGGGTACTCAGGTTGAAATTTAACGTCA‑ggac                             >  1:204224/1‑33 (MQ=12)        cTGGGTACTCAGGTCGAAATCTAACGCCAGG‑ACTGt                          >  1:985841/1‑36 (MQ=37)         tGGGTAGTCAGGTCGAGATCTAGCGCCAGG‑ACTGTa                         <  1:861215/36‑1 (MQ=255)                    gTCGAAATCTAACGCCAGG‑ACTGTACTGCACCCAtt              >  2:790420/1‑36 (MQ=37)                    gTCGAAATCTAACACCAGG‑CCTGTACTGCACCCAtt              >  1:833964/1‑36 (MQ=12)                      cGAAATCTAACGCCAGG‑ACTGAACTGCACCCAtttt            >  1:55673/1‑36 (MQ=21)                       aaaaTCTAACGCCAGG‑ACTGTACTGCACCCATTttg           <  1:943748/35‑1 (MQ=21)                         aaTCTAACGCCAGG‑ACTGTACTGCACCCATTttgtt         <  2:188578/36‑1 (MQ=18)                             tAACGCCA‑GTACTGTACTGCACCCATTTTGTTGGAc     <  2:318651/36‑1 (MQ=17)                              aaCGCCAGG‑ACTGTACTGCACCCATTTTGTTGGACg    >  2:799598/1‑36 (MQ=14)                              aaCGCCAGG‑ACTGTACTGCACCCATTTTGTTGGACg    >  1:339039/1‑36 (MQ=14)                                cGCCAGG‑ACTGTACTGCACCCATTTTGTTGGACGAt  >  1:261264/1‑36 (MQ=11)                                     |                                TATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA‑G‑ACTGTACTGCACCCATTTTGTTGGACGAT  >  NC_012967/3893515‑3893580

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