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% 199.8 / NA 66K207M (AAG→ATG) iadAisoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (30/36);  total (30/36)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGG  >  W3110S.gb/4563522‑4563654
                                                                   |                                                                 
ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:3430441/71‑1 (MQ=255)
ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:2008176/71‑1 (MQ=255)
 gTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                               <  1:4657627/71‑1 (MQ=255)
   cGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                             >  1:975977/1‑71 (MQ=255)
    ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGTCTGTAACGCCATTTTACt                                                            <  1:3835234/71‑1 (MQ=255)
     gcagcaGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTg                                                           <  1:2883089/71‑1 (MQ=255)
     gcagcaGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTg                                                           >  1:4310863/1‑71 (MQ=255)
     gcagcaGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTg                                                           >  1:5157758/1‑71 (MQ=255)
       agcagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGt                                                          >  1:4756294/1‑70 (MQ=255)
       agcagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTc                                                         >  1:1357193/1‑71 (MQ=255)
            cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTc                                                         <  1:1537302/66‑1 (MQ=255)
            cTTGATGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                    >  1:2483149/1‑71 (MQ=255)
              tGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:2509356/57‑1 (MQ=255)
              tGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                             >  1:4301366/1‑60 (MQ=255)
              tGCTGATCGGCAAATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                             >  1:5249447/1‑60 (MQ=255)
               gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                 >  1:5280630/1‑71 (MQ=255)
               gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                  <  1:1890754/70‑1 (MQ=255)
               gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                  <  1:3278744/70‑1 (MQ=255)
                 tGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCg                                                        >  1:5145468/1‑62 (MQ=255)
                   aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTg                                                           >  1:5032646/1‑57 (MQ=255)
                   aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                             >  1:4385363/1‑71 (MQ=255)
                   aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                             >  1:646252/1‑71 (MQ=255)
                   aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                             >  1:1213591/1‑71 (MQ=255)
                    tCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                             >  1:3158308/1‑70 (MQ=255)
                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           >  1:811347/1‑71 (MQ=255)
                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCAATGTGGAaca                                           >  1:3005712/1‑71 (MQ=255)
                      ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:3133407/49‑1 (MQ=255)
                      ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           <  1:3304253/70‑1 (MQ=255)
                       gCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:4526513/48‑1 (MQ=255)
                         acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:2941755/71‑1 (MQ=255)
                          caTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       >  1:1388904/1‑70 (MQ=255)
                          caTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       >  1:1394972/1‑70 (MQ=255)
                             cGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACg                                   >  1:2797666/1‑71 (MQ=255)
                              gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                                  <  1:2133707/71‑1 (MQ=255)
                               cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACg                                   <  1:2808602/69‑1 (MQ=255)
                               cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                 <  1:101433/71‑1 (MQ=255)
                               cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                 <  1:3509158/71‑1 (MQ=255)
                                 gTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                               <  1:1053830/71‑1 (MQ=255)
                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt                             <  1:3624935/70‑1 (MQ=255)
                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt                             <  1:5157124/70‑1 (MQ=255)
                                     tCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCa                                >  1:2299113/1‑66 (MQ=255)
                                     tCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                            <  1:3157831/70‑1 (MQ=255)
                                      ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                  <  1:2156882/47‑1 (MQ=255)
                                      ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                           <  1:4401335/70‑1 (MQ=255)
                                       cagcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTa                          <  1:2776586/70‑1 (MQ=255)
                                        agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           <  1:387585/52‑1 (MQ=255)
                                        agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAcc                        >  1:2406693/1‑71 (MQ=255)
                                           agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgcc                     <  1:3083485/71‑1 (MQ=255)
                                            ggTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgccg                    <  1:4697869/71‑1 (MQ=255)
                                             gTCATAAATAGTCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGa                   <  1:2753145/71‑1 (MQ=255)
                                             gTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                               <  1:4257992/43‑1 (MQ=255)
                                             gTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGa                   <  1:2696687/71‑1 (MQ=255)
                                               cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                               <  1:4559813/41‑1 (MQ=255)
                                               cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg                  >  1:5194684/1‑70 (MQ=255)
                                                aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGc                 <  1:4569205/70‑1 (MQ=255)
                                                aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGc                 >  1:4261508/1‑70 (MQ=255)
                                                aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCa                <  1:4121054/71‑1 (MQ=255)
                                                 tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaa               <  1:4092872/71‑1 (MQ=255)
                                                 tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaa               <  1:5265134/71‑1 (MQ=255)
                                                  aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaaa              >  1:15087/1‑71 (MQ=255)
                                                   aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTg                                      >  1:3226949/1‑46 (MQ=255)
                                                    aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc            >  1:2096062/1‑71 (MQ=255)
                                                    aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc            <  1:2045263/71‑1 (MQ=255)
                                                    aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc            <  1:1208744/71‑1 (MQ=255)
                                                      aGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAACAAACCGc          >  1:2398400/1‑71 (MQ=255)
                                                        gCTGTAACGCCATTTTACTGTCGCCCATGTGGAacact                                         <  1:4673831/38‑2 (MQ=38)
                                                         cTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                               >  1:5129192/1‑47 (MQ=255)
                                                           gTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGACAAcg     >  1:2580784/1‑71 (MQ=255)
                                                            tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                 <  1:4157268/42‑1 (MQ=255)
                                                            tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcg     <  1:4289977/70‑1 (MQ=255)
                                                              aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg  <  1:933043/71‑1 (MQ=255)
                                                              aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGACAAcgcg   <  1:1929918/70‑1 (MQ=255)
                                                                   |                                                                 
GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGG  >  W3110S.gb/4563522‑4563654

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: