Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 132,189 C→T 8.7% T161I (ACA→ATA) ‡ yadD → predicted transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE132,1890CT 8.7% 306.7 / 13.7 103T161I (ACA→ATA) ‡yadDpredicted transposase
Reads supporting (aligned to +/- strand):  ref base C (62/32);  new base T (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 = 8.30e-01

Only 100 of 109 total aligned reads displayed.

TTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCC  >  minE/132122‑132257
                                                                   |                                                                    
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: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)
               cTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCAt                                                    <  1:2266054/71‑1 (MQ=255)
                tGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGAt                                                               <  1:2069502/59‑1 (MQ=255)
                tGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGAt                                                               <  1:1568648/59‑1 (MQ=255)
                tGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGAt                                                               <  1:2905069/59‑1 (MQ=255)
                 ggCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACg                                                           >  1:1192265/1‑62 (MQ=255)
                 ggCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATg                                                   >  1:750181/1‑70 (MQ=255)
                 ggCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATg                                                   >  1:643475/1‑70 (MQ=255)
                  gcgcGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATg                                                   <  1:3744934/69‑1 (MQ=255)
                  gcgcGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                  >  1:922081/1‑70 (MQ=255)
                  gcgcGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                  <  1:1923041/70‑1 (MQ=255)
                  gcgcGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                  <  1:1973193/70‑1 (MQ=255)
                   cgcgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCacac                                                                    >  1:2933717/1‑51 (MQ=255)
                   cgcgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATc                                                      >  1:3571361/1‑65 (MQ=255)
                    gcgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCaa                                                >  1:1598612/1‑70 (MQ=255)
                    gcgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCaa                                                >  1:2757408/1‑70 (MQ=255)
                    gcgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAAc                                               <  1:746966/71‑1 (MQ=255)
                    gcgACGCGTCTATAACAGTCCTTTCCCGATGGTGGATATCACCATCACACCGGATGACGAAATCATGCaa                                                >  1:157050/1‑70 (MQ=255)
                     cgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACa                                              >  1:729607/1‑71 (MQ=255)
                     cgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACa                                              <  1:2993070/71‑1 (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)
                             cTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGa                                      >  1:2375227/1‑71 (MQ=255)
                              tataACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCg                                                                  >  1:2887487/1‑42 (MQ=255)
                              tataACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCg                                                                  >  1:1228138/1‑42 (MQ=255)
                                taACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                  >  1:3311360/1‑56 (MQ=255)
                                taACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                  >  1:1131343/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)
                                            cccGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCg                          <  1:1911822/68‑1 (MQ=255)
                                              cGCTGGTGGATATCACCATCACACCGGATGACGAAATCATg                                                   <  1:1888484/41‑1 (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:597024/1‑51 (MQ=255)
                                                       atatCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGa                                >  1:788052/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)
                                                          tCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACa          >  1:1589625/1‑70 (MQ=255)
                                                           cACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCg                          <  1:820672/53‑1 (MQ=255)
                                                            aCCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCatat       <  1:3378901/71‑1 (MQ=255)
                                                            aCCATCACACCGGATGACGAAATCATGCAACATcggcg                                        >  1:3106062/1‑38 (MQ=255)
                                                            aCCATCACACCGGATGACGAAATCATGCAACATcggcg                                        >  1:2493090/1‑38 (MQ=255)
                                                            aCCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACatat       <  1:3916233/71‑1 (MQ=255)
                                                            aCCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACatat       <  1:1866093/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:1094608/1‑71 (MQ=255)
                                                              cATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTc     >  1:3459369/1‑71 (MQ=255)
                                                              aaTCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTc     >  1:2295944/2‑71 (MQ=255)
                                                               aTCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                   >  1:3632017/1‑40 (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)
                                                                 cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACata        >  1:1445872/1‑65 (MQ=255)
                                                                  acacCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                        >  1:1816522/1‑48 (MQ=255)
                                                                   cacCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGcc  <  1:1085283/69‑1 (MQ=255)
                                                                   |                                                                    
TTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCC  >  minE/132122‑132257

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: