Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,631,410 T→C 22.2% E128E (GAA→GAG aegA ← fused predicted FeS binding subunit and predicted NAD/FAD‑binding subunit of oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,631,4100TC22.2% 64.1 / 15.5 36E128E (GAA→GAGaegAfused predicted FeS binding subunit and predicted NAD/FAD‑binding subunit of oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (16/12);  new base C (5/3);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CGGAGAGTGCGACGTCAGTGACCAGTTGCAGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCATTTATGCGCCGTGGCTTTTACTTTTCCTGCCGCGACGGGTGTCAGGA  >  minE/1631340‑1631480
                                                                      |                                                                      
cGGAGAGTGCGACGTCAGTGACCAGTTGCAGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGttt                                                                        <  1:3554292/71‑1 (MQ=255)
      gTGCGACGTCAGTGACCAGTTGCAGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGttt                                                                        <  1:423475/65‑1 (MQ=255)
      gTGCGACGTCAGTGACCAGTTGCAGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGttt                                                                        <  1:205845/65‑1 (MQ=255)
         cGACGTCAGTGACCAGTTGCAGCGCGTCCGGCGGGCAATTCTCAACACAGGCAGGACCGTTCTCGCGcccc                                                               >  1:3671550/1‑71 (MQ=255)
            cGTCAGTGACCAGTTGCAGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGTTCTCGCGCCCCGca                                                            >  1:2025305/1‑71 (MQ=255)
             gTCAGTGACCAGTTGCAGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGTTTTCGCGCCCCGca                                                            >  1:2251832/1‑70 (MQ=255)
                aGTGACCAGTTGCAGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGTTTTCGCGCCCCGcaca                                                          >  1:2908388/1‑69 (MQ=255)
                 gTGACCAGTTGCAGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGTTTTCGCGCCCCGCACAAAg                                                       <  1:1187495/71‑1 (MQ=255)
                  tGACCAGTTGCAGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGTTTTCGCGCCCCGCACAAAgg                                                      <  1:858635/71‑1 (MQ=255)
                     ccAGTTGCAGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCg                                                   >  1:3164279/1‑71 (MQ=255)
                      cAGTTGCAGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGTTTTc                                                                      >  1:3731809/1‑51 (MQ=255)
                      cAGTTGCAGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGTTTTc                                                                      >  1:2851924/1‑51 (MQ=255)
                        gTTGCAGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCAt                                                >  1:31327/1‑71 (MQ=255)
                           gCAGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCAttt                                              <  1:2029058/70‑1 (MQ=255)
                           gCAGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCAttt                                              <  1:2272105/70‑1 (MQ=255)
                           gCAGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGTTCTCGCGCCCCGCACAAAGGTCGCAttt                                              <  1:2868656/70‑1 (MQ=255)
                             aGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCATTTATg                                           <  1:586090/71‑1 (MQ=255)
                                 cgTCCGCCGGGCAATTCTCAACACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCATTTATGCGcc                                       <  1:142837/71‑1 (MQ=255)
                                  gTCCGCCGGGCAATTCTCAACACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCATTTATGCGCCg                                      >  1:3247097/1‑71 (MQ=255)
                                       ccgGGCAATTCTCAACACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCATTTATGCGCCGTGGCt                                 <  1:783107/71‑1 (MQ=255)
                                        cgGGCAATTCTCAACACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCAttt                                              >  1:3663783/1‑57 (MQ=255)
                                               ttCTCAACACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCATTTATGCGCCGTGGCTTTTACttt                         >  1:2739374/1‑71 (MQ=255)
                                               ttCTCAACACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCATTTATGCGCCGTGGCTTTTACttt                         >  1:639888/1‑71 (MQ=255)
                                               ttCTCAACACAGGCAGGACCGTTCTCGCGCCCCGCACAAAGGTCGCATTTATGCGCCGTGGCTTTTACttt                         >  1:3248370/1‑71 (MQ=255)
                                                tctcAACACAGGCAGGACCGTTTTCGCGCCCCGCACaaa                                                        >  1:2381198/1‑39 (MQ=255)
                                                tctcAACACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCATTTATGCGCCGTGGCTTTTACtttt                        >  1:1119072/1‑71 (MQ=255)
                                                 ctcAACACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCATTTATgc                                          >  1:1263081/1‑52 (MQ=255)
                                                 ctcAACACAGGCAGGACCGTTCTCGCGCCCCGCACAAAGGTCGCAt                                                >  1:1753559/1‑46 (MQ=255)
                                                    aaCACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCa                                                 >  1:1076904/1‑42 (MQ=255)
                                                    aaCACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCATTTATgc                                          >  1:3367759/1‑49 (MQ=255)
                                                    aaCACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCATTTATgc                                          >  1:2363017/1‑49 (MQ=255)
                                                    aaCACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCATTTATGCGCCGTGGCTTTTAc                            <  1:3918098/63‑1 (MQ=255)
                                                    aaCACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCATTTATGCGCCGTGGCTTTTACTTTTcc                      <  1:842149/69‑1 (MQ=255)
                                                    aaCACAGGCAGGACCGTTCTCGCGCCCCGCACAAAGGTCGCATTTATGCGCCGTGGCTTTTACTTTTcc                      <  1:1419129/69‑1 (MQ=255)
                                                      cacaGGCAGGACCGTTTTCGCGCCCCGCTCAAAGGTCGCATTTATGCGCCGTGGCTTTTa                             <  1:2675030/60‑1 (MQ=255)
                                                      cacaGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCATTTATGCGCCGTGGCTTTTa                             <  1:3700739/60‑1 (MQ=255)
                                                        caggcaggACCGTTCTCGCGCCCCGCACAAAGGTCGCATTTAt                                            <  1:552845/43‑1 (MQ=255)
                                                          ggcaggACCGTTTTCGCGCCCCGCACAAAGGTCGCATTTATGCGCCGTGGCTTTTACtt                          <  1:1415623/59‑1 (MQ=255)
                                                             aggACCGTTCTCGCGCCCCGCACAAAGGTCGCATTTATGCGCCGTGGCTTTTACTTTTCCt                     >  1:1692980/1‑61 (MQ=255)
                                                                    ttttCGCGCCCCGCACAAAGGTCGCATTTATGCGCCGTGGCt                                 >  1:3938523/1‑42 (MQ=255)
                                                                    ttttCGCGCCCCGCACAAAGGTCGCATTTATGCGCCGTGGCt                                 >  1:499777/1‑42 (MQ=255)
                                                                    ttttCGCGCCCCGCACAAAGGTCGCATTTATGCGCCGTGGCt                                 >  1:79932/1‑42 (MQ=255)
                                                                      ttCGCGCCCCGCACAAAGGTCGCATTTATGCGCCGTGGCTTTTACTTTTCCTGCCGCGACGGGTGTCAGGa  >  1:1537793/1‑71 (MQ=255)
                                                                      |                                                                      
CGGAGAGTGCGACGTCAGTGACCAGTTGCAGCGCGTCCGCCGGGCAATTCTCAACACAGGCAGGACCGTTTTCGCGCCCCGCACAAAGGTCGCATTTATGCGCCGTGGCTTTTACTTTTCCTGCCGCGACGGGTGTCAGGA  >  minE/1631340‑1631480

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: