Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,188,087	T→C	T110A (ACG→GCG)	yehA ←	putative fimbrial adhesin YehA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,188,087	0	T	C	100.0%	62.6 / NA	18	T110A (ACG→GCG)	yehA	putative fimbrial adhesin YehA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/10);  total (8/10)

TTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTT  >  NC_000913/2187864‑2188290                                                                                                                                                                                                                                |                                                                                                                                                                                                            ttAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGt                                                                                                                                                                                  >  1:189547/1‑251 (MQ=255)  tAAAGGTATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTc                                                                                                                                                                                 >  1:151982/1‑251 (MQ=255)       gAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGatt                                                                                                                                                                            <  1:269920/251‑1 (MQ=255)                            gCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCCCCGAATAAACACCGTAGAAATTGTACTGCCGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCACCCCCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGTTGTCACCGCAACTGAATAACCACTGGTGAAACTTGTCTGATTATTTGACATCGTGATtttgtt                                                                                                                                                       <  2:449079/251‑1 (MQ=255)                                   cATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCCCCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGAGTAACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACACCGTGCTGTTGTTTTTATTg                                                                                                                                                >  2:287942/1‑251 (MQ=255)                                                       caTGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAaatt                                                                                                                            >  1:9649/1‑251 (MQ=255)                                                                                                    aCTGCGGCACTACCGCGTGTCGATGAAGAAGCACTGGATCATGTCGCAACAACATTAACCATAACTGAACCACAGCCCCCCGCAGCAACATTCCCTTTTGCATTACTGGCCGCAGGTGTCACCACAACTGAATACCCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGATCCATTCAAACTGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCTGACgg                                                                               >  2:219218/1‑251 (MQ=255)                                                                                                       gCGGCACTCCCGCGGGTCGATGATGAAGCACTGGATAATGGCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGAGTATGTTACAACGTGATGTTGTTTTTATTGGTCCATTCAAAAGGTAAATTTAATTTTATGGTATGGTTATAAAGTCCAACAATAATATCCGACggg                                                                              <  1:300801/249‑1 (MQ=255)                                                                                                        cGGCCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCCCTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATCAAGTCCAACAATAATATCCGACGGGCTc                                                                           <  1:295482/251‑1 (MQ=255)                                                                                                                        ggatgatGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTAt                                                           <  2:31417/250‑1 (MQ=255)                                                                                                                                    aCTGGATAATGTCGCAACACCATTAATCATCACTGCACAGCCGCCCCCCGCAGAAACATTCCCTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGAtg                                                                                                                                                            >  2:491287/1‑142 (MQ=255)                                                                                                                                    aCTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGAtg                                                                                                                                                            <  1:491287/142‑1 (MQ=255)                                                                                                                                        gATAATGTCGCAACACCATTAATAATCACTGAACCGCCGACCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTATCACCGCAACTGACTAACCACTGGTGAAAATGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATg                                            >  1:180373/1‑250 (MQ=255)                                                                                                                                                  cAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGACGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATGTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTa                                 <  1:385245/251‑1 (MQ=255)                                                                                                                                                                        aCCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCa           <  1:287942/251‑1 (MQ=255)                                                                                                                                                                                 ccccGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGttt  <  2:189547/251‑1 (MQ=255)                                                                                                                                                                                 ccccGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGttt  <  2:113715/251‑1 (MQ=255)                                                                                                                                                                                 ccccGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGttt  >  2:466433/1‑251 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                            TTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTT  >  NC_000913/2187864‑2188290

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