Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 42,823 C→G 15.7% pseudogene (154/885 nt) caiB ← crotonobetainyl‑CoA:carnitineCoA‑transferase; b0038_1

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL60642,8230C→G15.7% 5.3 ‑4.6 0.043pseudogene (154/885 nt)caiBcrotonobetainyl‑CoA:carnitineCoA‑transferase; b0038_1
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (16/20):  total (23/20)
Fisher's exact test strand distribution p-value = 1.00e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CAGGCCAGAAAAGTAATCGGCGGTATACGGGAAGGCAGGCATTGGCTGGTCAACATCACCGTTCTGAATCA  >  REL606/42788‑42858
                                   |                                   
aaGGCCAGAAAAGTAATCGGCGGTATACGGGAaggc                                     <  1:5901433/35‑1 (MQ=255)
 aGGCCAGAAAAGTAATCGGCGGTATACGGGAaggca                                    >  1:376822/1‑36 (MQ=255)
 aGGCCAGAAAAGTAATCGGCGGGATAGGGGAaggca                                    >  1:2816458/1‑36 (MQ=37)
   gCCAGAAAAGTAATCGGCGGTATACGGGAaggcagg                                  <  1:3857520/36‑1 (MQ=255)
    ccAGAAAAGTAATCGGCGGTATACGGGAagggaggg                                 >  1:5600863/1‑35 (MQ=255)
    caaGAAAAGTAATCGGCGGTATACGGGAaggcaggc                                 <  1:3326536/34‑1 (MQ=255)
       gAAAAGTAATCGGCGGTATACGGGAAGGGAGGCAtt                              >  1:1512948/1‑36 (MQ=255)
       gAAAAGTAATCGGCGGTATACGGGAAGGCAGGCAtt                              <  1:263622/36‑1 (MQ=255)
       gAAAAGTAATCGGCGGGATACGGGAAGGGAGGCAtt                              >  1:2941238/1‑36 (MQ=38)
        aaaaGTAATCGGCGGTCTAAGGGGAGGCAGGGATTg                             >  1:4755970/1‑36 (MQ=19)
        aaaaGTAATCGGCGGTATACCGGAAGGGAGGCATTg                             >  1:7658490/1‑36 (MQ=38)
         aaaGTAATCGGCGGTATGCGGGGAGGGAGGGAtttg                            >  1:1747967/1‑34 (MQ=16)
          aaGTAATCGGCGGTATACGGGAAGGCAGGCATTGGc                           <  1:5587929/36‑1 (MQ=255)
           aGTAATCGGCGGTATACGGGAAGGCAGGCATTGGCt                          <  1:6999242/36‑1 (MQ=255)
           aGTAATCGGCGGTATACGGGAAGGCAGGCATTGGCt                          <  1:3636282/36‑1 (MQ=255)
            gTAATCGGCGGTATACGGGGAGGGAGGGGTTGGGTg                         >  1:4940165/1‑36 (MQ=255)
            gTAATCGGCGGTATACGGGGAGGGAGGGATTGGGTg                         >  1:3448629/1‑36 (MQ=255)
            gTAATCGGCGGTATACGGGAAGGCAGGCATTGGCTg                         <  1:7386452/36‑1 (MQ=255)
            gTAATCGGCGGTATACGGGAAGGCAGGCATTGGCTg                         <  1:5338708/36‑1 (MQ=255)
            gTAATCGGCGGTATAAGGGGAGGGAGGGATTGGGTg                         >  1:6745150/1‑36 (MQ=16)
               aTCGGCGGTATACGGGAAGGCAGGCATTGGCTGGTc                      <  1:5478306/36‑1 (MQ=255)
                  ggcggTATACGGGAAGGCAGGCATTGGCTGGTCAAc                   <  1:8021647/36‑1 (MQ=255)
                   gcggTATACGGGAAGGCAGGCATTGGCTGGTCCACa                  >  1:4401117/1‑36 (MQ=255)
                   gcggTATACGGGAAGGCAGGCATTGGCTGGTCAACa                  >  1:2456415/1‑36 (MQ=255)
                   gcggTATACGGGAAGGCAGGCATTGGCTGGTCAACa                  <  1:2240607/36‑1 (MQ=255)
                   gcggTATACGGGAAGGCAGGCATTGGCTGGTCAACa                  >  1:4790311/1‑36 (MQ=255)
                   gcggTATACGGGAAGGCAGGCATTGGCTGGTCAACa                  >  1:5125116/1‑36 (MQ=255)
                   gcggTATACGGGAAGGCAGGAGTTGGCTGGTCAAct                  >  1:5645838/1‑35 (MQ=37)
                    cggTATACGGGAAGGCAGGCATTGGCTGGTCAACAt                 <  1:4162452/36‑1 (MQ=255)
                     ggTATACGGGAAGGCAGGCATTGGCTGGTCAACATc                >  1:6659865/1‑36 (MQ=255)
                     ggTATACGGGAAGGCAGGCATTGGCTGGTCAACATc                <  1:859456/36‑1 (MQ=255)
                     ggTATACGGGAAGGCAGGCATTGGCTGGTCAACATc                <  1:864415/36‑1 (MQ=255)
                      gTATACGGGAAGGCAGGCATTGGCTGGTCAACATCa               <  1:313303/36‑1 (MQ=255)
                       tataCGGGAAGGCAGGCATTGGCTGGTCAACATCAc              <  1:8186756/36‑1 (MQ=255)
                        ataCGGGAAGGCAGGCATTGGCTGGTCAAGATCAcc             >  1:3298982/1‑36 (MQ=255)
                        ataCGGGAAGGCAGGCATTGGCTGGTCAACATCAcc             <  1:1178400/36‑1 (MQ=255)
                        ataCGGGAAGGCAGGCATTGGCTGGTCAACATCAcc             <  1:5480697/36‑1 (MQ=255)
                        ataCGGGAAGGCAGGCATTGGCTGGGCAACATCAcc             >  1:5794178/1‑36 (MQ=255)
                         taCGGGAAGGCAGGCATTGGCTGGTCAACATCACCg            >  1:2993049/1‑36 (MQ=255)
                         taCGGGAAGGCAGGCATTGGCTGGGCAACATCACCg            >  1:3019635/1‑36 (MQ=255)
                          aCGGGAAGGCAGGCATTGGCTGGTCAACATCACCGt           >  1:5092310/1‑36 (MQ=255)
                          aCGGGAAGGCAGGCATTGGCTGGTCAACATCACCGt           >  1:7333358/1‑36 (MQ=255)
                           cGGGACGACAGGCATGGGCTGGTCAACATCACCGtt          <  1:1931930/36‑1 (MQ=21)
                           cGGGAAGGCAGGCATTGGCTGGTCAACATCACCGtt          >  1:4293171/1‑36 (MQ=255)
                           cGGGAAGGCAGGCATTGGCTGGTCAACATCACCGtt          <  1:1746205/36‑1 (MQ=255)
                            gggAAGGCAGGCATTGGCTGGTCAACATCACCGTTc         <  1:7446200/36‑1 (MQ=255)
                            gggAAGGCAGGCATTGGCTGGTCAACATCACCGTTc         <  1:4867042/36‑1 (MQ=255)
                             ggAAGGCAGGCATTGGCTGGTCAACATCACCGTTCt        >  1:4142397/1‑36 (MQ=255)
                              gAAGGCAGGCATTGGCTGGTCAACATCACCGTTCTg       >  1:5205229/1‑36 (MQ=255)
                                 ggcaggcaTTGGCTGGTCAACATCACCGTTCTGAAt    <  1:7165329/36‑1 (MQ=255)
                                 ggcaggcaTTGGCTGGTCAACATCACCGTTCTGAAt    >  1:8178750/1‑36 (MQ=255)
                                 ggcaggcaTTGGCTGGACCACATCACCGTTCTGAAt    >  1:4608170/1‑36 (MQ=37)
                                  gcaggcaTTGGCTGGTGAAGATCACCGTTCTGCATc   >  1:1791092/1‑36 (MQ=25)
                                  gcaggcaTTGGCTGGTCAACATCACCGTTCTGAATc   >  1:5573516/1‑36 (MQ=255)
                                  gcaggcaTTGGCTGGTCAACATCACCGTTCTGAATc   <  1:4777875/36‑1 (MQ=255)
                                  gcaggcaTTGGCTGGTCAACATCACCGTTCTGAATc   >  1:3492789/1‑36 (MQ=255)
                                  gcaggcaTTGGCTGGTCAACATCAACGGTCTGAATc   >  1:774333/1‑36 (MQ=37)
                                   caggcaTTGGCTGGTCAACATCACCGTTCTGGATCa  >  1:1821176/1‑36 (MQ=255)
                                   caggcaTTGGCTGGTCAACATCACCGTTCTGAATCa  >  1:8459279/1‑36 (MQ=255)
                                   |                                   
CAGGCCAGAAAAGTAATCGGCGGTATACGGGAAGGCAGGCATTGGCTGGTCAACATCACCGTTCTGAATCA  >  REL606/42788‑42858

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 17 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG