Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 213,333 C→T 29.6% G92S (GGC→AGC)  gloB ← predicted hydroxyacylglutathione hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE213,3330CT29.6% 83.9 / 36.8 84G92S (GGC→AGC) gloBpredicted hydroxyacylglutathione hydrolase
Reads supporting (aligned to +/- strand):  ref base C (18/41);  new base T (14/11);  total (32/52)
Fisher's exact test for biased strand distribution p-value = 4.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TCCTAAAGTGTGACCCGGCGTAGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTG  >  minE/213268‑213396
                                                                 |                                                               
tCCTAAAGTGTGACCCGGCGTAGCAATTACACTAAATTCATGCCTCAAAACGAAGGCAGTTTCGCCATc                                                              >  1:1802469/1‑69 (MQ=255)
tCCTAAAGTGTGACCCGGCGTAGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATc                                                              >  1:1660221/1‑69 (MQ=255)
tCCTAAAGTGTGACCCGGCGTAGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATc                                                              >  1:791902/1‑69 (MQ=255)
tCCTAAAGTGTGACCCGGCGTAGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATc                                                              >  1:705095/1‑69 (MQ=255)
tCCTAAAGTGTGACCCGGCGTAGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATc                                                              >  1:701190/1‑69 (MQ=255)
tCCTAAAGTGTGACCCGGCGTAGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATc                                                              >  1:2862810/1‑69 (MQ=255)
    aaaGTGTGACCCGGCGTAGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCttt                                                           <  1:937583/68‑1 (MQ=255)
    aaaGTGTGACCCGGCGTAGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCttt                                                           <  1:925011/68‑1 (MQ=255)
    aaaGTGTGACCCGGCGTAGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCttt                                                           <  1:1715914/68‑1 (MQ=255)
       gtgtGACCCGGCGTAGCAATTACACTAAATTTATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACt                                                       <  1:561658/69‑1 (MQ=255)
       gtgtGACCCGGCGTAGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACt                                                       <  1:3058530/69‑1 (MQ=255)
              ccGGTGTGGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCttt                                                           >  1:3120036/1‑58 (MQ=255)
              ccGGTGTGGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCTTTGACTACCtgtg                                                >  1:2502182/1‑69 (MQ=255)
              ccGGTGTGGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCTTTGACTACCtgtg                                                >  1:2778198/1‑69 (MQ=255)
              ccGGTGTGGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCTTTGACTACCtgtg                                                >  1:202543/1‑69 (MQ=255)
              ccGGTGTGGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCTTTGACTACCtgtg                                                >  1:3007148/1‑69 (MQ=255)
              ccGGTGTGGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCTTTGACTACCtgtg                                                >  1:3318207/1‑69 (MQ=255)
              ccGGTGTGGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCTTTGACTACCtgtg                                                >  1:404950/1‑69 (MQ=255)
              ccGGTGTGGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCTTTGACTACCtgtg                                                >  1:633509/1‑69 (MQ=255)
              ccGGTGTGGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCTTTGACTACCtgtg                                                >  1:131182/1‑69 (MQ=255)
              ccGGTGTGGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCTTTGACTACCtgtg                                                >  1:694486/1‑69 (MQ=255)
              ccGGTGTGGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCTTTGACTACCtgtg                                                >  1:1013461/1‑69 (MQ=255)
              ccGGTGTGGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCTTTGACTACCtgt                                                 >  1:2441252/1‑68 (MQ=255)
              ccGGTGTGGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCTTTGACTACCtgt                                                 >  1:2296028/1‑68 (MQ=255)
              ccGGTGTGGCAATTACACTAAATTAATGCCCCAAAACGAAGGCAGTTTCGCTATCTTTGACTACCtgtg                                                >  1:3123746/1‑69 (MQ=255)
                         aTTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTa                                     <  1:163819/69‑1 (MQ=255)
                         aTTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTa                                     <  1:2624449/69‑1 (MQ=255)
                          ttACTCTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTAt                                    <  1:2307290/69‑1 (MQ=255)
                          ttACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTa                                     <  1:2686396/68‑1 (MQ=255)
                          ttACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTa                                     <  1:2852424/68‑1 (MQ=255)
                           tACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATc                                   <  1:2658454/69‑1 (MQ=255)
                           tACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATc                                   <  1:570206/69‑1 (MQ=255)
                            acacTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCt                                  >  1:2810361/1‑69 (MQ=255)
                            acacTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCt                                  >  1:2622411/1‑69 (MQ=255)
                            acacTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCt                                  >  1:981973/1‑69 (MQ=255)
                            acacTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCt                                  >  1:2417301/1‑69 (MQ=255)
                            acacTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCt                                  >  1:2129104/1‑69 (MQ=255)
                            acacTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCt                                  >  1:1445524/1‑69 (MQ=255)
                            acacTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCt                                  >  1:2093948/1‑69 (MQ=255)
                            acacTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCt                                  >  1:1847256/1‑69 (MQ=255)
                              acTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCttt                                                           <  1:2744601/42‑1 (MQ=255)
                              acTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCttt                                                           <  1:2863440/42‑1 (MQ=255)
                              acTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCttt                                                           <  1:162023/42‑1 (MQ=255)
                              acTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCttt                                                           <  1:1573079/42‑1 (MQ=255)
                              acTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCttt                                                           <  1:1549029/42‑1 (MQ=255)
                              acTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCttt                                                           <  1:2834152/42‑1 (MQ=255)
                              acTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCttt                                                           <  1:2829174/42‑1 (MQ=255)
                              acTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCttt                                                           <  1:380100/42‑1 (MQ=255)
                              acTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCttt                                                           <  1:1958585/42‑1 (MQ=255)
                              acTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCttt                                                           <  1:1363325/42‑1 (MQ=255)
                              acTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCTATCttt                                                           <  1:215516/42‑1 (MQ=255)
                                  aaTTCATGCCCCAAAACGAAGGCAGTTTCTCCATCTTTGACTACCTGtgtt                                              >  1:3207162/1‑51 (MQ=255)
                                  aaTTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGtc                            >  1:1121904/1‑69 (MQ=255)
                                  aaTTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCATATCTTGTGtc                            >  1:3243093/1‑69 (MQ=255)
                                    ttCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGtctc                          <  1:1657298/69‑1 (MQ=255)
                                     tCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTctct                         <  1:1279369/69‑1 (MQ=255)
                                     tCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTctct                         <  1:498211/69‑1 (MQ=255)
                                     tCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTctct                         <  1:337795/69‑1 (MQ=255)
                                     tCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTctct                         <  1:3309783/69‑1 (MQ=255)
                                     tCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTctct                         <  1:2845551/69‑1 (MQ=255)
                                      cATGCCCCATAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTctct                         <  1:1362724/68‑1 (MQ=255)
                                      cATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTctct                         <  1:1802339/68‑1 (MQ=255)
                                      cATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCtt                        <  1:3248519/69‑1 (MQ=255)
                                      cATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCtt                        <  1:1452715/69‑1 (MQ=255)
                                      cATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCtt                        <  1:3319180/69‑1 (MQ=255)
                                      cATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCtt                        <  1:1777227/69‑1 (MQ=255)
                                      cATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCtt                        <  1:1927640/69‑1 (MQ=255)
                                      cATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCtt                        <  1:116643/69‑1 (MQ=255)
                                      cATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCtt                        <  1:1207767/69‑1 (MQ=255)
                                      cATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCtt                        <  1:545269/69‑1 (MQ=255)
                                      cATGCCACAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCtt                        <  1:2441703/69‑1 (MQ=255)
                                            ccAAAACGAAGTCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGAc                  <  1:937701/69‑1 (MQ=255)
                                            ccAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGAc                  <  1:1147550/69‑1 (MQ=255)
                                            ccAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGAc                  <  1:1273822/69‑1 (MQ=255)
                                            ccAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGAc                  <  1:3054157/69‑1 (MQ=255)
                                            ccAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGAc                  <  1:1358231/69‑1 (MQ=255)
                                            ccAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGAc                  <  1:2110507/69‑1 (MQ=255)
                                            ccAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGAc                  <  1:37790/69‑1 (MQ=255)
                                            ccAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGAc                  <  1:3285527/69‑1 (MQ=255)
                                            ccAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGAc                  <  1:1769951/69‑1 (MQ=255)
                                              aaaaCGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGAc                  <  1:1189111/67‑1 (MQ=255)
                                                 aCGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTctct                         <  1:388843/57‑1 (MQ=255)
                                                     aGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCtt                        <  1:2603112/54‑1 (MQ=255)
                                                             tCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTtg  >  1:842517/1‑68 (MQ=255)
                                                                 |                                                               
TCCTAAAGTGTGACCCGGCGTAGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTG  >  minE/213268‑213396

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: