Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	563,150	A→G	N272S (AAC→AGC)	sfmH →	FimA homolog, function unknown

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	563,150	0	A	G	100.0%	179.2 / NA	49	N272S (AAC→AGC)	sfmH	FimA homolog, function unknown
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/27);  total (22/27)

TCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAACGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGT  >  NC_000913/563003‑563296                                                                                                                                                    |                                                                                                                                                   tCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCg                                                                                                                                                  >  1:600627/1‑150 (MQ=255)   aTGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGAt                                                                                                                                                >  2:554449/1‑150 (MQ=255)   aTGCCGCAAACCAAGAGCATCGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATa                                                                                                                                               <  1:519889/151‑1 (MQ=255)        gCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCg                                                                                                                                          <  2:776777/151‑1 (MQ=255)        gCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCg                                                                                                                                          >  1:304093/1‑151 (MQ=255)          aaaCCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGAt                                                                                                                                        <  1:499236/151‑1 (MQ=255)               aaGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGc                                                                                                                                   <  2:418793/151‑1 (MQ=255)               aaGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGc                                                                                                                                   >  2:372876/1‑151 (MQ=255)                       tGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGAt                                                                                                                           >  2:681645/1‑151 (MQ=255)                         cGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGAtct                                                                                                                         <  1:408730/151‑1 (MQ=255)                             cAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAAt                                                                                                                     <  1:309963/151‑1 (MQ=255)                                 tGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCg                                                                                                                 >  2:724879/1‑151 (MQ=255)                                  gTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGt                                                                                                                >  2:702319/1‑151 (MQ=255)                                              ctgcGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTcc                                                                                                     >  2:615041/1‑150 (MQ=255)                                                      gCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATg                                                                                            >  1:619543/1‑151 (MQ=255)                                                         tATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGc                                                                                           >  2:398154/1‑149 (MQ=255)                                                         tATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGAGGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCa                                                                                          <  1:710457/150‑1 (MQ=255)                                                            ttAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCt                                                                                       <  2:411719/150‑1 (MQ=255)                                                                 aaTGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGc                                                                                  >  2:439706/1‑150 (MQ=255)                                                                                     ccGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCgcgc                                                             <  1:111356/151‑1 (MQ=255)                                                                                         ttCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTgg                                                         >  1:30763/1‑151 (MQ=255)                                                                                         ttCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTg                                                          >  2:684024/1‑150 (MQ=255)                                                                                         ttCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACACCTAGCGATCTCAATAGCGTTATGCCATTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGTGCCTg                                                          >  2:372656/1‑150 (MQ=255)                                                                                              gTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTAt                                                    <  1:172123/151‑1 (MQ=255)                                                                                              gTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTAt                                                    <  1:396823/151‑1 (MQ=255)                                                                                                cAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATc                                                   <  1:328917/150‑1 (MQ=255)                                                                                                 aGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATc                                                   >  2:653109/1‑149 (MQ=255)                                                                                                  ggCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGt                                                <  2:739339/151‑1 (MQ=255)                                                                                                   gCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTa                                               <  2:401281/151‑1 (MQ=255)                                                                                                             cGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCa                                     <  1:753253/151‑1 (MQ=255)                                                                                                              ggACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCa                                     >  2:163587/1‑150 (MQ=255)                                                                                                                  aaTCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAac                                <  2:736044/151‑1 (MQ=255)                                                                                                                  aaTCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAac                                <  1:406235/151‑1 (MQ=255)                                                                                                                           ttAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGaa                        <  1:334733/150‑1 (MQ=255)                                                                                                                            tAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCaaaa                                  >  2:406705/1‑139 (MQ=255)                                                                                                                            tAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCaaaa                                  <  1:406656/139‑1 (MQ=255)                                                                                                                              ggTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGc                    >  2:241313/1‑151 (MQ=255)                                                                                                                                    aTTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGttt               >  2:530210/1‑150 (MQ=255)                                                                                                                                     ttGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGttt               >  2:346737/1‑149 (MQ=255)                                                                                                                                     ttGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGttt               >  1:642520/1‑149 (MQ=255)                                                                                                                                      tGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAATCGACCGAAGGGCCGTTTAg             <  1:577480/150‑1 (MQ=255)                                                                                                                                         cGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAgcgc          <  2:546410/150‑1 (MQ=255)                                                                                                                                          gCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAgcgcgc        >  1:66256/1‑151 (MQ=255)                                                                                                                                            cGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTg      <  1:278000/151‑1 (MQ=255)                                                                                                                                            cGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTg      <  1:96939/151‑1 (MQ=255)                                                                                                                                            cGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTg      <  1:570554/151‑1 (MQ=255)                                                                                                                                            cGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTg      <  2:144748/151‑1 (MQ=255)                                                                                                                                                cAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTgggg   <  1:700495/150‑1 (MQ=255)                                                                                                                                                cAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTgggg   <  2:275180/150‑1 (MQ=255)                                                                                                                                                 agagCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGt  >  1:632633/1‑150 (MQ=255)                                                                                                                                                 agagCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGt  >  1:690503/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAACGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGT  >  NC_000913/563003‑563296

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