Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,225,405 G→A 21.0% E423K (GAG→AAG)  yebU → predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,225,4050GA21.0% 125.6 / 28.5 57E423K (GAG→AAG) yebUpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (25/20);  new base A (4/8);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

AAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTC  >  minE/1225337‑1225473
                                                                    |                                                                    
aaGCGGTTATTGCTCTTGCCTCACCCTATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGgaggag                                                                    <  1:730935/71‑1 (MQ=255)
aaGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGgaggag                                                                    <  1:2733733/71‑1 (MQ=255)
    ggTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGCAGGAGTgg                                                                 >  1:2589597/1‑70 (MQ=255)
       tATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATc                                                             >  1:2039508/1‑71 (MQ=255)
        aTTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTg                                                                  >  1:982388/1‑65 (MQ=255)
        aTTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATcg                                                            >  1:314537/1‑71 (MQ=255)
        aTTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGAAGTGGTATcg                                                            >  1:1167567/1‑71 (MQ=255)
         ttGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTg                                                                  <  1:3255806/64‑1 (MQ=255)
          tGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCggagga                                                                     >  1:726360/1‑60 (MQ=255)
          tGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCggagga                                                                     >  1:3868287/1‑60 (MQ=255)
          tGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATcgcg                                                          >  1:703842/1‑71 (MQ=255)
           gCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCgg                                                         >  1:1509755/1‑71 (MQ=255)
                 gCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATcg                                                            >  1:1970779/1‑62 (MQ=255)
                 gCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATcg                                                            >  1:1144685/1‑62 (MQ=255)
                 gCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATcg                                                            >  1:2918469/1‑62 (MQ=255)
                 gCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGgcgc                                                     <  1:3937363/69‑1 (MQ=255)
                 gCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGgcgc                                                     <  1:578203/69‑1 (MQ=255)
                 gCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGcgcg                                                    <  1:2051049/70‑1 (MQ=255)
                 gCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGa                                                   >  1:3785461/1‑71 (MQ=255)
                  ccTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGAt                                                  <  1:950076/71‑1 (MQ=255)
                     cACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGtt                                               >  1:714986/1‑71 (MQ=255)
                      aCCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGTTATCGCGGGCGCGATGtt                                               <  1:3419134/70‑1 (MQ=255)
                      aCCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGtt                                               >  1:3394374/1‑70 (MQ=255)
                      aCCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGtt                                               <  1:3699510/70‑1 (MQ=255)
                       cccGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTa                                             <  1:2469223/71‑1 (MQ=255)
                        ccGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTa                                             <  1:1669722/70‑1 (MQ=255)
                        ccGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGAAGTGGTATCGCGGGcgcg                                                    <  1:333467/63‑1 (MQ=255)
                          gATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGAAGTGGTATCGCGGGcgcg                                                    <  1:2820768/61‑1 (MQ=255)
                           ataataTGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCg                                         <  1:3653605/71‑1 (MQ=255)
                            taataTGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCg                                         >  1:400624/1‑70 (MQ=255)
                            taataTGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGAAGTGGTATCGCGGGCGCGATGTTTACCCGc                                        <  1:654900/71‑1 (MQ=255)
                            taataTGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGAAGTGGTATCGCGGGCGCGATGTTTACCCGc                                        <  1:3764455/71‑1 (MQ=255)
                            taataTGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGAAGTGGTATCGCGGGCGCGATGTTTACCCGc                                        <  1:1828200/71‑1 (MQ=255)
                             aataTGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTTGTATCGCGGGCGCGATGTTTACCCGc                                        >  1:332923/1‑70 (MQ=255)
                              atatGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGgc                                                       >  1:3810813/1‑54 (MQ=255)
                              atatGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCa                                       <  1:1607565/70‑1 (MQ=255)
                              atatGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCa                                       <  1:2704457/70‑1 (MQ=255)
                               tatGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTa                                             <  1:3250986/63‑1 (MQ=255)
                               tatGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAg                                     <  1:2641074/71‑1 (MQ=255)
                               tatGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGAAGt                                                                   <  1:866870/41‑1 (MQ=255)
                                atGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGAAGt                                                                   <  1:1553111/40‑1 (MQ=255)
                                 tGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGAAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGc                                    >  1:3077547/1‑70 (MQ=255)
                                     tGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCgcgc                               <  1:419998/71‑1 (MQ=255)
                                       cATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCaa                                      >  1:3545878/1‑62 (MQ=255)
                                       cATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCaa                                      >  1:3125224/1‑62 (MQ=255)
                                       cATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCgcgc                               >  1:3136160/1‑69 (MQ=255)
                                       cATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGcc                              <  1:3221444/70‑1 (MQ=255)
                                          tCGAGCTGACACCGCAGGAAGCGGAGAAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTg                          >  1:3726788/1‑71 (MQ=255)
                                           cGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGt                                                >  1:919691/1‑48 (MQ=255)
                                           cGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGt                                                >  1:105921/1‑48 (MQ=255)
                                           cGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTg                          <  1:4067730/70‑1 (MQ=255)
                                           cGAGCAGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCa                             >  1:2834973/1‑67 (MQ=255)
                                               cTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGga                     >  1:2420626/1‑71 (MQ=255)
                                                  acacCGCAGGAAGCGGAGAAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGgatg                   <  1:3132206/70‑1 (MQ=255)
                                                       gCAGGAAGCGGAGAAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTa              >  1:2513493/1‑70 (MQ=255)
                                                        cAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGttt                                              >  1:1557821/1‑37 (MQ=255)
                                                         aGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTg           >  1:3253658/1‑71 (MQ=255)
                                                           gAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGatgatg                <  1:165314/64‑1 (MQ=255)
                                                           gAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGt         <  1:3136986/71‑1 (MQ=255)
                                                           gAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGt         <  1:2977688/71‑1 (MQ=255)
                                                               cGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACt     >  1:3145327/1‑71 (MQ=255)
                                                                  aggagTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTc  >  1:1016731/1‑71 (MQ=255)
                                                                    |                                                                    
AAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTC  >  minE/1225337‑1225473

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: