Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 132,232 C→G 22.2% L175L (CTC→CTG yadD → predicted transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE132,2320CG22.2% 191.5 / 47.7 98L175L (CTC→CTGyadDpredicted transposase
Reads supporting (aligned to +/- strand):  ref base C (32/44);  new base G (14/7);  total (46/52)
Fisher's exact test for biased strand distribution p-value = 5.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

Only 100 of 105 total aligned reads displayed.

CTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGA  >  minE/132162‑132299
                                                                      |                                                                   
 tttCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtctc                                                                     >  1:648428/1‑70 (MQ=255)
   tCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                                                                  >  1:5084483/1‑71 (MQ=255)
    cccGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCg                                                                    <  1:4883306/68‑1 (MQ=255)
     ccGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAac                                                                 <  1:4287679/70‑1 (MQ=255)
      cGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAac                                                                 <  1:3269773/69‑1 (MQ=255)
        cTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                                               >  1:5320433/1‑69 (MQ=255)
        cTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                                               >  1:1291540/1‑69 (MQ=255)
           gtggATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCa                                                          <  1:5028866/71‑1 (MQ=255)
              gATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaa                                                       <  1:2104096/71‑1 (MQ=255)
              gATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaa                                                       <  1:3186170/71‑1 (MQ=255)
              gATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaa                                                        <  1:5527922/70‑1 (MQ=255)
               atatCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaa                                                       <  1:2659004/70‑1 (MQ=255)
                  taaCCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                                               >  1:3366443/3‑59 (MQ=255)
                  tCACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCa                                                    >  1:4331714/1‑70 (MQ=255)
                  tCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                                               >  1:5194850/1‑59 (MQ=255)
                  tCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                                               >  1:3654363/1‑59 (MQ=255)
                  tCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                                               >  1:1715133/1‑59 (MQ=255)
                   cACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCata                                                  >  1:4456591/1‑71 (MQ=255)
                   cACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCat                                                   >  1:3151006/1‑70 (MQ=255)
                   cACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACata                                                  >  1:137935/1‑71 (MQ=255)
                    aCCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCata                                                  <  1:1682564/70‑1 (MQ=255)
                    aCCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACatat                                                 <  1:5411175/71‑1 (MQ=255)
                       aTCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCg                                              <  1:4109234/71‑1 (MQ=255)
                        tCACAGCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCg                                                                    <  1:793773/48‑1 (MQ=255)
                        tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                                               <  1:1364574/53‑1 (MQ=255)
                        tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                                             <  1:3733406/71‑1 (MQ=255)
                        tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                                             <  1:1990572/71‑1 (MQ=255)
                        tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                                             >  1:1724271/1‑71 (MQ=255)
                        tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                                             <  1:2811198/71‑1 (MQ=255)
                        tCACACCGGATGACGAAATCATGCAACATCGGAGGATTGCGATTCTCg                                                                    <  1:3171776/48‑1 (MQ=255)
                         cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                                                                  >  1:283599/1‑49 (MQ=255)
                         cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                                                                  >  1:1119891/1‑49 (MQ=255)
                         cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                                                                  >  1:2427586/1‑49 (MQ=255)
                         cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                                                                  >  1:1469273/1‑49 (MQ=255)
                         cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGa                                                                   >  1:824333/1‑48 (MQ=255)
                          acacCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGcc                                            <  1:5225852/70‑1 (MQ=255)
                          acacCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGcc                                            <  1:1456262/70‑1 (MQ=255)
                            acCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaa                                                       <  1:1591358/57‑1 (MQ=255)
                            acCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaa                                                       <  1:4922131/57‑1 (MQ=255)
                              cGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgcgc                                       >  1:5603944/1‑71 (MQ=255)
                              cGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgcgc                                       <  1:1982600/71‑1 (MQ=255)
                              cGGATGACGAAATCATGCAACATAGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgcg                                        >  1:4314735/1‑70 (MQ=255)
                               ggATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACa                                                    <  1:4411448/57‑1 (MQ=255)
                               ggATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACa                                                    <  1:5499548/57‑1 (MQ=255)
                                 aTGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaat                                                       <  1:1765041/52‑2 (MQ=255)
                                  tGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAac                                                                 <  1:5469016/41‑1 (MQ=255)
                                  tGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAac                                                                 <  1:3950972/41‑1 (MQ=255)
                                  tGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAac                                                                 <  1:1419189/41‑1 (MQ=255)
                                  tGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAac                                                                 <  1:2405323/41‑1 (MQ=255)
                                  tGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                                             <  1:5480886/61‑1 (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)
                                          tCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaaa                                                      <  1:5304759/44‑1 (MQ=255)
                                          tCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaaa                                                      <  1:2386059/44‑1 (MQ=255)
                                          tCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaaa                                                      <  1:3964826/44‑1 (MQ=255)
                                           cATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCata                                                  >  1:1096838/1‑47 (MQ=255)
                                           cATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCat                                                   >  1:3949494/1‑46 (MQ=255)
                                           cATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCat                                                   >  1:3643736/1‑46 (MQ=255)
                                           cATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCat                                                   >  1:3007904/1‑46 (MQ=255)
                                           cATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGat                                    >  1:2371657/1‑60 (MQ=255)
                                           cATGCAACATAGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGat                                    >  1:3248330/1‑60 (MQ=255)
                                            aTGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACa                                                    >  1:176267/1‑44 (MQ=255)
                                             tGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaa                                                        <  1:202621/39‑1 (MQ=255)
                                              gCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGcga                                           >  1:2947347/1‑49 (MQ=255)
                                               cAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgcg                                        >  1:1678405/1‑53 (MQ=255)
                                               cAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCt                      >  1:1755295/1‑71 (MQ=255)
                                               cAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCt                      <  1:4014833/71‑1 (MQ=255)
                                               cAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCt                      >  1:121453/1‑71 (MQ=255)
                                                aaCATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTa                                 >  1:2398642/1‑59 (MQ=255)
                                                aaCATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTa                                 >  1:2550432/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)
                                                     cggcggATTGCGATTCTCGAACTACTGCAAAAACatat                                                 >  1:2487730/1‑38 (MQ=255)
                                                     cggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAAt                               >  1:4727016/1‑56 (MQ=255)
                                                     cggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAAt                               >  1:4733468/1‑56 (MQ=255)
                                                     cggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAAt                               >  1:5388259/1‑56 (MQ=255)
                                                      ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa               >  1:5025523/1‑71 (MQ=255)
                                                      ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa               >  1:549698/1‑71 (MQ=255)
                                                      ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa               <  1:1932043/71‑1 (MQ=255)
                                                      ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa               <  1:1726524/71‑1 (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:3476834/1‑57 (MQ=255)
                                                                    ctcGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa               >  1:4439828/1‑57 (MQ=255)
                                                                      |                                                                   
CTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGA  >  minE/132162‑132299

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: