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%	182.1 / NA	50	N272S (AAC→AGC)	sfmH	FimA homolog, function unknown
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/28);  total (22/28)

GTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAACGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGG  >  NC_000913/563002‑563295                                                                                                                                                     |                                                                                                                                                  gtCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCg                                                                                                                                                 >  1:70653/1‑151 (MQ=255) gtCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCg                                                                                                                                                 <  1:902024/151‑1 (MQ=255) gtCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCg                                                                                                                                                 <  2:483299/151‑1 (MQ=255)   cATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGAt                                                                                                                                               <  2:753119/151‑1 (MQ=255)      gCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATAc                                                                                                                                             <  1:50525/150‑1 (MQ=255)           aaaCCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGAt                                                                                                                                       >  1:81069/1‑151 (MQ=255)            aaCCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATc                                                                                                                                      <  1:641132/151‑1 (MQ=255)               cAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACg                                                                                                                                   >  1:716689/1‑151 (MQ=255)                   agCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTa                                                                                                                               <  2:132634/151‑1 (MQ=255)                   agCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTa                                                                                                                               <  2:183934/151‑1 (MQ=255)                     cATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAAc                                                                                                                             <  2:878256/151‑1 (MQ=255)                          cGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGAtct                                                                                                                        >  1:520025/1‑151 (MQ=255)                              cAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAAt                                                                                                                    <  1:467330/151‑1 (MQ=255)                               aaaTGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATa                                                                                                                   <  2:460588/151‑1 (MQ=255)                                  tGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCg                                                                                                                >  1:247890/1‑151 (MQ=255)                                   gTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCg                                                                                                                <  1:392063/150‑1 (MQ=255)                                   gTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCg                                                                                                                <  1:22105/150‑1 (MQ=255)                                           gTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCtt                                                                                                       >  1:475756/1‑151 (MQ=255)                                           gTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCtt                                                                                                       <  1:243177/151‑1 (MQ=255)                                                 gcgcAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTc                                                                                                 >  2:621901/1‑151 (MQ=255)                                                          tATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGc                                                                                          >  1:3420/1‑149 (MQ=255)                                                             ttAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCt                                                                                      <  1:749969/150‑1 (MQ=255)                                                              tAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTg                                                                                     <  2:542239/150‑1 (MQ=255)                                                                 cAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGc                                                                                 >  1:146133/1‑151 (MQ=255)                                                                  aaTGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGc                                                                                 >  1:409867/1‑150 (MQ=255)                                                                       gTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCa                                                                            >  1:721250/1‑150 (MQ=255)                                                                               gCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCacac                                                                   <  2:361760/151‑1 (MQ=255)                                                                                  aGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTc                                                                <  1:750732/151‑1 (MQ=255)                                                                                   gTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTcg                                                               <  1:643901/151‑1 (MQ=255)                                                                                   gTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTcg                                                               <  1:441543/151‑1 (MQ=255)                                                                                   gTGCCGTTTCTGATCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTcg                                                               <  2:395025/151‑1 (MQ=255)                                                                                                 cAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATc                                                  <  1:621950/150‑1 (MQ=255)                                                                                                     cGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTAt                                             <  2:584059/151‑1 (MQ=255)                                                                                                            gtCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCgg                                       <  1:198997/150‑1 (MQ=255)                                                                                                              cGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCa                                    <  1:531988/151‑1 (MQ=255)                                                                                                               ggACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCa                                    >  2:629661/1‑150 (MQ=255)                                                                                                                   aaTCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCATCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAac                               >  1:142953/1‑151 (MQ=255)                                                                                                                   aaTCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAac                               >  2:239660/1‑151 (MQ=255)                                                                                                                   aaTCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAac                               <  2:317044/151‑1 (MQ=255)                                                                                                                     tCAGGATTTAGGTTTTATTGTCGCCGCTCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCGGGCCTATCAGTATTACCGGTCAACAaccg                             >  2:420158/1‑151 (MQ=255)                                                                                                                        ggATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAaccgacc                          <  2:98362/151‑1 (MQ=255)                                                                                                                         gATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAaccgaccg                         <  2:319044/151‑1 (MQ=255)                                                                                                                             tAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAgg                     >  1:472134/1‑151 (MQ=255)                                                                                                                              aGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAggg                    >  1:8608/1‑151 (MQ=255)                                                                                                                              aGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAggg                    >  1:253228/1‑151 (MQ=255)                                                                                                                                  tttCTTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGt                >  2:865588/1‑151 (MQ=255)                                                                                                                                  tttATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGt                >  1:244939/1‑151 (MQ=255)                                                                                                                                       tGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAgc           >  2:114328/1‑151 (MQ=255)                                                                                                                                             cGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTg     >  1:603134/1‑151 (MQ=255)                                                                                                                                                 cAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTgggg  <  2:253202/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAACGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGG  >  NC_000913/563002‑563295

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