Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 29,409 C→G 34.1% S305S (TCC→TCG carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE29,4090CG34.1% 39.1 / 32.1 44S305S (TCC→TCGcarBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base C (12/17);  new base G (6/9);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CAGTTTGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTCCCGATTGCTAAA  >  minE/29347‑29457
                                                              |                                                
cAGTTTGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGttct                                            >  1:1034269/1‑69 (MQ=255)
cAGTTTGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGttct                                            >  1:2977399/1‑69 (MQ=255)
cAGTTTGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGttc                                             >  1:2358939/1‑68 (MQ=255)
    ttGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCgg                                        <  1:834472/69‑1 (MQ=255)
    ttGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCgg                                        <  1:1831634/69‑1 (MQ=255)
           gAACCCGAAATACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGc                                 <  1:1572758/69‑1 (MQ=255)
           gAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGc                                 <  1:1172551/69‑1 (MQ=255)
           gAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGc                                 <  1:188589/69‑1 (MQ=255)
           gAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGc                                 <  1:2887483/69‑1 (MQ=255)
           gAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAAACCACGCGTGTCCCGTTCTTCGGCGCTGGc                                 <  1:3023341/69‑1 (MQ=255)
                cGAAGAACGGTCGCCTGATTGTTATCGAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGa                            <  1:1076309/69‑1 (MQ=255)
                cGAAGAACGGTCGCCTGATTGTTATCGAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGa                            <  1:928579/69‑1 (MQ=255)
                cGAAGAACGGTCGCCTGATTGTTATCGAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGa                            <  1:561452/69‑1 (MQ=255)
                cGAAGAACGGTCGCCTGATTGTTATCGAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGa                            <  1:459435/69‑1 (MQ=255)
                cGAAGAACGGTCGCCTGATTGTTATCGAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGa                            <  1:2086857/69‑1 (MQ=255)
                cGAAGAACGGTCGCCTGATTGTTATCGAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGa                            <  1:2037735/69‑1 (MQ=255)
                 gaagaaCGGTCGCCTGATTGTTATCGAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGa                            <  1:1515591/68‑1 (MQ=255)
                            gtcTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGtt                >  1:1123419/1‑69 (MQ=255)
                            gtcTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGtt                >  1:2292175/1‑69 (MQ=255)
                            gtcTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGtt                >  1:790950/1‑69 (MQ=255)
                            gtcTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGtt                >  1:3237272/1‑69 (MQ=255)
                            gtcTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGtt                >  1:1952156/1‑69 (MQ=255)
                              cTGATTGTTATCGAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGa                            <  1:2300043/55‑1 (MQ=255)
                              cTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTc              >  1:2548505/1‑69 (MQ=255)
                              cTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTc              >  1:1885181/1‑69 (MQ=255)
                               tGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTc              >  1:2737415/1‑68 (MQ=255)
                                 aTTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTCCCg           >  1:66657/1‑69 (MQ=255)
                                   tGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTccc            <  1:2029511/66‑1 (MQ=255)
                                   tGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTccc            <  1:729245/66‑1 (MQ=255)
                                   tGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTccc            <  1:2510448/66‑1 (MQ=255)
                                   tGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTccc            <  1:1996435/66‑1 (MQ=255)
                                   tGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTccc            <  1:2117219/66‑1 (MQ=255)
                                   tGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTCCCg           <  1:421408/67‑1 (MQ=255)
                                   tGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTCCCg           <  1:1840693/67‑1 (MQ=255)
                                   tGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTCCCg           <  1:1828050/67‑1 (MQ=255)
                                       aTCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGc                                 <  1:2407420/41‑1 (MQ=255)
                                          gagaTGAACCCGCGCGTGTCGCGTTCGTCGGCGCTGGCCTCGAAGGCCACCGGCTTCCCGATCGCGaa   <  1:2421518/65‑1 (MQ=255)
                                          gAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGAAAGCCACCGGCTTcc             >  1:257960/1‑58 (MQ=255)
                                          gAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGAAAGCCACCGGCTTCCCGATTGCCaaa  >  1:1917441/1‑69 (MQ=255)
                                          gAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGAAAGCCACCGGCTTCCCGATTGCCaaa  >  1:2833375/1‑69 (MQ=255)
                                          gAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGAAAGCCACCGGCTTCCCGATTGCCaaa  >  1:3089491/1‑69 (MQ=255)
                                          gAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGAAAGCCACCGGCTTCCCGATTGCCaaa  >  1:3278118/1‑69 (MQ=255)
                                          gAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGAAAGCCACCGGCTTCCCGATTGCCaaa  >  1:3304011/1‑69 (MQ=255)
                                          gAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTCCCGATTGCTaa   <  1:601925/68‑1 (MQ=255)
                                                              |                                                
CAGTTTGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTCCCGATTGCTAAA  >  minE/29347‑29457

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: