Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 131,899 C→T 20.7% L64L (CTC→CTT)  yadD → predicted transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE131,8990CT20.7% 125.4 / 21.3 86L64L (CTC→CTT) yadDpredicted transposase
Reads supporting (aligned to +/- strand):  ref base C (18/50);  new base T (4/14);  total (22/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAA  >  minE/131843‑131950
                                                        |                                                   
ttAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCa                                          <  1:2372601/68‑1 (MQ=255)
     gTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGc                                     <  1:1912286/68‑1 (MQ=255)
     gTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGc                                     <  1:51383/68‑1 (MQ=255)
     gTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGc                                     <  1:2147135/68‑1 (MQ=255)
     gTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGc                                     <  1:2301957/68‑1 (MQ=255)
     gTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGc                                     <  1:3085441/68‑1 (MQ=255)
     gTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGc                                     <  1:3065580/68‑1 (MQ=255)
     gTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGc                                     <  1:2483283/68‑1 (MQ=255)
      tCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAg                                   >  1:1049798/1‑69 (MQ=255)
      tCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAg                                   >  1:221524/1‑69 (MQ=255)
      tCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAg                                   >  1:820606/1‑69 (MQ=255)
      tCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAg                                   >  1:135358/1‑69 (MQ=255)
      tCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAg                                   >  1:1175934/1‑69 (MQ=255)
      tCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGCACGTGCTCTATTCCGTGCAAATGc                                     <  1:628739/67‑1 (MQ=255)
        ggggAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAggg                                 <  1:763496/69‑1 (MQ=255)
             gTTTTATTGAAGAAAGCCTGAAAGGGCACAGCACTGACGTGCTTTATTCCGTGCAAATGCAGGGCaacc                            >  1:3294066/1‑67 (MQ=255)
             gTTTTATTGAAGAAAGCCTGAAAGGGCACAGCACTGACGTGCTTTATTCCGTGCAAATGCAGGGCaacc                            >  1:2289745/1‑67 (MQ=255)
             gTTTTATTGAAGAAAGCCTGAAAGGGCACAGCACTGACGTGCTTTATTCCGTGCAAATGCAGGGCaacc                            >  1:876573/1‑67 (MQ=255)
             gTTTTATTGAAGAAAGCCTGAAAGGGCACAGCACTGACGTGCTTTATTCCGTGCAAATGCAGGGCaacc                            >  1:137475/1‑67 (MQ=255)
                tCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCg                         >  1:483496/1‑69 (MQ=255)
                    tGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                     <  1:1349575/69‑1 (MQ=255)
                    tGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                     <  1:827751/69‑1 (MQ=255)
                    tGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                     <  1:1364304/69‑1 (MQ=255)
                    tGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                     <  1:598433/69‑1 (MQ=255)
                    tGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                     <  1:2585326/69‑1 (MQ=255)
                    tGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                     <  1:2190606/69‑1 (MQ=255)
                    tGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                     <  1:504407/69‑1 (MQ=255)
                    tGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                     <  1:3168435/69‑1 (MQ=255)
                    tGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                     <  1:1642580/69‑1 (MQ=255)
                    tGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                     <  1:3185467/69‑1 (MQ=255)
                    tGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                     <  1:1736041/69‑1 (MQ=255)
                     gAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCATGGCAATCCCGGTTa                     <  1:2269822/68‑1 (MQ=255)
                     gAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                     <  1:797835/68‑1 (MQ=255)
                     gAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                     <  1:674053/68‑1 (MQ=255)
                      aaGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                     <  1:1977645/67‑1 (MQ=255)
                        gagtgCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTg                 <  1:3086438/69‑1 (MQ=255)
                        gagagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTg                 <  1:1317191/69‑1 (MQ=255)
                        gagagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTg                 <  1:3147318/69‑1 (MQ=255)
                        gagagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTg                 <  1:1738486/69‑1 (MQ=255)
                        gagagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTg                 <  1:2960995/69‑1 (MQ=255)
                        gagagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTg                 <  1:345222/69‑1 (MQ=255)
                        gagagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTg                 <  1:1244165/69‑1 (MQ=255)
                          gagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATg                                      >  1:585326/1‑46 (MQ=255)
                          gagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATg                                      >  1:85129/1‑46 (MQ=255)
                          gagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATg                                      >  1:2070788/1‑46 (MQ=255)
                          gagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATg                                      >  1:842067/1‑46 (MQ=255)
                          gagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATg                                      >  1:2151098/1‑46 (MQ=255)
                          gagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATg                                      >  1:3064009/1‑46 (MQ=255)
                          gagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATg                                      >  1:2201488/1‑46 (MQ=255)
                          gagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATg                                      >  1:2794878/1‑46 (MQ=255)
                          gagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATg                                      >  1:1585483/1‑46 (MQ=255)
                          gagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATg                                      >  1:471921/1‑46 (MQ=255)
                          gagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATg                                      >  1:2583443/1‑46 (MQ=255)
                          gagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCa               <  1:3192849/69‑1 (MQ=255)
                          gagCCTGAAAGGACAAAGCACGGACGTGCTCTATTCCGTGCAAATg                                      >  1:2061292/1‑46 (MQ=255)
                                   aGGACACATCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg        <  1:1107853/67‑1 (MQ=255)
                                   aGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg        <  1:1180712/67‑1 (MQ=255)
                                   aGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg        <  1:45981/67‑1 (MQ=255)
                                   aGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg        <  1:584540/67‑1 (MQ=255)
                                   aGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg        <  1:557912/67‑1 (MQ=255)
                                   aGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg        <  1:3243061/67‑1 (MQ=255)
                                   aGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg        <  1:1351791/67‑1 (MQ=255)
                                   aGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg        <  1:2727139/67‑1 (MQ=255)
                                   aGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg        <  1:2404060/67‑1 (MQ=255)
                                   aGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg        <  1:1562295/67‑1 (MQ=255)
                                   aGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg        <  1:217395/67‑1 (MQ=255)
                                   aGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg        <  1:889703/67‑1 (MQ=255)
                                   aGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg        <  1:2139469/67‑1 (MQ=255)
                                       cacaGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGaa  <  1:3423439/69‑1 (MQ=255)
                                       cacaGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGaa  <  1:1818036/69‑1 (MQ=255)
                                       cacaGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGaa  <  1:1746637/69‑1 (MQ=255)
                                              ctgACGTGCTTTATTCCGTGCAAATGCAGGGCAACCCCGGTTATTTGCa               <  1:901310/47‑1 (MQ=255)
                                              ctgACGTGCTTTATTCCGTGCAAATGCAGGGCAACCCCGGTTATTTGCa               <  1:991388/47‑1 (MQ=255)
                                              ctgACGTGCTTTATTCCGTGCAAATGCAGGGCAACCCCGGTTATTTGCa               <  1:46585/47‑1 (MQ=255)
                                              ctgACGTGCTTTATTCCGTGCAAATGCAGGGCAACCCCGGTTATTTGCa               <  1:462216/47‑1 (MQ=255)
                                              ctgACGTGCTTTATTCCGTGCAAATGCAGGGCAACCCCGGTTATTTGCa               <  1:3414582/47‑1 (MQ=255)
                                              ctgACGTGCTTTATTCCGTGCAAATGCAGGGCAACCCCGGTTATTTGCa               <  1:339742/47‑1 (MQ=255)
                                              ctgACGTGCTTTATTCCGTGCAAATGCAGGGCAACCCCGGTTATTTGCa               <  1:2505104/47‑1 (MQ=255)
                                              ctgACGTGCTTTATTCCGTGCAAATGCAGGGCAACCCCGGTTATTTGCa               <  1:2428571/47‑1 (MQ=255)
                                              ctgACGTGCTTTATTCCGTGCAAATGCAGGGCAACCCCGGTTATTTGCa               <  1:2130976/47‑1 (MQ=255)
                                              ctgACGTGCTTTATTCCGTGCAAATGCAGGGCAACCCCGGTTATTTGCa               <  1:191093/47‑1 (MQ=255)
                                              ctgACGTGCTTTATTCCGTGCAAATGCAGGGCAACCCCGGTTATTTGCa               <  1:1599567/47‑1 (MQ=255)
                                              ctgACGTGCTTTATTCCGTGCAAATGCAGGGCAACCCCGGTTATTTGCa               <  1:153853/47‑1 (MQ=255)
                                              ctgACGTGCTTTATTCCGTGCAAATGCAGGGCAACCCCGGTTATTTGCa               <  1:105521/47‑1 (MQ=255)
                                              cGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGATa                     <  1:77540/43‑1 (MQ=255)
                                               tgACGTGCTTTATTCCGTGCAAATGCAGGGCAACCCCGGTTATTTGct               <  1:2460423/47‑2 (MQ=255)
                                                        ctATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgttgt         >  1:536247/1‑45 (MQ=255)
                                                        ctATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgttgt         >  1:2963498/1‑45 (MQ=255)
                                                        ctATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgttgt         >  1:1769584/1‑45 (MQ=255)
                                                        |                                                   
TTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAA  >  minE/131843‑131950

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