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	100.0%	57.8 / NA	17	G169G (GGG→GGA)	uidA	beta‑glucuronidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/8);  total (9/8)

CCTGCCAGTCAACAGACGCGTGGTTACAGTCTTGCGCGACATGCGTCACCACGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATCCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCG  >  NC_000913/1695456‑1695790                                                                                                              |                                                                                                                                                                                                                                  catctcaactggctgACGCGTGGTTACAGTCTTGCGCGACATGCGTCACCACGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATa                                                                                      <  2:367452/237‑1 (MQ=255)         ctggctgACGCGTGGTTACAGTCTTGCGCGACATGCGTCACCACGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACgg                                                                              >  1:299802/7‑251 (MQ=255)            gctgACGCGTGGTTACAGTCTTGCGCGACATGCGTCACCACGGTGATATCGTCAACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTACAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGt                                                                           >  1:468853/4‑251 (MQ=255)                  gcgTGGTTACAGTCTTGCGCGACATGCGTCACCACGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGGTCAGTTCGTTGGTCACACAAACGGTGATACGTACGCTTTTCCCGGCAATAACATACGGCGGGACATc                                                                     <  1:173247/251‑1 (MQ=255)                         tACAGTCTTGCGCGACATGCGTCACCACGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGcc                                                                                                                                           <  1:13125/174‑1 (MQ=255)                         tACAGTCTTGCGCGACATGCGTCACCACGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGcc                                                                                                                                           >  2:13125/1‑174 (MQ=255)                          aCAGTCTTGCGCGACATGCGTCACCACGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCa                                                              >  1:328805/1‑250 (MQ=255)                                 tGCGCGCCATGCGTCACCACGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTGTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCAAATGGCg                                                       <  1:293862/250‑1 (MQ=255)                                        cATGCGTCACCACGGTGATATCGGCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCAAAGGGCGTATAgccg                                               <  2:178440/251‑1 (MQ=255)                                                          tatCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACt                             >  1:99407/1‑251 (MQ=255)                                                                 cACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATTGATTCCGGCATAGTTAAAGAAATCATTGAAGTAAGACTGCGTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGat                      <  2:455375/251‑1 (MQ=255)                                                                    ccAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGattat                   <  2:104186/251‑1 (MQ=255)                                                                          gtTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATTACGACATAGTTACAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGAccc             >  1:162030/1‑251 (MQ=255)                                                                                cGTGGTGTAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACtt       >  1:85634/1‑251 (MQ=255)                                                                                     tgtAGAGCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCg  >  1:265239/1‑251 (MQ=255)                                                                                          agCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCAt                                                                                                                                                                                                             >  2:110695/1‑43 (MQ=255)                                                                                          agCATTACGCTGCGATGGATTCCGGCATAGTTAAAGAAATCAt                                                                                                                                                                                                             <  1:110695/43‑1 (MQ=255)                                                                                                              |                                                                                                                                                                                                                                  CCTGCCAGTCAACAGACGCGTGGTTACAGTCTTGCGCGACATGCGTCACCACGGTGATATCGTCCACCCAGGTGTTCGGCGTGGTGTAGAGCATTACGCTGCGATGGATCCCGGCATAGTTAAAGAAATCATGGAAGTAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCG  >  NC_000913/1695456‑1695790

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