Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 132,202 A→G 8.8% E165E (GAA→GAG yadD → predicted transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE132,2020AG 8.8% 216.3 / 8.6 79E165E (GAA→GAGyadDpredicted transposase
Reads supporting (aligned to +/- strand):  ref base A (32/40);  new base G (5/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGT  >  minE/132132‑132272
                                                                      |                                                                      
cGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGaa                                                                        <  1:3473239/71‑1 (MQ=255)
           cgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAAc                                                              <  1:3049631/70‑1 (MQ=255)
                 gTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATcggcg                                                       <  1:3638659/71‑1 (MQ=255)
                  tCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCTTGCAACATcggcgg                                                      <  1:3914110/71‑1 (MQ=255)
                   cTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGaaa                                                                       <  1:1689273/53‑1 (MQ=255)
                    tataACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGAt                                                    >  1:4600577/1‑71 (MQ=255)
                     ataACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGAtt                                                   >  1:4143045/1‑71 (MQ=255)
                        aCAGTCCTTTCCCGCTGGTGGATATCACCATCACCCCGGATGACGAAATCATGCAACATcg                                                          >  1:2266476/1‑61 (MQ=255)
                        aCAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCaa                                                               >  1:3525250/1‑56 (MQ=255)
                        aCAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATcg                                                          >  1:541766/1‑61 (MQ=255)
                        aCAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATcg                                                          >  1:488148/1‑61 (MQ=255)
                        aCAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGc                                                 <  1:2161053/70‑1 (MQ=255)
                        aCAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGc                                                 <  1:3522708/70‑1 (MQ=255)
                         cAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGa                                               <  1:1986915/71‑1 (MQ=255)
                              cTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtctc                                          <  1:1800412/71‑1 (MQ=255)
                               tttCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtctc                                          >  1:648428/1‑70 (MQ=255)
                                ttCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCTACATCGGCGGATTGCGATTCTCg                                         >  1:2903375/1‑70 (MQ=255)
                                 tCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                                       >  1:5084483/1‑71 (MQ=255)
                                  cccGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCg                                         <  1:4883306/68‑1 (MQ=255)
                                   ccGCTGGTGGATATCACCATCACACCGGATGACGAAAt                                                                      <  1:4520174/38‑1 (MQ=255)
                                   ccGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAac                                      <  1:4287679/70‑1 (MQ=255)
                                    cGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAac                                      <  1:3269773/69‑1 (MQ=255)
                                      cTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAAc                                                              >  1:3978847/1‑43 (MQ=255)
                                      cTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                    >  1:5320433/1‑69 (MQ=255)
                                      cTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                    >  1:1291540/1‑69 (MQ=255)
                                         gtggATATCACCATCACACCGGATGACGAAATCATGCAACATcggcg                                                       <  1:152848/47‑1 (MQ=255)
                                         gtggATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCa                               <  1:5028866/71‑1 (MQ=255)
                                           ggATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaa                             <  1:879856/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)
                                             atatCACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTAtag                                                  <  1:3532925/48‑3 (MQ=255)
                                             atatCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaa                            <  1:2659004/70‑1 (MQ=255)
                                                taaCCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                    >  1:3366443/3‑59 (MQ=255)
                                                tCACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGct                                                >  1:1639548/1‑46 (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)
                                                tCACAATCACACCGGATGACGAAATCATGCAACATc                                                           >  1:1642589/1‑36 (MQ=255)
                                                 cACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCata                       >  1:4456591/1‑71 (MQ=255)
                                                 cACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCat                        >  1:5139991/1‑70 (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)
                                                      tCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                                  >  1:3622995/1‑39 (MQ=255)
                                                      tCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                                  >  1:3342140/1‑39 (MQ=255)
                                                      tCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                                  >  1:1589802/1‑39 (MQ=255)
                                                      tCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                                  >  1:2371299/1‑39 (MQ=255)
                                                      tCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                                  >  1:2733784/1‑39 (MQ=255)
                                                      tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                    <  1:1364574/53‑1 (MQ=255)
                                                      tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                  <  1:2811198/71‑1 (MQ=255)
                                                      tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                  <  1:3733406/71‑1 (MQ=255)
                                                      tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                  >  1:1724271/1‑71 (MQ=255)
                                                      tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                  <  1:1990572/71‑1 (MQ=255)
                                                      tCACACCGGATGACGAAATCATGCAACATCGGAGGATTGCGATTCTCg                                         <  1:3171776/48‑1 (MQ=255)
                                                       cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                                       >  1:1119891/1‑49 (MQ=255)
                                                       cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                                       >  1:2427586/1‑49 (MQ=255)
                                                       cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                                       >  1:283599/1‑49 (MQ=255)
                                                       cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                                       >  1:1469273/1‑49 (MQ=255)
                                                       cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGa                                        >  1:824333/1‑48 (MQ=255)
                                                        acacCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGcc                 <  1:1456262/70‑1 (MQ=255)
                                                        acacCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGcc                 <  1:5225852/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:5499548/57‑1 (MQ=255)
                                                             ggATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACa                         <  1:4411448/57‑1 (MQ=255)
                                                               aTGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaat                            <  1:1765041/52‑2 (MQ=255)
                                                                tGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAac                                      <  1:1419189/41‑1 (MQ=255)
                                                                tGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAac                                      <  1:3950972/41‑1 (MQ=255)
                                                                tGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAac                                      <  1:5469016/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)
                                                                      aaTCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGt  <  1:4346016/71‑1 (MQ=255)
                                                                      |                                                                      
CGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGT  >  minE/132132‑132272

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: