Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,563,587 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
*AP0090484,563,5870TA100.0% 178.2 / NA 59K207M (AAG→ATG) iadAisoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (30/29);  total (30/29)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGG  >  AP009048/4563520‑4563652
                                                                   |                                                                 
ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:199030/71‑1 (MQ=255)
ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:4266933/71‑1 (MQ=255)
 gtaggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                               <  1:1426843/68‑1 (MQ=255)
 gTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                               <  1:3684262/71‑1 (MQ=255)
  tCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                              <  1:2222115/71‑1 (MQ=255)
    ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                             <  1:2908896/70‑1 (MQ=255)
      cagcagCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGt                                                          >  1:4494414/1‑71 (MQ=255)
           gcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                                      >  1:2838317/1‑70 (MQ=255)
            cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGccc                                                     >  1:676889/1‑70 (MQ=255)
               gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                 >  1:865899/1‑71 (MQ=255)
                cTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:891339/55‑1 (MQ=255)
                   aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                               <  1:1971273/53‑1 (MQ=255)
                   aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                                >  1:2789754/1‑68 (MQ=255)
                   aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                                >  1:3148754/1‑68 (MQ=255)
                   aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                             <  1:2836812/71‑1 (MQ=255)
                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           >  1:3117901/1‑71 (MQ=255)
                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           >  1:327546/1‑71 (MQ=255)
                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           >  1:862796/1‑71 (MQ=255)
                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           >  1:3180743/1‑71 (MQ=255)
                      ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                               >  1:1262044/1‑66 (MQ=255)
                         acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       >  1:2058666/1‑71 (MQ=255)
                         acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:1752898/71‑1 (MQ=255)
                           aTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                             >  1:3635989/1‑63 (MQ=255)
                            tCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                                      >  1:692819/1‑53 (MQ=255)
                            tCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGTTGAc                                    <  1:2211458/71‑1 (MQ=255)
                              gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                                <  1:3443822/57‑1 (MQ=255)
                               cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                                  >  1:3243068/1‑70 (MQ=255)
                                aGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCa                                <  1:1672440/71‑1 (MQ=255)
                                  ttttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                             <  1:2136953/40‑1 (MQ=255)
                                  ttttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                             <  1:1018004/40‑1 (MQ=255)
                                  ttttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                                <  1:3263267/53‑1 (MQ=255)
                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACgg                                        >  1:339434/1‑59 (MQ=255)
                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                               >  1:97542/1‑68 (MQ=255)
                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                            >  1:2257499/1‑71 (MQ=255)
                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt                             <  1:1523504/70‑1 (MQ=255)
                                     tCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                            <  1:4046253/70‑1 (MQ=255)
                                      ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                           <  1:2883340/70‑1 (MQ=255)
                                       cagcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTa                          <  1:2487251/70‑1 (MQ=255)
                                       cagcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAc                         <  1:2377061/71‑1 (MQ=255)
                                          cagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgc                      >  1:4353286/1‑71 (MQ=255)
                                           agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgcc                     <  1:4013018/71‑1 (MQ=255)
                                            ggTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                    <  1:1349257/39‑1 (MQ=255)
                                               cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAt                                                   >  1:4178144/1‑37 (MQ=255)
                                               cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg                  >  1:793561/1‑70 (MQ=255)
                                                aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGc                 >  1:2930223/1‑70 (MQ=255)
                                                aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCa                <  1:1259729/71‑1 (MQ=255)
                                                  aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                               <  1:3269300/54‑1 (MQ=255)
                                                   aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAacac                                          >  1:2843663/1‑42 (MQ=255)
                                                   aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                                     >  1:1034177/1‑47 (MQ=255)
                                                   aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGCCGCCAGGTTTACCGCCGAGCAAAc             >  1:2794532/1‑71 (MQ=255)
                                                   aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaaa              <  1:1335121/70‑1 (MQ=255)
                                                   aaTAGGCGGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc             <  1:2484094/71‑1 (MQ=255)
                                                       ggCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCTCCGAGCAAACCGcc         <  1:1299799/71‑1 (MQ=255)
                                                       ggCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc         <  1:2623948/71‑1 (MQ=255)
                                                       ggCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc         <  1:867430/71‑1 (MQ=255)
                                                          tGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                           <  1:3437929/50‑1 (MQ=255)
                                                          tGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAc      >  1:4282352/1‑71 (MQ=255)
                                                           gTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc         <  1:3906858/67‑1 (MQ=255)
                                                            tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                                     >  1:3439726/1‑38 (MQ=255)
                                                            tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAgg                              <  1:4115695/45‑1 (MQ=255)
                                                            tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc    >  1:1582778/1‑71 (MQ=255)
                                                            tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc    >  1:2410313/1‑71 (MQ=255)
                                                            tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc    >  1:3409162/1‑71 (MQ=255)
                                                              aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc             <  1:2489404/60‑1 (MQ=255)
                                                              aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg  <  1:1160938/71‑1 (MQ=255)
                                                                   |                                                                 
GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGG  >  AP009048/4563520‑4563652

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