Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,481,302 G→C 46.9% A368A (GCG→GCC 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,3020GC46.9% 1.2 / 55.6 49A368A (GCG→GCCglpAsn‑glycerol‑3‑phosphate dehydrogenase (anaerobic), large subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (11/15);  new base C (16/7);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 8.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

ATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTGCCGATCTGGCACTGCC  >  minE/1481240‑1481358
                                                              |                                                        
atcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                                        <  1:1006317/65‑1 (MQ=255)
atcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                                        <  1:1147229/65‑1 (MQ=255)
atcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                                        <  1:2453045/65‑1 (MQ=255)
atcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                                        <  1:1901328/65‑1 (MQ=255)
atcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                                        <  1:251448/65‑1 (MQ=255)
atcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                                        <  1:3116875/65‑1 (MQ=255)
atcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                                        <  1:3031156/65‑1 (MQ=255)
 tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                                        >  1:2342331/1‑64 (MQ=255)
 tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                                        >  1:1058631/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)
 tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                                        >  1:3046202/1‑64 (MQ=255)
 tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                                        >  1:3271032/1‑64 (MQ=255)
 tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                                        >  1:501529/1‑64 (MQ=255)
 tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                                        >  1:781810/1‑64 (MQ=255)
 tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                                        >  1:871866/1‑64 (MQ=255)
 tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                                        >  1:981867/1‑64 (MQ=255)
 tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                                        >  1:1867345/1‑64 (MQ=255)
 tcaccatcaccGGTGGCAAACTGATGACCTATCGGCTGATGGCGGAATGGGCTACCGACGCCGt                                                        >  1:1616339/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:2168118/1‑69 (MQ=255)
                     cTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACgcgc                               >  1:821093/1‑69 (MQ=255)
                     cTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACgcgc                               >  1:635198/1‑69 (MQ=255)
                     cTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACgcgc                               >  1:1747186/1‑69 (MQ=255)
                     cTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACgcgc                               >  1:2842978/1‑69 (MQ=255)
                     cTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACgcgc                               >  1:2376815/1‑69 (MQ=255)
                     cTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACgcgc                               >  1:2476118/1‑69 (MQ=255)
                              tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTAc                       <  1:821843/68‑1 (MQ=255)
                              tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTAc                       <  1:658572/68‑1 (MQ=255)
                              tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg                      <  1:1795713/69‑1 (MQ=255)
                              tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg                      <  1:296483/69‑1 (MQ=255)
                              tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg                      <  1:2787748/69‑1 (MQ=255)
                              tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg                      <  1:3136174/69‑1 (MQ=255)
                              tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg                      <  1:900938/69‑1 (MQ=255)
                              tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg                      <  1:3408221/69‑1 (MQ=255)
                              tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg                      <  1:2643409/69‑1 (MQ=255)
                              tATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg                      <  1:808827/69‑1 (MQ=255)
                               aTCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg                      <  1:2808916/68‑1 (MQ=255)
                               aTCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACg                      <  1:184107/68‑1 (MQ=255)
                                  ggCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTg                  <  1:1139014/69‑1 (MQ=255)
                                  ggCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTg                  <  1:2018967/69‑1 (MQ=255)
                                  ggCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTg                  <  1:1325981/69‑1 (MQ=255)
                                                  ggCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTGCCGATCTGGCACTGcc  >  1:593250/1‑69 (MQ=255)
                                                  ggCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTGCCGATCTGGCACTGcc  >  1:1318359/1‑69 (MQ=255)
                                                  ggCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTGCCGATCTGGCACTGcc  >  1:3386506/1‑69 (MQ=255)
                                                              |                                                        
ATCACCATCACCGGTGGCAAACTGATGACCTATCGGCTGATGGCTGAATGGGCTACCGACGCGGTATGCCGCAAACTGGGCAACACGCGCCCCTGTACGACTGCCGATCTGGCACTGCC  >  minE/1481240‑1481358

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: