Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 132,199 C→T 10.3% D164D (GAC→GAT yadD → predicted transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE132,1990CT10.3% 271.8 / 16.9 97D164D (GAC→GATyadDpredicted transposase
Reads supporting (aligned to +/- strand):  ref base C (54/32);  new base T (7/3);  total (62/35)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 103 total aligned reads displayed.

GCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCG  >  minE/132130‑132263
                                                                     |                                                                
gCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGAc                                                                  >  1:1761879/1‑70 (MQ=255)
gCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGAc                                                                  >  1:3197789/1‑70 (MQ=255)
gCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGAc                                                                  >  1:3504281/1‑70 (MQ=255)
gCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGAc                                                                  >  1:142356/1‑70 (MQ=255)
 ccGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGa                                                                <  1:1618431/71‑1 (MQ=255)
 ccGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGa                                                                <  1:518685/71‑1 (MQ=255)
  cGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGa                                                                >  1:778636/1‑70 (MQ=255)
     aGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATc                                                            <  1:1731507/71‑1 (MQ=255)
      gCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCa                                                           >  1:2408109/1‑71 (MQ=255)
       cTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCAt                                                          <  1:2266054/71‑1 (MQ=255)
         ggCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACg                                                                 >  1:1192265/1‑62 (MQ=255)
         ggCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATg                                                         >  1:643475/1‑70 (MQ=255)
         ggCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATg                                                         >  1:750181/1‑70 (MQ=255)
          gcgcGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATg                                                         <  1:3744934/69‑1 (MQ=255)
          gcgcGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                        <  1:1973193/70‑1 (MQ=255)
          gcgcGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                        <  1:1923041/70‑1 (MQ=255)
          gcgcGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                        >  1:922081/1‑70 (MQ=255)
           cgcgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATc                                                            >  1:3571361/1‑65 (MQ=255)
            gcgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCaa                                                      >  1:2757408/1‑70 (MQ=255)
            gcgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCaa                                                      >  1:1598612/1‑70 (MQ=255)
            gcgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAAc                                                     <  1:746966/71‑1 (MQ=255)
            gcgACGCGTCTATAACAGTCCTTTCCCGATGGTGGATATCACCATCACACCGGATGACGAAATCATGCaa                                                      >  1:157050/1‑70 (MQ=255)
             cgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACa                                                    <  1:2993070/71‑1 (MQ=255)
             cgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACa                                                    >  1:729607/1‑71 (MQ=255)
             cgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGAAGAAATCATGCAACa                                                    >  1:1576723/1‑71 (MQ=255)
               aCGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATc                                                  >  1:3718495/1‑71 (MQ=255)
                  cgTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATctgc                                               <  1:1326904/71‑1 (MQ=255)
                    tCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCAt                                                          <  1:2143457/58‑1 (MQ=255)
                     cTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGa                                            >  1:2375227/1‑71 (MQ=255)
                        taACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                        >  1:1131343/1‑56 (MQ=255)
                        taACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                        >  1:3311360/1‑56 (MQ=255)
                         aaCAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGc                                        >  1:1834634/1‑71 (MQ=255)
                          aCAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGc                                        <  1:3546386/70‑1 (MQ=255)
                           cAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGa                                      >  1:2340633/1‑71 (MQ=255)
                            aGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGAt                                     >  1:2146896/1‑71 (MQ=255)
                              tCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtc                                   <  1:1865246/71‑1 (MQ=255)
                              tCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtc                                   <  1:2289256/71‑1 (MQ=255)
                                cTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtctc                                 <  1:1654283/71‑1 (MQ=255)
                                cTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtctc                                 >  1:1154686/1‑71 (MQ=255)
                                cTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtct                                  >  1:2360845/1‑70 (MQ=255)
                                 tttCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATcggc                                               <  1:2305727/56‑1 (MQ=255)
                                   tCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCg                                >  1:3672365/1‑69 (MQ=255)
                                    cccGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCg                                <  1:1911822/68‑1 (MQ=255)
                                      cGCTGGTGGATATCACCATCACACCGGATGACGAAATCATg                                                         <  1:1888484/41‑1 (MQ=255)
                                      cGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                           >  1:1559944/1‑71 (MQ=255)
                                       gCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAactac                          >  1:3282126/1‑71 (MQ=255)
                                        cTGGTGGATATCACCATCACACCGGATGACGAAATCATGCa                                                       >  1:1649488/1‑41 (MQ=255)
                                        atggtggATATCACCATCACACCGGATGACGAAATCATGCa                                                       >  1:3124053/2‑41 (MQ=255)
                                         tggtggATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                              >  1:2689060/1‑65 (MQ=255)
                                         tggtggATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTg                        >  1:1378704/1‑71 (MQ=255)
                                         tggtggATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTg                        >  1:3699296/1‑71 (MQ=255)
                                              gATATCACCATCGTACCTGACGATGAGATCATGCAACATcggcg                                              >  1:3257111/1‑44 (MQ=255)
                                              gATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaa                    <  1:2961241/70‑1 (MQ=255)
                                               atatCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGa                                      >  1:788052/1‑51 (MQ=255)
                                               atatCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGa                                      >  1:597024/1‑51 (MQ=255)
                                                tatCACCATCACACCGGATGACGAAATCATGCAAc                                                     >  1:113110/1‑35 (MQ=255)
                                                tatCACCATCACACCGGATGACGAAATCATGCAACATcggc                                               >  1:2884706/1‑41 (MQ=255)
                                                 atCACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCa                >  1:524861/1‑71 (MQ=255)
                                                 atCACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCa                >  1:312576/1‑71 (MQ=255)
                                                 atCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaa                    <  1:232113/67‑1 (MQ=255)
                                                  tCACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCa                >  1:3279009/1‑70 (MQ=255)
                                                  tCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACa                >  1:1589625/1‑70 (MQ=255)
                                                   cACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCg                                <  1:820672/53‑1 (MQ=255)
                                                    aCCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCatat             <  1:3378901/71‑1 (MQ=255)
                                                    aCCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTAGAact                            >  1:2996407/1‑56 (MQ=255)
                                                    aCCATCACACCGGATGACGAAATCATGCAACATcggcg                                              >  1:2493090/1‑38 (MQ=255)
                                                    aCCATCACACCGGATGACGAAATCATGCAACATcggcg                                              >  1:3106062/1‑38 (MQ=255)
                                                    aCCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACatat             <  1:1866093/71‑1 (MQ=255)
                                                    aCCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACatat             <  1:3916233/71‑1 (MQ=255)
                                                     ccATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCaaaa                   >  1:501805/1‑64 (MQ=255)
                                                     ccATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCatat             <  1:895428/70‑1 (MQ=255)
                                                     ccATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCatat             >  1:369444/1‑70 (MQ=255)
                                                     ccATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATAtt            >  1:3018776/1‑71 (MQ=255)
                                                      cATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTc           >  1:3459369/1‑71 (MQ=255)
                                                      cATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTc           >  1:1094608/1‑71 (MQ=255)
                                                      aaTCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTc           >  1:2295944/2‑71 (MQ=255)
                                                       aTCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                         >  1:3632017/1‑40 (MQ=255)
                                                       aTCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTc           >  1:1704122/1‑70 (MQ=255)
                                                        tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtct                                  >  1:2133987/1‑46 (MQ=255)
                                                        tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCg                                <  1:3474398/48‑1 (MQ=255)
                                                        tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACata              >  1:413424/1‑66 (MQ=255)
                                                        tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACa                <  1:3129806/64‑1 (MQ=255)
                                                        tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCg          >  1:2795204/1‑70 (MQ=255)
                                                        tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc         <  1:822055/71‑1 (MQ=255)
                                                         cgtacCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGc         <  1:353302/67‑1 (MQ=255)
                                                         cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACata              >  1:1445872/1‑65 (MQ=255)
                                                          acacCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                              >  1:1816522/1‑48 (MQ=255)
                                                          acacCGGATGACGAAATCATGCAACATCGGCGGATAGCGATTCTCg                                >  1:538351/1‑46 (MQ=255)
                                                           cacCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGcc        <  1:1085283/69‑1 (MQ=255)
                                                            acCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAactact                         <  1:2242153/51‑1 (MQ=255)
                                                            acCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgc     >  1:1215807/1‑71 (MQ=255)
                                                            acCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgc     >  1:3892751/1‑71 (MQ=255)
                                                              cGGATGACGAAATCATGCAACATCGGCGGATTGCg                                       >  1:192327/1‑35 (MQ=255)
                                                               ggATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGcgcg  <  1:1739920/71‑1 (MQ=255)
                                                                  tGACGAAATCATGCAACATCGGCGGATTGCGATtct                                  <  1:1927416/36‑1 (MQ=255)
                                                                  tGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgc     >  1:3322513/1‑65 (MQ=255)
                                                                   gACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaa                    <  1:2582517/49‑1 (MQ=255)
                                                                   gACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaa                    <  1:368264/49‑1 (MQ=255)
                                                                     cGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc         <  1:2036055/58‑1 (MQ=255)
                                                                     cGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc         <  1:3826494/58‑1 (MQ=255)
                                                                     |                                                                
GCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCG  >  minE/132130‑132263

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: