Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,548,234	A→G	I165V (ATC→GTC)	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,548,234	0	A	G	100.0%	202.6 / NA	56	I165V (ATC→GTC)	fimF	type 1 fimbriae minor subunit FimF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/33);  total (23/33)

GAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATC  >  NC_000913/4548099‑4548368                                                                                                                                        |                                                                                                                                       gAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCaa                                                                                                                                    >  1:315304/1‑140 (MQ=255)      gcaAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCa                                                                                                                               <  2:511769/140‑1 (MQ=255)      gcaAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCa                                                                                                                               <  2:293189/140‑1 (MQ=255)      gcaAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACACCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCa                                                                                                                               >  1:237591/1‑140 (MQ=255)        aaaaTCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACg                                                                                                                             <  2:553901/140‑1 (MQ=255)          aaTCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGc                                                                                                                           <  2:585667/140‑1 (MQ=255)          aaTCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGCGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGc                                                                                                                           <  2:111538/140‑1 (MQ=255)              aaaTACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTAcc                                                                                                                       <  1:298710/140‑1 (MQ=255)                   ccccTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCAc                                                                                                                  <  1:158241/140‑1 (MQ=255)                         aaTGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGa                                                                                                            >  1:427240/1‑140 (MQ=255)                           tGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAat                                                                                                          <  1:149654/140‑1 (MQ=255)                           tGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCATGAat                                                                                                          <  2:545417/140‑1 (MQ=255)                            gCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAata                                                                                                         <  2:466353/140‑1 (MQ=255)                              tCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATc                                                                                                       <  2:164774/140‑1 (MQ=255)                              tCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATc                                                                                                       <  2:628112/140‑1 (MQ=255)                               ccATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCa                                                                                                      <  1:85526/140‑1 (MQ=255)                               ccATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCa                                                                                                      <  2:315304/140‑1 (MQ=255)                               ccATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCa                                                                                                      <  1:613399/140‑1 (MQ=255)                                cATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAg                                                                                                     >  1:324749/1‑140 (MQ=255)                                 aTCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGt                                                                                                    >  1:246649/1‑140 (MQ=255)                                  tCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTa                                                                                                   >  2:434089/1‑140 (MQ=255)                                     tCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACt                                                                                                >  1:448620/1‑140 (MQ=255)                                       cgcgCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTgg                                                                                              >  2:254139/1‑140 (MQ=255)                                        gcgcTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGga                                                                                             <  1:611538/140‑1 (MQ=255)                                         cgcTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGgag                                                                                            <  2:616182/140‑1 (MQ=255)                                          gcTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGgaga                                                                                           <  2:259469/140‑1 (MQ=255)                                             ttCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGc                                                                                        <  2:437612/140‑1 (MQ=255)                                              tCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCt                                                                                       >  2:546689/1‑140 (MQ=255)                                                     cgacCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGaaa                                                                                <  2:693450/140‑1 (MQ=255)                                                        cccTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATgg                                                                             >  1:80048/1‑140 (MQ=255)                                                        cccTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATgg                                                                             >  1:555999/1‑140 (MQ=255)                                                        cccTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGC‑TGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGt                                                                            >  1:357304/1‑140 (MQ=255)                                                                ccGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACg                                                                     <  1:718009/140‑1 (MQ=255)                                                                ccGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACg                                                                     <  2:448620/140‑1 (MQ=255)                                                                ccGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACg                                                                     <  2:346694/140‑1 (MQ=255)                                                                 cGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGt                                                                    <  2:217391/140‑1 (MQ=255)                                                                   ggTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTgg                                                                  >  2:634208/1‑140 (MQ=255)                                                                   ggTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTgg                                                                  >  1:52255/1‑140 (MQ=255)                                                                     tAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGgt                                                                <  2:427240/140‑1 (MQ=255)                                                                               aCGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTggc                                                      <  1:254139/140‑1 (MQ=255)                                                                                 gCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTggcgg                                                    <  1:46637/140‑1 (MQ=255)                                                                                           tACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGgc                                          >  2:475809/1‑140 (MQ=255)                                                                                           tACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGgc                                          <  2:554971/140‑1 (MQ=255)                                                                                                       aTGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGc                              <  2:324749/140‑1 (MQ=255)                                                                                                        tGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGcg                             <  1:301428/140‑1 (MQ=255)                                                                                                        tGGCGACACAGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGcg                             <  2:690377/140‑1 (MQ=255)                                                                                                              cacaGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATa                       >  2:690722/1‑140 (MQ=255)                                                                                                              cacaGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATa                       >  2:769968/1‑140 (MQ=255)                                                                                                                 aGGTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAg                    >  1:282658/1‑140 (MQ=255)                                                                                                                   gTGCCTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGc                  >  1:684251/1‑140 (MQ=255)                                                                                                                      ccTGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCTACGATACAGGCAGc               >  2:490849/1‑140 (MQ=255)                                                                                                                        tGTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCg             <  2:52255/140‑1 (MQ=255)                                                                                                                         gTCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGa            >  1:739908/1‑140 (MQ=255)                                                                                                                          tCACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGAt           >  1:688033/1‑140 (MQ=255)                                                                                                                           cACTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATg          >  1:294703/1‑140 (MQ=255)                                                                                                                             cTGCGGGGCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTc        >  2:523462/1‑140 (MQ=255)                                                                                                                                   ggCATGTCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATc  <  1:634208/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                       GAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATC  >  NC_000913/4548099‑4548368

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