Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,168,136 C→G 22.1% R90G (CGA→GGA)  ‡ pabC → 4‑amino‑4‑deoxychorismate lyase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,168,1360C→G22.1% 3.1 ‑2.5 ‑0.018268pabC4‑amino‑4‑deoxychorismate lyase
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (5/9):  total (9/9)
Fisher's exact test strand distribution p-value = 8.24e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.31e-01

TACTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCAGCAACG  >  REL606/1168101‑1168171
                                   |                                   
tACTGAAAGTCGTGATCAGTCGCGGTAGTggcgggc                                     <  1:810967/36‑1 (MQ=255)
tACTGAAAGTCGTGATCAGTCGCGGTAGTggcgggc                                     <  1:3801357/36‑1 (MQ=255)
      aaGTCGTGATCAGTCGCGGTAGTGGCGGGGGGgggg                               >  1:1954491/1‑35 (MQ=37)
      aaGTCGTGATCAGTCGCGGTAGTGGCGGGCGGgggg                               >  1:1789389/1‑35 (MQ=255)
      aaGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGt                               <  1:33180/36‑1 (MQ=255)
       aGTCGTGATCAGTCGCGGTAGTGGCGGGGGGggggg                              >  1:189200/1‑34 (MQ=25)
       aGTCGTGATCAGTCGCGGTAGTGGCGGGGGGggggg                              >  1:3499161/1‑34 (MQ=25)
       aGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTa                              <  1:5016397/36‑1 (MQ=255)
              aTCAGTCGCGGTAGTGGCGGGCGAGGGTACAGcaca                       <  1:2215117/36‑1 (MQ=255)
               tCAGTCGCGGTAGTGGCGGGCGAGGGTACAGcactt                      <  1:2727851/36‑3 (MQ=255)
                    cgcgGTAGTGGCGGGCTGGTGTACAGCACATTGAAc                 <  1:667143/36‑1 (MQ=25)
                    cgcgGTAGTGGCGGGCGAGGGTACAGCACATTGAAc                 <  1:5130302/36‑1 (MQ=255)
                    cgcgGTAGTGGCGGGCGAGGGTACAGCACATTGAAc                 <  1:1737914/36‑1 (MQ=255)
                     gcgGTAGTGGCGGGGGGGGGGGGAGCACAATGaaaa                >  1:3816148/1‑34 (MQ=9)
                      cgGTAGTGGCGGGCGGGGGGACAGGACAATGAAAAg               >  1:261073/1‑36 (MQ=16)
                          aGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGa           <  1:477321/36‑1 (MQ=255)
                           gTGGCGGGCGTGGGTACAGCCCATTGAACAGCGGAc          >  1:1161171/1‑36 (MQ=37)
                            tGGCGGGCGAGGGTACAGCACATTGATCAGCGGAcc         <  1:3865600/36‑1 (MQ=255)
                            tGGCGGGCGAGGGTACAGCACATTGAACAGCGGAcc         >  1:5363078/1‑36 (MQ=255)
                            tGGCGGGCGAGGGTACAGCACATGGAACAGCGGAcc         >  1:4334752/1‑36 (MQ=255)
                              gcgggcgAGGGTACAGCACATTGAACAGCGGACCAg       <  1:4865052/36‑1 (MQ=255)
                              gcgggcgAGGGTACAGCACATTGAACAGCGGAACAg       >  1:2175603/1‑36 (MQ=255)
                                gggcgAGGGTACAGCACATTGAACAGCGGCCCAGcc     <  1:3714672/36‑2 (MQ=255)
                                gggcgAGGGTACAACACATTGAACAGCGGACCAGCa     >  1:3108442/1‑36 (MQ=255)
                                  gcgAGGGTACAGCACATTGAACAGCGGACCAGCAAc   >  1:549268/1‑36 (MQ=255)
                                  gcgAGGGTACAGCACATTGAACAGCGGACCAGCAAc   >  1:5791935/1‑36 (MQ=255)
                                   cgAGGGTACAGCACATTGAACAGCGGACCAGCAAcg  >  1:5844205/1‑36 (MQ=255)
                                   |                                   
TACTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCAGCAACG  >  REL606/1168101‑1168171

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 16 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG