Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,951,544 C→G 75.0% V14V (GTG→GTC visC ← predicted oxidoreductase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,951,5440CG75.0% 72.7 / 22.7 48V14V (GTG→GTCvisCpredicted oxidoreductase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base C (8/4);  new base G (20/16);  total (28/20)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CGCGCTGCTCCAGTACGGCAACGCGTAAGCCGCTCCCCTGTAAGCCACAGGCAACCGCCAGCCCCACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTATTCCT  >  minE/1951480‑1951596
                                                                |                                                    
cgcgCTGCTCCAGTACGGCAACGCGTAAGCCGCTCCCCTGTAAGCCACAGGCAACCGCCAGCCccac                                                    <  1:2397472/67‑1 (MQ=255)
     ttctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgccg                                             >  1:3425446/3‑69 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgccg                                             >  1:1962268/1‑69 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgccg                                             >  1:2504871/1‑69 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgccg                                             >  1:2479268/1‑69 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgccg                                             >  1:2773700/1‑69 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgccg                                             >  1:2873471/1‑69 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgccg                                             >  1:2158274/1‑69 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgccg                                             >  1:2123390/1‑69 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgccg                                             >  1:1970088/1‑69 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgccg                                             >  1:3233705/1‑69 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgccg                                             >  1:1385907/1‑69 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgccg                                             >  1:459086/1‑69 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgccg                                             >  1:510882/1‑69 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgccg                                             >  1:880396/1‑69 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgccg                                             >  1:1129/1‑69 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgccg                                             >  1:1010215/1‑69 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgcc                                              >  1:1930893/1‑68 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgcc                                              >  1:2610397/1‑68 (MQ=255)
     tgctCCAGTACGGCAACGCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgca                                              >  1:1926565/1‑67 (MQ=255)
     tgctCCAGTACGGCAACCCGTAAGCCGCTACCCTGTAAGCCGCAGGCAACCGCCAGCCCGACCATgccg                                             >  1:1283378/1‑69 (MQ=255)
              aCGGCAACGCGTAAGCCGCTCCCCTGTAAGCCACAGGCAACCGCCAGCCCCACCATGCCGCCGCcaaca                                    >  1:207576/1‑69 (MQ=255)
              aCGGCAACGCGTAAGCCGCTCCCCTGTAAGCCACAGGCAACCGCCAGCCCCACCATGCCGCCGCcaaca                                    >  1:1310379/1‑69 (MQ=255)
              aCGGCAACGCGTAAGCCGCTCCCCTGTAAGCCACAGGCAACCGCCAGCCCCACCATGCCGCCGCcaaca                                    >  1:3211730/1‑69 (MQ=255)
              aCGGCAACGCGTAAGCCGCTCCCCTGTAAGCCACAGGCAACCGCCAGCCCCACCATGCCGCCGCcaaca                                    >  1:2800938/1‑69 (MQ=255)
                              cGCTCCCCTGTAAGCCACAGGCAACCGCCAGCCCCACCATGCCGCCGCCAACAATGGCTACATCAacac                    >  1:1069926/1‑69 (MQ=255)
                                        tAAGCCACAGGCAACCGCCAGCCCCACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTg           >  1:797785/1‑68 (MQ=255)
                                        tAAGCCACAGGCAACCGCCAGCCCCACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTg           >  1:2280320/1‑68 (MQ=255)
                                            ccgcaGGCAACCGCCAGCCCGACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:2481321/66‑1 (MQ=255)
                                            ccgcaGGCAACCGCCAGCCCGACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:862679/66‑1 (MQ=255)
                                            ccgcaGGCAACCGCCAGCCCGACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:544717/66‑1 (MQ=255)
                                            ccgcaGGCAACCGCCAGCCCGACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:1164614/66‑1 (MQ=255)
                                            ccgcaGGCAACCGCCAGCCCGACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:481040/66‑1 (MQ=255)
                                            ccgcaGGCAACCGCCAGCCCGACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:1415614/66‑1 (MQ=255)
                                            ccgcaGGCAACCGCCAGCCCGACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:3356395/66‑1 (MQ=255)
                                            ccgcaGGCAACCGCCAGCCCGACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:330329/66‑1 (MQ=255)
                                            ccgcaGGCAACCGCCAGCCCGACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:3289891/66‑1 (MQ=255)
                                            ccgcaGGCAACCGCCAGCCCGACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:323384/66‑1 (MQ=255)
                                            ccgcaGGCAACCGCCAGCCCGACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:1762353/66‑1 (MQ=255)
                                            ccgcaGGCAACCGCCAGCCCGACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:2176457/66‑1 (MQ=255)
                                            ccgcaGGCAACCGCCAGCCCGACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:2311802/66‑1 (MQ=255)
                                            ccgcaGGCAACCGCCAGCCCGACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:2472760/66‑1 (MQ=255)
                                            ccgcaGGCAACCGCCAGCCCGACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:2539606/66‑1 (MQ=255)
                                            ccgcaGGCAACCGCCAGCCCGACCATGCAGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:1163205/66‑1 (MQ=255)
                                            ccACAGGCAACCGCCAGCCCCACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:316210/69‑1 (MQ=255)
                                            ccACAGGCAACCGCCAGCCCCACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAt      <  1:3049197/69‑1 (MQ=255)
                                             cacaGGCAACCGCCAGCCCCACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAtt     <  1:1291294/69‑1 (MQ=255)
                                              acaGGCAACCGCCAGCCCCACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTAtt     <  1:3180407/68‑1 (MQ=255)
                                                 ggCAACCGCCAGCCCCACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTATTCCt  >  1:3319133/1‑68 (MQ=255)
                                                                |                                                    
CGCGCTGCTCCAGTACGGCAACGCGTAAGCCGCTCCCCTGTAAGCCACAGGCAACCGCCAGCCCCACCATGCCGCCGCCAACAATGGCTACATCAACACTTTGCATTGTTTATTCCT  >  minE/1951480‑1951596

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: