Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 856,653 G→A 52.9% G239G (GGC→GGT trpB ← tryptophan synthase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE856,6530GA52.9% 1.3 / 59.4 68G239G (GGC→GGTtrpBtryptophan synthase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (9/23);  new base A (18/18);  total (27/41)
Fisher's exact test for biased strand distribution p-value = 8.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCC  >  minE/856587‑856719
                                                                  |                                                                  
aCCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATTAAATCAGCAAACATGcc                                                                  <  1:1833646/69‑1 (MQ=255)
aCCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGcc                                                                  >  1:1377666/1‑69 (MQ=255)
aCCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGcc                                                                  >  1:715721/1‑69 (MQ=255)
aCCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGcc                                                                  >  1:3275916/1‑69 (MQ=255)
     gACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATgg                                                             <  1:240635/69‑1 (MQ=255)
     gACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATgg                                                             <  1:3001422/69‑1 (MQ=255)
     gACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATgg                                                             <  1:2060699/69‑1 (MQ=255)
      accaccTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGc                                                            >  1:3254878/1‑69 (MQ=255)
       ccaccTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGcc                                                                  <  1:2331595/62‑1 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATTGCATTCGAaccacc                                                >  1:534816/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                                >  1:3331432/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                                >  1:2719830/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                                >  1:2380946/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                                >  1:2976422/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                                >  1:3283702/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                                >  1:3307522/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                                >  1:419462/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                                >  1:499496/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                                >  1:1527345/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                                >  1:592521/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                                >  1:1387417/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                                >  1:1308298/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                                >  1:1292320/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                                >  1:1248313/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                                >  1:8880/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAAccac                                                 >  1:303790/1‑66 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGAAAACATACCGATGGCATTCGAaccacc                                                >  1:995339/1‑69 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCa                                                           <  1:849728/53‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCa                                                           <  1:1549075/53‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                                           <  1:1598035/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                                           <  1:1156030/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                                           <  1:694425/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                                           <  1:615313/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                                           <  1:3397863/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                                           <  1:1617030/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                                           <  1:3140085/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                                           <  1:3010050/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                                           <  1:2976296/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                                           <  1:297942/69‑1 (MQ=255)
                        aaTCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                                           <  1:1714576/68‑1 (MQ=255)
                        aaTCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                                           <  1:2824501/68‑1 (MQ=255)
                                       ggTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGc                           >  1:552548/1‑69 (MQ=255)
                                          ttCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCa                          <  1:1881459/67‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:832166/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:732202/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:1453933/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:2975899/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:656838/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:586174/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:1503790/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:2123353/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:1161347/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:451216/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:431042/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:2761032/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:889188/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:941219/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:1633354/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:168317/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:3274273/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc                    <  1:3225144/68‑1 (MQ=255)
                                                   gaAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGTCAGGCg               <  1:436658/69‑1 (MQ=255)
                                                   gaAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCg               <  1:2302695/69‑1 (MQ=255)
                                                   gaAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCg               <  1:1463463/69‑1 (MQ=255)
                                                                 tGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTcc  >  1:369915/1‑68 (MQ=255)
                                                                 tGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTcc  >  1:2160312/1‑68 (MQ=255)
                                                                 tGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTcc  >  1:2805110/1‑68 (MQ=255)
                                                                 tGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTc   >  1:498670/1‑67 (MQ=255)
                                                                  |                                                                  
ACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCC  >  minE/856587‑856719

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: