Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 365,118 G→T 35.7% intergenic (‑133/+33) entD ← / ← fepA phosphopantetheinyltransferase component of enterobactin synthase multienzyme complex/iron‑enterobactin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE365,1180GT35.7% 87.2 / 60.8 84intergenic (‑133/+33)entD/fepAphosphopantetheinyltransferase component of enterobactin synthase multienzyme complex/iron‑enterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (29/25);  new base T (18/12);  total (47/37)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGT  >  minE/365057‑365183
                                                             |                                                                 
tatTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCa                                                            >  1:3302608/1‑69 (MQ=255)
    gAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACttt                                                        <  1:1466419/69‑1 (MQ=255)
    gAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACttt                                                        <  1:3394292/69‑1 (MQ=255)
    gAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACttt                                                        <  1:3120798/69‑1 (MQ=255)
               aTTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCt                                             >  1:2235634/1‑69 (MQ=255)
               aTTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCt                                             >  1:86094/1‑69 (MQ=255)
               aTTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCt                                             >  1:2623761/1‑69 (MQ=255)
                tttttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTg                                            >  1:2637766/1‑69 (MQ=255)
                          ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCa                                   <  1:2160782/68‑1 (MQ=255)
                          ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCa                                   <  1:1434323/68‑1 (MQ=255)
                           cGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGc                                    >  1:2827592/1‑66 (MQ=255)
                              aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                              <  1:967423/69‑1 (MQ=255)
                              aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                              <  1:2325290/69‑1 (MQ=255)
                              aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                              <  1:2666137/69‑1 (MQ=255)
                              aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                              <  1:2863787/69‑1 (MQ=255)
                              aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                              <  1:3381361/69‑1 (MQ=255)
                              aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                              <  1:927047/69‑1 (MQ=255)
                              aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                              <  1:740186/69‑1 (MQ=255)
                              aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                              <  1:1502301/69‑1 (MQ=255)
                              aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                              <  1:1529655/69‑1 (MQ=255)
                              aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                              <  1:1584568/69‑1 (MQ=255)
                              aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                              <  1:3368534/69‑1 (MQ=255)
                               tAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                              <  1:1029349/68‑1 (MQ=255)
                                   gCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTg                          <  1:78020/68‑1 (MQ=255)
                                   gCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTg                          <  1:2604239/68‑1 (MQ=255)
                                   gCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTg                          <  1:823174/68‑1 (MQ=255)
                                        cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                    >  1:3034848/1‑69 (MQ=255)
                                        cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                    >  1:2855850/1‑69 (MQ=255)
                                        cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                    >  1:409340/1‑69 (MQ=255)
                                        cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                    >  1:3171043/1‑69 (MQ=255)
                                        cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                    >  1:2634055/1‑69 (MQ=255)
                                        cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                    >  1:2591460/1‑69 (MQ=255)
                                        cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                    >  1:2396832/1‑69 (MQ=255)
                                        cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                    >  1:2542890/1‑69 (MQ=255)
                                        cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                    >  1:1661740/1‑69 (MQ=255)
                                        cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                    >  1:1718799/1‑69 (MQ=255)
                                        cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                    >  1:1942381/1‑69 (MQ=255)
                                        cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                    >  1:2092800/1‑69 (MQ=255)
                                               cATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTg                          <  1:207253/56‑1 (MQ=255)
                                                aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            <  1:2084405/69‑1 (MQ=255)
                                                aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            <  1:1756462/69‑1 (MQ=255)
                                                aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            <  1:2323848/69‑1 (MQ=255)
                                                aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            <  1:1507019/69‑1 (MQ=255)
                                                aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            <  1:2976494/69‑1 (MQ=255)
                                                aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            <  1:1490904/69‑1 (MQ=255)
                                                aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            <  1:2783450/69‑1 (MQ=255)
                                                aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            <  1:204085/69‑1 (MQ=255)
                                                aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            <  1:263629/69‑1 (MQ=255)
                                                aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            <  1:261508/69‑1 (MQ=255)
                                                aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            <  1:252310/69‑1 (MQ=255)
                                                aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            <  1:1167182/69‑1 (MQ=255)
                                                aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            <  1:2496100/69‑1 (MQ=255)
                                                aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            <  1:2430136/69‑1 (MQ=255)
                                                aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTAGGTGTTTACGCTCa            <  1:616472/69‑1 (MQ=255)
                                                  ccGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCata          <  1:3193085/69‑1 (MQ=255)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGCAGTGAGTATTAATGCTCATATACCacgt  >  1:1190122/1‑69 (MQ=21)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCAtata        >  1:1401212/1‑63 (MQ=21)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt  >  1:395254/1‑69 (MQ=25)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt  >  1:114176/1‑69 (MQ=25)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt  >  1:159115/1‑69 (MQ=25)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt  >  1:862452/1‑69 (MQ=25)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt  >  1:57234/1‑69 (MQ=25)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt  >  1:1256798/1‑69 (MQ=25)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt  >  1:1260868/1‑69 (MQ=25)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt  >  1:3095133/1‑69 (MQ=25)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt  >  1:2281444/1‑69 (MQ=25)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt  >  1:1881992/1‑69 (MQ=25)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt  >  1:195893/1‑69 (MQ=25)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt  >  1:3284071/1‑69 (MQ=25)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt  >  1:2608967/1‑69 (MQ=25)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt  >  1:209894/1‑69 (MQ=25)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt  >  1:2100706/1‑69 (MQ=25)
                                                          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt  >  1:2172339/1‑69 (MQ=25)
                                                          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            >  1:332060/1‑59 (MQ=255)
                                                          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            >  1:3188574/1‑59 (MQ=255)
                                                          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            >  1:1863334/1‑59 (MQ=255)
                                                          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa            >  1:536719/1‑59 (MQ=255)
                                                          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt  >  1:2620567/1‑69 (MQ=255)
                                                          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt  >  1:2268017/1‑69 (MQ=255)
                                                          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt  >  1:2417860/1‑69 (MQ=255)
                                                          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt  >  1:800278/1‑69 (MQ=255)
                                                          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt  >  1:1616896/1‑69 (MQ=255)
                                                          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt  >  1:3278883/1‑69 (MQ=255)
                                                          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacg   >  1:2517876/1‑68 (MQ=255)
                                                             |                                                                 
TATTGAATTGGCAAGATTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGT  >  minE/365057‑365183

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