Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 997,014 A→G 73.5% G90G (GGT→GGC ynfL ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE997,0140AG73.5% 30.3 / 7.6 30G90G (GGT→GGCynfLpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (4/4);  new base G (18/4);  total (22/8)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

AACGTATCGGACACCGCCCGAATAAAAGGAGCCGACGAAGTAAAACCAATGCGCAACTCCCCCGCTTCACCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCGTCATCCACCATAGACAGGATTTGCCGACTATC  >  minE/996946‑997079
                                                                    |                                                                 
aaCGTATCGGACACCGCCCGAATAAAAGGAGCCGACGAAGTAAAACCAATGCGCAACTCCCCCGCTTCa                                                                   <  1:3103492/69‑1 (MQ=255)
      tCGGACACCGCCCGAATAAAAGGAGCCGACGAAGTAAAACCAATGCGCAACTCCCCCGCTTCACCCTGa                                                             >  1:1137425/1‑69 (MQ=255)
           cacCGCCCGAATAAAAGGAGCCGACGAAGTAAAACCAATGCGCAACTCCCCCGCTTCACCCTGATGCAg                                                        <  1:1355776/69‑1 (MQ=255)
                     aTAAAAGGAGCCGACGAAGTAAAACCAATGAGCAACTCCCCCGCTTCACCCTGATGCAg                                                        <  1:1813447/59‑1 (MQ=255)
                        aaaGGAGCGGACGAAGTAAAACCAATGCGCAACTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGa                                           <  1:836566/69‑1 (MQ=255)
                        aaaGGAGCGGACGAAGTAAAACCAATGCGCAACTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGa                                           <  1:1337050/69‑1 (MQ=255)
                        aaaGGAGCGGACGAAGTAAAACCAATGCGCAACTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGa                                           <  1:3164572/69‑1 (MQ=255)
                                cgacgaAGTAAAACCAATGCGCAACTCCCCCGCTTCACCCTGATGCAGCCTTTCAGCGCGAGCGGCAGc                                   <  1:2057951/69‑1 (MQ=255)
                                                   cgcAACTCCCCCGCTTCACCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCGTCATCCACCATAGACAgg                >  1:989564/1‑69 (MQ=255)
                                                   cgcAACTCCCCCGCTTCACCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCGTCATCCACCATAGACAgg                <  1:3332855/69‑1 (MQ=255)
                                                   cgcAACTCCCCCGCTTCACCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCGTCATCCACCATAGACAgg                >  1:3069899/1‑69 (MQ=255)
                                                        cTCCCCCGCTTCGCCCTGATGCCGCCTTTCAGCGCGAGCGGCAGCATCATCCACCATAGACAGGAttt            >  1:422795/1‑68 (MQ=255)
                                                        cTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCATc                                >  1:828627/1‑48 (MQ=39)
                                                        cTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCATCATccacca                        >  1:2270727/1‑56 (MQ=255)
                                                        cTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCATCATCCACCATAg                     >  1:1818/1‑59 (MQ=255)
                                                        cTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCATCATCCACCATAGACAGGAttt            >  1:1609707/1‑68 (MQ=255)
                                                        cTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCATCATCCACCATAGACAGGAttt            >  1:1081684/1‑68 (MQ=255)
                                                        cTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCATCATCCACCATAGACAGGAttt            >  1:681276/1‑68 (MQ=255)
                                                        cTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCATCATCCACCATAGACAGGAttt            >  1:392083/1‑68 (MQ=255)
                                                        cTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCATCATCCACCATAGACAGGAttt            >  1:3117074/1‑68 (MQ=255)
                                                        cTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCATCATCCACCATAGACAGGAttt            >  1:3094159/1‑68 (MQ=255)
                                                        cTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCATCATCCACCATAGACAGGAttt            >  1:1633080/1‑68 (MQ=255)
                                                        cTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCATCATCCACCATAGACAGGAttt            >  1:2448391/1‑68 (MQ=255)
                                                        cTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCATCATCCACCATAGACAGGAttt            >  1:2242683/1‑68 (MQ=255)
                                                        cTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCATCATCCACCATAGACAGGAttt            >  1:2107816/1‑68 (MQ=255)
                                                        cTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCATCATCCACCATAGACAGGAttt            >  1:2059733/1‑68 (MQ=255)
                                                        cTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCATCATCCACCATAGACAGGAttt            >  1:1757023/1‑68 (MQ=255)
                                                        cTCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCATCATCCACCATAGACAGGAtt             >  1:1227081/1‑67 (MQ=255)
                                                         tCCCCCGCTTCGCCCTGATGCAGCCTTTCAGcgcgt                                           <  1:1392683/36‑2 (MQ=38)
                                                         tCCCCCGCTTCGCCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCATCATCCACCATAGACAGGAtt             >  1:1460218/1‑66 (MQ=255)
                                                                  tCACCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCGTCATCCACCATAGACAGTATTTGCCGACTATc  >  1:2140779/1‑68 (MQ=255)
                                                                    |                                                                 
AACGTATCGGACACCGCCCGAATAAAAGGAGCCGACGAAGTAAAACCAATGCGCAACTCCCCCGCTTCACCCTGATGCAGCCTTTCAGCGCGAGCGGCAGCGTCATCCACCATAGACAGGATTTGCCGACTATC  >  minE/996946‑997079

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: