Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 376,142 C→G 33.6% A315A (GCG→GCC fepD ← iron‑enterobactin transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE376,1420CG33.6% 106.6 / 81.0 104A315A (GCG→GCCfepDiron‑enterobactin transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (35/34);  new base G (18/17);  total (53/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

ACGTAAATCATGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGC  >  minE/376075‑376208
                                                                   |                                                                  
aCGTAAATCATGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACg                                                                   <  1:2199834/69‑1 (MQ=255)
aCGTAAATCATGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACg                                                                   <  1:781575/69‑1 (MQ=255)
     aaTCATGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACt                                                               >  1:967936/1‑68 (MQ=255)
     aaTCATGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTg                                                              >  1:294342/1‑69 (MQ=255)
     aaTCATGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTg                                                              >  1:433491/1‑69 (MQ=255)
     aaTCATGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTg                                                              >  1:362880/1‑69 (MQ=255)
     aaTCATGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTg                                                              >  1:1647561/1‑69 (MQ=255)
     aaTCATGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTg                                                              >  1:2711307/1‑69 (MQ=255)
     aaTCATGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTg                                                              >  1:113279/1‑69 (MQ=255)
     aaTCATGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTg                                                              >  1:1441435/1‑69 (MQ=255)
     aaTCATCCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTg                                                              >  1:5028/1‑69 (MQ=255)
              ccaccTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGa                                                      <  1:3269666/68‑1 (MQ=255)
              ccaccTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGa                                                      <  1:2963034/68‑1 (MQ=255)
              ccaccTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGa                                                      <  1:1171094/68‑1 (MQ=255)
               caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                    >  1:2302854/1‑69 (MQ=255)
               caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                    >  1:2586007/1‑69 (MQ=255)
               caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                    >  1:2066867/1‑69 (MQ=255)
               caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                    >  1:2817859/1‑69 (MQ=255)
               caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                    >  1:2139381/1‑69 (MQ=255)
               caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                    >  1:107525/1‑69 (MQ=255)
               caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                    >  1:291684/1‑69 (MQ=255)
               caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                    >  1:2963754/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:3383828/1‑69 (MQ=255)
               caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                    >  1:3396808/1‑69 (MQ=255)
               caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                    >  1:850105/1‑69 (MQ=255)
               caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                    >  1:205239/1‑69 (MQ=255)
               caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                    >  1:1154144/1‑69 (MQ=255)
               caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                    >  1:1886953/1‑69 (MQ=255)
               caccTCGCGTTTTACGTCGTACGAGGAAGATCAGCACTGGTGCGCCAATAAAGGCACTGACTACAGaaa                                                    >  1:1462957/1‑69 (MQ=255)
               caccTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGaaa                                                    >  1:1824301/1‑69 (MQ=255)
               caccTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGaaa                                                    >  1:263550/1‑69 (MQ=255)
               caccTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGaaa                                                    >  1:2950676/1‑69 (MQ=255)
                    cgcgTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGAc                                                            >  1:296462/1‑56 (MQ=255)
                    cgcgTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGAc                                                            >  1:3266408/1‑56 (MQ=255)
                    cgcgTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGAc                                                            >  1:2526047/1‑56 (MQ=255)
                    cgcgTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGAc                                                            >  1:1945026/1‑56 (MQ=255)
                    cgcgTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGAc                                                            >  1:3152547/1‑56 (MQ=255)
                    cgcgTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGAc                                                            >  1:1510098/1‑56 (MQ=255)
                                   aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                <  1:1788472/69‑1 (MQ=255)
                                   aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                <  1:359163/69‑1 (MQ=255)
                                   aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                <  1:3194258/69‑1 (MQ=255)
                                   aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                <  1:3056209/69‑1 (MQ=255)
                                   aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                <  1:3096427/69‑1 (MQ=255)
                                   aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                <  1:3341902/69‑1 (MQ=255)
                                   aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                <  1:576916/69‑1 (MQ=255)
                                   aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                <  1:101382/69‑1 (MQ=255)
                                   aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                <  1:2027245/69‑1 (MQ=255)
                                   aCGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                <  1:2638086/69‑1 (MQ=255)
                                    cGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                <  1:2004580/68‑1 (MQ=255)
                                    cGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGcaca                                <  1:1396376/68‑1 (MQ=255)
                                            aTCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCAcccgc                       >  1:1097145/1‑69 (MQ=255)
                                                     ggTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                  >  1:1812857/1‑65 (MQ=255)
                                                     ggTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                  >  1:49303/1‑65 (MQ=255)
                                                     ggTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                  >  1:1193724/1‑65 (MQ=255)
                                                     ggTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                  >  1:612932/1‑65 (MQ=255)
                                                     ggTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                  >  1:669083/1‑65 (MQ=255)
                                                        gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc            <  1:1214874/68‑1 (MQ=255)
                                                        gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc            <  1:486679/68‑1 (MQ=255)
                                                        gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc            <  1:3211615/68‑1 (MQ=255)
                                                        gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc            <  1:3129121/68‑1 (MQ=255)
                                                        gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc            <  1:3210638/68‑1 (MQ=255)
                                                        gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc            <  1:315084/68‑1 (MQ=255)
                                                        gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc            <  1:2349279/68‑1 (MQ=255)
                                                        gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc            <  1:154715/68‑1 (MQ=255)
                                                        gCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATc            <  1:1841068/68‑1 (MQ=255)
                                                              aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc     <  1:3048936/69‑1 (MQ=255)
                                                              aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc     <  1:892232/69‑1 (MQ=255)
                                                              aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc     <  1:856335/69‑1 (MQ=255)
                                                              aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc     <  1:668173/69‑1 (MQ=255)
                                                              aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc     <  1:2147403/69‑1 (MQ=255)
                                                              aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc     <  1:2613338/69‑1 (MQ=255)
                                                              aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc     <  1:383786/69‑1 (MQ=255)
                                                              aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc     <  1:369620/69‑1 (MQ=255)
                                                              aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc     <  1:2695835/69‑1 (MQ=255)
                                                              aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc     <  1:1301147/69‑1 (MQ=255)
                                                              aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc     <  1:1380971/69‑1 (MQ=255)
                                                              aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc     <  1:3125295/69‑1 (MQ=255)
                                                              aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc     <  1:2871248/69‑1 (MQ=255)
                                                              aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc     <  1:2899322/69‑1 (MQ=255)
                                                              aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc     <  1:3021503/69‑1 (MQ=255)
                                                              aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc     <  1:319903/69‑1 (MQ=255)
                                                              aTAAAGGCACTGACTACAGAAACGCGCAGTTCGCCGGGAACAATCACCCGCCCAATGATATCGGCAAAc     <  1:1513059/69‑1 (MQ=255)
                                                                 aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                  <  1:351661/53‑1 (MQ=255)
                                                                 aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                  <  1:3037533/53‑1 (MQ=255)
                                                                 aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                  <  1:2116871/53‑1 (MQ=255)
                                                                 aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                  <  1:2125142/53‑1 (MQ=255)
                                                                 aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                  <  1:505765/53‑1 (MQ=255)
                                                                 aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                  <  1:2428502/53‑1 (MQ=255)
                                                                 aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCgat                  <  1:3283536/53‑1 (MQ=255)
                                                                 aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc  >  1:589780/1‑69 (MQ=255)
                                                                 aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc  >  1:2457333/1‑69 (MQ=255)
                                                                 aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc  >  1:2140337/1‑69 (MQ=255)
                                                                 aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc  >  1:770802/1‑69 (MQ=255)
                                                                 aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc  >  1:2488340/1‑69 (MQ=255)
                                                                 aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc  >  1:1812190/1‑69 (MQ=255)
                                                                 aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc  >  1:1288717/1‑69 (MQ=255)
                                                                 aaCGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAAcagc  >  1:2068840/1‑69 (MQ=255)
                                                                   |                                                                  
ACGTAAATCATGCACCACCTCGCGTTTTACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGC  >  minE/376075‑376208

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: