Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,542,643	Δ1 bp	intergenic (+46/‑22)	b4316 → / → fimD	Pil/FimD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,542,640	0	G	.	100.0%	409.8 / NA	89	intergenic (+43/‑25)	b4316/fimD	Pil/FimD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (47/42);  total (47/42)

CTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTT  >  NC_000913/4542565‑4542711                                                                            |                                                                        cTTACCCCGAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCg                                                                         <  2:174327/76‑1 (MQ=255)   tACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCggg                                                                       >  1:1821803/1‑76 (MQ=255)    aCCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCgggg                                                                      >  2:151387/1‑76 (MQ=255)      ccccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCgggtga                                                                    <  1:2983933/76‑4 (MQ=255)        ccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGAtt                                                                 >  2:3258169/1‑76 (MQ=255)         cAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGAttt                                                                <  1:458409/76‑1 (MQ=255)         cAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGAttt                                                                <  2:3009300/76‑1 (MQ=255)         cAAAATGACGGGCGTAAGGGAATAACGCAGGGGGAATTTTTCGCCGGAATAAAAAGAATTGACTGCC‑GGGTGAttt                                                                <  2:3050769/76‑1 (MQ=255)          aaaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGAtttt                                                               <  2:4289535/76‑1 (MQ=255)           aaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTa                                                              >  2:1137875/1‑76 (MQ=255)            aaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTaa                                                             <  2:2800803/76‑1 (MQ=255)            aaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTaa                                                             <  1:2735616/76‑1 (MQ=255)            aaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTaa                                                             <  2:2806412/76‑1 (MQ=255)              tGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGc                                                           <  2:4456521/76‑1 (MQ=255)               gACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGcc                                                          >  2:257246/1‑76 (MQ=255)               gACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGcc                                                          <  2:1496439/76‑1 (MQ=255)               gACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGcc                                                          >  1:1127657/1‑76 (MQ=255)               gACGGGCGTAAGGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGcc                                                          <  2:4848429/76‑1 (MQ=255)                aCGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCg                                                         >  1:2425040/1‑76 (MQ=255)                 cGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCgg                                                        <  2:2425039/76‑1 (MQ=255)                 cGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCgg                                                        >  2:2232669/1‑76 (MQ=255)                  gggCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGAGATTTTAAGCCgga                                                       >  1:66727/1‑76 (MQ=255)                   ggCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCggag                                                      >  1:1791770/1‑76 (MQ=255)                   ggCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCggag                                                      <  2:3954104/76‑1 (MQ=255)                    gCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCggagg                                                     <  2:1357840/76‑1 (MQ=255)                      gTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaa                                                   >  1:1939101/1‑76 (MQ=255)                       tAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaat                                                  >  1:2098057/1‑76 (MQ=255)                        aaTGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaata                                                 >  2:3855986/1‑76 (MQ=255)                         aTGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaataa                                                <  1:4108133/76‑1 (MQ=255)                         aTGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaataa                                                >  1:4452867/1‑76 (MQ=255)                         aTGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaataa                                                <  2:1034893/76‑1 (MQ=255)                           ggAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATg                                              <  1:4128156/76‑1 (MQ=255)                           ggAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATg                                              <  1:1543537/76‑1 (MQ=255)                           ggAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATg                                              <  1:4466834/76‑1 (MQ=255)                            gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGt                                             <  1:4772650/76‑1 (MQ=255)                            gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGt                                             <  1:4737424/76‑1 (MQ=255)                               tAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCat                                          <  1:412747/76‑1 (MQ=255)                                aaCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCata                                         >  2:2178028/1‑76 (MQ=255)                                 aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCatat                                        <  2:3783364/76‑1 (MQ=255)                                 aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCatat                                        >  2:3984914/1‑76 (MQ=255)                                   gCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCt                                      <  2:281508/76‑1 (MQ=255)                                    cAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTg                                     <  2:2723594/76‑1 (MQ=255)                                     aGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGa                                    >  1:3517491/1‑76 (MQ=255)                                     aGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGa                                    >  2:1996939/1‑76 (MQ=255)                                     aGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGa                                    >  2:4814877/1‑76 (MQ=255)                                       ggggAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAt                                  >  1:3347067/1‑76 (MQ=255)                                        gggAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAtt                                 >  2:2874625/1‑76 (MQ=255)                                        gggAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAtt                                 <  1:2781235/76‑1 (MQ=255)                                         ggAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAttt                                >  1:1121459/1‑76 (MQ=255)                                          gAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTa                               >  2:4392431/1‑76 (MQ=255)                                           aaTTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTaa                              <  1:2538653/76‑1 (MQ=255)                                           aaTTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTaa                              >  2:2541334/1‑76 (MQ=255)                                            aTTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAACGTCATATCTGAATTTAAg                             >  2:4643696/1‑76 (MQ=255)                                                ttCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACtt                         >  2:2651934/1‑76 (MQ=255)                                                ttCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACtt                         <  2:3591796/76‑1 (MQ=255)                                                 tCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACttt                        <  2:229954/76‑1 (MQ=255)                                                 tCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACttt                        >  2:3503219/1‑76 (MQ=255)                                                 tCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACttt                        >  1:1395290/1‑76 (MQ=255)                                                 tCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACttt                        >  1:3708351/1‑76 (MQ=255)                                                  cGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTa                       >  2:3095761/1‑76 (MQ=255)                                                  cGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTa                       >  2:991329/1‑76 (MQ=255)                                                   gCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTAc                      >  1:1538218/1‑76 (MQ=255)                                                    ccTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTAcc                     >  2:738653/1‑76 (MQ=255)                                                    ccTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTAcc                     >  2:2429577/1‑76 (MQ=255)                                                    ccTCAATAAAAAGAATTGACTGCC‑GAGAGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTCAcc                     >  2:1298530/1‑76 (MQ=255)                                                     cTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCa                    >  1:1407155/1‑76 (MQ=255)                                                     cTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCa                    >  2:198335/1‑76 (MQ=255)                                                     cTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCa                    <  1:2094629/76‑1 (MQ=255)                                                      tGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAg                   >  2:2487830/1‑76 (MQ=255)                                                      tGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAg                   <  2:3862257/76‑1 (MQ=255)                                                      tGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAg                   >  2:1018292/1‑76 (MQ=255)                                                      tGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAg                   >  1:2390276/1‑76 (MQ=255)                                                        aaTAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCg                 <  2:174192/76‑1 (MQ=255)                                                          tAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGaa               <  1:2254507/76‑1 (MQ=255)                                                          tAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGaa               >  1:4205839/1‑76 (MQ=255)                                                          tAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGaa               <  1:4439762/76‑1 (MQ=255)                                                           aaaaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGaaa              >  1:2922834/1‑76 (MQ=255)                                                           aaaaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGaaa              >  2:1083439/1‑76 (MQ=255)                                                           aaaaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGaaa              <  2:4123978/76‑1 (MQ=255)                                                            aaaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAac             <  1:2861584/76‑1 (MQ=255)                                                                gAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAacacac         <  1:1279349/76‑1 (MQ=255)                                                                 aaTTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAAcacaca        <  2:2148270/76‑1 (MQ=255)                                                                 aaTTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAAcacaca        <  1:3618132/76‑1 (MQ=255)                                                                   ttGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTTATATCTGAATTTAAGACTTTACCAGCGAAACACACAAt      <  1:4552445/76‑1 (MQ=255)                                                                   ttGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAAt      >  1:3227382/1‑76 (MQ=255)                                                                    tGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAACAACACAATg     >  2:4742716/1‑76 (MQ=255)                                                                    tGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATg     <  2:1121458/76‑1 (MQ=255)                                                                     gACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGc    <  2:361179/76‑1 (MQ=255)                                                                     gACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGc    >  1:3663400/1‑76 (MQ=255)                                                                       cTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCtt  >  1:4474284/1‑76 (MQ=255)                                                                       cTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCtt  <  2:2735881/76‑1 (MQ=255)                                                                       cTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCtt  >  2:3988832/1‑76 (MQ=255)                                                                       cTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCtt  >  2:706352/1‑76 (MQ=255)                                                                            |                                                                        CTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTT  >  NC_000913/4542565‑4542711

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