Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 29,400 A→G 34.0% P302P (CCA→CCG carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE29,4000AG34.0% 28.7 / 22.2 44P302P (CCA→CCGcarBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base A (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 = 1.45e-01

GTTCAGTTTGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTCCCGATTGCTAAA  >  minE/29344‑29457
                                                        |                                                         
gTTCAGTTTGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCa                                                           >  1:1850509/1‑57 (MQ=255)
gTTCAGTTTGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCa                                                           >  1:2170772/1‑57 (MQ=255)
gTTCAGTTTGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCa                                                           >  1:2911709/1‑57 (MQ=255)
gTTCAGTTTGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCa                                                           >  1:437869/1‑57 (MQ=255)
gTTCAGTTTGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCa                                                           >  1:745514/1‑57 (MQ=255)
gTTCAGTTTGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCa                                                           >  1:1156396/1‑57 (MQ=255)
   cAGTTTGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGttct                                            >  1:2977399/1‑69 (MQ=255)
   cAGTTTGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGttct                                            >  1:1034269/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:561452/69‑1 (MQ=255)
                   cGAAGAACGGTCGCCTGATTGTTATCGAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGa                            <  1:1076309/69‑1 (MQ=255)
                   cGAAGAACGGTCGCCTGATTGTTATCGAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGa                            <  1:459435/69‑1 (MQ=255)
                   cGAAGAACGGTCGCCTGATTGTTATCGAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGa                            <  1:928579/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:2292175/1‑69 (MQ=255)
                               gtcTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGtt                >  1:1123419/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:2117219/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:2029511/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:3089491/1‑69 (MQ=255)
                                             gAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGAAAGCCACCGGCTTCCCGATTGCCaaa  >  1:3278118/1‑69 (MQ=255)
                                             gAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGAAAGCCACCGGCTTCCCGATTGCCaaa  >  1:3304011/1‑69 (MQ=255)
                                             gAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGAAAGCCACCGGCTTCCCGATTGCCaaa  >  1:1917441/1‑69 (MQ=255)
                                             gAAATGAACCCGCGTGTGTCGCGTTCTTCGGCGCTGGCGTCGAAAGCCACCGGCTTCCCGATTGCCaaa  >  1:2833375/1‑69 (MQ=255)
                                             gAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTCCCGATTGCTaa   <  1:601925/68‑1 (MQ=255)
                                                        |                                                         
GTTCAGTTTGCGGTGAACCCGAAAAACGGTCGTCTGATTGTTATCGAAATGAACCCACGCGTGTCCCGTTCTTCGGCGCTGGCGTCGAAAGCGACCGGTTTCCCGATTGCTAAA  >  minE/29344‑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: