Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,690,326 C→T 67.3% A557V (GCG→GTG)  yhgE → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,690,3260CT67.3% 52.9 / 36.1 71A557V (GCG→GTG) yhgEpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (10/13);  new base T (16/31);  total (27/44)
Fisher's exact test for biased strand distribution p-value = 5.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGAAAGAAAAATTACTCACAGCTGCGAAGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGAATGTGAATCTGAATGAGTGACAAAGCGC  >  minE/2690267‑2690389
                                                           |                                                               
aaCGAAAGAAAAATTACTCATGGCTGCGGAGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGa                                                        >  1:1743128/1‑69 (MQ=255)
aaCGAAAGAAAAATTACTCATGGCTGCGGAGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGa                                                        >  1:2591195/1‑69 (MQ=255)
aaCGAAAGAAAAATTACTCATGGCTGCGGAGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGa                                                        >  1:2666743/1‑69 (MQ=255)
aaCGAAAGAAAAATTACTCATGGCTGCGGAGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGa                                                        >  1:2438673/1‑69 (MQ=255)
aaCGAAAGAAAAATTACTCATGGCTGCGGAGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGa                                                        >  1:225336/1‑69 (MQ=255)
aaCGAAAGAAAAATTACTCATGGCTGCGGAGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGa                                                        >  1:2708315/1‑69 (MQ=255)
aaCGAAAGAAAAATTACTCATGGCTGCGGAGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGa                                                        >  1:2162083/1‑69 (MQ=255)
aaCGAAAGAAAAATTACTCATGGCTGCGGAGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGa                                                        >  1:1688312/1‑69 (MQ=255)
aaCGAAAGAAAAATTACTCATGGCTGCGGAGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGa                                                        >  1:3080951/1‑69 (MQ=255)
aaCGAAAGAAAAATTACTCATGGCTGCGGAGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGa                                                        >  1:3192143/1‑69 (MQ=255)
aaCGAAAGAAAAATTACTCATGGCTGCGGAGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGa                                                        >  1:387508/1‑69 (MQ=255)
aaCGAAAGAAAAATTACTCATGGCTGCGGAGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGa                                                        >  1:100227/1‑69 (MQ=255)
aaCGAAAGAAAAATTACTCATGGCTGCGGAGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGa                                                        >  1:1051600/1‑69 (MQ=255)
aaCGAAAGAAAAATTACTCATGGCTGCGGAGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGa                                                        >  1:86855/1‑69 (MQ=255)
aaCGAAAGAAAAATTACTCATGGCTGCGGAGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGa                                                        >  1:573825/1‑69 (MQ=255)
aaCGAAAGAAAAATTACTCATGGCTGCGGAGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGa                                                        >  1:1180709/1‑69 (MQ=255)
aaCGAAAGAAAAATTACTCATGGCTGCGGAGGATGGAAAACTGGGGACGAAGCCTAAAGAGTGGCGCGa                                                        >  1:2214219/1‑69 (MQ=255)
     aagaaaAATTACTCACAGCTGCGAAGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTg                                                   <  1:401242/69‑1 (MQ=255)
     aagaaaAATTACTCACAGCTGCGAAGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTg                                                   <  1:2673579/69‑1 (MQ=255)
     aagaaaAATTACTCACAGCTGCGAAGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTg                                                   <  1:2452524/69‑1 (MQ=255)
     aagaaaAATTACTCACAGCTGCGAAGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTg                                                   <  1:2167602/69‑1 (MQ=255)
     aagaaaAATTACTCACAGCTGCGAAGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTg                                                   <  1:1480468/69‑1 (MQ=255)
            aTTACTCACAGCTGCGAAGTATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGa                                            >  1:2914708/1‑69 (MQ=255)
                         gCGAAGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTg                               <  1:1642984/69‑1 (MQ=255)
                           ggaGGTTGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:2831672/66‑2 (MQ=255)
                           ggaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:2711532/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgt                              <  1:2025956/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:2828724/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:2716182/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:2713874/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:2831932/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:997544/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:637643/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:2902185/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:929079/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:56888/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:1919095/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:2573788/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:1631567/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:1474194/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:1394890/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:1636641/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:1685038/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:1108597/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:1769830/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:1957731/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:1102749/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:197442/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:203958/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:2080184/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:2223768/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:1057591/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:1023940/66‑2 (MQ=255)
                            gaGGATGGAAAACTGGGGACGAAGCCTAAACTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:3113384/66‑2 (MQ=255)
                            aaGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGa                              <  1:3119387/67‑1 (MQ=255)
                            aaGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGa                              <  1:2480558/67‑1 (MQ=255)
                            aaGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGa                              <  1:3266874/67‑1 (MQ=255)
                            aaGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGa                              <  1:3435245/67‑1 (MQ=255)
                            aaGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGa                              <  1:3080137/67‑1 (MQ=255)
                            aaGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGAAt                            >  1:951716/1‑69 (MQ=255)
                            aaGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGAAt                            >  1:1521398/1‑69 (MQ=255)
                            aaGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGAAt                            >  1:1240197/1‑69 (MQ=255)
                                   gAAAACTGGGGACGAAGCCTAAAGTGTGGCGCGATCTTACGGTGGATGGCGAAAGACTgg                              <  1:373792/60‑2 (MQ=255)
                                   gAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGAATGTGAAt                      >  1:2671237/1‑68 (MQ=255)
                                     aaaCTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGAATGTGAATCTg                   <  1:60526/69‑1 (MQ=255)
                                     aaaCTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGAATGTGAATCTg                   <  1:1568225/69‑1 (MQ=255)
                                                        aaGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGAATGTGAATCTGAATGAGTGACAAAgcgc  >  1:2946776/1‑67 (MQ=255)
                                                        aaGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGAATGTGAATCTGAATGAGTGACAAAgcgc  >  1:2849613/1‑67 (MQ=255)
                                                        aaGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGAATGTGAATCTGAATGAGTGACAAAgcgc  >  1:275756/1‑67 (MQ=255)
                                                        aaGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGAATGTGAATCTGAATGAGTGACAAAgcgc  >  1:3264741/1‑67 (MQ=255)
                                                        aaGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGAATGTGAATCTGAATGAGTGACAAAgcgc  >  1:3140667/1‑67 (MQ=255)
                                                           |                                                               
AACGAAAGAAAAATTACTCACAGCTGCGAAGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGAATGTGAATCTGAATGAGTGACAAAGCGC  >  minE/2690267‑2690389

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