Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,542 +G 8.7% coding (743/933 nt) thrB → homoserine kinase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,5421.→G 8.7% 4.3 ‑1.5 ‑0.623coding (743/933 nt)thrBhomoserine kinase
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (9/12):  total (11/12)
Fisher's exact test strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.47e-01

ACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC.GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCT  >  REL606/3507‑3577
                                    |                                   
aCGGTTACTGCCTGGCTTCcggcaggcgcggcaggc                                      <  1:1543474/36‑1 (MQ=255)
aCGGTTACTGCCTGGCTTCcggcaggcgcggcaggc                                      <  1:287278/36‑1 (MQ=255)
 cGGTTACTGCCTGGCGTCcggcaggcgcggcaggc.g                                    <  1:5336437/36‑1 (MQ=255)
  ggTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGG.gg                                   >  1:952617/1‑36 (MQ=255)
   gTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC.GGt                                  <  1:4417403/36‑1 (MQ=255)
    ttACTGCCTGGCTTCCGGCAGGCGCGGCAGGC.gggt                                 >  1:2468670/1‑34 (MQ=255)
     tACTGCCTGGCTTCCGGCAGGCGCGGCAGGC.GGGTg                                >  1:1383361/1‑36 (MQ=255)
       cTGCCTGGCTTCCGGCAGGCGCGGCAGGC.GGTTGCg                              <  1:894939/36‑1 (MQ=255)
       cTGCCTGGCTTCCGGCAGGCGCGGCAGGC.GGTTGCg                              <  1:6535831/36‑1 (MQ=255)
       cTGCCTGGCTTCCGGCAGGCGCGGCAGGC.GGTTGCg                              <  1:238778/36‑1 (MQ=255)
       cTGCCTGGCTTCCGGCAGGCGCGGCAGGC.GGGTggg                              >  1:3780058/1‑34 (MQ=37)
            tGGCTTCCGGCAGGCGCGGCAGGC.GGTTGCGGAAAt                         <  1:32313/36‑1 (MQ=255)
            tGGCTTCCGGCAGGCGCGGCAGGC.GGGTGGGGGAAt                         >  1:5846467/1‑36 (MQ=25)
             ggCTTCCGGCAGGCGCGGCAGGC.GGTTGCGGAAATc                        <  1:5451700/36‑1 (MQ=255)
             ggCTTCCGGCAGGCGCGGCAGGC.GGGTGGGGgaaac                        >  1:4280125/1‑34 (MQ=21)
             ggCTTCCGGCAGGCGCGGCAGGC.GGGTGCGGgaaaa                        >  1:539110/1‑34 (MQ=25)
              gCTTCCGGCAGGCGCGGCGGGC.GGGTGGGGGAATCg                       >  1:1297356/1‑36 (MQ=19)
              gCTTCCGGCAGGCGCGGCAGGCGGGTTGGGGaatag                        >  1:7942344/1‑33 (MQ=16)
               cTTCCGGCAGGCGCGGCAGGCGGGTTGGGGaaacgg                       >  1:2848470/1‑33 (MQ=16)
               cTTCCGGCAGGCGCGGCAGGC.GGTTGCGGAAATCgg                      <  1:3274279/36‑1 (MQ=255)
               cTTCCGGCAGGCGCGGCAGGC.GGTTGCGGAAATCgg                      <  1:1992498/36‑1 (MQ=255)
                ttCCGTCATGCGCGGCAGGC.GGTTGCGGAAATCGgc                     <  1:2380391/36‑1 (MQ=37)
                ttCCGGCAGGCGCGGCAGGC.GGTTGCGGAAATCGgc                     <  1:4483872/36‑1 (MQ=255)
                ttCCGGCAGGCGCGGCAGGC.GGTTGCGGAAATCGgc                     <  1:6424778/36‑1 (MQ=255)
                ttCCGGCAGGCGCGGCAGGC.GGGTGGGGGGATCggg                     >  1:4449223/1‑35 (MQ=16)
                ttCCGGCAGGCGCGGCAGGC.GGGTGCGGGAAACggg                     >  1:4781693/1‑35 (MQ=21)
                 tCCGGCAGGCGCGGCAGGC.GGTTGCGGGAAACGgcg                    >  1:8446295/1‑36 (MQ=37)
                 tCCGGCAGGCGCGGCAGGC.GGTTGCGGAAATCGgcg                    <  1:5505131/36‑1 (MQ=255)
                   cggcaggcgcggcaggc.gGTTGCGGAAATCGGcgcg                  <  1:892911/36‑1 (MQ=255)
                     gcaggcgcggcaggc.gGTTGGGGGAAATGGCGCggg                >  1:8017450/1‑35 (MQ=16)
                     gcaggcgcggcaggc.gGTTGCGGAAATCGGCGCGGt                <  1:6478288/36‑1 (MQ=255)
                      caggcgcggcaggc.gGTTGCGGAAATCGGCGCGGTa               <  1:5281201/36‑1 (MQ=255)
                       aggcgcggcaggc.gGTTGCGGGAATCGGGGCggggg              >  1:6301744/1‑33 (MQ=21)
                       aggcgcggcaggc.gGTTGCGGAAATCGGCGCGGTag              <  1:7052962/36‑1 (MQ=255)
                       aggcgcggcaggc.gGTTGCGGAAATCGGCGCGGTag              <  1:5670351/36‑1 (MQ=255)
                         gcgcggcaggc.gGTTTCGGAAATCGGCGGGGGaggg            >  1:700336/1‑34 (MQ=21)
                            cggcaggc.gGTTGCGGAAATCGGCGCGGTagcgagc         <  1:6317741/36‑1 (MQ=255)
                            cggcaggc.gGTTGCGGAAATCGGCGCGGTagcgagc         <  1:7189300/36‑1 (MQ=255)
                             ggcaggc.gGTTGCGGAAATCGGCGCGGTagcgagcg        <  1:5849613/36‑1 (MQ=255)
                              gcaggc.gGTTGCGGAAATCGGCGCGGTAGCGAGCgg       <  1:1905712/36‑1 (MQ=255)
                               caggc.gGTTGCGGAAATCGGCGCGGTAGCGAGCGGt      <  1:2265936/36‑1 (MQ=255)
                                aggc.gGTTGCGGAAATCGGCGCGGTAGCGAGCGGTa     >  1:5010869/1‑36 (MQ=255)
                                  gc.gGTTGCGGAAATCGGCGCGGTAGCGAGCGGTAtc   >  1:6167622/1‑36 (MQ=255)
                                  gc.gGTTGCGGAAATCGGCGCGGTAGCGAGCGGTAtc   <  1:6575446/36‑1 (MQ=255)
                                  gc.gGTTGCGGAAATCGGCGCGGTAGCGAGCGGTAtc   >  1:2124869/1‑36 (MQ=255)
                                  gc.gGTTGCGGAAATCGGCGCGGTAGCGAGCGGTAtc   >  1:8098713/1‑36 (MQ=255)
                                  gc.gGTTGCGGAAATCGGCGCGGTAGCGAGCGGTAtc   <  1:1613491/36‑1 (MQ=255)
                                   c.gGTTGCGGAAATCGGCGCGGTAGCGAGCGGTAtct  <  1:2004746/36‑1 (MQ=255)
                                    |                                   
ACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC.GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCT  >  REL606/3507‑3577

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG