Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 132,188 A→G 8.7% T161A (ACA→GCA) ‡ yadD → predicted transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE132,1880AG 8.7% 279.8 / 14.4 103T161A (ACA→GCA) ‡yadDpredicted transposase
Reads supporting (aligned to +/- strand):  ref base A (62/32);  new base G (7/2);  total (69/34)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 108 total aligned reads displayed.

TTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGC  >  minE/132122‑132256
                                                                  |                                                                    
ttttACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACAcc                                                                  >  1:328652/1‑71 (MQ=255)
 tttACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCg                                                                 >  1:3697165/1‑71 (MQ=255)
  ttACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCg                                                                 >  1:2328532/1‑70 (MQ=255)
  ttACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCg                                                                 >  1:1700667/1‑70 (MQ=255)
  ttACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCg                                                                 >  1:1483884/1‑70 (MQ=255)
   tACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCgg                                                                <  1:2188404/70‑1 (MQ=255)
    aCTCGCCTGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGa                                                               <  1:3727681/70‑1 (MQ=255)
      tCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGa                                                            >  1:3442030/1‑71 (MQ=255)
       cGCCGGAGCTGGCGCGACTCGTCTATAACAGTCCTTTCCCGCTGGTGGATATAACCATCACACCGGATg                                                             >  1:3597756/1‑69 (MQ=255)
        gCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGAc                                                           >  1:1761879/1‑70 (MQ=255)
        gCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGAc                                                           >  1:3197789/1‑70 (MQ=255)
        gCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGAc                                                           >  1:142356/1‑70 (MQ=255)
        gCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGAc                                                           >  1:3504281/1‑70 (MQ=255)
         ccGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGa                                                         <  1:518685/71‑1 (MQ=255)
         ccGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGa                                                         <  1:1618431/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)
                tGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCTGAt                                                              <  1:3322144/59‑1 (MQ=255)
                tGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGAt                                                              <  1:1568648/59‑1 (MQ=255)
                tGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGAt                                                              <  1:2905069/59‑1 (MQ=255)
                tGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGAt                                                              <  1:2069502/59‑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:922081/1‑70 (MQ=255)
                  gcgcGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                 <  1:1973193/70‑1 (MQ=255)
                  gcgcGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                 <  1:1923041/70‑1 (MQ=255)
                   cgcgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCacac                                                                   >  1:2933717/1‑51 (MQ=255)
                    gcgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCaa                                               >  1:1598612/1‑70 (MQ=255)
                    gcgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCaa                                               >  1:2757408/1‑70 (MQ=255)
                    gcgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAAc                                              <  1:746966/71‑1 (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)
                        cgcgTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGa                                                               >  1:3558635/1‑50 (MQ=255)
                         gcgTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGa                                                            >  1:2926153/1‑52 (MQ=255)
                          cgTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATctgc                                        <  1:1326904/71‑1 (MQ=255)
                            tCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCAt                                                   <  1:2143457/58‑1 (MQ=255)
                             cTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGa                                     >  1:2375227/1‑71 (MQ=255)
                              tataACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCg                                                                 >  1:1228138/1‑42 (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:2289256/71‑1 (MQ=255)
                                        cTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                  >  1:2128209/1‑63 (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:312576/1‑71 (MQ=255)
                                                         atCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaa             <  1:232113/67‑1 (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:3106062/1‑38 (MQ=255)
                                                            aCCATCACACCGGATGACGAAATCATGCAACATcggcg                                       >  1:2493090/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:369444/1‑70 (MQ=255)
                                                             ccATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCatat      <  1:895428/70‑1 (MQ=255)
                                                             ccATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATAtt     >  1:3018776/1‑71 (MQ=255)
                                                              cATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTc    >  1:3459369/1‑71 (MQ=255)
                                                              cATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTc    >  1:1094608/1‑71 (MQ=255)
                                                               aTCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                  >  1:3328833/1‑40 (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)
                                                                 cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACata       >  1:1445872/1‑65 (MQ=255)
                                                                  acacCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                       >  1:1816522/1‑48 (MQ=255)
                                                                  acacCGGATGACGAAATCATGCAACATCGGCGGATAGCGATTCTCg                         >  1:538351/1‑46 (MQ=255)
                                                                  |                                                                    
TTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGC  >  minE/132122‑132256

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: