Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,183,993 T→C 35.9% G1423G (GGT→GGC gltB → glutamate synthase, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,183,9930TC35.9% 46.3 / 44.6 75G1423G (GGT→GGCgltBglutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base T (36/12);  new base C (23/4);  total (59/16)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAG  >  minE/2183932‑2184056
                                                             |                                                               
gtggtATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCtt                                                          >  1:1462166/1‑69 (MQ=255)
 tggtATCGTCTGCATTCTGGGTAAAAGCGGCGTTAACTACGGTGCGGGCATGACCGGCGGTTTCGCTTa                                                         >  1:2004502/1‑69 (MQ=255)
 tggtATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTa                                                         >  1:668938/1‑69 (MQ=255)
      tCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGtt                                                     >  1:1377704/1‑68 (MQ=255)
      tCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTtc                                                    >  1:539861/1‑69 (MQ=255)
      tCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTtc                                                    >  1:1116383/1‑69 (MQ=255)
      tCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTtc                                                    >  1:18809/1‑69 (MQ=255)
        gTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTtctc                                                  >  1:1093396/1‑69 (MQ=255)
        gTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTtctc                                                  >  1:408952/1‑69 (MQ=255)
        gTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTtctc                                                  >  1:3277690/1‑69 (MQ=255)
        gTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTtctc                                                  >  1:2713262/1‑69 (MQ=255)
        gTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTtctc                                                  >  1:2559771/1‑69 (MQ=255)
        gTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTtctc                                                  >  1:2454631/1‑69 (MQ=255)
        gTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTtctc                                                  >  1:1950041/1‑69 (MQ=255)
        gTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTtctc                                                  >  1:151409/1‑69 (MQ=255)
        gTCTGCATTCTGGGTAAAACCGGAGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTtctc                                                  >  1:1594196/1‑69 (MQ=255)
           tGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGAt                                               <  1:11433/69‑1 (MQ=255)
           tGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGAt                                               <  1:63892/69‑1 (MQ=255)
                 cTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAgc                                         >  1:1389896/1‑69 (MQ=255)
                 cTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAgc                                         >  1:812544/1‑69 (MQ=255)
                 cTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAgc                                         >  1:528176/1‑69 (MQ=255)
                 cTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAg                                          >  1:2502503/1‑68 (MQ=255)
                 cTGGGTAAAACCGGCGTTAACTTCGGCGCGGGCATGACCGGTGGCTTCGCTTACGTTCTCGATGAAAgc                                         <  1:1082862/69‑1 (MQ=255)
                 cTGGGTAAAACCGGCGTTAACTTCGGCGCGGGCATGACCGGTGGCTTCGCTTACGTTCTCGATGAAAgc                                         <  1:3138120/69‑1 (MQ=255)
                 cTGGGTAAAACCGGCGTTAACTTCGGCGCGGGCATGACCGGTGGCTTCGCTTACGTTCTCGATGAAAgc                                         <  1:786047/69‑1 (MQ=255)
                  tGGGTAAAACCGGCGTTAACTTCGGCGCGGGCATGACCGGTGGCTTCGCTTACGTTCTCGATGAAAgc                                         <  1:1411166/68‑1 (MQ=255)
                       aaaaCCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGAt                                   >  1:1976468/1‑69 (MQ=255)
                       aaaaCCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGAt                                   >  1:1711929/1‑69 (MQ=255)
                                  tAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAAcg                        <  1:1573930/69‑1 (MQ=255)
                                   aaCTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAgc                                         >  1:2113963/1‑51 (MQ=255)
                                   aaCTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAAcgc                       >  1:2289564/1‑69 (MQ=255)
                                   aaCTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAAcgc                       >  1:2625508/1‑69 (MQ=255)
                                   aaCTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAAcgc                       >  1:467835/1‑69 (MQ=255)
                                   aaCTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAAcgc                       >  1:3127318/1‑69 (MQ=255)
                                   aaCTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAAcgc                       >  1:2252162/1‑69 (MQ=255)
                                   aaCTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAAcgc                       >  1:2153543/1‑69 (MQ=255)
                                   aaCTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAAcgc                       >  1:93903/1‑69 (MQ=255)
                                   aaCTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAAcgc                       >  1:944283/1‑69 (MQ=255)
                                   aaCTTCGGTGCGGGCATCACCGGCGGTTTCGCTTACGtt                                                     >  1:3176537/1‑39 (MQ=255)
                                        cGGCGCGGGCATGACCTGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTaa                  >  1:175930/1‑69 (MQ=255)
                                        cGGCGCGGGCATGACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTaa                  >  1:241118/1‑69 (MQ=255)
                                        cGGCGCGGGCATGACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTaa                  >  1:2457583/1‑69 (MQ=255)
                                        cGGCGCGGGCATGACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTaa                  >  1:681821/1‑69 (MQ=255)
                                        cGGCGCGGGCATGACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTaa                  >  1:1262471/1‑69 (MQ=255)
                                        cGGCGCGGGCATGACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTaa                  >  1:1903137/1‑69 (MQ=255)
                                        cGGCGCGGGCATGACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTaa                  >  1:3424001/1‑69 (MQ=255)
                                        cGGCGCGGGCATGACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTaa                  >  1:1192096/1‑69 (MQ=255)
                                        cGGCGCGGGCATGACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTaa                  >  1:3188512/1‑69 (MQ=255)
                                             cGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCgg             <  1:947451/69‑1 (MQ=255)
                                                    gACCGGTGGCTTCGCTTACGTTCTCGATGAAAGGGGCGATTTCCGCAAACGCGTTAACCCGGAACTgg       >  1:2753300/1‑68 (MQ=255)
                                                    gACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTa                   >  1:1289872/1‑56 (MQ=255)
                                                    gACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTgg       >  1:1567722/1‑68 (MQ=255)
                                                    gACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTgg       >  1:1548050/1‑68 (MQ=255)
                                                    gACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTgg       >  1:1904398/1‑68 (MQ=255)
                                                    gACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTgg       >  1:901198/1‑68 (MQ=255)
                                                    gACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTgg       >  1:1934879/1‑68 (MQ=255)
                                                    gACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTgg       >  1:791957/1‑68 (MQ=255)
                                                    gACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTgg       >  1:2325533/1‑68 (MQ=255)
                                                    gACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTgg       >  1:612398/1‑68 (MQ=255)
                                                    gACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTgg       >  1:2727397/1‑68 (MQ=255)
                                                    gACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTgg       >  1:2771311/1‑68 (MQ=255)
                                                    gACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTgg       >  1:33060/1‑68 (MQ=255)
                                                    gACCGGTGGCTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTg        >  1:360329/1‑67 (MQ=255)
                                                     aCCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCa                            >  1:390184/1‑46 (MQ=255)
                                                     aCCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGa            >  1:778827/1‑62 (MQ=255)
                                                     aCCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTgg       >  1:2639088/1‑67 (MQ=255)
                                                     aCCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTgg       >  1:3204979/1‑67 (MQ=255)
                                                        ggcggTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCTAg  <  1:111562/69‑1 (MQ=255)
                                                        ggcggTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAg  <  1:361075/69‑1 (MQ=255)
                                                        ggcggTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAg  <  1:2264328/69‑1 (MQ=255)
                                                        ggcggTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAg  <  1:747419/69‑1 (MQ=255)
                                                        ggcggTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAg  <  1:2091619/69‑1 (MQ=255)
                                                        ggcggTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAg  <  1:350340/69‑1 (MQ=255)
                                                        ggcggTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAg  <  1:829408/69‑1 (MQ=255)
                                                         gcggTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAg  <  1:2865398/68‑1 (MQ=255)
                                                             |                                                               
GTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAG  >  minE/2183932‑2184056

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: