Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,829,251	C→T	H69Y (CAT→TAT)	srlR →	DNA‑binding transcriptional repressor SrlR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,829,251	0	C	T	91.7%	29.7 / NA	12	H69Y (CAT→TAT)	srlR	DNA‑binding transcriptional repressor SrlR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (9/2);  minor base G (0/1);  total (9/3)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAACCTAAATCAGGTAATCACGCCCATGAAACCTCGTCAGCGTCAGGCCGCCATTCTGGAGTATCTGCAAAAGCAGGGTAAATGCTCGGTTGAAGAATTGGCGCAATACTTTGACACCACAGGCACAACCATTCGCAAAGATCTGGTCATTCTGGAACATGCCGGAACCGTCA‑TTCGTACTTATGGCGGAGTGGTGTTGAATAAAGAGGAATCCGATCCGCCTATCGATCATAAAACACTCATCAACACCCACAAGAAAGAGCTGATTGCAGAAGCTGCCGTTAGTTTTATCCATGATGGCGATTCGATCATTCTTGATGCTGGCAGTACCGTTTTGCAGATGGTTCCCCTGCTCTCGCGCTTTAATAACATCACGGTGATGACCAACAGCCTGCATATCGTCAATGCGCTATCCGAACTGG  >  NC_000913/2829021‑2829443                                                                                                                                                                                                                                        |                                                                                                                                                                                                 gAAACCTAAATCAGGTAATCACGCCCAGGAAACCTCGTCAGCGTCAGGCCGCCATTCTGGAGTATCTGCAAAAGCAGGGTAAATGCTCGGTTGAAGAATTGGCGCAATACTTTGACACCACAGGCACAACCATTCGCAAAGATCTGGTCATTCTGGAACATGCCGGAACCGTCA‑TTCGTACTTAGGGCGGAGTGGTGTTGAATAAAGAGGAATCCGATCCGCCTATCGATGATAAAACCCTCAGCAacac                                                                                                                                                                               <  1:43980/250‑1 (MQ=255)                                         gcgtcagGCCGCCATTCTGGAGTATATGCAAAAGCAGGGTAAATGCTCGGTTGAAGAATTGGCGCAATACTTTGACACCACAGGCACAACCATTCGCAAAGATCTGGTCATTCTGGAACATGCCGGAACCGTCA‑TTCGTACTTATGGCGGAGTGGTGTTGAATAAAGAGGAATCCGATCCGCCTATCGATTATAAAACACTCATCAACACCCACaagaaa                                                                                                                                                                     >  1:51384/1‑220 (MQ=255)                                                                                           ttGAAGAATTGGCGCAATACTTTGACACCACAGGCACAACCATTCGCAAAGATCTGGTCATTCTGGAACATGCCGGAACCGTCA‑TTCGTACTTATGGCGGAGTGGTGTTGAATAAAGAGGAATCCGATCCGCCTATCGATTATAAAACACTCATCAACACCCACAAGAAAGAGCTGATTGCAGAAGCTGCCGTTAGTTTTATCCATGATGGCGATTCGATCATTCTTGATGCTGGCAGTACCGTTTTGc                                                                                      >  1:66884/1‑249 (MQ=255)                                                                                                    tGGCGCAATACTTTGCCACCACAGGCACAACCATTCGCAAAGATCTGGTCATTCTGGAAAATGCCGGAACCGTCA‑ATCGTCCTTATGGCGGCGTGGTGTTGAATAAAGAGGAATCAGATACGCCTATCGATTCTACAACACTCATAAACACCcac                                                                                                                                                                           >  1:92147/1‑155 (MQ=255)                                                                                                     ggCGCAATAATTTGACACCACCGGCACAACCATTCGCAAAGATCTGGTCATTCTGAACCACGCAGGACCCGTCA‑TTCGTACTTATGGCGGAGTGGTGTTGAATAACGAGGAATCCGCTCCGCCTATCGATTATAAAACACTCATCAACTCCCACAAGAAAGCGCTGATTGCAGAAGCTGCCGTTCGTTTTATCCATGATGGCGATTCGATCATTCTTGATGCTGGCAGTACCGTTTTGCAGATGGTTcccc                                                                          >  1:164975/1‑251 (MQ=255)                                                                                                      gcgcAATACTTTGACACCACAGGCACAACCATTCGCAAAGATCTGGTCATTCTGGAACATGCCGGAACCGTCA‑TTCGTACTTATGGCGGAGTGGTGTTGAATAAAGAGGAATCCGATCCGCCTATCGATTATAAAACACTCATCAACACCCACAAGAAAGAGCTGATTGCAGAAGCTGCCGTTAGTTTTATCCATGATGGCGATTCGATCATTCTTGATGCTGGCaccgatgaaaacggcgacctcgggcc                                                                         >  1:115574/1‑226 (MQ=255)                                                                                                       cgcAATACTTTGACACCACAGGCACAACCATTCGCAAAGATCTGGTCATTCTGGAACATGCCGGAACCGTCA‑TTCGTACTTATGGCGGAGTGGTGTTGAATAAAGAGGAATCCGATCCGCCTATCGATTATAAAACACTCATCAACACCCACAAGAAAGAGCTGATTGCAGAAGCTGCCGTTAGTTTTATCCATGATGGCGATTCGATCATTCTTGATGCTGGCAGTACCGTTTTGCAGATGGTTCCCCTg                                                                        <  1:11226/251‑1 (MQ=255)                                                                                                         cAATACTTTGACACCACAGGCACAACCCTTCGCACAGATCTGGTCATTCTGGAACCTGCCGGAACCGTCA‑TTCGTACTTATGGCGGAGTGGGGTTGCATAAAGAGGAATCCGATCCGCCTATCGATTATAAAACACTCATCAACACCCACAAGAAAGAGCTGATTGCAGAAGCTGCCGTTAGTTTTATCCATGATGGCGATTCGATCATTCTTGATGCTGGCAGTACCGTTTTGCAGATGGTTCCCCTGct                                                                      >  1:57021/1‑251 (MQ=255)                                                                                                         cAATACTTTGACACCACAGGCACAACCATTCGCAAAGATCTGGTCATTCTGGACCATGCCGGAACCGTAA‑TTCGTACTTATGGAGGAGTGGTGTTGAATCCAGAGGAATCCGATCCGCCTATCGATTATAAAACACTCATCAACACCCTCAAGAAAGAGCTGATTGCAGAAGCTGCCGTTAGTTTTATCCATGATGGCGATTCGATCATTCTTGATg                                                                                                        >  1:165722/1‑217 (MQ=255)                                                                                                                           ngCACAACCATTCGCAAAGATCTGGTCATTCTGGAACATGCCGGAACCG‑CACTTCGTACTTATGGCGGAGTGGTGTTGAATACAGCGGAATACGATCCGACTCTCCATTATAAAAAACTCATCAAAAACCACAAGAAAGAGCTGATTGCAGAAGCTGCCGTTAGTTTTATCCATGATGGCGATTCGATCATTCTTTATGCTGGCAGTACCGTTTTg                                                                                       >  1:31482/2‑216 (MQ=255)                                                                                                                                                                         ncGTCA‑TTCGTACTTATGGCGGAGTGGTGTTGAATAAAGAGGAATCCGATCCGCCTATCGATTATAAAACACTCATCAACACCCACAAGAAAGAGCTGATTGCAGAAGCTGCCGTTAGTTTTATCCATGATGGCGATTCGATCATTCTTGATGCTGGCAGTACCGTTTTGCAGATGGTTCCCCTGCTCTCGCGCTTTAATAACATCACGGTGATGACCAACAGCCTGCATATCGTCAATGCGCTATCCGa       >  1:103665/2‑250 (MQ=255)                                                                                                                                                                             cA‑TTCGTACTTATGGCGGAGTGGTGTTGAATAAAGAGGAATCCGATCCGCCTATCGATTATAAAACACTCATCAACACCCACAAGAAAGAGCTGATTGCAGAAGCTTCCGTTAGTTTTATCCATGATGGCGATTCGATCATTCTTGATTCTGGCAGTACCTTTTTGCAGATTGTTCCCCTGCTCTCGCGCTTTAATAACATCGCGGTGATGACCAACAGCCTGCATATCGTCAATGCGCTATCCGAACTgg  <  1:27182/251‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                 GAAACCTAAATCAGGTAATCACGCCCATGAAACCTCGTCAGCGTCAGGCCGCCATTCTGGAGTATCTGCAAAAGCAGGGTAAATGCTCGGTTGAAGAATTGGCGCAATACTTTGACACCACAGGCACAACCATTCGCAAAGATCTGGTCATTCTGGAACATGCCGGAACCGTCA‑TTCGTACTTATGGCGGAGTGGTGTTGAATAAAGAGGAATCCGATCCGCCTATCGATCATAAAACACTCATCAACACCCACAAGAAAGAGCTGATTGCAGAAGCTGCCGTTAGTTTTATCCATGATGGCGATTCGATCATTCTTGATGCTGGCAGTACCGTTTTGCAGATGGTTCCCCTGCTCTCGCGCTTTAATAACATCACGGTGATGACCAACAGCCTGCATATCGTCAATGCGCTATCCGAACTGG  >  NC_000913/2829021‑2829443

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