Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	751,665	G→A	P169S (CCA→TCA)	ybgQ ←	putative fimbrial usher protein YbgQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	751,665	0	G	A	100.0%	28.9 / NA	11	P169S (CCA→TCA)	ybgQ	putative fimbrial usher protein YbgQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/4);  total (7/4)

TAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTAATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTCAGGTCGGGAACTGAAATCAATACACCGATCTATCTGTGGCAAGGACTGACGGATATCTGGTTTT  >  NC_000913/751485‑751891                                                                                                                                                                                     |                                                                                                                                                                                                                                   tAGAGCCTAATTGCGTTAATGGACGAAAAAGATAGGTGCGCTTTATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCAGTGCGTAAGCGCCATGTCCCGGCGTTAATTCCGGCGGTACCGTAGGGAATCAGGTTAGTGCTGCTGCTGCAAGCCTGTGAGCGGTAGCTGGTTGCAGAGAAGTTGTAATTCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAgg                                                                                                                                                                <  1:187913/249‑1 (MQ=255)              cGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTGCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTAATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGAGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAgggggg                                                                                                                                                            <  1:167579/240‑1 (MQ=255)                                                      gACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTAATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGAGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGCATACTAATAtt                                                                                                         >  1:190027/1‑251 (MQ=255)                                                           gTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTAATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGAGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGc                                                                                                                                         <  1:216487/214‑1 (MQ=255)                                                                                                                             gTTAATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGAGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGCGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTCAGGTCGGGAACTGAAATCAATACACCGatct                                   >  1:90629/1‑250 (MQ=255)                                                                                                                                                  ggCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGAGCGGTAGCTGCTGGCACACACGTTGTAATCCCTCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCCGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTCAGGTCGGGAACTGAAATCAATACACCGATCTATCTGTGGCAAGGACTGACgg              >  1:135339/1‑250 (MQ=255)                                                                                                                                                          ggTTAGTGCTGCTGCTGCCGTCCTGTGAGCGGTAGCTGCTGGCAAACAAGTTGTACTCAATCCGGACACCGGCACCACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACgc                                                                                                                                    >  1:39350/1‑124 (MQ=255)                                                                                                                                                              aGTGCTGCTGCTGCCGTCCTGTGAGCGGTAGCTGCTGGCAAACAAGTTGTAAGCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATAGTCAGTTGCTGATTGGCTTGATCGAAAGTGAAGAGCATTTCAGGTCGGGAAGTGAAATCAATACACCGAGCGATCGGTGGCAAGGACTGACGGATATCTGGtttn  <  1:28225/250‑2 (MQ=255)                                                                                                                                                                 gctgctgctgcCGTCCTGTGAGCGGTAGCTGCTGGCAAACCAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGCATACTAATATTCAGTTGc                                                                                                  >  1:229935/1‑151 (MQ=255)                                                                                                                                                                          gcCGTCCTGTGAGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAAAACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGACAGACAGGCTTGCGGAATACTAata                                                                                                           >  1:174401/1‑133 (MQ=255)                                                                                                                                                                             gTCCTGTGAGCGGTAGCTGCTGGCAAACAAGTTGTAATTCATCAGGACACCGGCAACACCTTCTTTACATGTAGAgggggg                                                                                                                                                            >  1:126686/1‑81 (MQ=255)                                                                                                                                                                                     |                                                                                                                                                                                                                                   TAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTAATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTCAGGTCGGGAACTGAAATCAATACACCGATCTATCTGTGGCAAGGACTGACGGATATCTGGTTTT  >  NC_000913/751485‑751891

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