Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 856,644 A→T 56.9% A242A (GCT→GCA trpB ← tryptophan synthase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE856,6440AT56.9% 22.3 / 66.6 65A242A (GCT→GCAtrpBtryptophan synthase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (5/23);  new base T (19/18);  total (24/41)
Fisher's exact test for biased strand distribution p-value = 8.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

ACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCG  >  minE/856587‑856706
                                                         |                                                              
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:3001422/69‑1 (MQ=255)
     gACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATgg                                                <  1:240635/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)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAAc                                                          >  1:2709059/1‑46 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATTGCATTCGAaccacc                                   >  1:534816/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                   >  1:2719830/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                   >  1:2380946/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                   >  1:3331432/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:592521/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                   >  1:3283702/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                   >  1:1527345/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                   >  1:2976422/1‑69 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAaccacc                                   >  1:8880/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)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGCAAACATACCGATGGCATTCGAAccac                                    >  1:303790/1‑66 (MQ=255)
                  cacaCCAATCAGGCCGACGTCGGTTTCGTTGATGAAATCTGAAAACATACCGATGGCATTCGAaccacc                                   >  1:995339/1‑69 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCa                                              <  1:1549075/53‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCa                                              <  1:849728/53‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                              <  1:3397863/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                              <  1:615313/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                              <  1:694425/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                              <  1:3140085/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                              <  1:297942/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                              <  1:1156030/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                              <  1:3010050/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                              <  1:2976296/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                              <  1:1617030/69‑1 (MQ=255)
                       cAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                              <  1:1598035/69‑1 (MQ=255)
                        aaTCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                              <  1:2824501/68‑1 (MQ=255)
                        aaTCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCaa                              <  1:1714576/68‑1 (MQ=255)
                                       ggTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGc              >  1:552548/1‑69 (MQ=255)
                                          ttCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCa             <  1:1881459/67‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc       <  1:1633354/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc       <  1:1161347/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc       <  1:656838/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc       <  1:1503790/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc       <  1:732202/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc       <  1:832166/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc       <  1:1453933/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc       <  1:889188/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc       <  1:941219/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc       <  1:586174/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc       <  1:168317/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc       <  1:451216/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc       <  1:431042/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc       <  1:2123353/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc       <  1:2761032/68‑1 (MQ=255)
                                               tgatgaAATCTGCAAACATACCGATGGCATTCGAACCACCGCCAACACAGGCGATAACGGCATCCGGc       <  1:2975899/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)
                                                         |                                                              
ACCGTGACCACCTGGCTCCACACCAATCAGGCCGACGTTGGTTTCATTGATGAAATCAGCAAACATGCCGATGGCATTCGAACCGCCGCCAACACAGGCGATAACGGCATCCGGCAGGCG  >  minE/856587‑856706

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: