Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,541,261	C→T	A102V (GCT→GTT)	fimB →	Type 1 fimbriae regulatory protein FimB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,541,261	0	C	T	97.0%	108.1 / NA	33	A102V (GCT→GTT)	fimB	Type 1 fimbriae regulatory protein FimB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/11);  minor base A (1/0);  total (22/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTAT  >  NC_000913/4541041‑4541484                                                                                                                                                                                                                             |                                                                                                                                                                                                                                ccTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTaa                                                                                                                                                                                                   <  2:136418/251‑1 (MQ=255)     aTGCAGCACGTAATTATTGTCTGACTTTGCTTCGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTCAgggg                                                                                                                                                                                               >  1:153019/1‑251 (MQ=255)           cacgacattattGCCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGCCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCATAGCGTATATATATCCATCGATTAAAAACAGGCTTTTCCACAACGCACCCGCTATTGAATAAAGACGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATccc                                                                                                                                                                                         >  2:185090/7‑250 (MQ=255)                           gACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGt                                                                                                                                                                         >  1:251808/1‑251 (MQ=255)                                  ctttgTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTAcc                                                                                                                                                                   >  1:353964/1‑250 (MQ=255)                                   tttgTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCat                                                                                                                                                                 >  1:258162/1‑251 (MQ=255)                                   tttgTTTTATCCATGGTTTCCGGGCGAGTGAAATTTGTCGCTTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCat                                                                                                                                                                 >  2:321614/1‑251 (MQ=255)                                         ttATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCAtattat                                                                                                                                                            <  1:295962/250‑1 (MQ=255)                                           aTTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTGTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTAtct                                                                                                                                                          <  1:295947/250‑1 (MQ=255)                                              cATGGTTTCCGGGCGAGTGAACTTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAACGTGTATATATATCGATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGAc                                                                                                                                                      >  2:36596/1‑251 (MQ=255)                                                         ggCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCggtggt                                                                                                                                            >  1:12364/1‑250 (MQ=255)                                                             aGTGAAATGTTTCGATTGAGGATTTCGGATATGGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACGTCCGGTGGTAATGc                                                                                                                                       <  2:254216/251‑1 (MQ=255)                                                               tGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAGAAAAAGGCTTTTCAACAACGCACCCACTATTGAATAAAGAAGTACAAGATTTATAAAACTGTTTGAGTATCCGTACTTCGTACACGCATGATGAGAGCAAGTGGGTATTTTTATCACGTAAGTGGAATCCGCTTTTGAGACAACAGGTTTACCATATTATCACGACTTCAGTTGGTTATGc                                                                                                                                       >  2:253506/1‑249 (MQ=255)                                                               tGAAATTTGTAGATTGAGGATGTGGGGTACTGATATTCAAGCAAAGTGTATATATATCAATAGATTACAACAAGGATTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGAATTTTTATCACGTAAGGGGATTCCGCTTTCTCGGCAACAGATTTACCATATTATCTCGACTTCCGGTGGTAATGCCg                                                                                                                                     >  2:412611/1‑251 (MQ=255)                                                                                ggATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGAtt                                                                                                                    >  2:104642/1‑251 (MQ=255)                                                                                 gATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTc                                                                                                                   >  2:145645/1‑251 (MQ=255)                                                                                       ggATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAACGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGCGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTCAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCAATGGAGATTCATCCGc                                                                                                             >  1:34242/1‑251 (MQ=255)                                                                                                       gCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTTAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTc                                                                                             <  2:258162/251‑1 (MQ=255)                                                                                                        cAAAGTGTATATATAGCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTTCTTCGTACCCGCATGTTGAGAGCGATTGGGGATTTTTATCTCGTAAGTGGAATCCGCTTTCTCGGCAACAGTTTTACCATAGTATCTCGACTTCCGGTGTTAATGCCGTGTTGTCGCTGGAGATTCATCCGCAAATGTTACGCCATTCg                                                                                            <  1:415347/251‑1 (MQ=255)                                                                                                          aaGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCgtg                                                                                          >  1:65873/1‑251 (MQ=255)                                                                                                             tgtATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTgg                                                                                       >  2:86701/1‑251 (MQ=255)                                                                                                                     tatCCATCGATTAAAAAAAGGCTTTTAAAGAACGAACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTAACCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCGCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCtt                                                                               >  2:299094/1‑251 (MQ=255)                                                                                                                         cATCGATTAAAAAAAGGATTTTCAACAACGCACCCGCTATTGAATACAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGc                                                                           >  1:325773/1‑251 (MQ=255)                                                                                                                                      aaggttattCAACAACGCACCCGCTATTGAATAAAGAAATACATGCTTGAAAAAACTGGTGTAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAAGCCTCTTTCTCGGCAACAGTTTTACCATATTATCTCGTCTTCCGGTGGTGATGCCTGGTTGTCACTGGAGATTCATCCGCGCATGTTACCCCCTTCGTGTGGTTTTGATTTTGCGAATATGGGAATa                                                              <  2:125243/244‑1 (MQ=255)                                                                                                                                       aGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAAt                                                               >  1:4759/1‑249 (MQ=255)                                                                                                                                              tCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAAAAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATAcgcg                                                      >  2:445450/1‑251 (MQ=255)                                                                                                                                                aacaacGCACCCGCTATTGTATAAAGAAGATCAGGCTTTAAAAAACTGGTTAAGTCACCGTACTTAGTACCCGCATGTTGAGAGCGAGTGGGTATTTT‑ATCACGTAAGGGGAATCCGATTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTTATGCCGGGGTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGAGGTTTTGCTTTGGCGAATATGGGAATAGATCCGCGACt                                                   >  2:120313/1‑251 (MQ=255)                                                                                                                                                       cACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATcc                                              >  1:431589/1‑250 (MQ=255)                                                                                                                                                                 tGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCt                                    <  1:130301/250‑1 (MQ=255)                                                                                                                                                                                     ttAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTc                <  1:86701/250‑1 (MQ=255)                                                                                                                                                                                         aaatcTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGTGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACt           <  1:461546/247‑1 (MQ=255)                                                                                                                                                                                                   gAGTATCCGTACTTCGTACCCGCATGTTTAGAGCTAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACTCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGtat  <  2:22080/250‑1 (MQ=255)                                                                                                                                                                                                   gAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGtat  >  2:213217/1‑250 (MQ=255)                                                                                                                                                                                                                             |                                                                                                                                                                                                                                CCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTAT  >  NC_000913/4541041‑4541484

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