Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,070,150	T→C	44.5%	R11G (AGA→GGA)	insC1 ←	IS2 repressor TnpA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,070,150	0	T	C	44.5%	‑1.0 / 6.9	9	R11G (AGA→GGA)	insC1	IS2 repressor TnpA
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (2/2);  total (4/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTT  >  NC_000913/2070129‑2070295                      |                                                                                                                                                  tttCCTGTGTGGTACGCCGTCCGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGGTTCAGGGGGCCAGTCTAGTGGc                     >  8:248806/1‑148 (MQ=17) tttCCTGTGTGGTACGCCGTCCGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGGTTCAGGGGGCCAGTCTAGTGGc                     >  8:254909/1‑148 (MQ=17)    ccTGTGTGGTCCGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGAATTAAGACTATCACTTATTTAAGTGATATTGGGTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGa                   <  5:391503/147‑1 (MQ=16)    ccTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGa                   <  3:174767/147‑1 (MQ=33)     cTGTGTGGTACGCCGTCCGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGGTTCAGGGGGCCAGTCTAGTGGCGAAg                 <  7:248806/148‑1 (MQ=25)     cTGTGTGGTACGCCGTCCGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGGTTCAGGGGGCCAGTCTAGTGGCGAAg                 <  7:254909/148‑1 (MQ=25)            gTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCAtcct           <  8:89388/147‑1 (MQ=40)             tACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAATTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCAtactcc         >  6:252392/1‑148 (MQ=255)                     cTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGttt  >  2:222089/1‑147 (MQ=255)                      |                                                                                                                                                  TTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTT  >  NC_000913/2070129‑2070295

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