Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,544,301	A→G	intergenic (+34/‑3)	fimI → / → fimC	putative fimbrial protein FimI/type 1 fimbriae periplasmic chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,544,301	0	A	G	97.7%	149.1 / ‑4.5	44	intergenic (+34/‑3)	fimI/fimC	putative fimbrial protein FimI/type 1 fimbriae periplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (26/17);  total (27/17)
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.91e-01

TCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACT  >  NC_000913/4544166‑4544432                                                                                                                                        |                                                                                                                                    tctACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGt                                                                                                                                 <  2:724042/140‑1 (MQ=255)   tACTTCGATACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTTTTCAGCAGATACTGCGATAACTGGAACAGGGCAGTGa                                                                                                                               >  1:338011/1‑140 (MQ=255)    aCTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGt                                                                                                                                 >  2:575044/1‑137 (MQ=255)    aCTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGt                                                                                                                                 <  1:575044/137‑1 (MQ=255)           tACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGGAGACAATGTGATAACACGAACAGGGCAGTGAGTAATaaa                                                                                                                       >  2:344617/1‑140 (MQ=255)           tACATTTCATCACCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATaaa                                                                                                                       >  1:713531/1‑140 (MQ=255)               tttCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACg                                                                                                                   >  1:689310/1‑140 (MQ=255)                ttCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGTCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGt                                                                                                                  <  2:36206/140‑1 (MQ=255)                ttCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGt                                                                                                                  >  1:704938/1‑140 (MQ=255)                ttCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGt                                                                                                                  >  2:508534/1‑140 (MQ=255)                ttCATCGCCAAATATAGTGCTATCGAGCGGCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGt                                                                                                                  <  1:344617/140‑1 (MQ=255)                  cATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCa                                                                                                                >  1:543982/1‑140 (MQ=255)                     cGCCAAATATCGTGCTACCGGGTGTTGGGATACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATg                                                                                                             <  1:291913/140‑1 (MQ=255)                       ccAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACGGGGAAGTGAGTAACAAAAACGTCAATGTa                                                                                                           >  1:330800/1‑140 (MQ=255)                       ccAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTa                                                                                                           >  1:504141/1‑140 (MQ=255)                       ccAAATATAGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGTAGTGAGTAAGAAAAACGTCAATGTa                                                                                                           >  1:306173/1‑140 (MQ=255)                                gTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCg                                                                                                  <  1:677355/140‑1 (MQ=255)                                  gCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCa                                                                                                <  1:381461/140‑1 (MQ=255)                                    tACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAgg                                                                                              <  2:713531/140‑1 (MQ=255)                                      ccGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGGCAATGTGAGGAAATCGCAGGaa                                                                                            >  1:331907/1‑140 (MQ=255)                                         ggCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATa                                                                                         >  2:22355/1‑140 (MQ=255)                                                    aCTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCtgctt                                                                              >  1:35525/1‑140 (MQ=255)                                                     cTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCtgcttg                                                                             <  2:504141/140‑1 (MQ=255)                                                            cATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAg                                                                      <  2:689310/140‑1 (MQ=255)                                                            cATCGCCAATGCCCAGGCCTGGGTCTCTTTAACTTGTCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAg                                                                      <  2:330800/140‑1 (MQ=255)                                                            cATCGCCAATGCCCAGGCCTGGCTCTCTTTAACCTCTCAGTAATTGTTTAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAg                                                                      <  2:306173/140‑1 (MQ=255)                                                               cGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTa                                                                   <  1:22355/140‑1 (MQ=255)                                                                                  ttCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAat                                                <  1:508534/140‑1 (MQ=255)                                                                                               aTCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAAGTAACCTGATGTTCATGTTAATGATGGTTGCCGGa                                   >  2:341780/1‑140 (MQ=255)                                                                                                    tAATTGTTCAGCAGACAATGTGATAGCAGAAACAGGGGAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACgcgc                              <  1:303347/140‑1 (MQ=255)                                                                                                      aTTGTTCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTg                            >  1:654601/1‑140 (MQ=255)                                                                                                          ttCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGc                        >  2:522678/1‑140 (MQ=255)                                                                                                          ttCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCAGATATTCGTGGCAATGGTGGTGGCCGGACGCGCTGAAGc                        >  2:190142/1‑140 (MQ=255)                                                                                                          ttCAGCAGATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGAAGGTATCCCGATCTTCTTGGCAATGTTGGTTGCCGGACGCGCTGAAGc                        >  2:271106/1‑140 (MQ=255)                                                                                                            cagcagATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCACGTATCCTGATGTCCATCGCAATGATGGTTGACGGACGCGCTGAAGCgg                      >  2:303347/1‑140 (MQ=255)                                                                                                              gcagATAATGTGATATCAGGAATAGAGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGa                    <  1:341780/140‑1 (MQ=255)                                                                                                              gcagATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGGATCCTGATGTTCATAGCGATGATGGTTGCCGTACGCGCTAAAGCGGGa                    >  2:348984/1‑140 (MQ=255)                                                                                                               cagATAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAg                   <  2:704938/140‑1 (MQ=255)                                                                                                                  aTAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTgg                >  1:9790/1‑140 (MQ=255)                                                                                                                  aTAATGTGATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTgg                >  1:448000/1‑140 (MQ=255)                                                                                                                         gATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAgg         >  2:168934/1‑140 (MQ=255)                                                                                                                         gATAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAgg         >  1:22748/1‑140 (MQ=255)                                                                                                                           tAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCgg                      <  2:676334/125‑1 (MQ=255)                                                                                                                           tAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCgg                      >  1:676334/1‑125 (MQ=255)                                                                                                                           tAACAGGAACAGGGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTg       >  1:70916/1‑140 (MQ=255)                                                                                                                            aacaggaacaggGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGGGGCCTTAGGTGc      >  2:723128/1‑140 (MQ=255)                                                                                                                            aacaggaacaggGCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGCTTGCCGGACGCCCTGAAGCGGAAGTGGCCTTAGGTGc      >  1:360081/1‑140 (MQ=255)                                                                                                                                 gaacagggCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACt  <  1:18215/139‑1 (MQ=255)                                                                                                                                 gaacagggCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACt  >  2:18215/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                    TCTACTTCGCTACATTTCATCGCCAAATATCGTGCTACCGGGCGTCGGGTTACTGGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACT  >  NC_000913/4544166‑4544432

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