Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,541,175	(A)7→6	coding (219/603 nt)	fimB →	regulator for fimA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,541,169	0	A	.	97.6%	169.2 / ‑3.0	42	coding (213/603 nt)	fimB	regulator for fimA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (30/11);  minor base G (1/0);  total (31/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 = 1.00e+00

GCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAA  >  NC_000913/4541040‑4541297                                                                                                                                  |                                                                                                                                 gCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTaaaa                                                                                                                               >  2:102326/1‑133 (MQ=255)   cTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTaaaaaa                                                                                                                             <  1:102326/133‑1 (MQ=255)                       tGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCAccc                                                                                                        >  1:318604/1‑133 (MQ=255)                         tCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGc                                                                                                      >  2:386160/1‑133 (MQ=255)                         tCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATAGATATCCATAGA‑TGAAATAAGGCTTTTCAACAACGCACCCGc                                                                                                      >  2:196615/1‑133 (MQ=255)                           tGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTa                                                                                                    >  2:44615/1‑133 (MQ=255)                             aCTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTAtt                                                                                                  >  1:122455/1‑133 (MQ=255)                                     ttgTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAaag                                                                                          <  2:88933/133‑1 (MQ=255)                                                  ggTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCGCCCGCTAGTGAGTATAGAAGTACAGGCttt                                                                             >  1:141837/1‑133 (MQ=255)                                                          ggCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACt                                                                     <  2:345928/133‑1 (MQ=255)                                                               gTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTg                                                                <  1:386160/133‑1 (MQ=255)                                                                 gAAATTTGCCGATTGTGGCTTTCGGATATTCATCTTAAGGCAAAGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAg                                                              <  1:196615/133‑1 (MQ=255)                                                                       tGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCg                                                        <  2:122455/133‑1 (MQ=255)                                                                       tGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCg                                                        <  2:62420/133‑1 (MQ=255)                                                                       tGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCg                                                        <  2:407378/133‑1 (MQ=255)                                                                         tCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAGTAAAGAAGTTCACGCTTTAAAAAACAGGTTAAGTATCCGTa                                                      >  2:173358/1‑133 (MQ=255)                                                                                aGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTa                                               >  2:48086/1‑133 (MQ=255)                                                                                        ggATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATg                                       <  1:48086/133‑1 (MQ=255)                                                                                             ttGATCTTAAGGCAAAGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTgag                                  >  1:475344/1‑133 (MQ=255)                                                                                               gATCTTAAGGCAAAGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGagag                                <  2:475344/133‑1 (MQ=255)                                                                                                            aGTGTATATATATCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTAttt                   >  1:61205/1‑133 (MQ=255)                                                                                                                tatatataTCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTAt               <  2:61205/133‑1 (MQ=255)                                                                                                                     atatCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGt          >  1:435906/1‑133 (MQ=255)                                                                                                                       atCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTaa        >  1:299978/1‑133 (MQ=255)                                                                                                                        tCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACg           <  1:105623/129‑1 (MQ=255)                                                                                                                        tCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACg           >  2:105623/1‑129 (MQ=255)                                                                                                                        tCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAg       >  2:357562/1‑133 (MQ=255)                                                                                                                        tCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAg       >  1:503561/1‑133 (MQ=255)                                                                                                                        tCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAg       >  1:79276/1‑133 (MQ=255)                                                                                                                        tCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAg       >  1:466051/1‑133 (MQ=255)                                                                                                                        tCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAg       >  2:47145/1‑133 (MQ=255)                                                                                                                        tCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAg       >  1:401708/1‑133 (MQ=255)                                                                                                                        tCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGCACCCGCATCCTGAGTGCGAGTGGGTGTTTTTCTCACGTAAg       >  2:240690/1‑133 (MQ=255)                                                                                                                        tCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGCACCCGCATACTGCGAGTGAGTGGGTGTTTTAATCACGTAAg       >  2:137189/1‑133 (MQ=255)                                                                                                                        tCCATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGATGAGAATCCGTACTTCGGACCCGCATACTGTTAGTGAGTGGGTATTTTAATCACGTAAg       >  2:196206/1‑133 (MQ=255)                                                                                                                         ccATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTGCCCGCATGGTGAGAGCGAGTGGGTAGTTTTGTCACGTAAgg      >  2:130520/1‑133 (MQ=255)                                                                                                                         ccATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAgg      >  1:312646/1‑133 (MQ=255)                                                                                                                          cATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTAGCCGCATGCGGAGAGCGAGTGGGTATGTTTACCACGTAAggg     >  2:186194/1‑133 (MQ=255)                                                                                                                          cATCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAggg     >  2:319385/1‑133 (MQ=255)                                                                                                                            tCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGGGAGAGAGTGGGTATTTGTATCTCGTAAGGGGa   >  2:227632/1‑133 (MQ=255)                                                                                                                            tCGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGa   >  1:286519/1‑133 (MQ=255)                                                                                                                             cGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGaa  >  1:112301/1‑133 (MQ=255)                                                                                                                             cGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGATATGTTAATCACGGAAGGGGaa  >  1:160507/1‑133 (MQ=255)                                                                                                                             cGATT‑AAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGACAGCGAGTAGGTATTTTGATCAGGTAAggggca  >  2:167371/1‑131 (MQ=255)                                                                                                                                  |                                                                                                                                 GCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAA  >  NC_000913/4541040‑4541297

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