Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,795,647 T→C 13.7% D40G (GAT→GGT)  celF ← cryptic phospho‑beta‑glucosidase, NAD(P)‑binding

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,795,6470T→C13.7% 2.1 ‑4.1 ‑0.036119celFcryptic phospho‑beta‑glucosidase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (19/12):  total (19/17)
Fisher's exact test strand distribution p-value = 1.64e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.23e-01

AATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGA  >  REL606/1795613‑1795682
                                  |                                   
aataatATCCAGTTTCGCTTTACCACCTTCGACATc                                    >  1:5218458/1‑36 (MQ=255)
aataatATCCAGTTTCGCTTTACCACCTTCGACATc                                    <  1:4655583/36‑1 (MQ=255)
 ataataTCCAGTTTCGCTTTACCACCTTCGACATcc                                   <  1:839608/36‑1 (MQ=255)
 ataataTCCAGTTTCGCTTTACCACCTTCGACATcc                                   >  1:2976653/1‑36 (MQ=255)
  taataTCCAGTTTCGCTTTACCACCTTCGACATcca                                  >  1:1470563/1‑36 (MQ=255)
  taataTCCAGTTTCGCTTTACCACCTTCGACATcca                                  >  1:5248772/1‑36 (MQ=255)
    atatCCAGTTTCGCTTTACCACCTTCGACATccacc                                >  1:4238444/1‑36 (MQ=255)
     tatCCAGTTTCTCTTTACCCCCTTCGACATccacca                               <  1:5353252/36‑1 (MQ=37)
       tCCAGTTTCGCTTTACCACCTTCGACATCCACCAGc                             <  1:2249482/36‑1 (MQ=255)
        ccAGTTTCGCTTTACCACCTTCGACATCCACCAGcc                            <  1:852836/36‑1 (MQ=255)
        ccAGTTTCGCTTTACCACCTTCGACATCCACCAGcc                            <  1:733799/36‑1 (MQ=255)
        ccAGTTTCGCTTTACCACCTTCGACATCCACCAGcc                            <  1:3863594/36‑1 (MQ=255)
        ccAGTGCCGCTTTACCCCCTTCGACATCCACCAGcc                            <  1:54067/36‑1 (MQ=25)
          aGTTCCCCTTTACCACCTTCGACATCCACCAGCCAt                          <  1:677688/36‑1 (MQ=37)
           gTTTCGCTTTACCACCTTCGACATCCACCAGCCATa                         >  1:1311573/1‑36 (MQ=255)
              tCGCTTTACCACCTTCGACATCCACCAGCCATAAtt                      >  1:1001242/1‑36 (MQ=255)
              tCGCTTTACCACCTTCGACATCCACCAGCCATAAtt                      >  1:5411477/1‑36 (MQ=255)
               cgttttACCTCCTCCGCGACCCCCCAGCCATAATTc                     <  1:5344418/33‑1 (MQ=2)
               cGCTTTACCACCTCCGACATCCCCCAGCCATAATTc                     <  1:5271415/36‑1 (MQ=37)
                gCTTTACCACCTTCGACATCCACCAGCCATAATTCg                    >  1:2902454/1‑36 (MQ=255)
                ccTTTCCCCCCTTCAACATCCACCAGCCATAATTCg                    <  1:4205300/35‑1 (MQ=21)
                 cTTTACCACCTTCGACATCCACCAGCCATAATTCGc                   >  1:4271299/1‑36 (MQ=255)
                 cTTTACCACCTTCGAAATCCACCAGCCATAATTCGc                   >  1:45739/1‑36 (MQ=255)
                  tttACCCCCTCCAACATCCCCCAGCCATAATTCGCt                  <  1:4541114/36‑1 (MQ=19)
                  tttACCACCTTCGACATCCACCAGCCATAATTCGCt                  >  1:4076744/1‑36 (MQ=255)
                  tttACCACCTTCGACATCCACCAGCCATAATTCGCt                  >  1:2223031/1‑36 (MQ=255)
                  ttaaccaccTTCCACACACACCACCCATAATTCGCt                  <  1:4207482/33‑1 (MQ=16)
                    tACCACCTTCGACATCCACCAGCCATAATTCGCTGa                >  1:3759325/1‑36 (MQ=255)
                     cccaccTACCACACCCCCCAGCCATAATTCGCTGAc               <  1:4157952/35‑1 (MQ=19)
                     accaccTTCGACATCCACCAGCCATAATTCGCTGAc               >  1:1189847/1‑36 (MQ=255)
                          ctccGACACCCCCCAGCCATAATTCGCTGACCGGCa          <  1:5133257/33‑1 (MQ=25)
                             ccaCACCCCCCAGCCATAATTCGCTGACCGGCAAtt       <  1:4674265/34‑1 (MQ=37)
                             cGACATCCACCAGCCATAATTCGCTGACCGGCAAtt       >  1:306606/1‑36 (MQ=255)
                               aCATCCACCAGCCATAATTCGCTGACCGGCAATTCg     >  1:1923730/1‑36 (MQ=255)
                                 aTCCACCAGCCATAATTCGCTGACCGGCAATTCGTg   >  1:1528296/1‑36 (MQ=255)
                                 aTCCACCAGACATCATTCGCTGACCGGCAAAACGTg   >  1:3823421/1‑36 (MQ=19)
                                  tCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGa  <  1:5311191/36‑1 (MQ=255)
                                  tCCACCAGCCATAATTCGCTGACCGGCAATTCGTGa  >  1:3942919/1‑36 (MQ=255)
                                  |                                   
AATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGA  >  REL606/1795613‑1795682

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