Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 498,918 A→C 26.6% R442R (CGT→CGG modF ← fused subunits of molybdate transporter and ATP‑binding components of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE498,9180AC26.6% 89.9 / 45.7 79R442R (CGT→CGGmodFfused subunits of molybdate transporter and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base A (42/16);  new base C (8/13);  total (50/29)
Fisher's exact test for biased strand distribution p-value = 7.97e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

TTCAGCGTGGTGCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAA  >  minE/498855‑498980
                                                               |                                                              
ttCAGCGTGGTGCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAAcggcg                                                           >  1:1143233/1‑69 (MQ=255)
ttCAGCGTGGTGCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAAcggcg                                                           >  1:837112/1‑69 (MQ=255)
ttCAGCGTGGTGCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAAcggcg                                                           >  1:737747/1‑69 (MQ=255)
ttCAGCGTGGTGCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAAcggcg                                                           >  1:1214106/1‑69 (MQ=255)
ttCAGCGTGGTGCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAAcggcg                                                           >  1:475872/1‑69 (MQ=255)
ttCAGCGTGGTGCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAAcggcg                                                           >  1:2112316/1‑69 (MQ=255)
ttCAGCGTGGTGCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAAcggcg                                                           >  1:1358367/1‑69 (MQ=255)
ttCAGCGTGGTGCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAAcggcg                                                           >  1:48595/1‑69 (MQ=255)
ttCAGCGTGGTGCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAAcggcg                                                           >  1:547212/1‑69 (MQ=255)
  cAGCGTGGTGCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAAcggcgg                                                          >  1:3283060/1‑68 (MQ=255)
  cAGCGTGGTGCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAAcggcgg                                                          >  1:1531024/1‑68 (MQ=255)
     cGTGGTGCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATa                                                       <  1:118772/68‑1 (MQ=255)
     cGTGGTGCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATa                                                       <  1:2394141/68‑1 (MQ=255)
         gtgCGAAGCAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTg                                                  >  1:1702263/1‑69 (MQ=255)
         gtgCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACGTCAACAAAACGGCGGATAAGCTg                                                  >  1:2134742/1‑69 (MQ=255)
         gtgCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTg                                                  >  1:3230751/1‑69 (MQ=255)
         gtgCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTg                                                  >  1:1584096/1‑69 (MQ=255)
         gtgCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTg                                                  >  1:1488479/1‑69 (MQ=255)
         gtgCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTg                                                  >  1:544426/1‑69 (MQ=255)
         gtgCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTg                                                  >  1:533729/1‑69 (MQ=255)
         gtgCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTg                                                  >  1:3004042/1‑69 (MQ=255)
         gtgCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTg                                                  >  1:386930/1‑69 (MQ=255)
         gtgCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTg                                                  >  1:482672/1‑69 (MQ=255)
                 caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:2238754/1‑69 (MQ=255)
                 caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:247304/1‑69 (MQ=255)
                 caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:1270720/1‑69 (MQ=255)
                 caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:123765/1‑69 (MQ=255)
                 caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:2646899/1‑69 (MQ=255)
                 caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:2748020/1‑69 (MQ=255)
                 caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:2792827/1‑69 (MQ=255)
                 caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:3315671/1‑69 (MQ=255)
                 caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:3282316/1‑69 (MQ=255)
                 caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:492365/1‑69 (MQ=255)
                            tGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGc                               >  1:1014085/1‑69 (MQ=255)
                                 ttCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTg                          <  1:410880/69‑1 (MQ=255)
                                      cTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGt                      <  1:2903530/68‑1 (MQ=255)
                                      cTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGt                      <  1:2859579/68‑1 (MQ=255)
                                      cTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGt                      <  1:1292758/68‑1 (MQ=255)
                                      cTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGt                      <  1:2624943/68‑1 (MQ=255)
                                      cTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGt                      <  1:2733672/68‑1 (MQ=255)
                                      cTTAGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGt                      <  1:1943547/68‑1 (MQ=255)
                                           cTAATCAGCACATCAAGAAACCGGCGGATAAGCTGTCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:1998065/68‑1 (MQ=255)
                                           cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:3430382/68‑1 (MQ=255)
                                           cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:1424623/68‑1 (MQ=255)
                                           cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:2651521/68‑1 (MQ=255)
                                           cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:2104729/68‑1 (MQ=255)
                                           cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:3330898/68‑1 (MQ=255)
                                           cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:1762098/68‑1 (MQ=255)
                                           cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:1610454/68‑1 (MQ=255)
                                           cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:123339/68‑1 (MQ=255)
                                           cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:1699083/68‑1 (MQ=255)
                                           cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:2224863/68‑1 (MQ=255)
                                             aaTCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:629693/66‑1 (MQ=255)
                                              aTCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCa                >  1:464077/1‑66 (MQ=255)
                                                cAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCgag           <  1:2606357/69‑1 (MQ=255)
                                                 aGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:2517250/62‑1 (MQ=255)
                                                       tCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGc                               >  1:2147243/1‑42 (MQ=255)
                                                       tCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGc    <  1:47104/69‑1 (MQ=255)
                                                       tCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGc    <  1:2508704/69‑1 (MQ=255)
                                                       tCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGc    >  1:2921784/1‑69 (MQ=255)
                                                       tCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGc    >  1:1183175/1‑69 (MQ=255)
                                                       tCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGc    <  1:1650101/69‑1 (MQ=255)
                                                       tCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGc    <  1:2502908/69‑1 (MQ=255)
                                                       tCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGc    >  1:2434371/1‑69 (MQ=255)
                                                       tCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGc    <  1:24242/69‑1 (MQ=255)
                                                       tCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGc    <  1:2418309/69‑1 (MQ=255)
                                                       tCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGc    >  1:2073091/1‑69 (MQ=255)
                                                         aacaaACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:938710/1‑69 (MQ=255)
                                                         aacaaACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:3152721/1‑69 (MQ=255)
                                                         aacaaACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:1171867/1‑69 (MQ=255)
                                                         aacaaACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:1794305/1‑69 (MQ=255)
                                                         aacaaACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:1868272/1‑69 (MQ=255)
                                                         aacaaACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:2814663/1‑69 (MQ=255)
                                                         aacaaACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:2988793/1‑69 (MQ=255)
                                                         aacaaACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:3055375/1‑69 (MQ=255)
                                                         aacaaAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:1882453/1‑69 (MQ=255)
                                                         aacaaAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:1934770/1‑69 (MQ=255)
                                                         aacaaAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:2799463/1‑69 (MQ=255)
                                                         aacaaAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCa   >  1:935521/1‑68 (MQ=255)
                                                               |                                                              
TTCAGCGTGGTGCGAAACAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAA  >  minE/498855‑498980

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: