Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,695,565	C→T	G169G (GGG→GGA)	uidA ←	beta‑glucuronidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,695,565	0	C	T	90.0%	23.7 / ‑2.9	10	G169G (GGG→GGA)	uidA	beta‑glucuronidase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (4/5);  total (5/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 = 2.41e-01

TGGCCACCACCTGCCAGTCAACAGACGCGTGGTTACAGTCTTGCGCGACATGCGTCACCACGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATCCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATG  >  NC_000913/1695447‑1695795                                                                                                                       |                                                                                                                                                                                                                                       tttagcattcatctcaactggctgACGCGTGGTGACAGTCTTGCGCGACATGCGTCACCACGGTGCTATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGGTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCCCACAAACGGTGATTCGTACACTTTTc                                                                                                    >  2:171635/24‑251 (MQ=255)    agcattcatctcaactggctgACGCGTGGTTACAGTCTTGCGCGACATGCGTCACCACGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCg                                                                                                 >  1:324478/21‑251 (MQ=255)    agcattcatctcaactggctgACGCGTGGTTACAGTCTTGCGCGACATGCGTCACCACGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAAGCATGGAAGTAAGACGGCTTTTTCTTGCCGTTTTCGTCTGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCg                                                                                                 <  1:159331/231‑1 (MQ=255)    agcattcatctcaactggctgACGCGTGGTTACAGTCCTGCGCGACATGCGTCACCACGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCTGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTccc                                                                                                  <  1:171635/230‑1 (MQ=255)     caccacATGCCAGTCAACAGACGCGTGGTTACAGTCTTGCGCGACATGCGTCACCACGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATCCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCgg                                                                                                >  2:135028/1‑251 (MQ=255)       attcatctcaactggctgACGCGTGGCTACAGTCTTGCGCGACATGCGTCACCACGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATCGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGc                                                                                               >  2:332435/18‑250 (MQ=255)            tctcaactggctgACGCGTGGTTACAGTCTTGCGCGACATGCGTCACCACGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATaa                                                                                          >  2:130973/13‑250 (MQ=255)                                                             cGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTcc                                        <  1:183380/251‑1 (MQ=255)                                                                                            ggtgTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGGTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTg          <  1:294815/250‑1 (MQ=255)                                                                                                   agCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGGCTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATg  <  2:168072/251‑1 (MQ=255)                                                                                                                       |                                                                                                                                                                                                                                       TGGCCACCACCTGCCAGTCAACAGACGCGTGGTTACAGTCTTGCGCGACATGCGTCACCACGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATCCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATG  >  NC_000913/1695447‑1695795

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: ‑