Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,892,348 G→A 33.3% intergenic (+46/‑94) ytfQ → / → ytfR predicted sugar transporter subunit/predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,892,3480GA33.3% 69.0 / 48.7 60intergenic (+46/‑94)ytfQ/ytfRpredicted sugar transporter subunit/predicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (21/19);  new base A (11/9);  total (32/28)
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

AAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTG  >  minE/2892279‑2892414
                                                                     |                                                                  
aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTcg                                                                    <  1:297379/70‑1 (MQ=255)
   tgcctgatgcgctacgcttatcaggcctacaattctcgcgccACTCGTAGGCCGGATAAGGCGTTCAcgcc                                                                >  1:212568‑M2/43‑71 (MQ=255)
   tgcctgatgcgctacgcttatcaggcctacaattctcgcgccACTCGTAGGCCGGATAAGGCGTTCAcgcc                                                                >  1:3638845‑M2/43‑71 (MQ=255)
   gcctgatgcgctacgcttatcaggcctacagattcttgagccATTCGTAGGCCGGATAAGGCGTTCacgc                                                                 >  1:434558‑M2/43‑70 (MQ=255)
    agaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcgc                                                                 >  1:2152367/1‑69 (MQ=255)
     ctgatgcgctacgcttatcaggcctacagattctcgcgccACTCGTAGGCCGGATAAGGCGTCCAcgccgc                                                              <  1:2768907‑M2/31‑1 (MQ=255)
     cctgatgcgctacgcttatcaggcctacaattctcgcgccACTCGTAGGCCGGATAAGGCGTTCAcgcc                                                                <  1:1742409‑M2/29‑1 (MQ=255)
       gatgcgctacgcttatcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                            >  1:3424724‑M2/39‑71 (MQ=255)
        atgcgctacgcttatcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                           >  1:1417802‑M2/38‑71 (MQ=255)
        atgcgctacgcttatcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                           >  1:1698281‑M2/38‑71 (MQ=255)
        atgcgctacgcttatcaggcctacagattctcgcgccACTCGTAGGCCGGATAAGGCGTCCACGCCGCATc                                                           <  1:3367228‑M2/34‑1 (MQ=255)
          tgcgctacgcttatcaggcctacaattctcgcgccACTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                         >  1:1593977‑M2/36‑71 (MQ=255)
           cgctacgcttatcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                         <  1:2959088‑M2/36‑1 (MQ=255)
           cgctacgcttatcaggcctacagattctcgcgccACTCGTAGGCCGGATAAGGCGTCCACGCCGCATCCgg                                                        <  1:2360872‑M2/37‑1 (MQ=255)
             tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGc                                                       <  1:104092/70‑1 (MQ=255)
             ctacgcttatcaggcctacagattctcgcgccACTCGTAGGCCGGATAAGGCGTCCACGCCGCATCCGgca                                                      <  1:1879157‑M2/39‑2 (MQ=255)
              tacgcttatcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCggcat                                                     >  1:2696707‑M2/32‑69 (MQ=255)
              tacgcttatcaggcctacagattctcgcgccACTCGTAGGCCGGATAAGGCGTCCACGCCGCATCCggcag                                                     >  1:1085164‑M2/32‑69 (MQ=255)
               tacgcttatcaggcctacaattctcgcgccACTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAt                                                    >  1:1709259‑M2/31‑71 (MQ=255)
               gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGa                                                     <  1:47971/70‑1 (MQ=255)
                ggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGa                                                     <  1:712316/69‑1 (MQ=255)
                cgcttatcaggcctacagattctcgcgccACTCGTAGGCCGGATAAGGCGTCCACGCCGCATCCggcag                                                     <  1:2136661‑M2/40‑3 (MQ=255)
                   ttatcaggcctacagattctcgcgccACTCGTAGGCCGGATAAGGCGTCCACGCCGCa                                                             >  1:3110272‑M2/27‑58 (MQ=255)
                    ttGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTg                                                >  1:2705674/1‑70 (MQ=255)
                     tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGcgccgc                                                              >  1:3850029/1‑55 (MQ=255)
                     tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGa                                                     >  1:1298471/1‑64 (MQ=255)
                     atcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTg                                                <  1:339089‑M2/45‑1 (MQ=255)
                     atcaggcctacagattcttgcgccATTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTg                                                <  1:585171‑M2/45‑1 (MQ=255)
                      gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCg                                                         >  1:3652576/1‑59 (MQ=18)
                      gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATcc                                                          <  1:1574127/58‑1 (MQ=255)
                      gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGt                                                 >  1:995263/1‑67 (MQ=255)
                      gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCa                                              <  1:3526392/70‑1 (MQ=255)
                      gAGGGGTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTcgcg                                                                  <  1:1045342/50‑1 (MQ=255)
                       aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACt                                            <  1:1563226/71‑1 (MQ=255)
                          gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGaa                                         <  1:3731858/71‑1 (MQ=255)
                           ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGcgcc                                                                >  1:2165178/1‑47 (MQ=255)
                           ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGcgcc                                                                >  1:3394226/1‑47 (MQ=255)
                           ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATg                                                   <  1:933582/60‑1 (MQ=255)
                              acagattcttgcgccATTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCggcaa                                                     >  1:41356‑M2/16‑53 (MQ=255)
                               tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAGGcat                                     >  1:1511977/1‑68 (MQ=255)
                                aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                   >  1:3206027/1‑71 (MQ=255)
                                 tGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                                    <  1:3604240/69‑1 (MQ=255)
                                  gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTg                                                <  1:3466736/56‑1 (MQ=255)
                                       aTGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACg                            >  1:3843917/1‑71 (MQ=255)
                                            gATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAc                       >  1:365581/1‑71 (MQ=255)
                                             aTTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                      >  1:72494/1‑71 (MQ=255)
                                              ttCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcg                     >  1:3310600/1‑71 (MQ=255)
                                                cGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgc                    <  1:3038261/70‑1 (MQ=255)
                                                 gTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcg                   <  1:1273212/70‑1 (MQ=255)
                                                    ggCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                >  1:3180584/1‑70 (MQ=255)
                                                    ggCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                <  1:2253388/70‑1 (MQ=255)
                                                     gCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa              <  1:275933/71‑1 (MQ=255)
                                                     gCCGGATAAGGCGCACGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa              <  1:2195048/71‑1 (MQ=255)
                                                       cGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt             >  1:501586/1‑70 (MQ=255)
                                                        ggATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa           <  1:124883/71‑1 (MQ=255)
                                                          aTAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg         <  1:2173154/71‑1 (MQ=255)
                                                          aTAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg         >  1:1935292/1‑71 (MQ=255)
                                                           tAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc        >  1:2525184/1‑71 (MQ=255)
                                                           tAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc        >  1:723037/1‑71 (MQ=255)
                                                               gcgcTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa     <  1:430234/70‑1 (MQ=255)
                                                                 gcTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCtt               <  1:450083/58‑1 (MQ=255)
                                                                 gcTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa     >  1:1043045/1‑68 (MQ=255)
                                                                 gcTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTg  >  1:1047539/1‑71 (MQ=255)
                                                                 gcTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTg  >  1:1121590/1‑71 (MQ=255)
                                                                     |                                                                  
AAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTG  >  minE/2892279‑2892414

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: