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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE132,1930GT 8.3% 266.7 / 6.9 84P162P (CCG→CCTyadDpredicted transposase
Reads supporting (aligned to +/- strand):  ref base G (37/40);  new base T (5/2);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

ACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGC  >  minE/132126‑132262
                                                                   |                                                                     
aCTCGCCGGAGCTGGCGCGACTCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGa                                                                     <  1:2029230/70‑1 (MQ=255)
 cTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGAt                                                                    <  1:4283319/70‑1 (MQ=255)
  gcGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGAt                                                                    <  1:2152463/68‑1 (MQ=255)
      cGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGaa                                                              <  1:3473239/71‑1 (MQ=255)
        gAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGAt                                                                    <  1:1125649/63‑1 (MQ=255)
                 cgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAAc                                                    <  1:3049631/70‑1 (MQ=255)
                     gcgTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACg                                                                <  1:5691943/54‑1 (MQ=255)
                     gcgTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACg                                                                <  1:5569931/54‑1 (MQ=255)
                     gcgTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACg                                                                <  1:1316281/54‑1 (MQ=255)
                       gTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATcggcg                                             <  1:3638659/71‑1 (MQ=255)
                        tCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCTTGCAACATcggcgg                                            <  1:3914110/71‑1 (MQ=255)
                         cTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGaaa                                                             <  1:1689273/53‑1 (MQ=255)
                          tataACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGAt                                                                    >  1:1353671/1‑45 (MQ=255)
                          tataACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGAt                                                                    >  1:2145813/1‑45 (MQ=255)
                          tataACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGAt                                                                    >  1:2034776/1‑45 (MQ=255)
                          tataACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGAt                                                                    >  1:4939152/1‑45 (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:488148/1‑61 (MQ=255)
                              aCAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATcg                                                >  1:541766/1‑61 (MQ=255)
                              aCAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGc                                       <  1:3522708/70‑1 (MQ=255)
                              aCAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGc                                       <  1:2161053/70‑1 (MQ=255)
                               cAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGa                                     <  1:1986915/71‑1 (MQ=255)
                                 gTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGa                                                               >  1:1980737/1‑43 (MQ=255)
                                 gTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGa                                                               >  1:5162473/1‑43 (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:3186170/71‑1 (MQ=255)
                                                  gATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaa                  <  1:2104096/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:1589802/1‑39 (MQ=255)
                                                            tCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                        >  1:2371299/1‑39 (MQ=255)
                                                            tCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                        >  1:2733784/1‑39 (MQ=255)
                                                            tCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                        >  1:3622995/1‑39 (MQ=255)
                                                            tCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                        >  1:3342140/1‑39 (MQ=255)
                                                            tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                          <  1:1364574/53‑1 (MQ=255)
                                                            tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc        <  1:1990572/71‑1 (MQ=255)
                                                            tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc        >  1:1724271/1‑71 (MQ=255)
                                                            tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc        <  1:3733406/71‑1 (MQ=255)
                                                            tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc        <  1:2811198/71‑1 (MQ=255)
                                                            tCACACCGGATGACGAAATCATGCAACATCGGAGGATTGCGATTCTCg                               <  1:3171776/48‑1 (MQ=255)
                                                             cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                             >  1:1469273/1‑49 (MQ=255)
                                                             cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                             >  1:1119891/1‑49 (MQ=255)
                                                             cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                             >  1:2427586/1‑49 (MQ=255)
                                                             cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                             >  1:283599/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:1982600/71‑1 (MQ=255)
                                                                  cGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgcgc  >  1:5603944/1‑71 (MQ=255)
                                                                  cGGATGACGAAATCATGCAACATAGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgcg   >  1:4314735/1‑70 (MQ=255)
                                                                   ggATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACa               <  1:5499548/57‑1 (MQ=255)
                                                                   ggATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACa               <  1:4411448/57‑1 (MQ=255)
                                                                   |                                                                     
ACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGC  >  minE/132126‑132262

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: