Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,014,607 T→C 52.4% V46V (GTA→GTG exbB ← membrane spanning protein in TonB‑ExbB‑ExbD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,014,6070TC52.4% 6.5 / 74.1 86V46V (GTA→GTGexbBmembrane spanning protein in TonB‑ExbB‑ExbD complex
Reads supporting (aligned to +/- strand):  ref base T (10/31);  new base C (13/32);  total (23/63)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

AACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAG  >  minE/2014546‑2014646
                                                             |                                       
aaCGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGctct                                  <  1:2832628/69‑1 (MQ=255)
aaCGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGctct                                  <  1:467787/69‑1 (MQ=255)
aaCGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGctct                                  <  1:2325246/69‑1 (MQ=255)
 aCGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCtt                                 >  1:804914/1‑69 (MQ=255)
         cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:916418/1‑66 (MQ=255)
         cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:501398/1‑66 (MQ=255)
         cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:2038047/1‑66 (MQ=255)
         cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:512382/1‑66 (MQ=255)
         cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:465832/1‑66 (MQ=255)
         cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:410608/1‑66 (MQ=255)
         cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:2194435/1‑66 (MQ=255)
         cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:1815074/1‑66 (MQ=255)
         cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:567286/1‑66 (MQ=255)
         cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:654356/1‑66 (MQ=255)
         cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:3212675/1‑66 (MQ=255)
         cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:2820659/1‑66 (MQ=255)
         cagccagcAGTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTga                            >  1:2556326/1‑66 (MQ=255)
           gccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTgaagaa                        <  1:3374512/68‑1 (MQ=255)
           gccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTgaagaa                        <  1:347442/68‑1 (MQ=255)
           gccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTgaagaa                        <  1:3266783/68‑1 (MQ=255)
           gccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTgaagaa                        <  1:3128304/68‑1 (MQ=255)
           gccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTgaagaa                        <  1:2083515/68‑1 (MQ=255)
           gccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTgaagaa                        <  1:11357/68‑1 (MQ=255)
           gccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTgaagaa                        <  1:835888/68‑1 (MQ=255)
           gccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATCGAAGAACTCTACGCTCTTACTgaagaa                        <  1:1909071/68‑1 (MQ=255)
            ccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:1934645/68‑1 (MQ=255)
            ccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:2811025/68‑1 (MQ=255)
            ccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:1115103/68‑1 (MQ=255)
            ccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:1837128/68‑1 (MQ=255)
            ccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:122761/68‑1 (MQ=255)
            ccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:2282383/68‑1 (MQ=255)
            ccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:1532221/68‑1 (MQ=255)
            ccagcAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:3341580/68‑1 (MQ=255)
             cagcagTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGAt                     <  1:2326508/69‑1 (MQ=255)
             cagcagTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGAt                     <  1:1763613/69‑1 (MQ=255)
             cagcagTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGAt                     <  1:1881443/69‑1 (MQ=255)
                  gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:734079/69‑1 (MQ=255)
                  gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:3227497/69‑1 (MQ=255)
                  gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:828943/69‑1 (MQ=255)
                  gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:1089154/69‑1 (MQ=255)
                  gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:520904/69‑1 (MQ=255)
                  gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:439675/69‑1 (MQ=255)
                  gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:2776336/69‑1 (MQ=255)
                  gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:2621651/69‑1 (MQ=255)
                  gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:2579543/69‑1 (MQ=255)
                  gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:2161221/69‑1 (MQ=255)
                  gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:2504122/69‑1 (MQ=255)
                  gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:1801152/69‑1 (MQ=255)
                  gTTGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:2372474/69‑1 (MQ=255)
                   ttGCTGCTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:2330944/68‑1 (MQ=255)
                    tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:1839333/69‑1 (MQ=255)
                    tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:3421818/69‑1 (MQ=255)
                    tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:3397026/69‑1 (MQ=255)
                    tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:3351657/69‑1 (MQ=255)
                    tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2106363/69‑1 (MQ=255)
                    tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2102210/69‑1 (MQ=255)
                    tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:1964518/69‑1 (MQ=255)
                    tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:619444/69‑1 (MQ=255)
                    tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2897311/69‑1 (MQ=255)
                    tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:1463468/69‑1 (MQ=255)
                    tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:1379231/69‑1 (MQ=255)
                    tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:1373434/69‑1 (MQ=255)
                    tgctgcTCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTAATGAAGAAGATTGCCCAg              <  1:3177254/69‑1 (MQ=255)
                     gctgctCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              >  1:1076748/1‑68 (MQ=255)
                     gctgctCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              >  1:41493/1‑68 (MQ=255)
                     gctgctCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              >  1:663262/1‑68 (MQ=255)
                     gctgctCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              >  1:2897169/1‑68 (MQ=255)
                     gctgctCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2206796/68‑1 (MQ=255)
                      ctgctCGCGCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:398981/65‑1 (MQ=255)
                         ctCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGaagaag                       <  1:2485577/55‑1 (MQ=255)
                           cgcgCTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGccc                <  1:2334793/60‑1 (MQ=255)
                             cgcTTAAGACGACGCTTCTGATTGAAGAACTCCAGGCTCTTACTGAAGAAGATTGCCCAg              <  1:1486925/60‑1 (MQ=255)
                              gcTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              <  1:3033336/59‑1 (MQ=255)
                              gcTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2719296/59‑1 (MQ=255)
                              gcTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2736329/59‑1 (MQ=255)
                              gcTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2507309/59‑1 (MQ=255)
                              gcTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              <  1:3200928/59‑1 (MQ=255)
                              gcTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2139532/59‑1 (MQ=255)
                               cTTAAGACGACGCTTCTGATTGAAGAACTCCACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2147980/58‑1 (MQ=255)
                                ttAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTa                                <  1:2147349/39‑1 (MQ=255)
                                 tAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAg              <  1:2718254/56‑1 (MQ=255)
                                 tAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAg  >  1:3125575/1‑68 (MQ=255)
                                 tAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAg  >  1:2925206/1‑68 (MQ=255)
                                 tAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAg  >  1:2340184/1‑68 (MQ=255)
                                 tAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAg  >  1:1431247/1‑68 (MQ=255)
                                 tAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAg  >  1:831126/1‑68 (MQ=255)
                                                             |                                       
AACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAG  >  minE/2014546‑2014646

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: