Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 365,150 A→T 34.5% intergenic (‑165/+1) 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,1500AT34.5% 112.6 / 87.1 110intergenic (‑165/+1)entD/fepAphosphopantetheinyltransferase component of enterobactin synthase multienzyme complex/iron‑enterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (30/42);  new base T (18/20);  total (48/62)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

CCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGC  >  minE/365083‑365216
                                                                   |                                                                  
ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCa                                                                    <  1:2160782/68‑1 (MQ=255)
ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCa                                                                    <  1:1434323/68‑1 (MQ=255)
    aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                                                               <  1:967423/69‑1 (MQ=255)
    aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                                                               <  1:2325290/69‑1 (MQ=255)
    aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                                                               <  1:927047/69‑1 (MQ=255)
    aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                                                               <  1:1584568/69‑1 (MQ=255)
    aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                                                               <  1:1529655/69‑1 (MQ=255)
    aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                                                               <  1:1502301/69‑1 (MQ=255)
    aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                                                               <  1:740186/69‑1 (MQ=255)
    aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                                                               <  1:2666137/69‑1 (MQ=255)
    aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                                                               <  1:3368534/69‑1 (MQ=255)
    aTAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                                                               <  1:3381361/69‑1 (MQ=255)
     tAAGGCGTTCACGCCGCATCCGGCATGAACTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGa                                                               <  1:1029349/68‑1 (MQ=255)
         gCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTg                                                           <  1:2604239/68‑1 (MQ=255)
         gCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTg                                                           <  1:78020/68‑1 (MQ=255)
         gCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTg                                                           <  1:823174/68‑1 (MQ=255)
              cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                                                     >  1:2634055/1‑69 (MQ=255)
              cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                                                     >  1:3034848/1‑69 (MQ=255)
              cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                                                     >  1:2591460/1‑69 (MQ=255)
              cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                                                     >  1:2092800/1‑69 (MQ=255)
              cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                                                     >  1:2855850/1‑69 (MQ=255)
              cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                                                     >  1:2542890/1‑69 (MQ=255)
              cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                                                     >  1:1718799/1‑69 (MQ=255)
              cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                                                     >  1:1661740/1‑69 (MQ=255)
              cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                                                     >  1:3171043/1‑69 (MQ=255)
              cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                                                     >  1:409340/1‑69 (MQ=255)
                     cATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTg                                                           <  1:207253/56‑1 (MQ=255)
                      aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             <  1:204085/69‑1 (MQ=255)
                      aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             <  1:2323848/69‑1 (MQ=255)
                      aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             <  1:1167182/69‑1 (MQ=255)
                      aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             <  1:2430136/69‑1 (MQ=255)
                      aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             <  1:2496100/69‑1 (MQ=255)
                      aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             <  1:252310/69‑1 (MQ=255)
                      aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             <  1:2084405/69‑1 (MQ=255)
                      aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             <  1:261508/69‑1 (MQ=255)
                      aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             <  1:263629/69‑1 (MQ=255)
                      aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             <  1:2783450/69‑1 (MQ=255)
                      aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             <  1:1507019/69‑1 (MQ=255)
                      aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             <  1:1756462/69‑1 (MQ=255)
                      aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             <  1:2976494/69‑1 (MQ=255)
                      aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTAGGTGTTTACGCTCa                                             <  1:616472/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:3284071/1‑69 (MQ=25)
                                aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                   >  1:862452/1‑69 (MQ=25)
                                aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                   >  1:3095133/1‑69 (MQ=25)
                                aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                   >  1:2281444/1‑69 (MQ=25)
                                aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                   >  1:1260868/1‑69 (MQ=25)
                                aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                   >  1:2172339/1‑69 (MQ=25)
                                aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                   >  1:159115/1‑69 (MQ=25)
                                aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                   >  1:209894/1‑69 (MQ=25)
                                aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                   >  1:1881992/1‑69 (MQ=25)
                                aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                   >  1:195893/1‑69 (MQ=25)
                                aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             >  1:3188574/1‑59 (MQ=255)
                                aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             >  1:536719/1‑59 (MQ=255)
                                aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             >  1:332060/1‑59 (MQ=255)
                                aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             >  1:1863334/1‑59 (MQ=255)
                                aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt                                   >  1:2620567/1‑69 (MQ=255)
                                aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt                                   >  1:800278/1‑69 (MQ=255)
                                aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt                                   >  1:2417860/1‑69 (MQ=255)
                                aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt                                   >  1:2268017/1‑69 (MQ=255)
                                aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt                                   >  1:3278883/1‑69 (MQ=255)
                                aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacg                                    >  1:2517876/1‑68 (MQ=255)
                                      gcgcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTcc                             >  1:11437/1‑69 (MQ=255)
                                      gcgcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTcc                             >  1:1173313/1‑69 (MQ=255)
                                      gcgcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTcc                             >  1:2325093/1‑69 (MQ=255)
                                      gcgcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTcc                             >  1:1749554/1‑69 (MQ=255)
                                      gcgcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTcc                             >  1:2136116/1‑69 (MQ=255)
                                         cACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGTGTGTTTACGCTCATATACCACGTACGTCCCgg                          <  1:900934/69‑1 (MQ=255)
                                         cACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                          <  1:1258097/69‑1 (MQ=255)
                                         cACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                          <  1:2875507/69‑1 (MQ=255)
                                         cACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                          <  1:1584762/69‑1 (MQ=255)
                                         cACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                          <  1:1236411/69‑1 (MQ=255)
                                         cACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                          <  1:894142/69‑1 (MQ=255)
                                         cACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                          <  1:1392747/69‑1 (MQ=255)
                                          aCTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             <  1:892423/49‑1 (MQ=255)
                                                tCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCACGTACGTCCCGGCTCGtt                    <  1:1401423/68‑1 (MQ=255)
                                                tCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCACGTACGTCCCGGCTCGtt                    <  1:762665/68‑1 (MQ=255)
                                                tCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCACGTACGTCCCGGCTCGtt                    <  1:1886658/68‑1 (MQ=255)
                                                tCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCACGTACGTCCCGGCTCGtt                    <  1:68819/68‑1 (MQ=255)
                                                tCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCACGTACGTCCCGGCTCGtt                    <  1:2083893/68‑1 (MQ=255)
                                                tCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCACGTACGTCCCGGCTCGtt                    <  1:882900/68‑1 (MQ=255)
                                                tCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCACGTACGTCCCGGCTCGtt                    <  1:1859740/68‑1 (MQ=255)
                                                tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                     <  1:1181239/67‑1 (MQ=255)
                                                tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                     <  1:1195690/67‑1 (MQ=255)
                                                tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                     <  1:2402551/67‑1 (MQ=255)
                                                tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                     <  1:1930080/67‑1 (MQ=255)
                                                tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                     <  1:1978392/67‑1 (MQ=255)
                                                tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                     <  1:73610/67‑1 (MQ=255)
                                                tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                     <  1:2935291/67‑1 (MQ=255)
                                                tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                     <  1:1296697/67‑1 (MQ=255)
                                                tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                     <  1:3051948/67‑1 (MQ=255)
                                                tCAACAATCTGACGTTAGCATCAGAAGTGGGTCTTTACGCTCATATACCACGTACGTCCCGGCTCGt                     <  1:3409283/67‑1 (MQ=255)
                                                  aacaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtat                  <  1:691065/68‑1 (MQ=255)
                                                  aacaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtat                  <  1:1483408/68‑1 (MQ=255)
                                                    caaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAgg               >  1:1500183/1‑69 (MQ=255)
                                                    caaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAgg               >  1:2801040/1‑69 (MQ=255)
                                                       tCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                             <  1:1611466/36‑1 (MQ=255)
                                                                   |                                                                  
CCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGC  >  minE/365083‑365216

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: