Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 132,259 G→C 13.6% Q184H (CAG→CAC) ‡ yadD → predicted transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE132,2590GC13.6% 220.1 / 22.5 88Q184H (CAG→CAC) ‡yadDpredicted transposase
Reads supporting (aligned to +/- strand):  ref base G (45/31);  new base C (5/7);  total (50/38)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 101 total aligned reads displayed.

CGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCGGAAGTCAGT  >  minE/132192‑132329
                                                                   |                                                                      
cGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgcgc                                                                     >  1:5603944/1‑71 (MQ=255)
cGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgcgc                                                                     <  1:1982600/71‑1 (MQ=255)
cGGATGACGAAATCATGCAACATAGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgcg                                                                      >  1:4314735/1‑70 (MQ=255)
    tGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACt                                                                 <  1:5105414/71‑1 (MQ=255)
      aCGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACtt                                                                >  1:4311932/1‑70 (MQ=255)
         agaTCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAAtc                                                            >  1:5264174/3‑70 (MQ=255)
          aaTCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGt                                                           <  1:4346016/71‑1 (MQ=255)
           aTCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACcgcg                                                                    <  1:1981471/61‑1 (MQ=255)
           aTCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAAtcgg                                                          <  1:3144933/71‑4 (MQ=255)
            tCATGCAACATCGGCGTATAGCTCTGCTTGAACTGATTCAAAAGCATATTCGCGACCGCGat                                                                  <  1:437409/62‑2 (MQ=255)
            tCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGat                                                                  <  1:3389478/62‑2 (MQ=255)
            tCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGat                                                                  <  1:5381492/62‑2 (MQ=255)
            tCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGAttt                                                                <  1:4614380/64‑1 (MQ=255)
            tCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAAtcgg                                                          <  1:4086207/70‑4 (MQ=255)
             cATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGat                                                                  >  1:2371657/1‑60 (MQ=255)
             cATGCAACATAGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGat                                                                  >  1:3248330/1‑60 (MQ=255)
                 cAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgcg                                                                      >  1:1678405/1‑53 (MQ=255)
                 cAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCt                                                    >  1:121453/1‑71 (MQ=255)
                 cAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCt                                                    >  1:1755295/1‑71 (MQ=255)
                 cAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCt                                                    <  1:4014833/71‑1 (MQ=255)
                  aaCATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTa                                                               >  1:2550432/1‑59 (MQ=255)
                  aaCATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTa                                                               >  1:2398642/1‑59 (MQ=255)
                   aCATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgcgc                                                                     <  1:825848/52‑1 (MQ=255)
                     aTCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTaa                                                              >  1:5389881/1‑57 (MQ=255)
                     aTCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTaa                                                              >  1:4756399/1‑57 (MQ=255)
                     aTCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAg                                                >  1:3152257/1‑71 (MQ=255)
                       cggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAAt                                                             >  1:4733468/1‑56 (MQ=255)
                       cggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAAt                                                             >  1:4727016/1‑56 (MQ=255)
                       cggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAAt                                                             >  1:5388259/1‑56 (MQ=255)
                        ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa                                             >  1:5025523/1‑71 (MQ=255)
                        ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa                                             <  1:1726524/71‑1 (MQ=255)
                        ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa                                             <  1:1932043/71‑1 (MQ=255)
                        ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa                                             >  1:549698/1‑71 (MQ=255)
                           ggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACt                                           <  1:1854585/70‑1 (MQ=255)
                                  gATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCtt                                                   <  1:1957148/55‑1 (MQ=255)
                                  gATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCtt                                                   <  1:1875035/55‑1 (MQ=255)
                                    ttCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCt                                  >  1:5450095/1‑70 (MQ=255)
                                    ttCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTg                                 <  1:2192961/71‑1 (MQ=255)
                                     tctcGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGa                                >  1:3072992/1‑71 (MQ=255)
                                     tctcGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGa                                >  1:1205416/1‑71 (MQ=255)
                                      ctcGAACTACTGCAAAAACATATTCGCCAGCGCGACTTaa                                                              >  1:1582883/1‑40 (MQ=255)
                                      ctcGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa                                             >  1:4439828/1‑57 (MQ=255)
                                      ctcGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa                                             >  1:3476834/1‑57 (MQ=255)
                                         gAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTAt                                                        >  1:1217928/1‑43 (MQ=255)
                                         gAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGAt                               <  1:498423/68‑1 (MQ=255)
                                         gAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATc                              <  1:2636217/69‑1 (MQ=255)
                                           actactGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATcgcc                           >  1:4945734/1‑68 (MQ=255)
                                           actactGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATcgac                           >  1:4301254/1‑70 (MQ=255)
                                           actactGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATcgac                           <  1:27477/70‑1 (MQ=255)
                                              actGCAAAAACATATTCGCCAGCGCGACTTTATGt                                                           <  1:2813211/35‑1 (MQ=255)
                                              actGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAg                       >  1:4110214/1‑71 (MQ=255)
                                               ctGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAAc                                            >  1:820289/1‑49 (MQ=255)
                                               ctGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAAc                                            >  1:2780334/1‑49 (MQ=255)
                                               ctGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATc                              <  1:1448313/63‑1 (MQ=255)
                                               ctGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGaa                        >  1:723550/1‑69 (MQ=255)
                                               ctGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAgg                      >  1:3767974/1‑71 (MQ=255)
                                                tGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGc                                               <  1:5630840/45‑1 (MQ=255)
                                                tGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGt                                        <  1:2358735/52‑1 (MQ=255)
                                                 gCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATcgacga                         >  1:4396399/1‑66 (MQ=255)
                                                 gCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATcgacga                         >  1:5393177/1‑66 (MQ=255)
                                                  cAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGAt                               >  1:671609/1‑59 (MQ=255)
                                                  cAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGAt                               >  1:3560665/1‑59 (MQ=255)
                                                  cAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGAt                               >  1:748620/1‑59 (MQ=255)
                                                  cAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGAt                               >  1:3749686/1‑59 (MQ=255)
                                                  cAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATcga                            >  1:692509/1‑62 (MQ=255)
                                                  cAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAggg                     <  1:4352719/69‑1 (MQ=255)
                                                   aaaaaCATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAg                       >  1:4778198/1‑66 (MQ=255)
                                                   aaaaaCATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTac                  <  1:3883655/71‑1 (MQ=255)
                                                   aaaaaCATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTa                   <  1:3088165/70‑1 (MQ=255)
                                                   aaaaaCATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTa                   <  1:3334884/70‑1 (MQ=255)
                                                   aaaaaCATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCAACGAAg                       >  1:3339152/1‑66 (MQ=255)
                                                     aaaCATATTCGCCAGCGCGACTTAATGTTATTGCtt                                                   <  1:1658124/36‑1 (MQ=255)
                                                     aaaCATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATcg                             <  1:279077/58‑1 (MQ=255)
                                                     aaaCATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTacac                >  1:4318894/1‑71 (MQ=255)
                                                       aCATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTacac                >  1:542714/1‑69 (MQ=255)
                                                       aCATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTacac                >  1:2959748/1‑69 (MQ=255)
                                                       aCATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTacac                >  1:1876175/1‑69 (MQ=255)
                                                       aCATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTa              <  1:3289267/71‑1 (MQ=255)
                                                          tatTCGCCAGCGCGACTTAATGTTATTGCTTGAGc                                               <  1:1829889/35‑1 (MQ=255)
                                                          tatTCGCCAGCGCGACTTAATGTTATTGCTTGAGc                                               <  1:1573545/35‑1 (MQ=255)
                                                          tatTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCg           <  1:1865938/71‑1 (MQ=255)
                                                            ttCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGc            >  1:4745196/1‑68 (MQ=255)
                                                            ttCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGc            >  1:3016066/1‑68 (MQ=255)
                                                            ttCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCgg          <  1:2048286/70‑1 (MQ=255)
                                                             tCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCgg          >  1:3781331/1‑69 (MQ=255)
                                                              cGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCGGAAg       >  1:1088716/1‑71 (MQ=255)
                                                              cGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCGGAAg       >  1:631586/1‑71 (MQ=255)
                                                               gCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAg                       >  1:661355/1‑54 (MQ=255)
                                                               gCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAg                       >  1:4285262/1‑54 (MQ=255)
                                                               gCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCGGAAGt      <  1:2896582/71‑1 (MQ=255)
                                                               gCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCGGAAGt      <  1:4878921/71‑1 (MQ=255)
                                                                 cAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCGGAAGTCa    <  1:5001049/71‑1 (MQ=255)
                                                                   gcgcGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTa                   <  1:2587922/54‑1 (MQ=255)
                                                                   gcgcGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTa                   <  1:2100752/54‑1 (MQ=255)
                                                                   gcgcGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCgg          <  1:3580560/63‑1 (MQ=255)
                                                                   gcgcGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCgg          <  1:3984407/63‑1 (MQ=255)
                                                                   gcgcGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCgg          <  1:3554239/63‑1 (MQ=255)
                                                                   gcgcGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCgg          <  1:1524346/63‑1 (MQ=255)
                                                                   gcgcGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCGGAAGTCagt  <  1:4992125/71‑1 (MQ=255)
                                                                   gcgcGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCGGAAGTCagt  <  1:1092622/71‑1 (MQ=255)
                                                                   |                                                                      
CGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGGGTACACTAGCGGAAGTCAGT  >  minE/132192‑132329

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: