Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 521,525 G→T 65.3% V515V (GTG→GTT uvrB → excinulease of nucleotide excision repair, DNA damage recognition component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE521,5250GT65.3% 32.3 / 38.4 49V515V (GTG→GTTuvrBexcinulease of nucleotide excision repair, DNA damage recognition component
Reads supporting (aligned to +/- strand):  ref base G (15/2);  new base T (20/12);  total (35/14)
Fisher's exact test for biased strand distribution p-value = 9.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCA  >  minE/521486‑521590
                                       |                                                                 
gATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGa                                         >  1:1314997/1‑66 (MQ=255)
gATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGa                                         >  1:820100/1‑66 (MQ=255)
gATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGa                                         >  1:2086640/1‑66 (MQ=255)
 aTCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACaa                                      >  1:1129919/1‑68 (MQ=255)
 aTCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACaa                                      >  1:84162/1‑68 (MQ=255)
 aTCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACaa                                      >  1:2915933/1‑68 (MQ=255)
 aTCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACaa                                      >  1:3377183/1‑68 (MQ=255)
 aTCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACaa                                      >  1:659217/1‑68 (MQ=255)
        tACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCTTGTGGCGATCCTCGACGCTGACAAAGAAGGCt                              <  1:2849758/69‑1 (MQ=255)
        tACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCt                              <  1:2553967/69‑1 (MQ=255)
               cgAAGGTCTGGATATGCCGGAAGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGc                        <  1:1896159/68‑1 (MQ=255)
               cgAAGGTCTGGATATGCCGGAAGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGc                        <  1:3429337/68‑1 (MQ=255)
               cgAAGGTCTGGATATGCCGGAAGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGc                        <  1:2394118/68‑1 (MQ=255)
               cgAAGGTCTGGATATGCCGGAAGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGc                        <  1:535281/68‑1 (MQ=255)
               cgAAGGTCTGGATATGCCGGAAGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGc                        <  1:2159006/68‑1 (MQ=255)
               cgAAGGTCTGGATATGCCGGAAGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGc                        <  1:664709/68‑1 (MQ=255)
               cgAAGGTCTGGATATGCCGGAAGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGc                        <  1:1871919/68‑1 (MQ=255)
               cgAAGGTCTGGATATGCCGGAAGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGc                        <  1:18297/68‑1 (MQ=255)
               cgAAGGTCTGGATATGCCGGAAGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGc                        <  1:3126705/68‑1 (MQ=255)
               cgAAGGTCTGGATATGCCGGAAGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGc                        <  1:3126163/68‑1 (MQ=255)
               cgAAGGTCTGGATATGCCGGAAGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGc                        <  1:3069952/68‑1 (MQ=255)
               cgAAGGTCTGGATATGCCGGAAGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGc                        <  1:3074237/68‑1 (MQ=255)
                            atGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAAcgttcg          >  1:2237834/1‑69 (MQ=255)
                            atGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAAcgttcg          >  1:3074592/1‑69 (MQ=255)
                            atGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAAcgttcg          >  1:516354/1‑69 (MQ=255)
                            atGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAAcgttcg          >  1:1842557/1‑69 (MQ=255)
                            atGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAAcgttcg          >  1:1821867/1‑69 (MQ=255)
                            atGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAAcgttcg          >  1:912194/1‑69 (MQ=255)
                            atGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAAcgttc           >  1:1622626/1‑68 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:462329/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:2754085/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:3281145/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:743323/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:97019/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:3003588/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:2995875/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:2902810/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:2901120/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:2673552/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:2591527/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:231209/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:2239169/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:2018021/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:1940492/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:1713563/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:1689149/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:1541879/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:1220986/1‑69 (MQ=255)
                                    aGTTTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCa  >  1:1213938/1‑69 (MQ=255)
                                       |                                                                 
GATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACAAAGAAGGCTTCCTGCGTTCCGAACGTTCGTTGATCCA  >  minE/521486‑521590

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: