Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 365,168 C→T 26.4% V742I (GTA→ATA) ‡ fepA ← iron‑enterobactin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE365,1680CT26.4% 116.9 / 36.9 98V742I (GTA→ATA) ‡fepAiron‑enterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (22/50);  new base T (18/8);  total (40/58)
Fisher's exact test for biased strand distribution p-value = 9.45e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

ATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGC  >  minE/365105‑365234
                                                               |                                                                  
aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                               <  1:204085/69‑1 (MQ=255)
aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                               <  1:1490904/69‑1 (MQ=255)
aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                               <  1:2323848/69‑1 (MQ=255)
aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                               <  1:1507019/69‑1 (MQ=255)
aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                               <  1:2084405/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:1756462/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:1167182/69‑1 (MQ=255)
aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                               <  1:2976494/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:2281444/1‑69 (MQ=25)
          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                                     >  1:395254/1‑69 (MQ=25)
          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                                     >  1:2172339/1‑69 (MQ=25)
          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                                     >  1:2100706/1‑69 (MQ=25)
          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                                     >  1:209894/1‑69 (MQ=25)
          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                                     >  1:114176/1‑69 (MQ=25)
          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                                     >  1:57234/1‑69 (MQ=25)
          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                                     >  1:3284071/1‑69 (MQ=25)
          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                                     >  1:2608967/1‑69 (MQ=25)
          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                                     >  1:3095133/1‑69 (MQ=25)
          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                                     >  1:195893/1‑69 (MQ=25)
          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                                     >  1:1881992/1‑69 (MQ=25)
          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                                     >  1:1256798/1‑69 (MQ=25)
          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                                     >  1:1260868/1‑69 (MQ=25)
          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                                     >  1:862452/1‑69 (MQ=25)
          aaCTACGCGCACTTTGTCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCacgt                                                     >  1:159115/1‑69 (MQ=25)
          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                               >  1:1863334/1‑59 (MQ=255)
          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                               >  1:332060/1‑59 (MQ=255)
          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                               >  1:3188574/1‑59 (MQ=255)
          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                               >  1:536719/1‑59 (MQ=255)
          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt                                                     >  1:1616896/1‑69 (MQ=255)
          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt                                                     >  1:3278883/1‑69 (MQ=255)
          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt                                                     >  1:800278/1‑69 (MQ=255)
          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt                                                     >  1:2268017/1‑69 (MQ=255)
          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt                                                     >  1:2620567/1‑69 (MQ=255)
          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgt                                                     >  1:2417860/1‑69 (MQ=255)
          aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacg                                                      >  1:2517876/1‑68 (MQ=255)
                gcgcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTcc                                               >  1:1749554/1‑69 (MQ=255)
                gcgcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTcc                                               >  1:2136116/1‑69 (MQ=255)
                gcgcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTcc                                               >  1:1173313/1‑69 (MQ=255)
                gcgcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTcc                                               >  1:2325093/1‑69 (MQ=255)
                gcgcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTcc                                               >  1:11437/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:1392747/69‑1 (MQ=255)
                   cACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                                            <  1:1236411/69‑1 (MQ=255)
                   cACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                                            <  1:1584762/69‑1 (MQ=255)
                   cACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                                            <  1:894142/69‑1 (MQ=255)
                    aCTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                               <  1:892423/49‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCACGTACGTCCCGGCTCGtt                                      <  1:1886658/68‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCACGTACGTCCCGGCTCGtt                                      <  1:1401423/68‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCACGTACGTCCCGGCTCGtt                                      <  1:882900/68‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCACGTACGTCCCGGCTCGtt                                      <  1:809667/68‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCACGTACGTCCCGGCTCGtt                                      <  1:2083893/68‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCACGTACGTCCCGGCTCGtt                                      <  1:1859740/68‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCACGTACGTCCCGGCTCGtt                                      <  1:68819/68‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCTTCAGAAGTGAGTATTAATGCTCATATACCACGTACGTCCCGGCTCGtt                                      <  1:762665/68‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:73610/67‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:3051948/67‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:2402551/67‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:2935291/67‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:1181239/67‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:1195690/67‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:1296697/67‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:1930080/67‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:1978392/67‑1 (MQ=255)
                          tCAACAATCTGACGTTAGCATCAGAAGTGGGTCTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:3409283/67‑1 (MQ=255)
                            aacaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtat                                    <  1:1483408/68‑1 (MQ=255)
                            aacaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtat                                    <  1:691065/68‑1 (MQ=255)
                              caaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAgg                                 >  1:2801040/1‑69 (MQ=255)
                              caaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAgg                                 >  1:1500183/1‑69 (MQ=255)
                                 tCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                               <  1:1611466/36‑1 (MQ=255)
                                            cATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGc                    <  1:969651/68‑1 (MQ=255)
                                                       ggTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:2090025/38‑1 (MQ=255)
                                                       ggTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:2050489/38‑1 (MQ=255)
                                                       ggTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:2027965/38‑1 (MQ=255)
                                                       ggTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:678745/38‑1 (MQ=255)
                                                       ggTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:1017309/38‑1 (MQ=255)
                                                       ggTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:3343355/38‑1 (MQ=255)
                                                       ggTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:1911197/38‑1 (MQ=255)
                                                       ggTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:1867394/38‑1 (MQ=255)
                                                       ggTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:2499216/38‑1 (MQ=255)
                                                       ggTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:2567703/38‑1 (MQ=255)
                                                       ggTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:2732662/38‑1 (MQ=255)
                                                       ggTGTTTACGCTCATATACCACGTACGTCCCGGCTCGt                                       <  1:2910499/38‑1 (MQ=255)
                                                          gtTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGtt     >  1:362947/1‑69 (MQ=255)
                                                          gtTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGtt     >  1:3283090/1‑69 (MQ=255)
                                                          gtTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGtt     >  1:2588711/1‑69 (MQ=255)
                                                             tACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGc  >  1:2974628/1‑69 (MQ=255)
                                                               |                                                                  
ATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGC  >  minE/365105‑365234

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: