Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,546,601	C→T	I502I (ATC→ATT)	fimD →	type I fimbriae usher protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,546,601	0	C	T	93.5%	88.1 / ‑1.3	31	I502I (ATC→ATT)	fimD	type I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (17/12);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

CGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAA  >  NC_000913/4546356‑4546844                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    cGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGaaa                                                                                                                                                                                                                                                 <  1:38992/250‑1 (MQ=255)       gttAATTTCGGTATAGGGAAAAACATGGGGGTAATGGGCGGTCTGTCTGTGGATAGGACGCAGGCTAATTCCACACTTCCCTATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAATCACAgg                                                                                                                                                                                                                                           <  2:80811/249‑1 (MQ=255)                   aTCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTc                                                                                                                                                                                                                               <  1:60324/250‑1 (MQ=255)                             aCATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGCTACAACATACAGTCGAATGCATGGCACCAACATTGAAACACAGTCAGGAGTTATTCAGGTTAAGCCg                                                                                                                                                                                                                    >  2:87625/1‑251 (MQ=255)                                     gCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATGTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAc                                                                                                                                                                                                            <  1:12477/251‑1 (MQ=255)                                              gcTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAtt                                                                                                                                                                                                    >  1:157030/1‑250 (MQ=255)                                               ctctGACTGAGGATAGGAAGCAGGATAATTACAAACTTCCCGCTGACAGTAAGCATGACGGACAAACGGAGCGTTTTCTATATAAAAAATCCCTCGATGAATCAGGCACGAATATTAAGTTAGTGGGATACAGTTATTCGACCAGCGGATATTTTAATTTCGCTGGTAAAACATACAGTGGAATGAAGGGCTACAAAATTGAAACACAGGACGGAGTTATTCAGTATAAGACGAAATTAAACGACAATTa                                                                                                                                                                                                   >  2:147552/1‑250 (MQ=255)                                                   gtctgtGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCTAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACaa                                                                                                                                                                                                <  2:185307/249‑1 (MQ=255)                                                       ggggATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAAGCGGTGCGTTTTCTCTATAACAAAGCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCg                                                                                                                                                                                           <  2:157030/248‑1 (MQ=255)                                                         ggATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAAGCGCACAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCtt                                                                                                                                                                                        <  2:156089/251‑1 (MQ=255)                                                                   cAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcg                                                                                                                                                                              >  1:157074/1‑251 (MQ=255)                                                                                           gACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTAc                                                                                                                                                      >  1:66204/1‑251 (MQ=255)                                                                                             cAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTc                                                                                                                                                    >  2:70442/1‑251 (MQ=255)                                                                                                         cGGACAATCGGTGCGTTTTCTCTATCACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGgc                                                                                                                                        >  1:207653/1‑251 (MQ=255)                                                                                                                    tGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTAGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGAAGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATAACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAac                                                                                                                              >  2:124214/1‑250 (MQ=255)                                                                                                                              ctATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTAtt                                                                                                                      <  1:234197/248‑1 (MQ=255)                                                                                                                                   aCAAATCGCTCAATGAATCATGCCAGATTATTCAGTTAGTGGGTTACCGTGATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATTCAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGTGCACATCAACACTGTATTTGAGTGGt                                                                                                              <  2:207653/251‑1 (MQ=255)                                                                                                                                    cAAATCGCTCAATGAATCAGGCACGCATATTCAGTTAGTGGGTAACAGTTATTCGACCAGCGGATATTTTAATTTCGCTGATCCAACATACAGTCGAATGAATGGCTACAACATTGAACCACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACACGGGAAATTACAACTCACCGTTAATCAGCAACTCGGGCGCCCATCAACACTGTATTTGAGTGGt                                                                                                              >  2:42054/1‑250 (MQ=255)                                                                                                                                    cAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGt                                                                                                              >  2:176388/1‑250 (MQ=255)                                                                                                                                           cTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGAACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATc                                                                                                       >  1:290275/1‑250 (MQ=255)                                                                                                                                                                  tCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATg                                                                               >  1:259163/1‑251 (MQ=255)                                                                                                                                                                                     tATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGc                                                            >  1:224458/1‑251 (MQ=255)                                                                                                                                                                                             cAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGTCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGAt                                                                          <  2:237227/229‑1 (MQ=255)                                                                                                                                                                                             cAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGAt                                                                          >  1:237227/1‑229 (MQ=255)                                                                                                                                                                                                              tCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAGTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGa                                    <  1:37455/250‑1 (MQ=255)                                                                                                                                                                                                                gCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGatat                                 >  1:208982/1‑251 (MQ=255)                                                                                                                                                                                                                            tACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCt                     >  1:67469/1‑251 (MQ=255)                                                                                                                                                                                                                                   gaatgaatGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCtat              >  2:100466/1‑251 (MQ=255)                                                                                                                                                                                                                                       gaatgGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGcc          <  1:124214/251‑1 (MQ=255)                                                                                                                                                                                                                                           ggCTACAACATTGAACCACAGAACGGAGTTATTCAGGTTAAGCCGAAATTCACCTCCTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGACATCAAACTTATTGGGGACCGAGTCATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTGCGAAGATATCAACTGGACGCTCAGCTATAGACTGAc      >  2:177707/1‑251 (MQ=255)                                                                                                                                                                                                                                               acaacaTTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAGTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGaaa  <  2:67469/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    CGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAA  >  NC_000913/4546356‑4546844

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: ‑