Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 498,939 C→G 37.5% L435L (CTG→CTC 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,9390CG37.5% 37.7 / 48.1 56L435L (CTG→CTCmodFfused subunits of molybdate transporter and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base C (21/14);  new base G (8/13);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

CAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAA  >  minE/498872‑498980
                                                                   |                                         
caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:2238754/1‑69 (MQ=255)
caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:3282316/1‑69 (MQ=255)
caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:3315671/1‑69 (MQ=255)
caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:2792827/1‑69 (MQ=255)
caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:2748020/1‑69 (MQ=255)
caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:2646899/1‑69 (MQ=255)
caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:247304/1‑69 (MQ=255)
caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:492365/1‑69 (MQ=255)
caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:123765/1‑69 (MQ=255)
caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCa                                          >  1:1270720/1‑69 (MQ=255)
           tGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGc                               >  1:1014085/1‑69 (MQ=255)
                ttCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTg                          <  1:410880/69‑1 (MQ=255)
                     cTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGt                      <  1:1292758/68‑1 (MQ=255)
                     cTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGt                      <  1:2903530/68‑1 (MQ=255)
                     cTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGt                      <  1:2859579/68‑1 (MQ=255)
                     cTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGt                      <  1:2733672/68‑1 (MQ=255)
                     cTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGt                      <  1:2624943/68‑1 (MQ=255)
                     cTTAGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGt                      <  1:1943547/68‑1 (MQ=255)
                          cTAATCAGCACATCAAGAAACCGGCGGATAAGCTGTCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:1998065/68‑1 (MQ=255)
                          cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:123339/68‑1 (MQ=255)
                          cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:1424623/68‑1 (MQ=255)
                          cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:2651521/68‑1 (MQ=255)
                          cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:3430382/68‑1 (MQ=255)
                          cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:1762098/68‑1 (MQ=255)
                          cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:1610454/68‑1 (MQ=255)
                          cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:1699083/68‑1 (MQ=255)
                          cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:3330898/68‑1 (MQ=255)
                          cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:2224863/68‑1 (MQ=255)
                          cTAATCAGCACATCAACAAACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTc                 <  1:2104729/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:2921784/1‑69 (MQ=255)
                                      tCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGc    >  1:1183175/1‑69 (MQ=255)
                                      tCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGc    <  1:1650101/69‑1 (MQ=255)
                                      tCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGc    <  1:47104/69‑1 (MQ=255)
                                      tCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGc    <  1:2508704/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:3055375/1‑69 (MQ=255)
                                        aacaaACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:2988793/1‑69 (MQ=255)
                                        aacaaACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:2814663/1‑69 (MQ=255)
                                        aacaaACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:1868272/1‑69 (MQ=255)
                                        aacaaACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:1794305/1‑69 (MQ=255)
                                        aacaaACCGGCGGATAAGCTGGCGATTGAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:1171867/1‑69 (MQ=255)
                                        aacaaAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:2799463/1‑69 (MQ=255)
                                        aacaaAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:1934770/1‑69 (MQ=255)
                                        aacaaAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCaa  >  1:1882453/1‑69 (MQ=255)
                                        aacaaAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCa   >  1:935521/1‑68 (MQ=255)
                                                                   |                                         
CAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAA  >  minE/498872‑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: