Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,037,993 G→A 100% T33T (ACG→ACA rsxG → predicted oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,037,9930GA81.0% 67.2 / 8.2 47T33T (ACG→ACArsxGpredicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (8/1);  new base A (25/13);  total (33/14)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACGTTGGCGCTATTTGCAGCGGGTTCAACAGGGTTAACTGCGGCCATCAACCAGATGACCAAAACGACGATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTATTTGATCAGGTGCTGCCAGCCGAACGC  >  minE/1037925‑1038056
                                                                    |                                                               
aCGTTGGCGCTATTTGCAGCGGGTTCAACAGGGTTAACTGCGGCCATCAACCAGATGACCAAAacgacg                                                                 <  1:666109/69‑1 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:61331/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:3103620/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:3254630/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:3416440/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:2852020/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:2817946/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:42754/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:2623759/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:59200/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:2370970/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:2210758/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:2124537/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:2084606/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:1963829/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:1105424/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:1326670/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:1097706/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:1649369/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:166911/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:1745119/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:1865512/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAg                                   >  1:1870541/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTATGCAACAAAAg                                   >  1:312301/1‑69 (MQ=255)
                              gggTTAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGAACAGGCAAGTCTGCAACAAAAg                                   >  1:2564582/1‑69 (MQ=255)
                              gggATAACTGCGGCCATCAACCAGATGACCAAAACGACAATTGCTGaa                                                        >  1:397746/1‑48 (MQ=39)
                                gaaaaCTGCGGCCATCAACCAGATGACCAAAACGACGATTGCTGAACAGGCCAGTCTGCAACAAAAGGc                                 <  1:2010565/66‑1 (MQ=255)
                                gTTAACTGCGGCCATCAACCAGATGACCAAAACGACGATTGCTGAACAGGCCAGTCTGCAACAAAAGGc                                 >  1:576856/1‑69 (MQ=255)
                                gTTAACTGCGGCCATCAACCAGATGACCAAAACGACGATTGCTGAACAGGCCAGTCTGCAACAAAAGGc                                 >  1:507536/1‑69 (MQ=255)
                                gTTAACTGCGGCCATCAACCAGATGACCAAAACGACGATTGCTGAACAGGCCAGTCTGCAACAAAAGGc                                 >  1:486568/1‑69 (MQ=255)
                                gTTAACTGCGGCCATCAACCAGATGACCAAAACGACGATTGCTGAACAGGCCAGTCTGCAACAAAAGGc                                 >  1:67862/1‑69 (MQ=255)
                                gTTAACTGCGGCCATCAACCAGATGACCAAAACGACGATTGCTGAACAGGCCAGTCTGCAACAAAAGGc                                 >  1:1270822/1‑69 (MQ=255)
                                gTTAACTGCGGCATCAACCAGATGACCAAAACGACGATTGCTGAACAGGCCAGTCTGCAACAAAAGGc                                 >  1:2985400/1‑68 (MQ=255)
                                      tGCGGCCATCAACCAGATGACCAAAACGACGATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTAt                            >  1:751120/1‑68 (MQ=255)
                                      tGCGGCCATCAACCAGATGACCAAAACGACGATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTAt                            >  1:247467/1‑68 (MQ=255)
                                                aaCCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTATTTGATCAGGtg                 <  1:3094241/69‑1 (MQ=255)
                                                aaCCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTATTTGATCAGGtg                 <  1:3132462/69‑1 (MQ=255)
                                                aaCCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTATTTGATCAGGtg                 <  1:305352/69‑1 (MQ=255)
                                                aaCCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTATTTGATCAGGtg                 <  1:2919736/69‑1 (MQ=255)
                                                aaCCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTATTTGATCAGGtg                 <  1:2801806/69‑1 (MQ=255)
                                                aaCCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTATTTGATCAGGtg                 <  1:2683846/69‑1 (MQ=255)
                                                aaCCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTATTTGATCAGGtg                 <  1:2597091/69‑1 (MQ=255)
                                                aaCCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTATTTGATCAGGtg                 <  1:2555934/69‑1 (MQ=255)
                                                aaCCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTATTTGATCAGGtg                 <  1:2225102/69‑1 (MQ=255)
                                                aaCCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTATTTGATCAGGtg                 <  1:1678502/69‑1 (MQ=255)
                                                aaCCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTATTTGATCAGGtg                 <  1:1394220/69‑1 (MQ=255)
                                                aaCCAGATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTATTTGATCAGGtg                 <  1:134234/69‑1 (MQ=255)
                                                     gATGACCAAAACGACAATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTATTTGATCAGGtg                 <  1:633675/64‑1 (MQ=255)
                                                               acgacgATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTATTTGATCAGGTGCTGCCAGCCGAACGc  >  1:3289623/1‑69 (MQ=255)
                                                               acgacgATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTATTTGATCAGGTGCTGCCAGCCGAACGc  >  1:2950479/1‑69 (MQ=255)
                                                                    |                                                               
ACGTTGGCGCTATTTGCAGCGGGTTCAACAGGGTTAACTGCGGCCATCAACCAGATGACCAAAACGACGATTGCTGAACAGGCCAGTCTGCAACAAAAGGCGTTATTTGATCAGGTGCTGCCAGCCGAACGC  >  minE/1037925‑1038056

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: