Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,563,589 T→A 100% K207M (AAG→ATG)  iadA ← isoaspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,563,5890TA100.0% 169.4 / NA 58K207M (AAG→ATG) iadAisoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (31/27);  total (31/27)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGA  >  W3110S.gb/4563522‑4563656
                                                                   |                                                                   
tgTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                  <  1:154319/70‑1 (MQ=255)
ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                  <  1:430884/71‑1 (MQ=255)
   cGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                               >  1:2699073/1‑71 (MQ=255)
     gcagcaGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTg                                                             >  1:3517532/1‑71 (MQ=255)
         cagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGc                                                         >  1:3030493/1‑71 (MQ=255)
           gcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGccc                                                       <  1:3632234/71‑1 (MQ=255)
            cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                      >  1:207069/1‑71 (MQ=255)
               gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                    <  1:145110/70‑1 (MQ=255)
               gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                    >  1:2624245/1‑70 (MQ=255)
               gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                    <  1:3066923/70‑1 (MQ=255)
                   aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                                        <  1:4402721/62‑1 (MQ=255)
                   aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                                 <  1:562490/69‑1 (MQ=255)
                   aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                                 <  1:3115212/69‑1 (MQ=255)
                      ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCTATTTTACTGTCGCCCATGTGGAacac                                            <  1:4580904/71‑1 (MQ=255)
                      ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                  <  1:4404042/49‑1 (MQ=255)
                      ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAacac                                            <  1:945926/71‑1 (MQ=255)
                       gCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACg                                           <  1:3633641/71‑1 (MQ=255)
                       gCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACg                                           >  1:2961780/1‑71 (MQ=255)
                             cGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCg                                                          >  1:4123070/1‑50 (MQ=255)
                             cGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCg                                                          >  1:3875312/1‑50 (MQ=255)
                             cGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                             >  1:306051/1‑63 (MQ=255)
                             cGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                             >  1:2156251/1‑63 (MQ=255)
                             cGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                             >  1:2680938/1‑63 (MQ=255)
                              gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGt                                                            >  1:3236306/1‑47 (MQ=255)
                              gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                                    >  1:2398974/1‑71 (MQ=255)
                               cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                   <  1:5031313/71‑1 (MQ=255)
                                  ttttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                   >  1:4152829/1‑52 (MQ=255)
                                  ttttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGAc                                      <  1:4133977/65‑1 (MQ=255)
                                   tttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                                       >  1:1946368/1‑63 (MQ=255)
                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                                        >  1:3514647/1‑45 (MQ=255)
                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                              >  1:4237505/1‑71 (MQ=255)
                                      ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                             <  1:1097557/70‑1 (MQ=255)
                                      ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                             <  1:573533/70‑1 (MQ=255)
                                      ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                             <  1:1266680/70‑1 (MQ=255)
                                        agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTg                                        >  1:1210982/1‑57 (MQ=255)
                                        agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAcc                          >  1:1370927/1‑71 (MQ=255)
                                        agcagGTCATAAATAGGCTGTAACGCCATTCTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAcc                          >  1:4334806/1‑71 (MQ=255)
                                           agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgcc                       <  1:5042374/71‑1 (MQ=255)
                                            ggTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                      <  1:1387413/39‑1 (MQ=255)
                                             gTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                         <  1:2340970/51‑1 (MQ=255)
                                             gTCATAAATAGGCTGTAACGCCATTTTACTCTCGCCCATGTGGAACACGGt                                         <  1:4011457/51‑1 (MQ=255)
                                               cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                   <  1:1682570/39‑1 (MQ=255)
                                               cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg                    >  1:1625908/1‑70 (MQ=255)
                                               cATAAATAGGCTGTAACGCCATTTTACTGTCCCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg                    >  1:2452200/1‑70 (MQ=255)
                                                aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                               <  1:4534403/42‑1 (MQ=255)
                                                  aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaaa                >  1:701390/1‑71 (MQ=255)
                                                  aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaaa                >  1:1947642/1‑71 (MQ=255)
                                                   aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc               >  1:1109093/1‑71 (MQ=255)
                                                   aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc               >  1:2704648/1‑71 (MQ=255)
                                                    aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                             <  1:2123373/40‑1 (MQ=255)
                                                    aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAgg                                >  1:258011/1‑53 (MQ=255)
                                                         cTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCaa         <  1:4158865/71‑1 (MQ=255)
                                                         cTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCaa         <  1:3439452/71‑1 (MQ=255)
                                                         cTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCaa         <  1:863818/71‑1 (MQ=255)
                                                          tGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCaa         >  1:403441/1‑70 (MQ=255)
                                                            tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc      <  1:3953239/71‑1 (MQ=255)
                                                              aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                                 >  1:3556288/1‑42 (MQ=255)
                                                              aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg    <  1:1652313/71‑1 (MQ=255)
                                                               cGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc              <  1:438694/60‑1 (MQ=255)
                                                                gCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCaa         >  1:2709802/1‑64 (MQ=255)
                                                                gCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGa  >  1:3158576/1‑71 (MQ=255)
                                                                   |                                                                   
GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGA  >  W3110S.gb/4563522‑4563656

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