Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	150,947	T→C	T3A (ACT→GCT)	yadC ←	putative fimbrial‑like adhesin protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	150,947	0	T	C	100.0%	225.7 / NA	63	T3A (ACT→GCT)	yadC	putative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (32/31);  total (32/31)

AGGTCATCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGTCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCAGGATCGA  >  NC_000913/150812‑151073                                                                                                                                        |                                                                                                                               aGGTCATCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCtt                                                                                                                             >  1:858971/1‑139 (MQ=255)      aTCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCAttt                                                                                                                        <  2:154419/139‑1 (MQ=255)       tCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCAtttt                                                                                                                       <  2:400771/139‑1 (MQ=255)                 aGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATg                                                                                                             <  1:918709/139‑1 (MQ=255)                       ggggCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCtt                                                                                                       <  2:545446/139‑1 (MQ=255)                       ggggCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCtt                                                                                                       <  1:879992/139‑1 (MQ=255)                         ggCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                     <  2:70793/139‑1 (MQ=255)                         ggCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                     <  1:1037830/139‑1 (MQ=255)                          gCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTAt                                                                                                    >  2:312521/1‑139 (MQ=255)                             tCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGt                                                                                                 >  2:421257/1‑139 (MQ=255)                                    ttATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGtt                                                                                          >  1:427831/1‑139 (MQ=255)                                       tCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgc                                                                                       >  1:823404/1‑139 (MQ=255)                                        cAACTCGTGCATTATTTCCCCCTATAAAACCTGTATATGCAGTGGCTGTACTACAACAATAATACAGTAGTAAAAAAAGTATGTACCTGAAGATAGCCTGCATTTTTACTCGAATGTTTCTTTATTGTTATGGTTTgcg                                                                                      <  2:28015/139‑1 (MQ=255)                                         aaCTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcg                                                                                      >  2:586941/1‑138 (MQ=255)                                         aaCTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcg                                                                                      <  1:586941/138‑1 (MQ=255)                                               tGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCAt                                                                                                                          <  2:883873/96‑1 (MQ=255)                                               tGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCAt                                                                                                                          >  1:883873/1‑96 (MQ=255)                                                gCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGttttt                                                                              <  1:312521/139‑1 (MQ=255)                                                gCATTATTTCCCACAATAAAACTAGAATATCCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGACTGTACCTGAAGCTAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGttttt                                                                              >  2:91855/1‑139 (MQ=255)                                                                tAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCgtagt                                                              <  2:823404/139‑1 (MQ=255)                                                                        gTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTcc                                                                                                                 >  2:119946/1‑80 (MQ=255)                                                                        gTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTcc                                                                                                                 <  1:119946/80‑1 (MQ=255)                                                                        gTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATcc                                                      >  1:1038254/1‑139 (MQ=255)                                                                          atatGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCAc                                                    <  2:73730/139‑1 (MQ=255)                                                                                  tgactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATc                                                       <  2:46495/128‑1 (MQ=255)                                                                                  tgactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATc                                                       >  1:46495/1‑128 (MQ=255)                                                                                  tgactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACAt                                            <  2:427831/139‑1 (MQ=255)                                                                                   gactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATc                                                       >  1:333887/1‑127 (MQ=255)                                                                                   gactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATc                                                       <  2:333887/127‑1 (MQ=255)                                                                                    actgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGtt                                                                                                >  1:631981/1‑85 (MQ=255)                                                                                    actgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGtt                                                                                                <  2:631981/85‑1 (MQ=255)                                                                                    actgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcg                                                                                      <  2:189649/95‑1 (MQ=255)                                                                                    actgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcg                                                                                      >  1:189649/1‑95 (MQ=255)                                                                                        tATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCtt                                                                                                       >  1:890560/1‑74 (MQ=255)                                                                                        tATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCtt                                                                                                       <  2:890560/74‑1 (MQ=255)                                                                                        tATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                     >  2:88146/1‑76 (MQ=255)                                                                                        tATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                     <  1:88146/76‑1 (MQ=255)                                                                                        tATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTtttgt                                                                            <  1:579224/101‑1 (MQ=255)                                                                                        tATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTtttgt                                                                            >  2:579224/1‑101 (MQ=255)                                                                                         aTTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCAt                                     >  2:932929/1‑139 (MQ=255)                                                                                                  aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGAcc                            >  1:309353/1‑139 (MQ=255)                                                                                                  aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGAcc                            >  2:324060/1‑139 (MQ=255)                                                                                                           tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCg                                                                  >  1:73929/1‑92 (MQ=255)                                                                                                           tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCg                                                                  <  2:73929/92‑1 (MQ=255)                                                                                                           tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAac                   >  1:331616/1‑139 (MQ=255)                                                                                                           tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAac                   >  2:152173/1‑139 (MQ=255)                                                                                                            gCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAaca                  >  2:961974/1‑139 (MQ=255)                                                                                                             cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCtt                                                                                                       >  2:811500/1‑53 (MQ=255)                                                                                                             cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCtt                                                                                                       <  1:811500/53‑1 (MQ=255)                                                                                                             cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGtt                                                                                                >  2:1075296/1‑60 (MQ=255)                                                                                                             cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGtt                                                                                                <  1:1075296/60‑1 (MQ=255)                                                                                                             cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGt                                                                                           <  1:867497/65‑1 (MQ=255)                                                                                                             cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGt                                                                                           >  2:867497/1‑65 (MQ=255)                                                                                                             cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAAcaca                                             <  2:502033/111‑1 (MQ=255)                                                                                                             cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAAcaca                                             >  1:502033/1‑111 (MQ=255)                                                                                                             cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAacac                 >  1:232333/1‑139 (MQ=255)                                                                                                              tAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcg                                                                                      <  2:269918/69‑1 (MQ=255)                                                                                                              tAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcg                                                                                      >  1:269918/1‑69 (MQ=255)                                                                                                              tAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACt                <  1:932929/139‑1 (MQ=255)                                                                                                                     aGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCa         >  1:928409/1‑139 (MQ=255)                                                                                                                            aCCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCAGGATCGa  <  1:324060/139‑1 (MQ=255)                                                                                                                                 aaGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcg                                                                                      <  2:22657/50‑1 (MQ=255)                                                                                                                                 aaGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcg                                                                                      >  1:22657/1‑50 (MQ=255)                                                                                                                                        |                                                                                                                               AGGTCATCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGTCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCAGGATCGA  >  NC_000913/150812‑151073

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