Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,186,297	G→A	R37R (CGG→CGA)	insC5 →	IS2 insertion element repressor InsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,186,297	0	G	A	88.9%	15.0 / ‑5.8	9	R37R (CGG→CGA)	insC5	IS2 insertion element repressor InsA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (3/5);  minor base C (1/0);  total (4/5)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCT‑CCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAAT  >  NC_000913/3186057‑3186343                                                                                                                                                                                                                                                  |                                               cGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCT‑CCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGACAACATGgtg                                      <  2:443352/251‑1 (MQ=255)         ttAAAAATATTTATTGTTCATTATCCGGTATTCGACCGGCACCCTGAATAT‑CCTGACAACCAATATCACTTAAATAAGTGATAGTCTTAATCATAGTTATTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGGAGAAACGCAGAAGGCGTACCACACAGGCAAAGATCCCAATTGTTCAGCAGAGCTTGGAACCGGGGATGACGGTCTCCATCGATGACCGCCAACATGGAGTAGcagcc                              >  2:323986/1‑251 (MQ=14)           aaaaaTATTTATTGTTCATTATCCGTTATTAGACTAGCCCCCTGAATCT‑CCAGACATCCAATAGCACTTAAATAAGTGATAGTCTTAATTCTAGTTGTTAGACTAGTCATTGGAGAACAGATGATTGAGGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATTTCTATTGTTCAGCAGAGCTTTGAACCTGGGATGACGGTCTCCCTCGTTGCCCGACAACATGGTGTAGcagcca                             <  1:95383/250‑1 (MQ=14)              aaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCT‑CCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGACAACATGGTGTAGCAgccagcca                         <  1:168030/251‑1 (MQ=255)               ctatttgttGATCATTATCAGTTATTAGACAGGCCCAATGAAACTACCGCACAAC‑AATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAAAACAGGTGATTGATGCCTTAGGGCCGGAGAAACGCAGACGGCGTCCAACACAGGATAATATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGACAAAATGGTGTAGcagccagcg                          >  2:396174/2‑248 (MQ=11)                     attGTTCATTATCCGTTATTCGACTGGCCCCCTGCATCT‑CCAGACAAACAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCAGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGACAACATGGTGTAGCAGCCAGCCAGTTAttt                  >  1:505270/1‑251 (MQ=255)                       tGTTCATTATCCGTTATTAGACTGTCACCCTGAATCT‑CCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAGACTAGTTTTTAGACTAGTCATTTGAGAACAGATGATTGATGTCTTAGTGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCGTTTGTTCAGCAGATCTTTTAACCGGGGATGACGGTCTCCCTCGTTTCCCGACAACAGAGTGCAACAGCCAGCCAGtaagttc                 <  2:279566/250‑7 (MQ=11)                          tCATTATCCGTTATTAGACTGGCCCCCTGAATCT‑CCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATGGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGGTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGACAACATTGTGTAGCAGCCAGCCAGTTATTTctct              <  2:587338/250‑1 (MQ=255)                                     tattaGACTGGCCCCCTGAATCT‑CCAGACAACCCATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGACAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAAt  >  1:488959/1‑251 (MQ=11)                                                                                                                                                                                                                                                  |                                               CGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCT‑CCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAAT  >  NC_000913/3186057‑3186343

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