Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,343,422	C→T	A250T (GCT→ACT)	yfaL ←	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,343,422	0	C	T	100.0%	174.1 / NA	58	A250T (GCT→ACT)	yfaL	adhesin
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/30);  total (28/30)

AGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGC  >  NC_000913/2343287‑2343556                                                                                                                                        |                                                                                                                                       aGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGt                                                                                                                                     <  1:436103/139‑1 (MQ=255)     gagaTACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACg                                                                                                                                 <  2:685519/139‑1 (MQ=255)        aTACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTa                                                                                                                              <  1:157323/139‑1 (MQ=255)                     cccGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCg                                                                                                                 >  1:557264/1‑139 (MQ=255)                     cccGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCg                                                                                                                 >  1:407346/1‑139 (MQ=255)                      ccgccgCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGc                                                                                                                >  2:635278/1‑139 (MQ=255)                       cgccgcAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGcc                                                                                                               <  2:557264/139‑1 (MQ=255)                             agagGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTACTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGTTAGGGAGTCGTCACGCTATTGTTATATTTCGCCCCGCCa                                                                                                         <  2:271591/139‑1 (MQ=255)                               aggaaggaCCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTa                                                                                                                              >  1:454415/1‑116 (MQ=255)                               aggaaggaCCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTa                                                                                                                              <  2:454415/116‑1 (MQ=255)                                         aTCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAAAAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTAAATTTCGCCCCGCCATAACCTTCGGCa                                                                                             <  1:173903/139‑1 (MQ=255)                                                  gCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTAttt                                                                                    <  2:407346/139‑1 (MQ=255)                                                     tgctgcGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACg                                                                                                                                 <  1:263904/91‑1 (MQ=255)                                                     tgctgcGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACg                                                                                                                                 >  2:263904/1‑91 (MQ=255)                                                     tgctgcGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGc                                                                                              <  2:555766/126‑1 (MQ=255)                                                     tgctgcGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGc                                                                                              >  1:555766/1‑126 (MQ=255)                                                      gctgcGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTa                                                                                                                              >  1:529534/1‑93 (MQ=255)                                                      gctgcGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTa                                                                                                                              <  2:529534/93‑1 (MQ=255)                                                      gctgcGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTaa                                                                                >  1:946253/1‑139 (MQ=255)                                                           gcTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCa                                                                           <  2:946253/139‑1 (MQ=255)                                                                    ctcATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTAttt                                                                  >  1:850603/1‑139 (MQ=255)                                                                       aTCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAgtg                                                                        >  2:769753/1‑130 (MQ=255)                                                                       aTCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATAAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAgtg                                                                        <  1:769753/130‑1 (MQ=255)                                                                         cGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCa                                                                           >  1:55835/1‑125 (MQ=255)                                                                         cGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCa                                                                           <  2:55835/125‑1 (MQ=255)                                                                            cTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCa                                                                                             <  2:951790/104‑1 (MQ=255)                                                                            cTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCa                                                                                             >  1:951790/1‑104 (MQ=255)                                                                              gacacGCCTCCGTTCTGGTTTTAGCTGTCATAAAAAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGtat                                                        <  1:261771/138‑1 (MQ=255)                                                                              aaCCCGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGtat                                                        <  2:51018/139‑1 (MQ=255)                                                                               acacGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGtata                                                       >  1:828384/1‑139 (MQ=255)                                                                                     cTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCa                                                                                                                                   >  1:923638/1‑57 (MQ=255)                                                                                     cTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCa                                                                                                                                   <  2:923638/57‑1 (MQ=255)                                                                                         gTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTAtt                                                                   <  2:717478/117‑1 (MQ=255)                                                                                         gTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTAtt                                                                   >  1:717478/1‑117 (MQ=255)                                                                                          ttCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGtt                                                                      >  2:695869/1‑113 (MQ=255)                                                                                          ttCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGac                                            >  2:261771/1‑138 (MQ=255)                                                                                          ttCTGGCTGTAGCTGTCATCAAAAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGATATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGtt                                                                      <  1:695869/113‑1 (MQ=255)                                                                                               gctgtagctgtaATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGAtgcttg                                       <  1:970700/139‑1 (MQ=255)                                                                                                  gTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTa                                                                                                                              >  1:248525/1‑49 (MQ=255)                                                                                                  gTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTa                                                                                                                              <  2:248525/49‑1 (MQ=255)                                                                                                   tAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGt                                   >  1:495345/1‑139 (MQ=255)                                                                                                    aGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTa                                                                 >  2:1001685/1‑108 (MQ=255)                                                                                                    aGCTGTCATCAACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTa                                                                 <  1:1001685/108‑1 (MQ=255)                                                                                                         tcatcaACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGACATAACCTTCGGCAGTGTTa                                                                                       >  1:43595/1‑81 (MQ=255)                                                                                                         tcatcaACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGACATAACCTTCGGCAGTGTTa                                                                                       <  2:43595/81‑1 (MQ=255)                                                                                                           atcaACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATACATCGCCCCGCCa                                                                                                         >  1:434420/1‑61 (MQ=255)                                                                                                           atcaACAGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATACATCGCCCCGCCa                                                                                                         <  2:434420/61‑1 (MQ=255)                                                                                                               aaaGAAATATCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCtattattat                       <  1:632441/137‑1 (MQ=255)                                                                                                                     atatCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACt                                                               >  2:989701/1‑93 (MQ=255)                                                                                                                     atatCAATAAGATAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACt                                                               <  1:989701/93‑1 (MQ=255)                                                                                                                                aTAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGt                                                                                  >  1:452354/1‑63 (MQ=255)                                                                                                                                aTAGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGt                                                                                  <  2:452354/63‑1 (MQ=255)                                                                                                                                  aGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCgg                                                                                               >  1:396458/1‑48 (MQ=255)                                                                                                                                  aGGGAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCgg                                                                                               <  2:396458/48‑1 (MQ=255)                                                                                                                                    ggAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGAtgcttgct                                     <  2:593115/104‑1 (MQ=255)                                                                                                                                    ggAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGAtgcttgct                                     >  1:593115/1‑104 (MQ=255)                                                                                                                                    ggAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCcgc  >  1:5671/1‑139 (MQ=255)                                                                                                                                    ggAGTCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCcgc  >  1:340030/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                       AGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGC  >  NC_000913/2343287‑2343556

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