Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,364,664	G→A	G620R (GGA→AGA)	yhcD →	putative fimbrial usher protein YhcD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,364,664	0	G	A	100.0%	103.3 / NA	35	G620R (GGA→AGA)	yhcD	putative fimbrial usher protein YhcD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/22);  total (13/22)

TGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGGGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGT  >  NC_000913/3364422‑3364903                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                tGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGAcc                                                                                                                                                                                                                                         >  1:28625/1‑251 (MQ=255)                     cATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCaa                                                                                                                                                                                                                     <  2:24705/251‑1 (MQ=255)                               aCTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTc                                                                                                                                                                                                           <  2:213390/251‑1 (MQ=255)                                  aTGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGc                                                                                                                                                                                                         <  2:45415/250‑1 (MQ=255)                                        cTGGTATCATGGGAGTCGTCTCAGCTGAGCTGTCCTACTCCGTGACTGCTGCACGAGATCACGCCGAACGTGAAACGAGCTTGAAAGACAGTATCAGTTCCAATCCCCCTTACACCCAAGTCAGCCTTGCCGCAGGAGCACGCGGCAGCGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGAACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGc                                                                                                                                                                                                   >  2:261582/1‑250 (MQ=255)                                           gtgtAATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGAAGGAGCAAGCGGCAGTGATAGCAGTACTTACAACGGTACGATGTCAGGTGGGATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGACATCGCaaa                                                                                                                                                                                                >  1:142624/1‑250 (MQ=255)                                             gtgaTGGGAGTCGTCTCAGATGAGATGTACTACTATGTGAATGCTTAACTAGATCACGACGAACGTGAAACGAGCTTGAACGTCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGGACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCAACGTGGACTATCAATGACACTTTCGCCATCGCAAAAAt                                                                                                                                                                                             <  1:153145/248‑1 (MQ=255)                                                aTGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATgg                                                                                                                                                                                           >  2:89412/1‑250 (MQ=255)                                                aTGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATgg                                                                                                                                                                                           >  2:478/1‑250 (MQ=255)                                                                            actaTGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTTATAGCCGTACTTACAACTGTACGATGTCATGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCGCCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTa                                                                                                                                                                <  1:20373/249‑1 (MQ=255)                                                                            actaTGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGTTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACTCTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACa                                                                                                                                                              <  1:108007/251‑1 (MQ=255)                                                                               aTGTGAATGCTGAACGAGATCACGACGAACGTGAAACGATCTTGAACGGCAGTATCAGTTCCAATCACCATTACTCCCAAGTCAGCCTTGCCTCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCTTACATGATCAGAGAGTGACCATTTCACCGTGGACTATCAATGTCACTGTCGCCATCGCAAAAATGGATAACAAGATTGCAGGTGTCATAATTACAtct                                                                                                                                                           <  2:100785/251‑1 (MQ=255)                                                                                   gAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTcagg                                                                                                                                                       <  1:89412/251‑1 (MQ=255)                                                                                    aaTGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTcaggc                                                                                                                                                      >  1:20920/1‑251 (MQ=255)                                                                                                   cACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCAGTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGAc                                                                                                                                       <  2:22628/251‑1 (MQ=255)                                                                                                      gacgaACGTGAAACGAGCTTGAACGGCAGTAACAGTTCCAATCGCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAAGGACACTTTCGCCATCGCAAAAATGGATAACAATAGTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAgc                                                                                                                                    <  2:17488/251‑2 (MQ=255)                                                                                                            cGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCg                                                                                                                              >  2:71064/1‑251 (MQ=255)                                                                                                                        ttGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGt                                                                                                                  >  1:22390/1‑251 (MQ=255)                                                                                                                        ttGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCACAGTATGGACAGATTTTCGGGGTAATGCCGt                                                                                                                  >  1:151177/1‑251 (MQ=255)                                                                                                                           aaCGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCAt                                                                                                               <  2:62292/251‑1 (MQ=255)                                                                                                                            aCGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTGCGGGGTAATGCCGTCAtt                                                                                                              <  2:115275/251‑1 (MQ=255)                                                                                                                                            cAATCTCCATGACACCCAAGTCAGCCTGGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATGTTCGGGGTAATTCCGTCATTCCATCAATCCAGCCGt                                                                                              <  2:29378/251‑1 (MQ=255)                                                                                                                                                       acacCCAAGTCAGCCATGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGGCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTGCACCGTGGACTATCAATGACACTGTCGCCATCGCAAAAAGTGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATGTGCGGGGTAATGCCGAAATTCCATCAATCCAGCCGTGGCGAACATCa                                                                                   <  1:247626/251‑1 (MQ=255)                                                                                                                                                          cccAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAgg                                                                                 >  1:52349/1‑250 (MQ=255)                                                                                                                                                              aGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGAGCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTTCATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAgatg                                                                            <  2:259709/251‑4 (MQ=255)                                                                                                                                                                tCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGGCCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAAGTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGcc                                                                                                                    <  1:98817/209‑1 (MQ=255)                                                                                                                                                                tCAGCCTCACCGCAGGAGCACGCGGCAGTGATAGCCGTACTTACAACGGTACGATCTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGTCTATCAGTGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCCGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATgac                                                                                                                    >  2:98817/1‑207 (MQ=255)                                                                                                                                                                                         cAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGAGCAGAGAGGGACCTTTTCACCGTGGACTATCAATGACACGTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGGTTTTCGGGGTAATGCCGGCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCaaaaa                                                 <  1:261582/251‑1 (MQ=255)                                                                                                                                                                                          aGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGCCAGAATTACATCTCAGGCAGGCCCAGTAGGGACAGAGTTTCGGGGTAAGGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATCCCGCCAGCTTGCCaaaaaa                                                <  2:28625/251‑1 (MQ=255)                                                                                                                                                                                             gATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTGTAGATCGATACCGCCAGCTTGCCAAAAAATGt                                             <  1:71064/251‑1 (MQ=255)                                                                                                                                                                                                         tACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTa                                  <  1:54072/250‑1 (MQ=255)                                                                                                                                                                                                                  acgctgtaaGGTGGTATTGCCGTACATGCTCAGAGAGGGACCTTTTCCCCGTGGACTATCAATGACCCTTTCGCCCTCGCACAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGCCAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACaa                         >  2:237993/9‑250 (MQ=255)                                                                                                                                                                                                                          aGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATgttc                 >  1:9748/1‑247 (MQ=255)                                                                                                                                                                                                                              gggatggCCGTACATGATCAGAGAGTTACCTTTTCACCGTGGACTATCAATGACACTTTCTCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCTATACCGCCAGCTTGCCAAAAAATGGCGATATCGGTAACGGCACAAAAATGATCAAAcc            <  2:193087/245‑2 (MQ=255)                                                                                                                                                                                                                                         aCATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCgtggt  <  1:223566/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                TGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGGGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGT  >  NC_000913/3364422‑3364903

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