Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,549,889	T→C	intergenic (+179/+64)	fimH → / ← gntP	type 1 fimbriae D‑mannose specific adhesin/fructuronate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,549,889	0	T	C	100.0%	90.4 / NA	26	intergenic (+179/+64)	fimH/gntP	type 1 fimbriae D‑mannose specific adhesin/fructuronate transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/10);  total (16/10)

AGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCA  >  NC_000913/4549756‑4550017                                                                                                                                      |                                                                                                                                 aGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGc                                                                                                                            >  2:224707/1‑140 (MQ=255) aGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGc                                                                                                                            >  1:127448/1‑140 (MQ=255)  gCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGca                                                                                                                           <  1:224707/140‑1 (MQ=255)        aaCCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGCACAtttt                                                                                                                     <  2:197030/140‑1 (MQ=255)                   tAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTCGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTGtc                                                                                                          >  2:278285/1‑140 (MQ=255)                               ggCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGc                                                                                              <  2:264451/140‑1 (MQ=255)                               ggCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGc                                                                                              <  2:127448/140‑1 (MQ=255)                                 ccAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCt                                                                                            <  1:82772/140‑1 (MQ=255)                                 ccAGATAAGACGCGCCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTg                                                                                                            <  1:126485/124‑1 (MQ=255)                                 ccAGATAAGACGCGCCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTg                                                                                                            >  2:126485/1‑124 (MQ=255)                                     aTAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGc                                                                                              <  2:112889/134‑1 (MQ=255)                                     aTAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGc                                                                                              >  1:112889/1‑134 (MQ=255)                                     aTAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCa                                                                                        >  1:347970/1‑140 (MQ=255)                                            gcgTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAg                                                                                 >  1:37002/1‑140 (MQ=255)                                                   cGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTgg                                                                          >  1:264600/1‑140 (MQ=255)                                                     tCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgtg                                                                        <  2:264600/140‑1 (MQ=255)                                                                   aaCAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTAc                                                          <  2:37002/140‑1 (MQ=255)                                                                     cAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCa                                                        <  2:347970/140‑1 (MQ=255)                                                                                           ttCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACaataat                                  >  2:125246/1‑140 (MQ=255)                                                                                                          cccTCTCCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGa                   >  1:171895/1‑140 (MQ=255)                                                                                                               tcCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGa              >  2:158456/1‑140 (MQ=255)                                                                                                               tcCCCGGAGGGGCGAGGGGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGa              >  1:25461/1‑140 (MQ=255)                                                                                                                     gAGGGGCGAGGCGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGcc        >  1:175982/1‑139 (MQ=255)                                                                                                                         ggCGAGGGGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcag    >  1:324043/1‑140 (MQ=255)                                                                                                                           cGAGGGGACTGCCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagca  >  2:348761/1‑140 (MQ=255)                                                                                                                           cGAGGGGACTGCCCGGACACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagca  >  2:77576/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                 AGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCA  >  NC_000913/4549756‑4550017

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