Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,547,852	C→A	G37G (GGC→GGA)	fimF →	type 1 fimbriae minor subunit FimF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,547,852	0	C	A	100.0%	151.4 / NA	46	G37G (GGC→GGA)	fimF	type 1 fimbriae minor subunit FimF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (24/22);  total (24/22)

GGCGTGATGAGAAACAAACCTTTTTATCTTCTGTGCGCTTTTTTGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGCTGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTGTCACCCTGTGGTAATGCCGT  >  NC_000913/4547736‑4547986                                                                                                                     |                                                                                                                                       ggCGTGATGAGAAACAAACCTTTTTATCTTCTGTGCGCTTTTTTGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAAc                                                                                                                 <  1:877246/140‑1 (MQ=255) ggCGTGATGAGAAACAAACCTTTTTATCTTCTGTGCGCTTTTTTGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAAc                                                                                                                 >  2:12817/1‑140 (MQ=255) ggCGTGATGAGAAACAAACCTTTTTATCTTCTGTGCGCTTTTTTGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAAc                                                                                                                 <  2:779622/140‑1 (MQ=255)   cGTGATGAGAAACAAACCTTTTTATCTTCTGTGCGCTTTTTTGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCa                                                                                                               >  2:113628/1‑140 (MQ=255)    gTGATGAGAAACAAACCTTTTTATCTTCTGTGCGCTTTTTTGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCaa                                                                                                              >  1:183204/1‑140 (MQ=255)       atgaGAAACAAACCTTTTTATCTTCTGTGCGCTTTTTTGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAAttt                                                                                                           <  1:418619/140‑1 (MQ=255)        tgaGAAACAAACCTTTTTATCTTCTGTGCGCTTTTTTGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAAtttt                                                                                                          >  1:723216/1‑140 (MQ=255)        tgaGAAACAAACCTTTTTATCTTCTGTGCGCTTTTTTGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAAtttt                                                                                                          >  2:777712/1‑140 (MQ=255)          agaAACAAACCTTTTTATCTTCTGTGCGCTTTTTTGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTAc                                                                                                        >  2:396062/1‑140 (MQ=255)                  acCTTTTTATCTTCTGTGCGCTTTTTTGTGGCTGGCGGTGAGTCACGCTTTTGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGGGGCCGCTGAATCAACCAATTTTACTGTTGATc                                                                                                >  2:383783/1‑140 (MQ=255)                              tctGTGCGCTTTTTTGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAAcg                                                                                    >  1:183015/1‑140 (MQ=255)                                 gtgCGCTTTTTTGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACgcgg                                                                                 >  2:666923/1‑140 (MQ=255)                                  tgCGCTTTTTTGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGTATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACgcggc                                                                                >  2:768145/1‑140 (MQ=255)                                    cgcTTTTTTGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGa                                                                              >  2:16719/1‑140 (MQ=255)                                      cTTTTTTGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAg                                                                            <  1:16719/140‑1 (MQ=255)                                      cTTTTTTGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAg                                                                            <  2:394265/140‑1 (MQ=255)                                         ttttGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCaa                                                                         >  1:721895/1‑140 (MQ=255)                                          tttGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAAt                                                                        <  1:396062/140‑1 (MQ=255)                                                   ggcGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAacaaca                                                               >  2:914833/1‑140 (MQ=255)                                                   ggcGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAacaaca                                                               >  1:597968/1‑140 (MQ=255)                                                           gTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGcgcg                                                       <  2:629420/140‑1 (MQ=255)                                                             cACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGcgac                                                     <  2:597968/140‑1 (MQ=255)                                                                gCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCgacgac                                                  <  1:364527/140‑1 (MQ=255)                                                                gCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCgacgac                                                  <  1:777712/140‑1 (MQ=255)                                                                     ggCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATTACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCtg                                             <  1:634709/140‑1 (MQ=255)                                                                            gATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTcc                                      >  1:435999/1‑140 (MQ=255)                                                                                gCACGATTACTATCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCAttt                                  <  1:117078/140‑1 (MQ=255)                                                                                           aTCCGCGGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCt                       >  2:573053/1‑140 (MQ=255)                                                                                              cgcgGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTGTc                    <  2:336343/140‑1 (MQ=255)                                                                                               gcgGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTGTCa                   <  1:573053/140‑1 (MQ=255)                                                                                               gcgGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTGTCa                   <  1:560638/140‑1 (MQ=255)                                                                                                cgGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTGTCAc                  >  1:932571/1‑140 (MQ=255)                                                                                                cgGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTGTCAc                  >  1:205133/1‑140 (MQ=255)                                                                                                cgGCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTGTCAc                  >  1:135209/1‑140 (MQ=255)                                                                                                 ggCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTg                      >  1:308937/1‑135 (MQ=255)                                                                                                 ggCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTg                      <  2:308937/135‑1 (MQ=255)                                                                                                 ggCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTGTCAcc                 >  1:419676/1‑140 (MQ=255)                                                                                                  gCTATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTGTCAccc                <  2:183204/140‑1 (MQ=255)                                                                                                    tATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTGTCACCCtg              <  1:12817/140‑1 (MQ=255)                                                                                                    tATGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTGTCACCCtg              <  1:914833/140‑1 (MQ=255)                                                                                                      tGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTAttt                           >  1:890739/1‑125 (MQ=255)                                                                                                      tGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTAttt                           <  2:890739/125‑1 (MQ=255)                                                                                                      tGTCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTGTCACCCtgtg            >  1:295529/1‑140 (MQ=255)                                                                                                        tCAGGGATAACGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTGTCACCCTGTGGt          <  2:721895/140‑1 (MQ=255)                                                                                                                aaCGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTGTCACCCTGTGGTAATGCCGt  >  2:216066/1‑140 (MQ=255)                                                                                                                aaCGGATGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTGTCACCCTGTGGTAATGCCGt  <  1:325077/140‑1 (MQ=255)                                                                                                                     |                                                                                                                                       GGCGTGATGAGAAACAAACCTTTTTATCTTCTGTGCGCTTTTTTGTGGCTGGCGGTGAGTCACGCTTTGGCTGCGGATAGCACGATTACTATCCGCGGCTATGTCAGGGATAACGGCTGTAGTGTGGCCGCTGAATCAACCAATTTTACTGTTGATCTGATGGAAAACGCGGCGAAGCAATTTAACAACATTGGCGCGACGACTCCTGTTGTTCCATTTCGTATTTTGCTGTCACCCTGTGGTAATGCCGT  >  NC_000913/4547736‑4547986

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