Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,187,903	C→T	R171K (AGA→AAA)	yehA ←	putative fimbrial adhesin YehA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,187,903	0	C	T	100.0%	86.5 / NA	25	R171K (AGA→AAA)	yehA	putative fimbrial adhesin YehA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/12);  total (13/12)

TTCCGTGATGTTTGATTTTGTTGACCTAATAAATTTTTACAACGCAGAATAATATCCCCTTGTTTACTATTTATTGTTGCTTTGTTACCTGTGGCTGGTAATTGTGAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTT  >  NC_000913/2187678‑2188138                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                            ttCCGTGATGTTTGATTTTGTTGACCTAATAAATTTTTACAACGCAGAATAATATCCCCTTGTTTACTATTTATTGTTGCTTTGTTACCTGTGGCTGGTAATTGTGAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATccc                                                                                                                                                                                                                     <  2:167862/250‑1 (MQ=255)  tcagTGATGTTTGATTTTGTTGACCTAATAAATTTTTACAACGCAGAATAATATCCCCTTGTTTACTATTTATTGTTGCTTTGTTACCTGTGGCTGGTAATTGTGAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCAt                                                                                                                                                                                                                   >  2:131472/4‑251 (MQ=255)              gATTTTGTTGACCTAATAAATTTTTACAACGCAGAATAATATCCCCTTGTTTACTATTTATTGTTGCTTTGTTACCTGTGGCTGGTAATTGTGAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCt                                                                                                                                                                                                       >  1:352880/1‑251 (MQ=255)                      tgACCTAATAAATTTTTACAACGCAGAATAATATCCCCTTGTTTACTATTTATTGTTGCTTTGTTACCTGTGGCTGGTAATTGTGAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATaaaa                                                                                                                                                                                                >  2:268117/1‑249 (MQ=255)                                       aCAACGCAGAATAATATCCCCTTGTTTACTATTTATTGTTGCTTTGTTACCTGTGGCTGGTAATTGTGAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACt                                                                                                                                                                               <  1:365714/250‑1 (MQ=255)                                               gAATAATATCCCCTTGTTTACTATTTATTGTTGCTTTGTTACCTGTGGCTGGTAATTGTGAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTc                                                                                                                                                                      >  1:47416/1‑251 (MQ=255)                                                          ccTTGTTTACTATTTATTGTTGCTTTGTTACCTGTGGCTGGTAATTGTGAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCga                                                                                                                                                           <  2:121646/251‑1 (MQ=255)                                                                ttACTATTTATTGTTGCTTTGTTACCTGTGGCTGGTAATTGTGAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGatgatg                                                                                                                                                      <  1:300032/250‑1 (MQ=255)                                                                  tcTATTTATTGTTGCTTTGTTACCTGTGGCTGGTAATTGTGAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGTAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGaa                                                                                                                                                    <  1:227857/249‑1 (MQ=255)                                                                               gCTTTGTTACCTGTGGCTGGTAATTGTGAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCg                                                                                                                                                                            <  1:188557/213‑1 (MQ=255)                                                                               gCTTTGTTACCTGTGGCTGGTAATTGTGAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCg                                                                                                                                                                            >  2:188557/1‑213 (MQ=255)                                                                                    gTTACCTGTGGCTGGTAATTGTGAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa                                                                                                                                  >  1:80844/1‑250 (MQ=255)                                                                                          tgtgGCTGGTAATTGTGAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTacac                                                                                                                                                                                                                              <  1:284807/152‑1 (MQ=255)                                                                                          tgtgGCTGGTAATTGAGAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTacac                                                                                                                                                                                                                              >  2:284807/1‑152 (MQ=255)                                                                                                        gtgAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGa                                                                                                              >  1:127124/1‑250 (MQ=255)                                                                                                          gAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAAcc                                                                                                           <  2:127124/251‑1 (MQ=255)                                                                                                                               gTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACAtt                                                                                      >  1:360936/1‑251 (MQ=255)                                                                                                                                                       cAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAgg                                                              <  1:240843/251‑1 (MQ=255)                                                                                                                                                               ggTTATCGGGTTGTATGACAATGTCCTATTAAAGGAATAGACACCTGCATTTTTTGCCATTGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTCCCACCTCATAACCACATTATAAATTGTACTGCGCCTCTCCCGCGTGTCGATGATGAAGCACTGGCTAATGTCGCAACACCATTAATCATCTCTGTACCGCCGCCCCCCGCAGAAACATTCGCTTTTGCAGTACTGGCCGCatttttcacc                                                       <  1:373981/250‑10 (MQ=255)                                                                                                                                                                                                 gAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTATCTATTTACACATGCATCCCATGACTTGCCCCCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACCCCATTAATCATCGCTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTgg                           <  2:414243/244‑1 (MQ=255)                                                                                                                                                                                                 gAATAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTgg                           >  1:414243/1‑244 (MQ=255)                                                                                                                                                                                                  aaTAGACACCTTCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGCATAACCACTGGTGAAACTGGTCTGatt                    >  2:124679/1‑250 (MQ=255)                                                                                                                                                                                                   aTAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGCCTGattat                  >  2:241338/1‑251 (MQ=255)                                                                                                                                                                                                   aTAGACACCTGCATTTTGTGCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACAAAGTAGAAATTGTACTGCGGCACTCCCGCGAGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAAACGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACAACTGGTGAAACTGGTCTGattat                  >  1:441730/1‑251 (MQ=255)                                                                                                                                                                                                                    gtgCCAATGCATTTTTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGAtgtt  <  2:47416/250‑1 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                            TTCCGTGATGTTTGATTTTGTTGACCTAATAAATTTTTACAACGCAGAATAATATCCCCTTGTTTACTATTTATTGTTGCTTTGTTACCTGTGGCTGGTAATTGTGAAACGGGAATACTGTCTAAAGTAATTAATAAATCGTCCGGTTTGCAGGTTGTGGTTATCGGGTTGTATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTT  >  NC_000913/2187678‑2188138

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