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

AAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAAC  >  NC_000913/2187866‑2188318                                                                                                                                                                                                                              |                                                                                                                                                                                                                                        aaaGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCt                                                                                                                                                                                                            <  1:40635/251‑1 (MQ=255)       aTAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGattat                                                                                                                                                                                                      <  2:113683/251‑1 (MQ=255)         aGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTAttt                                                                                                                                                                                                    <  2:249042/251‑1 (MQ=255)                  cATTTTGTGCCACTGCATTAATGGAGCTATTTACACATGCATCCCATGAATTGACACCTAATAAACACAGTAGAAATTGTACTGCGACACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTAACTTTTGCATTACTGGCCGCCGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTg                                                                                                                                                                                           >  2:209492/1‑251 (MQ=255)                    ttttGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGAt                                                                                                                                                                                         >  1:279910/1‑251 (MQ=255)                       tgtgCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGAtgt                                                                                                                                                                                       <  1:73841/250‑1 (MQ=255)                                 cATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTg                                                                                                                                                                            <  2:195718/251‑1 (MQ=255)                                        gTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCAtt                                                                                                                                                                     <  1:180375/251‑1 (MQ=255)                                            cTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAACCATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCaa                                                                                                                                                                  >  2:230814/1‑250 (MQ=255)                                                                                   cacaGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCcaac                                                                                                                          >  1:238439/1‑251 (MQ=255)                                                                                                            tCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTAt                                                                                                                                          >  1:165457/1‑210 (MQ=255)                                                                                                            tCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTAt                                                                                                                                          <  2:165457/210‑1 (MQ=255)                                                                                                             cccGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAgggg                                                                                                >  1:285466/1‑251 (MQ=255)                                                                                                                    gtCGATGATGAAGCACTGTATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAg                                                                                          <  2:326150/250‑1 (MQ=255)                                                                                                                         tgatgaAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTg                                                                                    >  1:36930/1‑251 (MQ=255)                                                                                                                                 cACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCGCATAGGGGTGATAGTATTTGTACTTACg                                                                            <  2:351685/251‑1 (MQ=255)                                                                                                                                  aCTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGc                                                                           >  1:127340/1‑251 (MQ=255)                                                                                                                                       aTAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGAc                                                                      >  2:123607/1‑251 (MQ=255)                                                                                                                                                                                         aCATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGc                    <  1:182602/251‑1 (MQ=255)                                                                                                                                                                                           aTTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCa                   >  1:284310/1‑250 (MQ=255)                                                                                                                                                                                                    tGCATTACTGGCCGCAGCTGTCACCGCCACTGACTCACCACTGGTGAAACTGGTCTGATTCTTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAACTGTAAATTTAATTTTATGGTATCGTTATAACGTCCAACAATAATATCCGACGGGCTCATAGGGGCGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATCTCGCCATTTGTGTTTGAGAAAGTGCAAAATCCCCCt         >  1:106686/1‑251 (MQ=255)                                                                                                                                                                                                             ggCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGaacaac  <  2:36930/249‑1 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                                        AAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAAC  >  NC_000913/2187866‑2188318

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