Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,003,243	G→T	G119W (GGG→TGG)	ycbU →	putative fimbrial protein YcbU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,003,243	0	G	T	92.0%	69.4 / NA	25	G119W (GGG→TGG)	ycbU	putative fimbrial protein YcbU
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (8/15);  total (10/15)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGATTAACGTG  >  NC_000913/1002999‑1003485                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   ggTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTg                                                                                                                                                                                                                                              <  2:192033/251‑1 (MQ=255)  gTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGc                                                                                                                                                                                                                                             <  2:98309/251‑1 (MQ=255)              ggtgATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAAGACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATTGCACCGTGGATTGCCGTGCAAATTCt                                                                                                                                                                                                                                 <  2:292306/250‑1 (MQ=255)                             aaaaGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTc                                                                                                                                                                                                                   >  1:185313/1‑250 (MQ=255)                               aaGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGCTATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGCCACTGCTCAAACTAACTGACACAGGCGGAAGCGGCGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGc                                                                                                                                                                                                                >  1:169995/1‑251 (MQ=255)                                aGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGc                                                                                                                                                                                                                >  1:343771/1‑250 (MQ=255)                                          aaaaTATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCCACGGGGATTGCCGTGCAAATTCTTTATGCGGAAAATCAGCAAAAAATcccc                                                                                                                                                                                                     >  2:291230/1‑250 (MQ=255)                                                   atatTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATc                                                                                                                                                                                              >  1:108897/1‑249 (MQ=255)                                                       tGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTcc                                                                                                                                                                                        >  1:282122/1‑251 (MQ=255)                                                                  atatCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACg                                                                                                                                                                             <  1:281235/251‑1 (MQ=255)                                                                                             ccGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTc                                                                                                                                                  <  2:24842/251‑1 (MQ=255)                                                                                                 aaGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCGTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAAGCAGGTCCAGCATCTTACGCCCGTAAATGCCGGGGATAACACACTCaaa                                                                                                                                               <  2:335546/250‑1 (MQ=255)                                                                                                     aCTGCTGCGGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCa                                                                                                                                          <  2:164172/251‑1 (MQ=255)                                                                                                                       atCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTTATAAGTCcac                                                                                                                        >  1:64727/1‑251 (MQ=255)                                                                                                                               cACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCGGGGGATAACACACTCAAAGATCAACTTCGTTATAAGTCCACAAAGGCg                                                                                                                 <  2:160554/250‑1 (MQ=255)                                                                                                                                   ggTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACCGGCGGAAGCGGGGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTTCTAAGTCCACAAAGGCGGGAGc                                                                                                            >  2:161284/1‑251 (MQ=255)                                                                                                                                         gCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTTATAAGTCCACAAAGGCGGGAGCAACggg                                                                                                      >  2:48555/1‑251 (MQ=255)                                                                                                                                                        ggCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACGAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGAcggc                                                                                       <  2:288750/251‑1 (MQ=255)                                                                                                                                                             ggaagtgcTTTTAAGCGGAGAGGCTGACACCCTTGCGCAGGAACGGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTGAAAAGCCGAGGATAACAGACTCAAATATCAACTGCGTTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTtc                                                                                  <  2:352443/244‑1 (MQ=255)                                                                                                                                                             ggACGTCCTTGTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTtc                                                                                  <  2:352466/251‑1 (MQ=255)                                                                                                                                                                                       acacCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGTTTTACCAGTGa                                                        <  2:147851/251‑1 (MQ=255)                                                                                                                                                                                          cccTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAgg                                                     <  2:344429/251‑1 (MQ=255)                                                                                                                                                                                            cgtGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTGATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAg                                                   <  1:319336/249‑1 (MQ=255)                                                                                                                                                                                                    gACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCGACTTCGTTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTg                                           <  2:64727/251‑1 (MQ=255)                                                                                                                                                                                                                                             tGGCAACGTGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGATTAACGTg  >  1:82654/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   GGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGATTAACGTG  >  NC_000913/1002999‑1003485

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