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%	76.0 / NA	25	T110A (ACG→GCG)	yehA	putative fimbrial adhesin YehA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/10);  total (15/10)

TTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGA‑CTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAAC  >  NC_000913/2187846‑2188315                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     ttGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGA‑CTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGCTGAAGCACTGGATAATGTCGCAAAACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTgtat                                                                                                                                                                                                                              >  2:346003/1‑247 (MQ=255)               cAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGA‑CTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCAGTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGaaa                                                                                                                                                                                                                <  1:324503/250‑1 (MQ=255)                           aTAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGA‑CTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAAGGTCGCAACACCATTAATCATCACTGAACCGCCGCCACCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGattat                                                                                                                                                                                                   <  1:424030/251‑1 (MQ=255)                                            gtgCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGA‑CTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCACCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTAttttcgttgcaggatgtc                                                                                                                                                                                  >  1:237988/1‑236 (MQ=255)                                                cAATGCATTTCTGTAGCTATTTATACATGCATCCCATGA‑CTTGCCACCTAATAAACCCAGTAGAAATTGTACTGCGGCTCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGGCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTTATGTTGtttt                                                                                                                                                                              <  2:375101/251‑1 (MQ=255)                                                cAATGCATTTCTGTAGCTATTTACACATGCATCCCATGA‑CTTGCCACCTAATAAACACAGTGGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCACCCGCAGAAACATTCACTTCTGCATTACTGGCCGCAGGTGTCACCGCACCTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTACTGTTGtttt                                                                                                                                                                              >  2:1108/1‑251 (MQ=255)                                                       tttCTGTAGCTATTTACACATGCATCCCATGA‑CTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGt                                                                                                                                                                       >  1:322305/1‑251 (MQ=255)                                                                 tATTTACACATGCATCCCATGA‑CTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCaaa                                                                                                                                                              >  2:240137/1‑250 (MQ=255)                                                                     tACACATGCATCCCATGA‑CTTGCCACCTAATAAACACAGTCGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGt                                                                                                                                                          >  1:405063/1‑250 (MQ=255)                                                                         caTGCATCCCATGATGTTGC‑ACCTAATAAACACAGTAGAAATGGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATATGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAaatt                                                                                                                                                     <  1:12535/251‑1 (MQ=255)                                                                                        cTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTAt                                                                                                                                       >  1:455020/1‑251 (MQ=255)                                                                                                                    tGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCg                                                                                                            <  1:19353/250‑1 (MQ=255)                                                                                                                       aCTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACgg                                                                                                        >  1:192398/1‑251 (MQ=255)                                                                                                                              ccttccctcttgtagatgatGAAGCACTGGTAAATGTCGCAACCCCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGTTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATa                                                                                                 <  1:207796/237‑1 (MQ=255)                                                                                                                                             atgaggAAGCACTGACTAATGTCGAAACACCATTAAACATCACTGAACCGACGCCCCCAGCAGAAACATTAACTTTTGCATTACTGGCCGCACGTGTCACCGCGACTGCATAACCACTGGTAAAACTGGTCTGCTTATTTGTCAACGTGATGTTGTTTTTATTGGTCTATTCATATTGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAAAAATAATATCCGACGGCCTCATGGGGGTGagatacttt                                                                                  >  2:36694/1‑243 (MQ=255)                                                                                                                                                                                              gccgccCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGGGAAACTGGTCTGATTATTTGACAACGTTATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCAttt                                 <  2:455020/251‑1 (MQ=255)                                                                                                                                                                                                    ccccGCAGAAACATTCCCTTTTGCATTACTTGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGttt                           >  2:50416/1‑251 (MQ=255)                                                                                                                                                                                                    ccccGCAGAAACATTCACTTTTGCATTACTGTCCGCAGGTGTCACCGCAACTTAATAACCACTGGTGAAACTGGTCTGATTATTTGACACCGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTCCGCATGACGAATCATTAAAATTATGCTCAAATGTCGCCATTTGTtttt                           >  2:454815/1‑251 (MQ=255)                                                                                                                                                                                                        gCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGt                       <  1:240137/251‑1 (MQ=255)                                                                                                                                                                                                          aGAAACATTCACTTTTGCATTACTGGCCGCAGGAGTCACCGCAACTGAATAACCACTGGAGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTACAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGCCGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTa                      >  2:177141/1‑250 (MQ=255)                                                                                                                                                                                                            aaaCATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGt                   <  1:33550/251‑1 (MQ=255)                                                                                                                                                                                                             aaCATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTAAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTg                  >  1:243970/1‑251 (MQ=255)                                                                                                                                                                                                                    aCTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATc           >  1:347860/1‑251 (MQ=255)                                                                                                                                                                                                                         tGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATcccact      >  1:100929/1‑248 (MQ=255)                                                                                                                                                                                                                              tACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTTAGTAAGTGCAAAATCCCCCTGaac  <  1:1108/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     TTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGA‑CTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAAC  >  NC_000913/2187846‑2188315

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