Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 382,496 G→A 18.7% V510V (GTG→GTA entE → 2,3‑dihydroxybenzoate‑AMP ligase component of enterobactin synthase multienzyme complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE382,4960GA18.7% 114.2 / 15.3 74V510V (GTG→GTAentE2,3‑dihydroxybenzoate‑AMP ligase component of enterobactin synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base A (12/2);  total (43/31)
Fisher's exact test for biased strand distribution p-value = 3.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TGCGTCGTTTCCTGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCAC  >  minE/382432‑382563
                                                                |                                                                   
tGCGTCGTTTCCTGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGAtt                                                                 >  1:3419727/1‑69 (MQ=255)
            tGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTg                                                      <  1:582096/68‑1 (MQ=255)
            tGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTg                                                      <  1:1828531/68‑1 (MQ=255)
            tGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTg                                                      <  1:205592/68‑1 (MQ=255)
            tGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTg                                                      <  1:1686224/68‑1 (MQ=255)
            tGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTg                                                      <  1:1650801/68‑1 (MQ=255)
            tGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTg                                                      <  1:2376704/68‑1 (MQ=255)
            tGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTg                                                      <  1:687515/68‑1 (MQ=255)
            tGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTg                                                      <  1:2846813/68‑1 (MQ=255)
            tGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTg                                                      <  1:1530997/68‑1 (MQ=255)
            tGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTg                                                      <  1:1360199/68‑1 (MQ=255)
            tGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTg                                                      <  1:3030573/68‑1 (MQ=255)
            tGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTg                                                      <  1:568573/68‑1 (MQ=255)
            tGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTg                                                      <  1:515194/68‑1 (MQ=255)
            tGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTg                                                      <  1:3351184/68‑1 (MQ=255)
            tGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTg                                                      <  1:3395660/68‑1 (MQ=255)
             gCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGAc                                                    >  1:757221/1‑69 (MQ=255)
             gCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGAc                                                    >  1:948721/1‑69 (MQ=255)
                       ggTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCgg                                           >  1:359081/1‑68 (MQ=255)
                       ggTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCgg                                           >  1:253218/1‑68 (MQ=255)
                       ggTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCgg                                           >  1:1567294/1‑68 (MQ=255)
                       ggTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCgg                                           >  1:1571604/1‑68 (MQ=255)
                       ggTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCgg                                           >  1:943995/1‑68 (MQ=255)
                       ggTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCgg                                           >  1:1349152/1‑68 (MQ=255)
                       ggTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCgg                                           >  1:1608253/1‑68 (MQ=255)
                             gCCGATTTTAAATTACCGGATCGCGTGGAGTGTCTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTc                                    >  1:1193745/1‑69 (MQ=255)
                             gCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGt                                     >  1:2647285/1‑68 (MQ=255)
                             gCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTc                                    >  1:819056/1‑69 (MQ=255)
                             gCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTc                                    >  1:2810980/1‑69 (MQ=255)
                             gCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTc                                    >  1:2806247/1‑69 (MQ=255)
                             gCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTc                                    >  1:1790820/1‑69 (MQ=255)
                             gCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTc                                    >  1:1335758/1‑69 (MQ=255)
                             gCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTc                                    >  1:2527664/1‑69 (MQ=255)
                             gCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTc                                    >  1:73108/1‑69 (MQ=255)
                             gCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTc                                    >  1:1065807/1‑69 (MQ=255)
                             gCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTc                                    >  1:1913449/1‑69 (MQ=255)
                             gCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCAATTCCGCTGACGGCGGTCGGGAAAGTc                                    >  1:1151655/1‑69 (MQ=255)
                               cGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCg                                   >  1:617563/1‑68 (MQ=255)
                                gAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGa                                  >  1:248888/1‑68 (MQ=255)
                                gAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGa                                  >  1:3310136/1‑68 (MQ=255)
                                     tAAATTACCGGATCTCGTGGAGTGTGTAGATTCACTTCCGCTGACGGCGGTCGGGAAAg                                      >  1:2177537/1‑59 (MQ=255)
                                     tAAATTACCGGATCGCGTGGAGTGTGTAGATTCACTTCCGCTGACGGCGGTCGGGAAAg                                      >  1:2734532/1‑59 (MQ=255)
                                     tAAATTACCGGATCGCGTGGAGTGTGTAGATTCACTTCCGCTGACGGCGGTCGGGAAAg                                      >  1:3207518/1‑59 (MQ=255)
                                     tAAATTACCGGATCGCGTGGAGTGTGTAGATTCACTTCCGCTGACGGCGGTCGGGAAAg                                      >  1:429025/1‑59 (MQ=255)
                                     tAAATTACCGGATCGCGTGGAGTGTGTAGATTCACTTCCGCTGACGGCGGTCGGGAAAg                                      >  1:3220404/1‑59 (MQ=255)
                                     tAAATTACCGGATCGCGTGGAGTGTGTAGATTCACTTCCGCTGACGGCGGTCGGGAAAg                                      >  1:3022501/1‑59 (MQ=255)
                                     tAAATTACCGGATCGCGTGGAGTGTGTAGATTCACTTCCGCTGACGGCGGTCGGGAAAg                                      >  1:2965484/1‑59 (MQ=255)
                                     tAAATTACCGGATCGCGTGGAGTGTGTAGATTCACTTCCGCTGACGGCGGTCGGGAAAg                                      >  1:2791185/1‑59 (MQ=255)
                                     tAAATTACCGGATCGCGTGGAGTGTGTAGATTCACTTCCGCTGACGGCGGTCGGGAAAg                                      >  1:1011928/1‑59 (MQ=255)
                                     tAAATTACCGGATCGCGTGGAGTGTGTAGATTCACTTCCGCTGACGGCGGTCGGGAAAg                                      >  1:2548771/1‑59 (MQ=255)
                                     tAAATTACCGGATCGCGTGGAGTGTGTAGATTCACTTCCGCTGACGGCGGTCGGGAAAg                                      >  1:2200580/1‑59 (MQ=255)
                                     tAAATTACCGGATCGCGTGGAGTGTGTAGATTCACTTCCGCTGACGGCGGTCGGGAAAg                                      >  1:2113715/1‑59 (MQ=255)
                                        aTTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAAc                          <  1:1641925/68‑1 (MQ=255)
                                         ttACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACa                         <  1:2567209/68‑1 (MQ=255)
                                          tACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACTGCGGTCGGGAAAGTCGATAAAAAACAAt                       <  1:168651/69‑1 (MQ=255)
                                          tACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAAt                       <  1:671609/69‑1 (MQ=255)
                                          tACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAAt                       <  1:726656/69‑1 (MQ=255)
                                          tACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAAt                       <  1:754226/69‑1 (MQ=255)
                                          tACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAAt                       <  1:2706824/69‑1 (MQ=255)
                                          tACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAAt                       <  1:1165653/69‑1 (MQ=255)
                                          tACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAAt                       <  1:2947011/69‑1 (MQ=255)
                                          tACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAAt                       <  1:1596863/69‑1 (MQ=255)
                                          tACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAAt                       <  1:2493750/69‑1 (MQ=255)
                                          tACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAAt                       <  1:2409270/69‑1 (MQ=255)
                                          tACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGGCGATAAAAAACAAt                       <  1:1020045/69‑1 (MQ=255)
                                           aCCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAAtt                      >  1:386715/1‑69 (MQ=255)
                                           aCCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAAtt                      >  1:3210204/1‑69 (MQ=255)
                                           aCCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAAtt                      >  1:636424/1‑69 (MQ=255)
                                           aCCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAAtt                      >  1:640427/1‑69 (MQ=255)
                                           aCCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAAtt                      >  1:2349083/1‑69 (MQ=255)
                                           aCCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAAtt                      >  1:2024621/1‑69 (MQ=255)
                                                aTCGCTTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTc                 <  1:309946/69‑1 (MQ=255)
                                                       ggAGTGTGTAGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGtggct          <  1:409274/69‑1 (MQ=255)
                                                       ggAGTGTGTAGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGtggct          <  1:1759090/69‑1 (MQ=255)
                                                            gtgtgGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGt     <  1:2952128/69‑1 (MQ=255)
                                                               tgGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCAc  >  1:2618941/1‑69 (MQ=255)
                                                               tgGATTAACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCAc  >  1:3052915/1‑69 (MQ=255)
                                                                |                                                                   
TGCGTCGTTTCCTGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCGCGTGGAGTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCAC  >  minE/382432‑382563

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