Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 376,151 C→T 37.0% V312V (GTG→GTA fepD ← iron‑enterobactin transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE376,1510CT37.0% 60.6 / 59.7 97V312V (GTG→GTAfepDiron‑enterobactin transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (27/34);  new base T (18/18);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

CCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGC  >  minE/376089‑376216
                                                              |                                                                 
ccaccTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGa                                                              <  1:1171094/68‑1 (MQ=255)
ccaccTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGa                                                              <  1:2963034/68‑1 (MQ=255)
ccaccTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGa                                                              <  1:3269666/68‑1 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:107525/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:2302854/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:2139381/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:2586007/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:2817859/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:2066867/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:205239/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:291684/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:2963754/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:1886953/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:3011159/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:3138732/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:3251989/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:1462957/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:1154144/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:3383828/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:3396808/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                            >  1:850105/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGaaa                                                            >  1:1824301/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGaaa                                                            >  1:2950676/1‑69 (MQ=255)
 caccTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGaaa                                                            >  1:263550/1‑69 (MQ=255)
                     aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                        <  1:99644/69‑1 (MQ=255)
                     aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                        <  1:3194258/69‑1 (MQ=255)
                     aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                        <  1:101382/69‑1 (MQ=255)
                     aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                        <  1:2638086/69‑1 (MQ=255)
                     aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                        <  1:576916/69‑1 (MQ=255)
                     aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                        <  1:359163/69‑1 (MQ=255)
                     aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                        <  1:3341902/69‑1 (MQ=255)
                     aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                        <  1:3056209/69‑1 (MQ=255)
                     aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                        <  1:3096427/69‑1 (MQ=255)
                     aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                        <  1:1788472/69‑1 (MQ=255)
                     aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                        <  1:2027245/69‑1 (MQ=255)
                      cGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                        <  1:1396376/68‑1 (MQ=255)
                      cGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                        <  1:2004580/68‑1 (MQ=255)
                              aTCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCAcccgc                               >  1:1097145/1‑69 (MQ=255)
                                       ggTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                          >  1:669083/1‑65 (MQ=255)
                                       ggTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                          >  1:612932/1‑65 (MQ=255)
                                       ggTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                          >  1:49303/1‑65 (MQ=255)
                                       ggTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                          >  1:1193724/1‑65 (MQ=255)
                                       ggTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                          >  1:1812857/1‑65 (MQ=255)
                                          gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc                    <  1:1214874/68‑1 (MQ=255)
                                          gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc                    <  1:3211615/68‑1 (MQ=255)
                                          gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc                    <  1:3210638/68‑1 (MQ=255)
                                          gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc                    <  1:3129121/68‑1 (MQ=255)
                                          gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc                    <  1:154715/68‑1 (MQ=255)
                                          gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc                    <  1:2349279/68‑1 (MQ=255)
                                          gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc                    <  1:486679/68‑1 (MQ=255)
                                          gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc                    <  1:315084/68‑1 (MQ=255)
                                          gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc                    <  1:1841068/68‑1 (MQ=255)
                                                aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:1301147/69‑1 (MQ=255)
                                                aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:668173/69‑1 (MQ=255)
                                                aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:3125295/69‑1 (MQ=255)
                                                aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:1380971/69‑1 (MQ=255)
                                                aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:2147403/69‑1 (MQ=255)
                                                aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:319903/69‑1 (MQ=255)
                                                aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:1513059/69‑1 (MQ=255)
                                                aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:892232/69‑1 (MQ=255)
                                                aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:3048936/69‑1 (MQ=255)
                                                aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:369620/69‑1 (MQ=255)
                                                aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:3021503/69‑1 (MQ=255)
                                                aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:856335/69‑1 (MQ=255)
                                                aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:383786/69‑1 (MQ=255)
                                                aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:2613338/69‑1 (MQ=255)
                                                aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:2899322/69‑1 (MQ=255)
                                                aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:2871248/69‑1 (MQ=255)
                                                aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:2695835/69‑1 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                          <  1:351661/53‑1 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                          <  1:505765/53‑1 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                          <  1:2116871/53‑1 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                          <  1:3283536/53‑1 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                          <  1:2428502/53‑1 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                          <  1:3037533/53‑1 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                          <  1:2125142/53‑1 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc          >  1:1993307/1‑69 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc          >  1:2140337/1‑69 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc          >  1:770802/1‑69 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc          >  1:589780/1‑69 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc          >  1:2068840/1‑69 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc          >  1:2872964/1‑69 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc          >  1:3281181/1‑69 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc          >  1:1288717/1‑69 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc          >  1:2488340/1‑69 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc          >  1:1812190/1‑69 (MQ=255)
                                                   aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc          >  1:2457333/1‑69 (MQ=255)
                                                         cTGACTACAGAAACTCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc             <  1:212267/60‑1 (MQ=255)
                                                          tGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAgcaggg   <  1:1912680/69‑1 (MQ=255)
                                                          tGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAgcaggg   <  1:793666/69‑1 (MQ=255)
                                                           gACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAgcagggc  >  1:3162614/1‑69 (MQ=255)
                                                           gACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAgcagggc  >  1:592438/1‑69 (MQ=255)
                                                           gACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAgcagggc  >  1:1343097/1‑69 (MQ=255)
                                                           gACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAgcagggc  >  1:3289831/1‑69 (MQ=255)
                                                           gACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAgcagggc  >  1:2968442/1‑69 (MQ=255)
                                                           gACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAgcagggc  >  1:2399767/1‑69 (MQ=255)
                                                           gACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAgcagggc  >  1:2276426/1‑69 (MQ=255)
                                                              |                                                                 
CCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGC  >  minE/376089‑376216

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: