Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,546,637	Δ1 bp	coding (1542/2637 nt)	fimD →	fimbrial usher outer membrane porin protein; FimCD chaperone‑usher

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,546,635	0	A	.	100.0%	252.9 / NA	52	coding (1540/2637 nt)	fimD	fimbrial usher outer membrane porin protein; FimCD chaperone‑usher
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (35/17);  total (35/17)

ACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTC  >  NC_000913/4546486‑4546769                                                                                                                                                      |                                                                                                                                       aCAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGaa                                                                                                                                       <  1:717844/151‑1 (MQ=255)    aaTCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGaattc                                                                                                                                    >  2:638055/1‑148 (MQ=255)     aTCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGaattc                                                                                                                                    >  2:307946/1‑147 (MQ=255)       cGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCAcc                                                                                                                                >  2:603984/1‑151 (MQ=255)        gCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCg                                                                                                                               >  2:682316/1‑151 (MQ=255)        gCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCg                                                                                                                               >  1:680806/1‑151 (MQ=255)        gCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCg                                                                                                                               >  1:584394/1‑151 (MQ=255)           cAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACt                                                                                                                            <  2:94539/151‑1 (MQ=255)           cAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACt                                                                                                                            <  1:821879/151‑1 (MQ=255)               gAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTa                                                                                                                        <  2:65210/151‑1 (MQ=255)                   cAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACaa                                                                                                                     <  1:209282/150‑1 (MQ=255)                       cACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCg                                                                                                                <  2:452728/151‑1 (MQ=255)                        aCGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGc                                                                                                               >  2:78112/1‑151 (MQ=255)                        aCGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGc                                                                                                               >  1:347057/1‑151 (MQ=255)                              atTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTtata                                                                                                          <  2:120036/150‑1 (MQ=255)                              atTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTtata                                                                                                          >  2:351699/1‑150 (MQ=255)                              atTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTtata                                                                                                          >  2:454660/1‑150 (MQ=255)                              atTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTtata                                                                                                          >  1:283488/1‑150 (MQ=255)                              atTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATaa                                                                                                         <  1:682344/151‑1 (MQ=255)                                tCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAAc                                                                                                        >  2:646275/1‑150 (MQ=255)                                 cAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACaa                                                                                                      <  1:202703/151‑1 (MQ=255)                                 cAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACa                                                                                                       >  2:129541/1‑150 (MQ=255)                                         gggTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCggg                                                                                               <  1:663589/150‑1 (MQ=255)                                              aCCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAAtt                                                                                          <  2:287519/150‑1 (MQ=255)                                                  ttATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAc                                                                                     >  1:8972/1‑151 (MQ=255)                                                    aTTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTc                                                                                   <  2:61063/151‑1 (MQ=255)                                                          ccAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGtt                                                                             >  2:592044/1‑151 (MQ=255)                                                           cAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTa                                                                            <  1:499727/151‑1 (MQ=255)                                                                    ttttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCaa                                                                    <  1:833325/150‑1 (MQ=255)                                                                             cGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGgcgc                                                           >  2:198072/1‑150 (MQ=255)                                                                             cGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGgcgc                                                           >  1:879494/1‑150 (MQ=255)                                                                              gCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGcac                                                         >  1:434208/1‑151 (MQ=255)                                                                               cTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGcac                                                         >  2:22833/1‑150 (MQ=255)                                                                                tGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACAt                                                       >  1:107349/1‑151 (MQ=255)                                                                                    acaacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCaa                                                    >  1:38878/1‑150 (MQ=255)                                                                                      aacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAacac                                                 <  2:510315/151‑1 (MQ=255)                                                                                      aacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAacac                                                 >  2:450316/1‑151 (MQ=255)                                                                                           aCAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTAt                                            >  2:419472/1‑151 (MQ=255)                                                                                             aGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTAttt                                          >  1:385161/1‑151 (MQ=255)                                                                                                     gaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAAc                                                                                                    >  1:773345/1‑85 (MQ=255)                                                                                                     gaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAAc                                                                                                    <  2:773554/85‑1 (MQ=255)                                                                                                        tgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGcc                               >  2:269032/1‑151 (MQ=255)                                                                                                        tgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGc                                >  1:424852/1‑150 (MQ=255)                                                                                                            tACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATc                            >  1:516047/1‑150 (MQ=255)                                                                                                               aacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAAc                        >  1:773998/1‑151 (MQ=255)                                                                                                                acaTCGAAACACAGGACGGAGTTATTCATGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACt                       >  2:359472/1‑151 (MQ=255)                                                                                                                acaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACt                       >  2:108561/1‑151 (MQ=255)                                                                                                                   tCGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTAt                    >  1:62120/1‑151 (MQ=255)                                                                                                                    cGAAACACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTAtt                   >  1:456812/1‑151 (MQ=255)                                                                                                                       aaCACAGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTggg                <  1:227456/151‑1 (MQ=255)                                                                                                                          acaGGACGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGa              >  1:397216/1‑150 (MQ=255)                                                                                                                                cGGAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTa       <  1:134494/151‑1 (MQ=255)                                                                                                                                 ggAGTTATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTa       >  2:555534/1‑150 (MQ=255)                                                                                                                                     ttATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTc  >  1:647563/1‑151 (MQ=255)                                                                                                                                     ttATTCAGGTTAAGCCG‑AATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTc  >  1:459771/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                       ACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTC  >  NC_000913/4546486‑4546769

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