Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,721,165 A→G L139S (TTG→TCG)  yhfK ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,721,1650AG100.0% 285.8 / NA 78L139S (TTG→TCG) yhfKconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (47/31);  total (47/31)

GCGCCAGATCCAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCAACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATGGCGGCGAGCAGCGATGCAGG  >  minE/2721098‑2721232
                                                                   |                                                                   
gcgcCAGATCCAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGAc                                                                   <  1:182620/70‑1 (MQ=255)
gcgcCAGATCCAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACg                                                                  >  1:1156749/1‑71 (MQ=255)
gcgcCAGATCCAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACg                                                                  >  1:3911820/1‑71 (MQ=255)
gcgcCAGATCCAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACg                                                                  >  1:3378157/1‑71 (MQ=255)
 cgcCAGATCCAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACgg                                                                 <  1:78569/71‑1 (MQ=255)
   ccAGATCCAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGtt                                                               <  1:2635650/71‑1 (MQ=255)
   ccAGATCCAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGtt                                                               <  1:1496619/71‑1 (MQ=255)
     aGATCCAGAACCAAAACCAGTTAAACAATCCGTACCAGAGGGTGCCCAACGCATAGATGAGCGACGGTTcc                                                             >  1:641594/1‑71 (MQ=255)
     aGATCCAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTcc                                                             >  1:2655598/1‑71 (MQ=255)
     aGATCCAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTcc                                                             <  1:2930965/71‑1 (MQ=255)
      gATCCAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTcc                                                             >  1:3086376/1‑70 (MQ=255)
       aTCCAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTccag                                                           >  1:3042476/1‑69 (MQ=255)
       aTCCAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTcc                                                             >  1:222509/1‑69 (MQ=255)
         ccAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGAc                                                                   >  1:3226738/1‑61 (MQ=255)
         ccAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGa                                                         >  1:2772031/1‑71 (MQ=255)
         ccAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGa                                                         >  1:2272595/1‑71 (MQ=255)
         ccAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGa                                                         >  1:3514163/1‑71 (MQ=255)
          cAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGAc                                                        >  1:885697/1‑71 (MQ=255)
          cAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGAc                                                        <  1:701991/71‑1 (MQ=255)
          cAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTGCCCAGa                                                         <  1:3831936/70‑1 (MQ=255)
           aGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGtt                                                               >  1:3365057/1‑63 (MQ=255)
             aaCCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCgg                                                     <  1:1456761/71‑1 (MQ=255)
              aCCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCgg                                                     >  1:3594042/1‑70 (MQ=255)
              aCCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGc                                                    <  1:1306954/71‑1 (MQ=255)
              aCCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGc                                                    <  1:2491684/71‑1 (MQ=255)
                  aaaCCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTcc                                                             <  1:3166948/58‑1 (MQ=255)
                  aaaCCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTcc                                                             <  1:1719518/58‑1 (MQ=255)
                  aaaCCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGt                                                <  1:2585032/71‑1 (MQ=255)
                  aaaCCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGt                                                <  1:1622127/71‑1 (MQ=255)
                   aaCCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTa                                               >  1:3155916/1‑71 (MQ=255)
                    aCCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTAt                                              >  1:3117428/1‑71 (MQ=255)
                     ccAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACAGGCATGTATc                                             >  1:3863188/1‑71 (MQ=255)
                      cAGTTAAACAATCCGTACCAGAGAGGGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATc                                             <  1:2976387/70‑1 (MQ=255)
                        gTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCg                                          >  1:2610567/1‑71 (MQ=255)
                           aaaCAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACg                                                                  >  1:2201477/1‑44 (MQ=255)
                           aaaCAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGcc                                        >  1:277540/1‑70 (MQ=255)
                           aaaCAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGcc                                        >  1:2835431/1‑70 (MQ=255)
                           aaaCAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGcc                                        >  1:19177/1‑70 (MQ=255)
                             aCAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCg                                                      >  1:4019413/1‑54 (MQ=255)
                             aCAATCCGTACCAGAGAGTGCCCAACGCATAGAAGAGCGACGGTTCCCAGACCGGCa                                                   >  1:1457135/1‑57 (MQ=255)
                              cAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAAc                                    >  1:1709844/1‑71 (MQ=255)
                               aaTCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGAcc                                                       >  1:3959660/1‑51 (MQ=255)
                               aaTCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACt                                   <  1:61107/71‑1 (MQ=255)
                                aTCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCaaa                                     >  1:3948889/1‑68 (MQ=255)
                                aTCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTg                                  <  1:46851/71‑1 (MQ=255)
                                 tCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTAt                                              >  1:2939889/1‑58 (MQ=255)
                                 tCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTAt                                              >  1:2339905/1‑58 (MQ=255)
                                 tCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACt                                   >  1:2486235/1‑69 (MQ=255)
                                 tCCGTACCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACt                                   >  1:2481678/1‑69 (MQ=255)
                                      aCCAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAggg                              >  1:667549/1‑69 (MQ=255)
                                       ccAGAGAGTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGt                             <  1:34032/69‑1 (MQ=255)
                                         agagagTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTggg                                >  1:3450639/1‑62 (MQ=255)
                                         agagagTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAgg                               >  1:532898/1‑65 (MQ=255)
                                         agagagTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTaaa                          >  1:3784845/1‑70 (MQ=255)
                                           agagTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAAt                       >  1:3430012/1‑71 (MQ=255)
                                             agTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATgg                     <  1:2861917/71‑1 (MQ=255)
                                             agTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATgg                     <  1:217747/71‑1 (MQ=255)
                                             agTGCCCAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATgg                     <  1:8174/71‑1 (MQ=255)
                                                 cccAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTaaa                          >  1:2104196/1‑62 (MQ=255)
                                                  ccAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATGgcggcg                <  1:3649905/71‑1 (MQ=255)
                                                   cAACGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATggcggc                 >  1:3446706/1‑69 (MQ=255)
                                                    aaCGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATGGCGGCGag              >  1:3676446/1‑71 (MQ=255)
                                                    aaCGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATGGCGGCGag              >  1:1411707/1‑71 (MQ=255)
                                                    aaCGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATGGCGGCGag              >  1:1416322/1‑71 (MQ=255)
                                                     aCGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAAc                                    <  1:2600276/48‑1 (MQ=255)
                                                     aCGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATggc                    >  1:2848544/1‑64 (MQ=255)
                                                     aCGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATGGCGGCGagc             <  1:2644201/71‑1 (MQ=255)
                                                     aCGCATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGACAAAc                                    <  1:3574417/48‑1 (MQ=255)
                                                      cGCTTAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATGGCGGCGag              >  1:971906/1‑69 (MQ=255)
                                                        cATAGATGAGCGCCGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATGGCGGCGagcagc          >  1:980721/1‑71 (MQ=255)
                                                        cATAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATGGCGGCGagcagc          >  1:921199/1‑71 (MQ=255)
                                                         aTAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATGGCGGCGagcagc          <  1:422568/70‑1 (MQ=255)
                                                         aTAGATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATGGCGGCGagcagc          >  1:3159738/1‑70 (MQ=255)
                                                            gATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCa                                       <  1:1784993/38‑1 (MQ=255)
                                                            gATGAGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATGGCGGCGAGCAGCGAt       <  1:3388922/70‑1 (MQ=255)
                                                                aGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATGGCGGCGAGCAGCGATGCAgg  <  1:2838180/71‑1 (MQ=255)
                                                                aGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATGGCGGCGAGCAGCGATGCAgg  <  1:1732115/71‑1 (MQ=255)
                                                                aGCGACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATGGCGGCGAGCAGCGATGCAgg  <  1:1201352/71‑1 (MQ=255)
                                                                   |                                                                   
GCGCCAGATCCAGAACCAAAACCAGTTAAACAATCCGTACCAGAGAGTGCCCAACGCATAGATGAGCAACGGTTCCCAGACCGGCATGTATCCCGCCAAACTGAGGGTAAAAATGGCGGCGAGCAGCGATGCAGG  >  minE/2721098‑2721232

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: