Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,679,505 C→T 24.3% V395V (GTG→GTA feoB ← fused ferrous iron transporter, protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,679,5050CT24.3% 121.1 / 31.1 95V395V (GTG→GTAfeoBfused ferrous iron transporter, protein B
Reads supporting (aligned to +/- strand):  ref base C (46/25);  new base T (14/9);  total (61/34)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CAAGCGTACGTGCACCCATTACCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCA  >  minE/2679447‑2679572
                                                          |                                                                   
cAAGCGTACGTGCGCCCATTACCGACGGTACGTTACAACCGAAACCAACGATCAGCGGTACAAAGGAtt                                                           >  1:1973684/1‑69 (MQ=255)
cAAGCGTACGTGCGCCCATTACCGACGGTACGTTACAACCGAAACCAACGATCAGCGGTACAAAGGAtt                                                           >  1:968406/1‑69 (MQ=255)
cAAGCGTACGTGCGCCCATTACCGACGGTACGTTACAACCGAAACCAACGATCAGCGGTACAAAGGAtt                                                           >  1:882639/1‑69 (MQ=255)
cAAGCGTACGTGCGCCCATTACCGACGGTACGTTACAACCGAAACCAACGATCAGCGGTACAAAGGAtt                                                           >  1:659985/1‑69 (MQ=255)
cAAGCGTACGTGCGCCCATTACCGACGGTACGTTACAACCGAAACCAACGATCAGCGGTACAAAGGAtt                                                           >  1:3380305/1‑69 (MQ=255)
cAAGCGTACGTGCGCCCATTACCGACGGTACGTTACAACCGAAACCAACGATCAGCGGTACAAAGGAtt                                                           >  1:3119520/1‑69 (MQ=255)
cAAGCGTACGTGCGCCCATTACCGACGGTACGTTACAACCGAAACCAACGATCAGCGGTACAAAGGAtt                                                           >  1:3100704/1‑69 (MQ=255)
cAAGCGTACGTGCGCCCATTACCGACGGTACGTTACAACCGAAACCAACGATCAGCGGTACAAAGGAtt                                                           >  1:2007895/1‑69 (MQ=255)
cAAGCGTACGTGCGCCCATTACCGACGGTACGTTACAACCGAAACCAACGATCAGCGGTACAAAGGAtt                                                           >  1:19146/1‑69 (MQ=255)
cAAGCGTACGTGCGCCCATTACCGACGGTACGTTACAACCGAAACCAACGATCAGCGGTACAAAGGAtt                                                           >  1:1832489/1‑69 (MQ=255)
cAAGCGTACGTGCGCCCATTACCGACGGTACGTTACAACCGAAACCAACGATCAGCGGTACAAAGGAtt                                                           >  1:1804878/1‑69 (MQ=255)
cAAGCGTACGTGCGCCCATTACCGACGGTACGTTACAACCGAAACCAACGATCAGCGGTACAAAGGAtt                                                           >  1:1379463/1‑69 (MQ=255)
cAAGCGTACGTGCGCCCATTACCGACGGTACGTTACAACCGAAACCAACGATCAGCGGTACAAAGGAt                                                            >  1:1595018/1‑68 (MQ=255)
cAAGCGTACGTGCGCCCATTACCGACGGTACGTTAAAACCGAAACCAACGATCAGCGGTACAAAGGAtt                                                           >  1:1961603/1‑69 (MQ=255)
cAAGCGTACGTGCACCCATTACCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGAtt                                                           >  1:1366718/1‑69 (MQ=255)
        cGTGCACCCATTACCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCg                                                     >  1:2143856/1‑67 (MQ=255)
        cGTGCACCCATTACCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCg                                                     >  1:3101247/1‑67 (MQ=255)
                 aTTACCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAgc                                          <  1:519535/69‑1 (MQ=255)
                 aTTACCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAgc                                          <  1:2051312/69‑1 (MQ=255)
                 aTTACCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAgc                                          <  1:728862/69‑1 (MQ=255)
                 aTTACCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAgc                                          <  1:736528/69‑1 (MQ=255)
                 aTTACCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAgc                                          <  1:1169645/69‑1 (MQ=255)
                 aTTACCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAgc                                          <  1:459755/69‑1 (MQ=255)
                 aTTACCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAgc                                          <  1:391181/69‑1 (MQ=255)
                 aTTACCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAgc                                          <  1:1489397/69‑1 (MQ=255)
                    aCCGACGGTTCGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              <  1:94198/46‑1 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              <  1:1855014/46‑1 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              <  1:1115275/46‑1 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              >  1:2091201/1‑46 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              >  1:2334366/1‑46 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              >  1:1321962/1‑46 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              >  1:2714540/1‑46 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              >  1:294670/1‑46 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              >  1:1010248/1‑46 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              >  1:1289124/1‑46 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              >  1:1945242/1‑46 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              <  1:3352941/46‑1 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              >  1:841719/1‑46 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              >  1:1581779/1‑46 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              <  1:185107/46‑1 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              >  1:1275237/1‑46 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              <  1:3089009/46‑1 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              <  1:1788892/46‑1 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              >  1:1448041/1‑46 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              >  1:1457602/1‑46 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              <  1:1662887/46‑1 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              >  1:156903/1‑46 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAgg                                                              <  1:1579503/46‑1 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGa                                                             <  1:2999731/47‑1 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGa                                                             <  1:1266023/47‑1 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGa                                                             <  1:1637202/47‑1 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGcc                                       >  1:3276023/1‑69 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGcc                                       >  1:3094540/1‑69 (MQ=255)
                    aCCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGcc                                       >  1:358773/1‑69 (MQ=255)
                     ccGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAgc                                          <  1:659597/65‑1 (MQ=255)
                       gACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTg                                     <  1:1382578/68‑1 (MQ=255)
                       gACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTg                                     <  1:1590566/68‑1 (MQ=255)
                         cGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGc                                    >  1:1432522/1‑67 (MQ=255)
                         cGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGc                                    >  1:2438178/1‑67 (MQ=255)
                            tACGTTACAACCGAAACCGACGATCAGCGGCACAAAGTATTTCCCCGGCAAGCCCAGCGCCTGCATCAg                               <  1:306055/69‑1 (MQ=255)
                             aCGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAg                               >  1:543263/1‑68 (MQ=255)
                             aCGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGa                              >  1:2783462/1‑69 (MQ=255)
                             aCGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGa                              >  1:1013834/1‑69 (MQ=255)
                             aCGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGa                              >  1:227077/1‑69 (MQ=255)
                             aCGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGa                              >  1:110275/1‑69 (MQ=255)
                             aCGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGa                              >  1:1648824/1‑69 (MQ=255)
                             aCGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGa                              >  1:1664385/1‑69 (MQ=255)
                             aCGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGa                              >  1:930303/1‑69 (MQ=255)
                             aCGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGa                              >  1:1020590/1‑69 (MQ=255)
                             aCGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGa                              >  1:2125815/1‑69 (MQ=255)
                             aCGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGa                              >  1:1904316/1‑69 (MQ=255)
                                      ccgaaaccgaCGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCAt                      >  1:1250114/1‑68 (MQ=255)
                                           accgaCGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACaaa                >  1:1282202/1‑69 (MQ=255)
                                           accgaCGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACaaa                >  1:2850590/1‑69 (MQ=255)
                                              gacgaTCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGAc                             >  1:1768203/1‑53 (MQ=255)
                                              gacgaTCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCAt                      >  1:685695/1‑60 (MQ=255)
                                              gacgaTCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAAcgc             >  1:654545/1‑69 (MQ=255)
                                              gacgaTCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAAcgc             >  1:1320265/1‑69 (MQ=255)
                                              gacgaTCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAAcgc             >  1:1705285/1‑69 (MQ=255)
                                               acgaTCAGCGGTACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAAcgcc            <  1:1618622/69‑1 (MQ=255)
                                               acgaTCAGCGGTACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAAcgcc            <  1:50228/69‑1 (MQ=255)
                                               acgaTCAGCGGTACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAAcgcc            <  1:1218797/69‑1 (MQ=255)
                                               acgaTCAGCGGTACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAAcgcc            <  1:3012104/69‑1 (MQ=255)
                                               acgaTCAGCGGTACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAAcgcc            <  1:2270436/69‑1 (MQ=255)
                                               acgaTCAGCGGTACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAAcgcc            <  1:2103976/69‑1 (MQ=255)
                                               acgaTCAGCGGTACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAAcgcc            <  1:2060230/69‑1 (MQ=255)
                                               acgaTCAGCGGTACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAAcgcc            <  1:1719121/69‑1 (MQ=255)
                                               acgaTCAGCGGTACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAAcgcc            <  1:1876180/69‑1 (MQ=255)
                                                    cAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCAcg       >  1:131912/1‑69 (MQ=255)
                                                    cAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCAcg       >  1:320863/1‑69 (MQ=255)
                                                    cAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCAcg       >  1:1989196/1‑69 (MQ=255)
                                                    cAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCAcg       >  1:3123770/1‑69 (MQ=255)
                                                    cAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCAc        >  1:1941940/1‑68 (MQ=255)
                                                    cAGCGGAACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCAcg       >  1:437678/1‑69 (MQ=255)
                                                         gCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCa  <  1:1867274/69‑1 (MQ=255)
                                                          |                                                                   
CAAGCGTACGTGCACCCATTACCGACGGTACGTTACAACCGAAACCGACGATCAGCGGCACAAAGGATTTCCCCGGCAAGCCCAGCGCCTGCATCAGACGGTCCATCACAAACGCCGCACGCGCCA  >  minE/2679447‑2679572

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: