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

GTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAACGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCT  >  NC_000913/563002‑563291                                                                                                                                                     |                                                                                                                                              gtCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCg                                                                                                                                             >  2:631259/1‑151 (MQ=255)      gCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACg                                                                                                                                        >  2:742055/1‑151 (MQ=255)        cGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGcc                                                                                                                                      >  2:355551/1‑151 (MQ=255)           aaaCCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGAt                                                                                                                                   >  2:792662/1‑151 (MQ=255)            aaCCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATc                                                                                                                                  <  1:156891/151‑1 (MQ=255)              ccAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCAc                                                                                                                                >  2:332469/1‑151 (MQ=255)              ccAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCAc                                                                                                                                <  2:486723/151‑1 (MQ=255)              ccAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCAc                                                                                                                                >  2:140917/1‑151 (MQ=255)                 agagCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGcc                                                                                                                             >  2:331241/1‑151 (MQ=255)                   agCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCt                                                                                                                            >  1:554312/1‑150 (MQ=255)                    gCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTaa                                                                                                                          <  1:108252/151‑1 (MQ=255)                     cATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAAc                                                                                                                         <  2:520438/151‑1 (MQ=255)                      aTTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACg                                                                                                                        <  1:148637/151‑1 (MQ=255)                         gCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGAtc                                                                                                                     >  1:147617/1‑151 (MQ=255)                                   gTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGt                                                                                                           >  2:469668/1‑151 (MQ=255)                                   gTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGGTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCg                                                                                                            >  2:34400/1‑150 (MQ=255)                                     aCGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGtt                                                                                                          >  1:838915/1‑150 (MQ=255)                                     aCGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTa                                                                                                         <  1:327054/151‑1 (MQ=255)                                     aCGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTa                                                                                                         <  2:319196/151‑1 (MQ=255)                                          tGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGGGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATTTCAATATCGTTATTCCt                                                                                                    >  1:213022/1‑151 (MQ=255)                                           gTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCtt                                                                                                   >  2:576261/1‑151 (MQ=255)                                           gTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCtt                                                                                                   <  1:674461/151‑1 (MQ=255)                                           gTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCt                                                                                                    <  1:31170/150‑1 (MQ=255)                                            ttGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCttt                                                                                                  >  2:21690/1‑151 (MQ=255)                                             tgctgcGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCttt                                                                                                  <  2:468145/150‑1 (MQ=255)                                              gctgcGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTcc                                                                                                <  1:752083/151‑1 (MQ=255)                                              gctgcGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTcc                                                                                                <  1:432538/151‑1 (MQ=255)                                              gctgcGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTc                                                                                                 <  2:291855/150‑1 (MQ=255)                                                      ggCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGAt                                                                                        <  1:649378/151‑1 (MQ=255)                                                           aTTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGc                                                                                   >  1:78913/1‑151 (MQ=255)                                                               aaCAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGTGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGTGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCg                                                                               <  2:525542/151‑1 (MQ=255)                                                               aaCAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGc                                                                                <  2:302790/150‑1 (MQ=255)                                                               aaCAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCg                                                                               <  1:513582/151‑1 (MQ=255)                                                                  aaTGCGTCTGGAAGCCAGTGCCGTTTCTGGCCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCATTTCCGTCTGGATGCAGCTGCGGc                                                                             >  1:98372/1‑150 (MQ=255)                                                                        tCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAAt                                                                      >  1:724285/1‑151 (MQ=255)                                                                          tGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGt                                                                    >  2:321319/1‑151 (MQ=255)                                                                                  aGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTc                                                            >  2:480682/1‑151 (MQ=255)                                                                                  aGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTc                                                            <  1:21668/151‑1 (MQ=255)                                                                                     gCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTcgc                                                          <  1:735999/150‑1 (MQ=255)                                                                                        gTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCt                                                      <  1:766170/151‑1 (MQ=255)                                                                                        gTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCt                                                      <  1:530438/151‑1 (MQ=255)                                                                                                   ggCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGTAGCCAATGTCACACTTCGCGCCTGGCCTATCAGt                                           >  2:397808/1‑151 (MQ=255)                                                                                                    gCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGt                                           >  1:792611/1‑150 (MQ=255)                                                                                                      gATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTAtt                                        <  1:754546/151‑1 (MQ=255)                                                                                                           tgtCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTAGCGCCTGGCCTATCAGTAttgcag                                    >  1:787068/1‑146 (MQ=255)                                                                                                               ggACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCa                                >  2:840813/1‑150 (MQ=255)                                                                                                                        ggATTTAGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAaccgacc                      >  2:84934/1‑151 (MQ=255)                                                                                                                              aGGTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAggg                >  1:684430/1‑151 (MQ=255)                                                                                                                               ggTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAaccgacc                      >  2:379947/1‑144 (MQ=255)                                                                                                                               ggTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAaccgacc                      <  1:379923/144‑1 (MQ=255)                                                                                                                               ggTTTTATTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGc               <  1:187346/151‑1 (MQ=255)                                                                                                                                 ttttATTGTCGCCGATCAGAGCGATACGCCGATCAAGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTAGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCg             >  1:681642/1‑151 (MQ=255)                                                                                                                                     aTTGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTa         <  2:396460/151‑1 (MQ=255)                                                                                                                                      ttGTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGttt          >  1:460219/1‑149 (MQ=255)                                                                                                                                        gTCGCCGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAgcg      >  2:87807/1‑151 (MQ=255)                                                                                                                                          cGCCGATCAGAGCGAGACGCCGATCACGCCTAACGATCGCAATAGCGTGAGTCCTTTCCGTCGGGATGCAGCTGCAGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAGGGGCCGTTTAgcgc     <  1:446717/150‑1 (MQ=255)                                                                                                                                            ccGATCAGAGCGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCt  <  1:589124/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              GTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAATCAGGATTTAGGTTTTATTGTCGCCGATCAGAACGATACGCCGATCACGCCTAACGATCTCAATAGCGTTATTCCTTTCCGTCTGGATGCAGCTGCGGCAGCCAATGTCACACTTCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCT  >  NC_000913/563002‑563291

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