Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,971,610 A→G 38.6% E142E (GAA→GAG serB → 3‑phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,971,6100AG38.6% 17.5 / 26.9 44E142E (GAA→GAGserB3‑phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base A (7/20);  new base G (11/6);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 1.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

AGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACG  >  minE/2971552‑2971676
                                                          |                                                                   
aGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGAAc                                                           >  1:1309882/1‑69 (MQ=255)
aGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGAAc                                                           >  1:923903/1‑69 (MQ=255)
aGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGAAc                                                           >  1:557603/1‑69 (MQ=255)
aGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGAAc                                                           >  1:1402357/1‑69 (MQ=255)
aGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGAAc                                                           >  1:1406093/1‑69 (MQ=255)
aGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGAAc                                                           >  1:1413233/1‑69 (MQ=255)
aGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGAAc                                                           >  1:319889/1‑69 (MQ=255)
aGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGAAc                                                           >  1:2534017/1‑69 (MQ=255)
aGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGAAc                                                           >  1:2333798/1‑69 (MQ=255)
aGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGAAc                                                           >  1:1991907/1‑69 (MQ=255)
aGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGAAc                                                           >  1:1574719/1‑69 (MQ=255)
        tGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGAt                                                   >  1:1348819/1‑68 (MQ=255)
        tGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATg                                                  >  1:85573/1‑69 (MQ=255)
        tGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATg                                                  >  1:1385206/1‑69 (MQ=255)
        tGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATg                                                  >  1:2550391/1‑69 (MQ=255)
           aTTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGt                                                                   >  1:971247/1‑50 (MQ=255)
           aTTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGt                                                                   <  1:3161993/50‑1 (MQ=255)
           aTTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGt                                                                   <  1:3154265/50‑1 (MQ=255)
           aTTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGt                                                                   >  1:3080361/1‑50 (MQ=255)
           aTTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGt                                                                   <  1:1502814/50‑1 (MQ=255)
           aTTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGt                                                                   >  1:2905747/1‑50 (MQ=255)
           aTTGATGAAATTGCAAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGt                                                                   <  1:2509530/50‑1 (MQ=255)
              gatgaAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACg                                                         <  1:1977921/56‑1 (MQ=255)
              gatgaAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACg                                                         <  1:1383247/56‑1 (MQ=255)
              gatgaAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACg                                                         <  1:3043656/56‑1 (MQ=255)
              gatgaAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACg                                                         <  1:254769/56‑1 (MQ=255)
              gatgaAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACg                                                         <  1:23569/56‑1 (MQ=255)
              gatgaAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACg                                                         <  1:1533911/56‑1 (MQ=255)
               atgaAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACg                                                         <  1:2353793/55‑1 (MQ=255)
                     ttGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGAACGGGCGATGCGCGGCGAACTCg                                     <  1:1756949/69‑1 (MQ=255)
                     ttGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGAACGGGCGATGCGCGGCGAACTCg                                     <  1:1553995/69‑1 (MQ=255)
                     ttGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGAACGGGCGATGCGCGGCGAACTCg                                     <  1:2062490/69‑1 (MQ=255)
                     ttGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGAACGGGCGATGCGCGGCGAACTCg                                     <  1:3345406/69‑1 (MQ=255)
                     ttGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGAACGGGCGATGCGCGGCGAACTCg                                     <  1:1428777/69‑1 (MQ=255)
                     ttGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCg                                     <  1:772586/69‑1 (MQ=255)
                     ttGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCg                                     <  1:908634/69‑1 (MQ=255)
                          aaaCTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGAtttt                                <  1:3060695/69‑1 (MQ=255)
                          aaaCTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGAtttt                                <  1:595819/69‑1 (MQ=255)
                           aaCTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGAtttt                                <  1:2273172/68‑1 (MQ=255)
                              tGGCCGGAAGGGGCGAGATGGTGGCGGAAGTAACCGAACg                                                         <  1:3260603/40‑1 (MQ=255)
                                           cGAGATGGTGGCGGAAGTAACCGAACGGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCa                <  1:32140/69‑1 (MQ=255)
                                           cGAGATGGTGGCGGAAGTAACCGAACGGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCa                <  1:2558982/69‑1 (MQ=255)
                                           cGAGATGGTGGCGGAAGTAACCGAACGGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCATGCGCa                <  1:2382001/69‑1 (MQ=39)
                                               aTGGTGGCGGAGGTAACCGAACGGGCGATGCGCGGCGAACTCg                                     <  1:359770/43‑1 (MQ=255)
                                                          aGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACg  >  1:3047221/1‑67 (MQ=255)
                                                          |                                                                   
AGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACG  >  minE/2971552‑2971676

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