Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	558,618	A→G	N136S (AAT→AGT)	sfmA →	putative fimbrial protein SfmA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	558,618	0	A	G	100.0%	202.2 / NA	58	N136S (AAT→AGT)	sfmA	putative fimbrial protein SfmA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/25);  total (33/25)

TGGTCAGGCTGATGCAATCAATGATAATTTATTGGCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGA  >  NC_000913/558493‑558750                                                                                                                              |                                                                                                                                     tGGTCAGGCTGATGCAATCAATGATAATTTATTGGCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTATCAAACCTGATGGGAGTAGCTTCTcaacca                                                                                                                        >  2:292377/1‑140 (MQ=255)   gTCAGGCTGATGCAATCAATGATAATTTATTGGCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTcaaccaac                                                                                                                      <  2:390164/140‑1 (MQ=255)   gTCAGGCTGATGCAATCAATGATAATTTATTGGCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTcaaccaac                                                                                                                      >  2:678433/1‑140 (MQ=255)   gTCAGGCTGATGCAATCAATGATAATTTATTGGCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTcaaccaac                                                                                                                      <  2:778389/140‑1 (MQ=255)    tCAGGCTGATGCAATCAATGATAATTTATTGGCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCaaccaacc                                                                                                                     <  1:436298/140‑1 (MQ=255)        gCTGATGCAATCAATGATAATTTATTGGCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGaa                                                                                                                 >  1:140233/1‑140 (MQ=255)          tGATGCAATCAATGATAATTTATTGGCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACt                                                                                                               <  1:299405/140‑1 (MQ=255)             tGCAATCAATGATAATTTATTGGCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGa                                                                                                            >  1:33112/1‑140 (MQ=255)             tGCAATCAATGATAATTTATTGGCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGa                                                                                                            <  1:344586/140‑1 (MQ=255)             tGCAATCAATGATAATTTATTGGCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGa                                                                                                            >  1:598337/1‑140 (MQ=255)                       gATAATTTATTGGCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGAc                                                                                                  <  2:140233/140‑1 (MQ=255)                                   gCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCa                                                                                      <  2:536226/140‑1 (MQ=255)                                       ttGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGtt                                                                                           >  1:96530/1‑131 (MQ=255)                                       ttGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGtt                                                                                           <  2:96530/131‑1 (MQ=255)                                       ttGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCAtttt                                                                                  >  2:111538/1‑140 (MQ=255)                                        tGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCAttttt                                                                                 >  1:269231/1‑140 (MQ=255)                                        tGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCAttttt                                                                                 >  1:813493/1‑140 (MQ=255)                                            tcgaGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGACGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGc                                                                             <  2:661781/137‑1 (MQ=255)                                            tCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTATCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGc                                                                             <  1:292377/140‑1 (MQ=255)                                                gTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGt                                                                         <  2:610892/140‑1 (MQ=255)                                                 tACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGtt                                                                        <  2:557488/140‑1 (MQ=255)                                                 tACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGtt                                                                        >  2:625525/1‑140 (MQ=255)                                                 tACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGtt                                                                        >  1:627412/1‑140 (MQ=255)                                                      aTACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATaaa                                                                   <  2:269231/140‑1 (MQ=255)                                                       tACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAg                                                                  >  2:150434/1‑140 (MQ=255)                                                           acaacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCa                                                                            <  1:479321/126‑1 (MQ=255)                                                           acaacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCa                                                                            >  2:479321/1‑126 (MQ=255)                                                           acaacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCAc                                                              <  2:813493/140‑1 (MQ=255)                                                           acaacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCAc                                                              >  1:197987/1‑140 (MQ=255)                                                               caGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGt                                                          >  2:601723/1‑140 (MQ=255)                                                               caGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGt                                                          >  2:705067/1‑140 (MQ=255)                                                                 gCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTAc                                                        >  1:141617/1‑140 (MQ=255)                                                                  cAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACa                                                       >  2:387190/1‑140 (MQ=255)                                                                  cAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACa                                                       >  2:584852/1‑140 (MQ=255)                                                                    aCGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATTGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                     >  1:815670/1‑140 (MQ=255)                                                                    aCGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                     <  1:387190/140‑1 (MQ=255)                                                                         tgtCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCAt                                                >  2:182477/1‑140 (MQ=255)                                                                         tgtCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCAt                                                >  1:659013/1‑140 (MQ=255)                                                                          gtCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACGTGATGGGAGTAGCTTCTCGACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATc                                               >  2:739553/1‑140 (MQ=255)                                                                             ggTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATcagc                                            >  1:670156/1‑140 (MQ=255)                                                                                 ttGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAggg                                        <  2:670156/140‑1 (MQ=255)                                                                                 ttGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAggg                                        <  2:61059/140‑1 (MQ=255)                                                                                         cTTGATAATACATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCaaa                                >  1:794229/1‑140 (MQ=255)                                                                                                  aCATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGAcgc                       <  2:92258/140‑1 (MQ=255)                                                                                                   cATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGAcgcg                      >  2:218272/1‑140 (MQ=255)                                                                                                   cATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGAcgcg                      <  2:659013/140‑1 (MQ=255)                                                                                                   cATCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGAcgcg                      >  1:788480/1‑140 (MQ=255)                                                                                                    aTCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGa                     <  1:652754/140‑1 (MQ=255)                                                                                                    aTCCGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGa                     >  2:473997/1‑140 (MQ=255)                                                                                                      ccGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACt                   >  1:292507/1‑140 (MQ=255)                                                                                                      ccGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACt                   >  1:522779/1‑140 (MQ=255)                                                                                                       cGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAATTGCATCAGCAGGGCAAGCAAATGCTGACGCGACtt                  >  1:203926/1‑140 (MQ=255)                                                                                                       cGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACtt                  >  2:444723/1‑140 (MQ=255)                                                                                                       cGCAATTCTCAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACtt                  <  2:197987/140‑1 (MQ=255)                                                                                                              ctcAAACCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTtattat           <  2:141617/140‑1 (MQ=255)                                                                                                                   aCCTGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGat      <  2:636304/140‑1 (MQ=255)                                                                                                                      tGATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATg   >  1:648597/1‑140 (MQ=255)                                                                                                                       gATGGGAGTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGa  <  2:794229/140‑1 (MQ=255)                                                                                                                              |                                                                                                                                     TGGTCAGGCTGATGCAATCAATGATAATTTATTGGCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGA  >  NC_000913/558493‑558750

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