Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,451,573	T→C	G4G (GGA→GGG)	yfcR ←	putative fimbrial‑like adhesin protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,451,573	0	T	C	100.0%	174.1 / NA	49	G4G (GGA→GGG)	yfcR	putative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/24);  total (25/24)

AGGGTGGCGGAGCATCGACCAACACGCGCAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACTCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGTCAGCACCGGTACGCTA  >  NC_000913/2451446‑2451695                                                                                                                                |                                                                                                                           aGGGTGGCGGAGCATCGACCAACACGCGCAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGtcat                                                                                                                 <  2:1480402/139‑1 (MQ=255)   ggTGGCGGAGCATCGACCAACACGCGCAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGtcatca                                                                                                               >  2:71645/1‑139 (MQ=255)          gAGCATCGACCAACACGCGCAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACGGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTg                                                                                                        <  2:1273208/139‑1 (MQ=255)          gAGCATCGACCAACACGCGCAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAGCTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTg                                                                                                        <  2:25362/139‑1 (MQ=255)          gAGCATCGACCAACACGCGCAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTg                                                                                                        <  2:1435969/139‑1 (MQ=255)                  aCCAACACGCGCAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTg                                                                                                        <  2:576744/131‑1 (MQ=255)                  aCCAACACGCGCAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTg                                                                                                        >  1:576744/1‑131 (MQ=255)                            gcAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATcc                                                                                      >  2:873045/1‑139 (MQ=255)                                aaTTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCt                                                                                  >  2:1461732/1‑139 (MQ=255)                                          aCGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATc                                                                                       >  1:1011189/1‑124 (MQ=255)                                          aCGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATc                                                                                       <  2:1011189/124‑1 (MQ=255)                                           cGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACcccccc                                                                                                                       >  2:1398105/1‑91 (MQ=255)                                           cGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACcccccc                                                                                                                       <  1:1398105/91‑1 (MQ=255)                                                 tactacCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACGAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTg                                                                                                        <  1:1266519/100‑1 (MQ=255)                                                 tactacCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTg                                                                                                        >  2:1266519/1‑100 (MQ=255)                                                  actacCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGc                                                                <  1:873045/139‑1 (MQ=255)                                                        aGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGgctct                                                                                                 <  2:598227/100‑1 (MQ=255)                                                        aGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGgctct                                                                                                 >  1:598227/1‑100 (MQ=255)                                                        aGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGga                                                          >  1:854157/1‑139 (MQ=255)                                                              aGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCAt                                                                   >  1:1151285/1‑124 (MQ=255)                                                              aGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCAt                                                                   <  2:1151285/124‑1 (MQ=255)                                                                    tAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCa                                              >  1:1357104/1‑139 (MQ=255)                                                                          aCAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCa                                                                              <  2:1185651/101‑1 (MQ=255)                                                                          aCAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCa                                                                              >  1:1185651/1‑101 (MQ=255)                                                                                 cTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCaaa                                 >  1:577721/1‑139 (MQ=255)                                                                                  tGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAAc                                                     >  2:1125452/1‑118 (MQ=255)                                                                                  tGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAAc                                                     <  1:1125452/118‑1 (MQ=255)                                                                                  tGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAg                                >  2:1525772/1‑139 (MQ=255)                                                                                          agcagcCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGAc                        >  2:1137163/1‑139 (MQ=255)                                                                                           gcagcCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACg                       <  1:71645/139‑1 (MQ=255)                                                                                              gcCCTGCGCCCACTACAAACTTTTGACTCCTAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGc                                                                <  2:381995/95‑1 (MQ=255)                                                                                              gcCCTGCGCCCACTACAAACTTTTGACTCATAAcccccc                                                                                                                        >  1:305692/1‑39 (MQ=255)                                                                                              gcCCTGCGCCCACTACAAACTTTTGACTCATAAcccccc                                                                                                                        <  2:305692/39‑1 (MQ=255)                                                                                              gcCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGc                                                                >  1:381995/1‑95 (MQ=255)                                                                                              gcCCTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTaa                    >  1:1452457/1‑139 (MQ=255)                                                                                                ccTGCGCCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGt                  <  2:577721/139‑1 (MQ=255)                                                                                                   gcgcCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGTCAg               >  2:568186/1‑139 (MQ=255)                                                                                                    cgcCCACTACAAACTTTTGACTCATAACCCCCCCCGtcatca                                                                                                               <  1:899881/42‑1 (MQ=255)                                                                                                    cgcCCACTACAAACTTTTGACTCATAACCCCCCCCGtcatca                                                                                                               >  2:899881/1‑42 (MQ=255)                                                                                                    cgcCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCAt                                                                                        <  1:1307323/65‑1 (MQ=255)                                                                                                    cgcCCACTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCAt                                                                                        >  2:1307323/1‑65 (MQ=255)                                                                                                         aCTACAAACTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGTCAGCACCgg         <  1:1137163/139‑1 (MQ=255)                                                                                                                aCTTTTGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGTCAGCACCGGTACGCTa  <  2:1175931/139‑1 (MQ=255)                                                                                                                    ttGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGTGCCAAAGTCGTTGACGTAAGTCAg               <  1:1201890/122‑1 (MQ=255)                                                                                                                    ttGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGTGCCAAAGTCGTTGACGTAAGTCAg               >  2:1201890/1‑122 (MQ=255)                                                                                                                     tGACTCATAACCCCCCCCGTCATCAGCCTTTGCGgct                                                                                                   <  1:459105/37‑1 (MQ=38)                                                                                                                     tGACTCATAACCCCCCCCGTCATCAGCCTTTGCGgct                                                                                                   >  2:459105/1‑37 (MQ=38)                                                                                                                     tGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGgaag                                                        <  1:926420/80‑1 (MQ=255)                                                                                                                     tGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGgaac                                                        >  2:926420/1‑79 (MQ=255)                                                                                                                     tGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGt                            >  2:538420/1‑108 (MQ=255)                                                                                                                     tGACTCATAACCCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGt                            <  1:538420/108‑1 (MQ=255)                                                                                                                           aTAACCCCCCCCGTCATCAGCCTTCGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAAc                                                     >  2:1036299/1‑77 (MQ=255)                                                                                                                           aTAACCCCCCCCGTCATCAGCCTTCGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAAc                                                     <  1:1036299/77‑1 (MQ=255)                                                                                                                                |                                                                                                                           AGGGTGGCGGAGCATCGACCAACACGCGCAGAATTAAATCCACGGGTTTACTACCAGCCATAGCAGATAAACTACAGACACTGCATACCAGCAGCCCTGCGCCCACTACAAACTTTTGACTCATAACTCCCCCCGTCATCAGCCTTTGCGGCTCTGCTCTTCATCCACCTGGCAGCTGTTGCCATTGCACTGGAAGAACAGCGGCATCCGTGCGCCAAAGTCGTTGACGTAAGTCAGCACCGGTACGCTA  >  NC_000913/2451446‑2451695

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