Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,481,284 T→G 34.2% A362A (GCT→GCG glpA → sn‑glycerol‑3‑phosphate dehydrogenase (anaerobic), large subunit, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,481,2840TG34.2% 63.2 / 61.5 79A362A (GCT→GCGglpAsn‑glycerol‑3‑phosphate dehydrogenase (anaerobic), large subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (17/35);  new base G (20/7);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 7.48e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-03

CGCGATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTG  >  minE/1481219‑1481342
                                                                 |                                                          
ggcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:1085324/66‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:1045210/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:2368819/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:2203096/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:21939/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:2593614/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:2055344/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:2035998/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:2007900/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:1989093/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:3120416/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:3167952/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:1703422/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:3302229/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:160457/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:775701/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:717609/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:624980/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:122497/67‑1 (MQ=255)
cgcgATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTg                                                           <  1:556618/67‑1 (MQ=255)
    aTGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATgg                                                      >  1:3425078/1‑68 (MQ=255)
    aTGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATgg                                                      >  1:674748/1‑68 (MQ=255)
    aTGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATgg                                                      >  1:2160957/1‑68 (MQ=255)
    aTGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATgg                                                      >  1:2971005/1‑68 (MQ=255)
    aTGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATgg                                                      >  1:362531/1‑68 (MQ=255)
    aTGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATgg                                                      >  1:1938256/1‑68 (MQ=255)
    aTGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATgg                                                      >  1:353402/1‑68 (MQ=255)
    aTGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATgg                                                      >  1:3157910/1‑68 (MQ=255)
    aTGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATgg                                                      >  1:1410876/1‑68 (MQ=255)
                     atcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        <  1:3031156/65‑1 (MQ=255)
                     atcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        <  1:3116875/65‑1 (MQ=255)
                     atcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        <  1:1006317/65‑1 (MQ=255)
                     atcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        <  1:2453045/65‑1 (MQ=255)
                     atcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        <  1:251448/65‑1 (MQ=255)
                     atcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        <  1:1147229/65‑1 (MQ=255)
                     atcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        <  1:1901328/65‑1 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGc                                                    >  1:860095/1‑52 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGAcgccg                                         >  1:2239019/1‑61 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGAcgccg                                         >  1:1172525/1‑61 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGAcgccg                                         >  1:3147149/1‑61 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        >  1:2342331/1‑64 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        >  1:981867/1‑64 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        >  1:871866/1‑64 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        >  1:781810/1‑64 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        >  1:1058631/1‑64 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        >  1:501529/1‑64 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        >  1:1616339/1‑64 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        >  1:3271032/1‑64 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        >  1:1867345/1‑64 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        >  1:3046202/1‑64 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        >  1:2345436/1‑64 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        >  1:271791/1‑64 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        >  1:2850545/1‑64 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGTATga                                    >  1:2243668/1‑67 (MQ=255)
                      tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGAGGAATGGGCTACCGACGCCGt                                        >  1:1512343/1‑64 (MQ=255)
                      tcacaatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                        >  1:1265165/1‑64 (MQ=255)
                                          cTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGCCGCGGTATGCCGCAAACTGGGCAACACgcgc               >  1:1294954/1‑69 (MQ=255)
                                          cTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACgcgc               >  1:2842978/1‑69 (MQ=255)
                                          cTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACgcgc               >  1:2376815/1‑69 (MQ=255)
                                          cTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACgcgc               >  1:821093/1‑69 (MQ=255)
                                          cTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACgcgc               >  1:2476118/1‑69 (MQ=255)
                                          cTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACgcgc               >  1:2168118/1‑69 (MQ=255)
                                          cTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACgcgc               >  1:1747186/1‑69 (MQ=255)
                                          cTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACgcgc               >  1:635198/1‑69 (MQ=255)
                                                   tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTAc       <  1:658572/68‑1 (MQ=255)
                                                   tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTAc       <  1:821843/68‑1 (MQ=255)
                                                   tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg      <  1:3136174/69‑1 (MQ=255)
                                                   tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg      <  1:2643409/69‑1 (MQ=255)
                                                   tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg      <  1:2787748/69‑1 (MQ=255)
                                                   tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg      <  1:808827/69‑1 (MQ=255)
                                                   tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg      <  1:1795713/69‑1 (MQ=255)
                                                   tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg      <  1:3408221/69‑1 (MQ=255)
                                                   tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg      <  1:900938/69‑1 (MQ=255)
                                                   tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg      <  1:296483/69‑1 (MQ=255)
                                                    aTCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg      <  1:184107/68‑1 (MQ=255)
                                                    aTCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg      <  1:2808916/68‑1 (MQ=255)
                                                       ggCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTg  <  1:2018967/69‑1 (MQ=255)
                                                       ggCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTg  <  1:1139014/69‑1 (MQ=255)
                                                       ggCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTg  <  1:1325981/69‑1 (MQ=255)
                                                                 |                                                          
CGCGATGGTCTGGACGGATTTATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTG  >  minE/1481219‑1481342

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