Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,026,281 A→C 6.6% R197R (CGT→CGG) ‡ malI ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,026,2810AC 6.6% 218.4 / 11.3 91R197R (CGT→CGG) ‡malIDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (37/48);  new base C (4/2);  total (41/50)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGC  >  minE/1026216‑1026347
                                                                 |                                                                  
cAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGt                                                               >  1:1864089/1‑71 (MQ=255)
cAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGt                                                               >  1:487951/1‑71 (MQ=255)
      cTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAAt                                                           <  1:1355247/69‑1 (MQ=255)
        aTGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAg                                                        >  1:1974277/1‑70 (MQ=255)
         tGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAgg                                                       >  1:407766/1‑70 (MQ=255)
          gAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCa                                                                    >  1:2277226/1‑56 (MQ=255)
           aaaCGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGt                                                               >  1:3535286/1‑60 (MQ=255)
           aaaCGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAAt                                                           >  1:1074950/1‑64 (MQ=255)
           aaaCGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGaa                                                     <  1:731033/70‑1 (MQ=255)
           aaaCGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAAc                                                    >  1:3959190/1‑71 (MQ=255)
           aaaCGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAAc                                                    <  1:664021/71‑1 (MQ=255)
           aaaCGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAAc                                                    <  1:3346185/71‑1 (MQ=255)
             actgCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACtt                                                  >  1:3860929/4‑71 (MQ=255)
              cGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACggg                                                                >  1:800629/1‑56 (MQ=255)
              cGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACttt                                                 >  1:1507395/1‑71 (MQ=255)
              cGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACttt                                                 >  1:2320379/1‑71 (MQ=255)
              cGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCAAACCCGTTCTGCACGGGTTAATGAGGAACttt                                                 >  1:545630/1‑71 (MQ=255)
               ggCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTg                                                <  1:4351545/71‑1 (MQ=255)
               ggCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTg                                                <  1:3133852/71‑1 (MQ=255)
               ggCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTg                                                <  1:1143929/71‑1 (MQ=255)
                gCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGc                                               <  1:4486738/71‑1 (MQ=255)
                 cAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGGTAATGAGGAACTTTGcc                                              >  1:3935471/1‑71 (MQ=255)
                    gCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCa                                                                    <  1:3057346/46‑1 (MQ=255)
                    gCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTc                                           <  1:3675724/71‑1 (MQ=255)
                    gCCAAATTTTAGTAGAGTCGCACAATAGCCCCCCACCCGTTCTGCCAGGGTTAAt                                                           <  1:2481311/55‑1 (MQ=255)
                    gCCAAATTTTAGTAGAGTCGCACAATAGCCCCCCACCCGTTCTGCCAGGGTTAAt                                                           <  1:2486913/55‑1 (MQ=255)
                     ccAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTcc                                          >  1:2620407/1‑71 (MQ=255)
                      cAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTccc                                         <  1:4147251/71‑1 (MQ=255)
                      cAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTcc                                          >  1:404403/1‑70 (MQ=255)
                        aaTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAgg                                                       <  1:1467789/55‑1 (MQ=255)
                        aaTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAgg                                                       <  1:320617/55‑1 (MQ=255)
                         aTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGATGAACTTTGCCCTCccagc                                      <  1:2925210/71‑1 (MQ=255)
                         aTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCccagc                                      <  1:3876115/71‑1 (MQ=255)
                         aTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCccagc                                      <  1:3498971/71‑1 (MQ=255)
                                tagAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGa                               <  1:1224232/71‑1 (MQ=255)
                                tagAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGa                               <  1:3192031/71‑1 (MQ=255)
                                  gagTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATc                             >  1:3726274/1‑71 (MQ=255)
                                      tGCACAATAGCCCCCAAGCCGTTCTGCACGGGTTAATGAGGaa                                                     <  1:2491385/43‑1 (MQ=255)
                                       gCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                         <  1:4476218/70‑1 (MQ=255)
                                       gCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                         <  1:3480297/70‑1 (MQ=255)
                                       gCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                         >  1:1746715/1‑70 (MQ=255)
                                       gCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                         <  1:3031291/70‑1 (MQ=255)
                                       gCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                         <  1:2430290/70‑1 (MQ=255)
                                       gCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                         >  1:1547410/1‑70 (MQ=255)
                                       gCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                         >  1:3836022/1‑70 (MQ=255)
                                        cacaATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                         <  1:4178244/69‑1 (MQ=255)
                                         acaATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGAt                      <  1:2698070/71‑1 (MQ=255)
                                          caATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTc                                           <  1:3155229/49‑1 (MQ=255)
                                          caATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATg                     >  1:4210805/1‑71 (MQ=255)
                                           aaTAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                         <  1:2956926/66‑1 (MQ=255)
                                            aTAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGcc                   <  1:3759363/71‑1 (MQ=255)
                                             tAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGa                       >  1:4620860/1‑66 (MQ=255)
                                              aGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCa                 >  1:2711161/1‑71 (MQ=255)
                                               gCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCAt                >  1:1743139/1‑71 (MQ=255)
                                               gCCCCCAACCCGTTCTGCACGGGTAAATGAGGAACTTTGccc                                             >  1:2004074/1‑42 (MQ=255)
                                                  cccAACCCGTTCTGCACGGGTTAATGAGGAACttt                                                 <  1:757539/35‑1 (MQ=255)
                                                  cccAACCCGTTCTGCACGGGTTAATGAGGAACttt                                                 <  1:1413542/35‑1 (MQ=255)
                                                  cccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                         <  1:4602983/59‑1 (MQ=255)
                                                  cccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                         <  1:2309540/59‑1 (MQ=255)
                                                  cccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTg                        <  1:4460691/60‑1 (MQ=255)
                                                  cccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCAtttc             <  1:2293929/71‑3 (MQ=255)
                                                   ccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                         <  1:529317/58‑1 (MQ=255)
                                                   ccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCg            <  1:3883145/71‑1 (MQ=255)
                                                   ccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCg            <  1:3430805/71‑1 (MQ=255)
                                                   ccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCg            <  1:2748201/71‑1 (MQ=255)
                                                    ccaCCCGTTCTGCCAGGGTTAATGAGGAACTTTGTCCTCCCAGCCAGGCGATACGCTGATGCCCATTGCg            >  1:3033247/3‑70 (MQ=255)
                                                    cAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATc                             <  1:1416749/53‑1 (MQ=255)
                                                    cAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCg            >  1:2843065/1‑70 (MQ=255)
                                                     aaCCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                         <  1:2210166/56‑1 (MQ=255)
                                                     aaCCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                         <  1:6938/56‑1 (MQ=255)
                                                     aaCCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCa                 >  1:4238012/1‑64 (MQ=255)
                                                     aaCCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCg            >  1:1451187/1‑69 (MQ=255)
                                                     aaCCCGTTCTGCACGGGTTAATGAGGAACTTTGCACTCCCAGCCAGGCGAt                              >  1:3673000/1‑51 (MQ=255)
                                                      aCCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCccagc                                      <  1:4091481/42‑1 (MQ=255)
                                                       cccGTTCTGCCAGGGTTAATGAGGAACTTTGTCCTCCCAGCCAGGCGATACGCTGATGCCCATTGCGAAt         >  1:892211/1‑70 (MQ=255)
                                                       cccGTTCTGCCAGGGTTAATGAGGAACTTTGTCCTCCCAGCCAGGCGATACGCTGATGCCCATTGCGAAt         >  1:2814333/1‑70 (MQ=255)
                                                       cccGTTCTGCCAGGGTTAATGAGGAACTTTGTCCTCCCAGCCAGGCGATACGCTGATGCCCATTGCGAAt         >  1:2800439/1‑70 (MQ=255)
                                                       cccGTTCTGCACGGGTTAATGAGGAACTTTGGCCTCccagccag                                   <  1:292527/44‑1 (MQ=255)
                                                       cccGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCccagccat                                   <  1:1157882/44‑2 (MQ=255)
                                                       cccGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGAt                              >  1:455471/1‑49 (MQ=255)
                                                        ccGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAAt         >  1:1819795/1‑69 (MQ=255)
                                                        ccGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAAt         >  1:3622447/1‑69 (MQ=255)
                                                        ccGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATga       <  1:2103489/71‑1 (MQ=255)
                                                         cGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGc                    <  1:1242914/57‑1 (MQ=255)
                                                         cGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAAt         >  1:4363192/1‑68 (MQ=255)
                                                           ttCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                         <  1:2134950/50‑1 (MQ=255)
                                                           ttCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGcc                   <  1:1232607/56‑1 (MQ=255)
                                                           ttCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCAt                >  1:3504996/1‑59 (MQ=255)
                                                           ttCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGACCAt                >  1:2591618/1‑59 (MQ=255)
                                                           ttCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCAAGGCGATCCGCt                         <  1:3488346/50‑1 (MQ=255)
                                                             cTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATc                             >  1:3011095/1‑44 (MQ=255)
                                                             cTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATgag      >  1:3230069/1‑67 (MQ=255)
                                                             cTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGc  >  1:1800747/1‑71 (MQ=255)
                                                                 |                                                                  
CAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGC  >  minE/1026216‑1026347

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: