Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,395,675 T→C 7.8% F127F (TTT→TTC setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,395,6750TC 7.8% 174.8 / 5.5 64F127F (TTT→TTCsetBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base T (25/34);  new base C (3/2);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

CTACTTTGTTTTGCTATTCGTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTACGCGCT  >  minE/1395606‑1395744
                                                                     |                                                                     
ctactTTGTTTTGCTATTCGTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTttg                                                                      >  1:2390772/1‑71 (MQ=255)
ctactTTGTTTTGCTATTCGTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTttg                                                                      <  1:677794/71‑1 (MQ=255)
   ctTTGTTTTGCTATTCGTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTttgccc                                                                   >  1:3017967/1‑71 (MQ=255)
   ctTTGTTTTGCTATTCGTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTttgcc                                                                    <  1:3909760/70‑1 (MQ=255)
      tGTTTTGCTATTCGTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTttgcccttg                                                                >  1:3815234/1‑71 (MQ=255)
      tGTTTTGCTATTCGTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTttgccctt                                                                 >  1:2447194/1‑70 (MQ=255)
      tGTTTTGCTATTCGTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTttgccctt                                                                 <  1:606811/70‑1 (MQ=255)
           tGCTATTCGTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCg                                                            <  1:1879064/70‑1 (MQ=255)
              tATTCGTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGaa                                                        >  1:2355869/1‑71 (MQ=255)
               aTTCGTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTttgccc                                                                   <  1:530827/59‑1 (MQ=255)
               aTTCGTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGaa                                                        >  1:3877655/1‑70 (MQ=255)
               aTTCGTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGaa                                                        >  1:3142771/1‑70 (MQ=255)
                ttCGTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACa                                                      >  1:2269993/1‑71 (MQ=255)
                  cGTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATg                                                    >  1:56878/1‑71 (MQ=255)
                   gTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTttgcc                                                                    >  1:2839136/1‑54 (MQ=255)
                   gTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACAt                                                     <  1:2236280/69‑1 (MQ=255)
                    cgggCGTACTTCTGAGTAGTTTTGCTTCCACCGCAAACCCGCAAATGTTCGCCCTCGCCCGTGAACACGcc                                                  >  1:19748/3‑71 (MQ=255)
                            tttcttAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGcc                                                  <  1:2619620/63‑1 (MQ=255)
                            tttcttAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAAc                                          >  1:2334121/1‑71 (MQ=255)
                            cttctgAGTAGTTTTGCTTCCACCGCAAACCCGCAAATGTtcgccct                                                                  >  1:1451895/2‑47 (MQ=255)
                             ttcttAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAAc                                                       <  1:2137535/57‑1 (MQ=255)
                             ttcttAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAAcc                                         >  1:466338/1‑71 (MQ=255)
                             ttcttAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAAcc                                         >  1:3868049/1‑71 (MQ=255)
                             ttctgAGTAGTTTTGCTTCCACCGCAAACCCGCAAATGTTCGCCCTCGCCCGTGaa                                                        <  1:3101074/56‑1 (MQ=255)
                              tcttAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAAcc                                         <  1:2629621/70‑1 (MQ=255)
                              tcttAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAAc                                          <  1:1861699/69‑1 (MQ=255)
                              tcttAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCg                                        >  1:2452205/1‑71 (MQ=255)
                               cttAGCAGCTTTGGCTGGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCg                                        <  1:3744246/70‑1 (MQ=255)
                               cttAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTttgcccttgccc                                                             <  1:3381556/49‑1 (MQ=255)
                               cttAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGa                                                         <  1:3887590/53‑1 (MQ=255)
                                   gcagcTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCg                                        >  1:815844/1‑66 (MQ=255)
                                   gcagcTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGt                                   <  1:2099785/71‑1 (MQ=255)
                                   gcagcTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGt                                   <  1:1834735/71‑1 (MQ=255)
                                   gcagcTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGt                                   <  1:3876625/71‑1 (MQ=255)
                                      gcTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACa                                                      <  1:2257204/49‑1 (MQ=255)
                                        tttGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACaa                                             <  1:2550579/56‑1 (MQ=255)
                                        tttGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGt                                   <  1:80215/66‑1 (MQ=255)
                                        tttGGCTCGACCGCTAACCCG‑‑AATGTTTGCCCTTGCCCGTGAACATGCCGACaaa                                            <  1:3818484/55‑1 (MQ=255)
                                         ttGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGaa                                                        <  1:1013837/44‑1 (MQ=255)
                                          tGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGc                              <  1:398573/69‑1 (MQ=255)
                                          tGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCg                             <  1:3145439/70‑1 (MQ=255)
                                           ggCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGt                           >  1:3757543/1‑71 (MQ=255)
                                            gCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGTCGTGAGGCGGt                           <  1:1883105/70‑1 (MQ=255)
                                              tCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACaaa                                            >  1:3579251/1‑51 (MQ=255)
                                              tCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGa                         <  1:51195/70‑1 (MQ=255)
                                               cGACCGCTAACCCGGAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTg                          >  1:3791775/1‑68 (MQ=255)
                                               cGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGAt                        >  1:54471/1‑70 (MQ=255)
                                               cGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATg                       >  1:1513304/1‑71 (MQ=255)
                                                gACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGAt                        <  1:3629887/69‑1 (MQ=255)
                                                gACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGt                      >  1:3622715/1‑71 (MQ=255)
                                                 aCCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGtt                     <  1:1174679/71‑1 (MQ=255)
                                                 aCCGCAAACCCGCAAATGTTCGCCCTCGCCCGTGAACACGCCGACAGAACAGGCCGTGAGACGGTCATGtt                     <  1:3363924/71‑1 (MQ=255)
                                                   cGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGAc                                               <  1:1007913/43‑1 (MQ=255)
                                                   cGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCa                   <  1:1119772/71‑1 (MQ=255)
                                                   cGCAAACCCGCAAATGTTCGCCCTCGCCCGTGAACACGCCGACAGAACAGGCCGTGAGACGGTCATGtt                     >  1:897252/1‑69 (MQ=255)
                                                    gCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGa                                      <  1:3576171/51‑1 (MQ=255)
                                                    gCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAg                  <  1:1033888/71‑1 (MQ=255)
                                                    gCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAg                  <  1:1626097/71‑1 (MQ=255)
                                                    gCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAg                  >  1:3007713/1‑71 (MQ=255)
                                                      tAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGct                >  1:2686959/1‑71 (MQ=255)
                                                      aaaCCCGCAAATGTTCGCCCTCGCCCGTGAACACGCCGACAGAACAGGCCGTGAGACGGTCATGt                      >  1:2405965/2‑65 (MQ=255)
                                                       aaCCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGctc               >  1:1816886/1‑71 (MQ=255)
                                                        aCCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGctc               <  1:2101159/70‑1 (MQ=255)
                                                         cccGCAAATGTTTGCCCTTGCCCGTGAACATTCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCtt             <  1:3527938/71‑1 (MQ=255)
                                                         cccGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCtt             <  1:2093670/71‑1 (MQ=255)
                                                          ccGAAAATGTTTGCCCTTGCCCGTGAACATGCCGACaaa                                            <  1:3598227/39‑1 (MQ=255)
                                                            gCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAATACCGGAc                                     <  1:2289944/44‑1 (MQ=255)
                                                            gCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAACCCGGACGTGAGGCGGTGATGTTCAGCTCttttt          >  1:2842291/1‑71 (MQ=255)
                                                            gCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAg                  >  1:1523258/1‑63 (MQ=255)
                                                             cAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGtt                     >  1:3234989/1‑59 (MQ=255)
                                                             cAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGt                      >  1:1455035/1‑58 (MQ=255)
                                                             cAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCtttttt         <  1:212846/71‑1 (MQ=255)
                                                              aaaTGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTa        >  1:2160573/1‑71 (MQ=255)
                                                              aaaTGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTa        >  1:2685030/1‑71 (MQ=255)
                                                               aaTGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTAc       <  1:2510478/71‑1 (MQ=255)
                                                                   tttGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTACgcgc   >  1:25089/1‑71 (MQ=255)
                                                                     tgcccttgcccGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTACGCGCt  <  1:2173089/70‑1 (MQ=255)
                                                                     tgcccttgcccGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTACGCGCt  <  1:2013326/70‑1 (MQ=255)
                                                                     |                                                                     
CTACTTTGTTTTGCTATTCGTTGGCGTCTTTCTTAGCAGCTTTGGCTCGACCGCTAACCCGCAAATGTTTGCCCTTGCCCGTGAACATGCCGACAAAACCGGACGTGAGGCGGTGATGTTCAGCTCTTTTTTACGCGCT  >  minE/1395606‑1395744

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: