Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,551,927 A→C 8.3% K178N (AAA→AAC yibD → predicted glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,551,9270AC 8.3% 149.4 / 9.4 60K178N (AAA→AACyibDpredicted glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base C (2/3);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

GGGCTTTCTTCGCGCCGTTGGACTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGACAACAGAATTCATGTTTAA  >  minE/2551859‑2551995
                                                                    |                                                                    
gggCTTTCTTCGCGCCGTTGGACTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTaaaa                                                                     <  1:2443717/70‑1 (MQ=255)
 ggCTTTCTTCGCGCCGTTGGACTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAAt                                                                   <  1:3233518/71‑1 (MQ=255)
  gCTTTCTTCGCGCCGTTGGACTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATa                                                                  <  1:2310621/71‑1 (MQ=255)
   cTTTCTTCGCGCCGTTGGACTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATaa                                                                 >  1:892137/1‑71 (MQ=255)
     ttcttcGCGCCGTTGGACTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATAACa                                                               >  1:155422/1‑71 (MQ=255)
        ttcGCGCCGTTGGACTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATaa                                                                 >  1:3218329/1‑66 (MQ=255)
        ttcGCGCCGTTGGACTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATa                                                                  >  1:2098230/1‑65 (MQ=255)
          ttcgcCGTTGGACTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATaa                                                                 >  1:1582568/3‑64 (MQ=255)
            cgcCGTTGGACTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTTTTAACAATa                                                                  <  1:653166/61‑1 (MQ=255)
            cgcCGTTGGACTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTTTTAACAATa                                                                  <  1:2093323/61‑1 (MQ=255)
            cgcCGTTGGACTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTTTTAACAATa                                                                  <  1:1856344/61‑1 (MQ=255)
             gcCGTTGGACTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATAACATTAAAttt                                                       <  1:3443236/71‑1 (MQ=255)
                 ttGGACTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTg                                                   <  1:1307972/71‑1 (MQ=255)
                 ttGGACTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTg                                                   >  1:1767405/1‑71 (MQ=255)
                   ggACTCACGTTGTCTGGATGGGGGTTTATCGCCGTAATGTTATTTTTAACAATAACAtc                                                             >  1:2725236/1‑58 (MQ=25)
                     aCTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAAt                                                                   >  1:3209822/1‑51 (MQ=255)
                     aCTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCgg                                               <  1:355752/71‑1 (MQ=255)
                       tCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTa                                                      <  1:1532797/62‑1 (MQ=255)
                       tCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGa                                              >  1:487251/1‑70 (MQ=255)
                       tCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGAt                                             >  1:3426747/1‑71 (MQ=255)
                         aCGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTAttt                                                   >  1:1567990/1‑62 (MQ=255)
                         aCGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGAtt                                            <  1:3059257/70‑1 (MQ=255)
                          cGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGATTAc                                          <  1:3479223/71‑1 (MQ=255)
                               tCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGATTAcatcat                                     <  1:2947823/71‑1 (MQ=255)
                                       ggggTTTATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGatat                             >  1:1451598/1‑71 (MQ=255)
                                          gTTTATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGt                          >  1:2888363/1‑71 (MQ=255)
                                           tttATCGCCGTGATGTTATTGTTAAAAATAACATTAAAtt                                                        <  1:1170453/40‑1 (MQ=255)
                                           tttATCGCCGTGATGTTATTGTTAAAAATAACATTAAAtt                                                        <  1:960500/40‑1 (MQ=255)
                                           tttATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGAtt                                            <  1:1092978/52‑1 (MQ=255)
                                            ttATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTGc                                                  <  1:3440323/45‑1 (MQ=255)
                                            ttATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCt                        >  1:1055664/1‑71 (MQ=255)
                                            ttATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCt                        >  1:1866762/1‑71 (MQ=255)
                                            ttATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCt                        >  1:1084617/1‑71 (MQ=255)
                                                 gCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGac                    <  1:2831132/70‑1 (MQ=255)
                                                 gCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGac                    <  1:2553805/70‑1 (MQ=255)
                                                  ccGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGATTACatcatc                                    >  1:1543071/1‑53 (MQ=255)
                                                  ccGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGATTACatcatc                                    >  1:640033/1‑53 (MQ=255)
                                                    gTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGacaa                  >  1:457525/1‑69 (MQ=255)
                                                    gTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGacaa                  >  1:954684/1‑69 (MQ=255)
                                                       aTGTTATTGTTAAAAATAACATTAAATTTATTGCCg                                                <  1:412018/36‑1 (MQ=255)
                                                       aTGTTATTGTTAAAAATAACATTAAATTTATTGCCg                                                <  1:1866034/36‑1 (MQ=255)
                                                       aTGTTATTGTTAAAAATAACATTAAATTTATTGCCg                                                <  1:232747/36‑1 (MQ=255)
                                                       aTGTTATTGTTAAAAATAACATTAAATTTATTGCCGGATTAcatc                                       >  1:2434814/1‑45 (MQ=255)
                                                       aTGTTATTGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGACAACAGaa             <  1:3429235/71‑1 (MQ=255)
                                                         gTTATTGTTAAAAATAACATTAAATTTATTGCCGGa                                              >  1:977866/1‑36 (MQ=255)
                                                          ttATTGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGACAACAGAAtt           <  1:899090/70‑1 (MQ=255)
                                                           tATTGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGACAACAGAAtt           <  1:2093576/69‑1 (MQ=255)
                                                            aTTGTTAAAAATAACATTAAATTTATTGCCGGATTAca                                         >  1:924675/1‑38 (MQ=255)
                                                            aTTGTTAAAAATAACATTAAATTTATTGCCGGATTAca                                         >  1:1145619/1‑38 (MQ=255)
                                                            aTTGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTgg                      <  1:582846/57‑1 (MQ=255)
                                                            aTTGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTgg                      <  1:2555893/57‑1 (MQ=255)
                                                            aTTGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGACAACAGAATTCAt        >  1:831674/1‑71 (MQ=255)
                                                             ttGTTAAAAATAACATTAAATTTATTGCCGGATTa                                           >  1:2269931/1‑35 (MQ=255)
                                                             ttGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGACAACAGAATTCATg       >  1:1091237/1‑71 (MQ=255)
                                                              tGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGACAACAGAATTCATGt      >  1:1088609/1‑71 (MQ=255)
                                                                ttAACAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGACAACAg               >  1:2894289/1‑60 (MQ=255)
                                                                ttAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGACAACAGAATTc          <  1:440086/65‑1 (MQ=255)
                                                                ttAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGACAACAGAATTc          <  1:3533278/65‑1 (MQ=255)
                                                                ttAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGACAACAGAATTCATGttt    <  1:3419540/71‑1 (MQ=255)
                                                                ttAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGACAACAGAATTCATGttt    <  1:2011032/71‑1 (MQ=255)
                                                                ttAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGACAACAGAATTCATGt      <  1:1940903/69‑1 (MQ=255)
                                                                  aaaaaTAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGACAACAGAATTCATGTTTaa  >  1:3575314/1‑71 (MQ=255)
                                                                  aaaaaTAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGACAACAGAATTCATGTTTaa  >  1:1039124/1‑71 (MQ=255)
                                                                  aaaaaTAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGACAACAGAATTCATGTTTa   <  1:2557407/70‑1 (MQ=255)
                                                                    |                                                                    
GGGCTTTCTTCGCGCCGTTGGACTCACGTTGTCTGGATGGGGGTTTATCGCCGTGATGTTATTGTTAAAAATAACATTAAATTTATTGCCGGATTACATCATCAGGATATTGTCTGGACAACAGAATTCATGTTTAA  >  minE/2551859‑2551995

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