Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 132,193 G→T 10.4% P162P (CCG→CCT yadD → predicted transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE132,1930GT10.4% 330.6 / 14.7 105P162P (CCG→CCTyadDpredicted transposase
Reads supporting (aligned to +/- strand):  ref base G (62/32);  new base T (7/4);  total (69/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

Only 100 of 113 total aligned reads displayed.

TTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCG  >  minE/132123‑132263
                                                                      |                                                                      
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: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)
               tGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCTGAt                                                                     <  1:3322144/59‑1 (MQ=255)
               tGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGAt                                                                     <  1:2905069/59‑1 (MQ=255)
               tGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGAt                                                                     <  1:2069502/59‑1 (MQ=255)
               tGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGAt                                                                     <  1:1568648/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:1923041/70‑1 (MQ=255)
                 gcgcGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                        <  1:1973193/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)
                   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)
                       cgcgTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGa                                                                      >  1:3558635/1‑50 (MQ=255)
                           tCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCAt                                                          <  1:2143457/58‑1 (MQ=255)
                            cTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGa                                            >  1:2375227/1‑71 (MQ=255)
                             tataACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCg                                                                        >  1:2887487/1‑42 (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:1865246/71‑1 (MQ=255)
                                     tCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtc                                   <  1:2289256/71‑1 (MQ=255)
                                       cTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                         >  1:2128209/1‑63 (MQ=255)
                                       cTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtctc                                 >  1:1154686/1‑71 (MQ=255)
                                       cTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtctc                                 <  1:1654283/71‑1 (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:3699296/1‑71 (MQ=255)
                                                tggtggATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTg                        >  1:1378704/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)
                                                        atCACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCa                >  1:524861/1‑71 (MQ=255)
                                                        atCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaa                    <  1:232113/67‑1 (MQ=255)
                                                         tCACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCa                >  1:3279009/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)
                                                            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:1094608/1‑71 (MQ=255)
                                                             aaTCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTc           >  1:2295944/2‑71 (MQ=255)
                                                              aTCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                         >  1:3328833/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)
                                                                      ggATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGcgcg  <  1:1739920/71‑1 (MQ=255)
                                                                      |                                                                      
TTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCG  >  minE/132123‑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: