Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	562,201	C→T	G835G (GGC→GGT)	sfmD →	putative fimbrial usher protein SfmD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	562,201	0	C	T	100.0%	66.9 / NA	24	G835G (GGC→GGT)	sfmD	putative fimbrial usher protein SfmD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/12);  total (12/12)

ATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCT  >  NC_000913/561960‑562427                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                   aTGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCAATGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTc                                                                                                                                                                                                                            >  1:214142/1‑250 (MQ=255)     tGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGTATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTTGGTGATGACGGGCAAGTTTATTTAAGTTGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGc                                                                                                                                                                                                                        <  1:257283/250‑1 (MQ=255)            aaTATTCGCAGCTTTGCGCCGGCGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGTCGTGCGGGCGCTCATTACCTGTAACCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATGTAAGTGTTGCGCCATTGTCTGGTGAAGTACTGGTTCAGTGGGGAGACGTTTTGAACTCACGCTGCATTgc                                                                                                                                                                                                                <  2:409880/251‑1 (MQ=255)            aaTATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGAGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGAGAACTCAAGCTGCATTgc                                                                                                                                                                                                                >  2:108109/1‑251 (MQ=255)                   gcagcGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTAt                                                                                                                                                                                                         <  1:168716/251‑1 (MQ=255)                                                    gTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTc                                                                                                                                                                        >  1:44202/1‑251 (MQ=255)                                                          aTTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACGGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGtt                                                                                                                                                                  <  2:387076/251‑1 (MQ=255)                                                                                 cGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATcc                                                                                                                                           <  1:99384/251‑1 (MQ=255)                                                                                 cGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATcc                                                                                                                                           >  1:107696/1‑251 (MQ=255)                                                                                  gggCGCTCAGTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGAGCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAAGTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCt                                                                                                                                          <  1:51792/251‑1 (MQ=255)                                                                                         cATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCa                                                                                                                                    <  1:93950/250‑1 (MQ=255)                                                                                                   aCCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGGATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTAGGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAAt                                                                                                                          <  2:257687/250‑1 (MQ=255)                                                                                                           cGGAAAAACGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATa                                                                                                                  >  1:108105/1‑250 (MQ=255)                                                                                                                aaCCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATa                                                                                                             >  1:411170/1‑250 (MQ=255)                                                                                                                  ccGGTGACGATTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGAGGGTGATGACGGGCAAGTTTATTTCAGAGGAACGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCCGCATTGCTCACTATGTATTGCCGAAGCAAAGCTGACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGAAAACATCCTAGCTCATAAAGGAAATTATCACTAAGATAATc                                                                                                          >  2:17156/1‑251 (MQ=255)                                                                                                                                      ggTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGAt                                                                                      >  1:364723/1‑251 (MQ=255)                                                                                                                                         aCGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGc                                                                                   <  1:229908/251‑1 (MQ=255)                                                                                                                                                 aCAGTACCGGAATAACCAGTATGGTGGGCGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTtgt                                                                           >  2:121459/1‑251 (MQ=255)                                                                                                                                                  cAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGGTGTTCGCCATTGTCTTGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAATCCGTCACTTTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGTAAATTATCAAGAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTtgtg                                                                          <  1:309692/251‑1 (MQ=255)                                                                                                                                                                       ggtgggtgATGACGGGCAAGTTTATTTAAGTGGTGAGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGc                                                     >  1:221679/1‑251 (MQ=255)                                                                                                                                                                                        aaGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGCTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAAt                                     >  1:454973/1‑250 (MQ=255)                                                                                                                                                                                             tatttaAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGTCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTAc                                <  2:108105/250‑1 (MQ=255)                                                                                                                                                                                                                    tGGTGAATTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACt         >  1:364322/1‑250 (MQ=255)                                                                                                                                                                                                                          aTTACTGGTTCAGTGGGGAGACGGTGCGAACTCACGCTTCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCTGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGAGAATCTGCAGATTATTATTGGCGATGGCAGGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTTCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGAtct  <  1:269256/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                   ATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCT  >  NC_000913/561960‑562427

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