Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,478,949 A→G 42.8% V326V (GTT→GTC glpT ← sn‑glycerol‑3‑phosphate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,478,9490AG42.8% 3.5 / 18.0 28V326V (GTT→GTCglpTsn‑glycerol‑3‑phosphate transporter
Reads supporting (aligned to +/- strand):  ref base A (13/3);  new base G (8/4);  total (21/7)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTCATCCAGTAAACGATAGTCGCGATGGTCACCAGTGTCATAAAGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGACTTTATCCGACATCCAGCCGCACAGCAGAGT  >  minE/1478900‑1479011
                                                 |                                                              
ggTTCATCCAGTAAACGATAGTCGCGATGGTCACCAGTGTCATAAAGAAGACGCCGGTTGCCCCACGGt                                             <  1:1190026/69‑1 (MQ=255)
ggTTCATCCAGTAAACGATAGTCGCGATGGTCACCAGTGTCATAAAGAAGACGCCGGTTGCCCCACGGt                                             <  1:1458410/69‑1 (MQ=255)
ggTTCATCCAGTAAACGATAGTCGCGATGGTCACCAGTGTCATAAAGAAGACGCCGGTTGCCCCACGGt                                             <  1:672290/69‑1 (MQ=255)
     aTCCAGTAAACGATAGTCGCGATGGTCACCAGTGTCATAAAGAAGACGCCGGTTGCCCCACGGt                                             <  1:1515840/64‑1 (MQ=255)
             aaCGATAGTCGCGATGGTCACCAGTGTCATAAAGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGAc                                >  1:1794214/1‑69 (MQ=255)
             aaCGATAGTCGCGATGGTCACCAGTGTCATAAAGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGAc                                >  1:3083280/1‑69 (MQ=255)
             aaCGATAGTCGCGATGGTCACCAGTGTCATAAAGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGAc                                >  1:2695186/1‑69 (MQ=255)
             aaCGATAGTCGCGATGGTCACCAGTGTCATAAAGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGAc                                >  1:735796/1‑69 (MQ=255)
             aaCGATAGTCGCGATGGTCACCAGTGTCATAAAGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGAc                                >  1:945449/1‑69 (MQ=255)
               ggATAGTCGCGATGGTCACCAGTGTCATAAAGAAGACGCCGGTTGCCCCACGGTTACAACGAAAGACtt                              >  1:246290/2‑69 (MQ=255)
               cGATAGTCGCGATGGTCACCAGTGTCATAAAGAAGACGCCGGTTGccc                                                   >  1:1522590/1‑48 (MQ=255)
               cGATAGTCGCGATGGTCACCAGTGTCATAAAGAAGACGCCGGTTGCCCCACGGTTACCACGaaa                                   >  1:1734312/1‑64 (MQ=255)
               cGATAGTCGCGATGGTCACCAGTGTCATAAAGAAGACGCCGGTTGCCCCACGGTTACCACGAAAGACtt                              >  1:1674711/1‑69 (MQ=255)
               cGATAGTCGCGATGGTCACCAGTGTCATAAAGAAGACGCCGGTTGCCCCACGGTTACCACGAAAGACtt                              >  1:692577/1‑69 (MQ=255)
               cGATAGTCGCGATGGTCACCAGTGTCATAAAGAAGACGCCGGTTGCCCCACGGTTACCACGAAAGACtt                              >  1:563778/1‑69 (MQ=255)
               cGATAGTCGCGATGGTCACCAGTGTCATAAAGAAGACGCCGGTTGCCCCACGGTTACCACGAAAGACtt                              >  1:3101852/1‑69 (MQ=255)
               cGATAGTCGCGATGGTCACCAGTGTCATAAAGAAGACGCCGGTTGCCCCACGGTTACCACGAAAGACtt                              >  1:3197476/1‑69 (MQ=255)
                  tAGTCGCGATGGTCACCAGTGTCATAAAGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGACTTTAt                           >  1:2664932/1‑69 (MQ=255)
                  tAGTCGCGATGGTCACCAGTGTCATAAAGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGACTTTAt                           >  1:596511/1‑69 (MQ=255)
                  tAGTCGCGATGGTCACCAGTGTCATAAAGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGACTTTAt                           >  1:3224017/1‑69 (MQ=255)
                  tAGTCGCGATGGTCACCAGTGTCATAAAGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGACTTTAt                           >  1:1043853/1‑69 (MQ=255)
                  tAGTCGCGATGGTCACCAGTGTCATAAAGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGACTTTAt                           >  1:2505229/1‑69 (MQ=255)
                                         aTAAAGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGACTTTATCCGACATCCAGCCGCACAGCaga    >  1:259817/1‑69 (MQ=255)
                                         aTAAAGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGACTTTATCCGACATCCAGCCGCACAGCaga    >  1:2025565/1‑69 (MQ=255)
                                         aTAAAGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGACTTTATCCGACATCCAGCCGCACAGCaga    >  1:1386960/1‑69 (MQ=255)
                                           aaaGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGACTTTATCCGACATCCAGCCGCACAGCAGAGt  <  1:978088/69‑1 (MQ=255)
                                            aaGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGACTTTATCCGACATCCAGCCGCACAGCAGAGt  <  1:747628/68‑1 (MQ=255)
                                            aaGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGACTTTATCCGACATCCAGCCGCACAGCAGAGt  <  1:1313634/68‑1 (MQ=255)
                                                 |                                                              
GGTTCATCCAGTAAACGATAGTCGCGATGGTCACCAGTGTCATAAAGAAAACGCCGGTTGCCCCACGGTTGCCACGGAAGACTTTATCCGACATCCAGCCGCACAGCAGAGT  >  minE/1478900‑1479011

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: