Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 30,714 G→A 62.8% T740T (ACG→ACA carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE30,7140GA62.8% 20.6 / 24.3 43T740T (ACG→ACAcarBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base G (9/7);  new base A (20/7);  total (29/14)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATG  >  minE/30657‑30766
                                                         |                                                    
cggTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGc                                             <  1:567940/67‑1 (MQ=255)
cggTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGc                                             <  1:1995324/67‑1 (MQ=255)
cggTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGc                                             <  1:1806892/67‑1 (MQ=255)
cggTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGc                                             <  1:951719/67‑1 (MQ=255)
cggTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGc                                             <  1:543686/67‑1 (MQ=255)
cggTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGc                                             <  1:593919/67‑1 (MQ=255)
cggTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGc                                             <  1:2300317/67‑1 (MQ=255)
    cGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCt                                       >  1:1302967/1‑69 (MQ=255)
                  cGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATgcg                              >  1:358174/1‑64 (MQ=255)
                  cGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAgt                         <  1:1661371/69‑1 (MQ=255)
                  cGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAgt                         <  1:1281427/69‑1 (MQ=255)
                  cGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAgt                         <  1:2183257/69‑1 (MQ=255)
                  cGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAgt                         >  1:1850484/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:2379544/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:2736969/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:2966717/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:743646/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:3126262/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:3300237/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:492192/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:664437/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:692319/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:712643/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:1160017/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:1428752/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:1549972/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:1662074/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:1737311/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:1817094/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:2108064/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:798180/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGCGGTCAGCGTATCTAACGATGCGACAGTGTTg                     >  1:702354/1‑69 (MQ=255)
                      tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACAGAGGTCAGCGTATCTAACGATGCGCCAGTGTTg                     >  1:1960383/1‑69 (MQ=255)
                        tGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGc                    >  1:1963027/1‑68 (MQ=255)
                        tGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGc                    >  1:2176852/1‑68 (MQ=255)
                        tGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGc                    >  1:2842126/1‑68 (MQ=255)
                        tGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGc                    >  1:1259552/1‑68 (MQ=255)
                        tGACGAAGCTGACCTGCGACGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGc                    >  1:97544/1‑68 (MQ=255)
                                     cTGCGTCGCTACTTCCATACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTc      <  1:2756013/69‑1 (MQ=255)
                                     cTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTc      <  1:2337473/69‑1 (MQ=255)
                                     cTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTc      <  1:1785787/69‑1 (MQ=255)
                                     cTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTc      <  1:1251050/69‑1 (MQ=255)
                                          tCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATg  >  1:2344058/1‑68 (MQ=255)
                                                         |                                                    
CGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATG  >  minE/30657‑30766

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: