Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,357,820 C→T 30.2% T434T (ACG→ACA yehB ← predicted outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,357,8200CT30.2% 89.8 / 41.1 93T434T (ACG→ACAyehBpredicted outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (28/36);  new base T (5/23);  total (34/59)
Fisher's exact test for biased strand distribution p-value = 1.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGTGATCGTTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATC  >  minE/1357757‑1357882
                                                               |                                                              
aaCGTGATCGTTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCg                                                            >  1:1071688/1‑68 (MQ=255)
aaCGTGATCGTTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGa                                                           >  1:3379587/1‑69 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCTACGtt                                                       <  1:439054/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:2885034/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:1616919/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:3026506/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:1673263/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:3020235/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:2906455/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:1761163/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:1054792/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:1924974/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:2858691/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:2475275/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:2180084/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:2375158/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:2353713/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:605533/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:3403184/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:1415762/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:1300621/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:1491436/69‑1 (MQ=255)
    tGATCATTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCCCCTGACCAAAACGTGTCGACGtt                                                       <  1:1163253/69‑1 (MQ=255)
      aTCGTTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTg                                                     >  1:2715500/1‑69 (MQ=255)
        cattAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:1469234/63‑1 (MQ=255)
        cGTTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTgg                                                    <  1:2178858/68‑1 (MQ=255)
        cGTTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTgg                                                    <  1:1993335/68‑1 (MQ=255)
           tAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTc                                                 >  1:1430237/1‑68 (MQ=255)
           tAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTc                                                 >  1:1816303/1‑68 (MQ=255)
                   cGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATttg                                        <  1:2068739/69‑1 (MQ=255)
                     gTAATCACGCGACGAATATCGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATttgtt                                      <  1:2419176/69‑1 (MQ=255)
                     gTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATttgt                                       <  1:46135/68‑1 (MQ=255)
                     gTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATttgt                                       <  1:1863712/68‑1 (MQ=255)
                       aaTCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGTTTGGCTCACAAATTtgttgt                                    >  1:293433/1‑69 (MQ=255)
                       aaTCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGTTTGGCTCACAAATTtgttgt                                    >  1:1421275/1‑69 (MQ=255)
                       aaTCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGTTTGGCTCACAAATTtgttgt                                    >  1:30618/1‑69 (MQ=255)
                       aaTCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGTTTGGCTCACAAATTtgttgt                                    >  1:2193697/1‑69 (MQ=255)
                       aaTCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGTTTGGCTCACAAATTtgttgt                                    >  1:1579173/1‑69 (MQ=255)
                       aaTCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTtgttgt                                    >  1:485846/1‑69 (MQ=255)
                       aaTCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTtgttgt                                    >  1:240731/1‑69 (MQ=255)
                             gcgACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCaa                              <  1:1346464/69‑1 (MQ=255)
                                acgaATAACGCCAGGCTGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAAttt                           <  1:160846/69‑1 (MQ=255)
                                acgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAAttt                           <  1:422624/69‑1 (MQ=255)
                                acgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAAttt                           <  1:2652034/69‑1 (MQ=255)
                                acgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAAttt                           <  1:493852/69‑1 (MQ=255)
                                acgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAAttt                           <  1:313292/69‑1 (MQ=255)
                                acgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAAttt                           <  1:2330312/69‑1 (MQ=255)
                                acgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAAttt                           <  1:1002443/69‑1 (MQ=255)
                                acgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAAttt                           <  1:2098635/69‑1 (MQ=255)
                                acgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAAttt                           <  1:1367025/69‑1 (MQ=255)
                                acgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAAttt                           <  1:2037176/69‑1 (MQ=255)
                                acgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAAttt                           <  1:1631766/69‑1 (MQ=255)
                                acgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAAttt                           <  1:1986581/69‑1 (MQ=255)
                                 cgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTg                          <  1:933242/69‑1 (MQ=255)
                                 cgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTg                          <  1:122589/69‑1 (MQ=255)
                                 cgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTg                          <  1:1539028/69‑1 (MQ=255)
                                 cgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTg                          <  1:3075172/69‑1 (MQ=255)
                                 cgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTg                          <  1:1694290/69‑1 (MQ=255)
                                  gaATAACGCCAGGCCGCCAGACCAAAACGTGTCGACGtt                                                       <  1:512232/39‑1 (MQ=255)
                                   aaTAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGAt                        <  1:3254465/69‑1 (MQ=255)
                                      aaCGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAAc                     <  1:2397634/69‑1 (MQ=255)
                                      aaCGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAAc                     <  1:1466985/69‑1 (MQ=255)
                                      aaCGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAAc                     <  1:3221242/69‑1 (MQ=255)
                                      aaCGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAAc                     <  1:2703848/69‑1 (MQ=255)
                                      aaCGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAAc                     <  1:2987712/69‑1 (MQ=255)
                                      aaCGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAAc                     <  1:2221905/69‑1 (MQ=255)
                                      aaCGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAAc                     <  1:1746585/69‑1 (MQ=255)
                                                     gACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAaca      <  1:489881/69‑1 (MQ=255)
                                                     gACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAaca      <  1:173544/69‑1 (MQ=255)
                                                     gACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAaca      <  1:3231808/69‑1 (MQ=255)
                                                       ccAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAAttt                           <  1:1079544/46‑1 (MQ=255)
                                                       ccAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAAcaca    >  1:1345872/1‑69 (MQ=255)
                                                       ccAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAAcaca    >  1:198826/1‑69 (MQ=255)
                                                        cAAAACGCGTCGACGTTTGGCTGACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACAt   >  1:3123083/1‑69 (MQ=255)
                                                        cAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACAt   >  1:1872647/1‑69 (MQ=255)
                                                        cAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACAt   >  1:1247409/1‑69 (MQ=255)
                                                        cAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACAt   >  1:1358169/1‑69 (MQ=255)
                                                        cAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACAt   >  1:1661957/1‑69 (MQ=255)
                                                         aaaaCGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATc  >  1:3391514/1‑69 (MQ=255)
                                                         aaaaCGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATc  >  1:2453151/1‑69 (MQ=255)
                                                         aaaaCGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATc  >  1:3421395/1‑69 (MQ=255)
                                                         aaaaCGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATc  >  1:2060020/1‑69 (MQ=255)
                                                         aaaaCGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATc  >  1:3376248/1‑69 (MQ=255)
                                                         aaaaCGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATc  >  1:461159/1‑69 (MQ=255)
                                                         aaaaCGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATc  >  1:318138/1‑69 (MQ=255)
                                                         aaaaCGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATc  >  1:2599940/1‑69 (MQ=255)
                                                         aaaaCGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATc  >  1:3099408/1‑69 (MQ=255)
                                                         aaaaCGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATc  >  1:2001878/1‑69 (MQ=255)
                                                         aaaaCGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATc  >  1:1649325/1‑69 (MQ=255)
                                                         aaaaCGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATc  >  1:53283/1‑69 (MQ=255)
                                                         aaaaCGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATc  >  1:872676/1‑69 (MQ=255)
                                                         aaaaCGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATc  >  1:2967595/1‑69 (MQ=255)
                                                         aaaaCGAGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATc  >  1:436641/1‑69 (MQ=255)
                                                               |                                                              
AACGTGATCGTTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATC  >  minE/1357757‑1357882

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: