Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,446,572	C→T	G238S (GGT→AGT)	feaR ←	DNA‑binding transcriptional activator FeaR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,446,572	0	C	T	81.8%	12.9 / ‑0.0	12	G238S (GGT→AGT)	feaR	DNA‑binding transcriptional activator FeaR
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (6/3);  total (6/6)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

ATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCG  >  NC_000913/1446332‑1446760                                                                                                                                                                                                                                                 |                                                                                                                                                                                             atCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTGGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTTCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACTTTTATCGGc                                                                                                                                                                                     <  1:232585/250‑1 (MQ=255)          aaCTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACTTTTATCGGCAAACATTCg                                                                                                                                                                            <  1:239872/250‑1 (MQ=255)              aaTTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAATCGTTGCTTAAATACCTTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACTTTTATCGGCAAACATTCGGTAc                                                                                                                                                                        <  1:24250/250‑1 (MQ=255)                                                    ggAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATACTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACTTACTGACATACCTGTCTCTGTGGCGATCCAc                                                                                                                                 <  1:274563/251‑1 (MQ=255)                                                                                                                              gAAAAGCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCGGCACAAAAATCGAGACGACGGTTACGAATAGATTGCGCGACGACCAAACGTTTAGCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCAc                                                                                                                                 <  1:324327/177‑1 (MQ=255)                                                                                                                                     cccAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACTTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTg                                                 >  1:282682/1‑250 (MQ=255)                                                                                                                                                                                   gCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACTTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCg  >  1:248822/1‑251 (MQ=255)                                                                                                                                                                                          gCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACTTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTNTTCNNNNntannnnnn                                                                                                >  2:195309/1‑144 (MQ=255)                                                                                                                                                                                                 cacaAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACTTTTATCGGCAAACATTCGGAACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTACGGGCGTAATATCTCTTCGCGAATATTATCGTCTATTAACGTAACCACTTTTTGAAACTGACGTTCacgac                                                        >  1:325367/1‑183 (MQ=255)                                                                                                                                                                                                                         gAATATATTGCGCGACTACCAAACTTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAatat                                                                                                                   >  1:318913/1‑100 (MQ=255)                                                                                                                                                                                                                                    ncgACTACCAAACTTTTCTCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGaca          >  1:261130/2‑194 (MQ=255)                                                                                                                                                                                                                                          aCCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGc                                    <  1:269907/162‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                             ATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCG  >  NC_000913/1446332‑1446760

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