Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 43,125 A→G 21.6% pseudogene (184/420 nt) caiB ← crotonobetainyl‑CoA:carnitineCoA‑transferase; b0038_2

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL60643,1250A→G21.6% 8.1 ‑3.2 0.045pseudogene (184/420 nt)caiBcrotonobetainyl‑CoA:carnitineCoA‑transferase; b0038_2
Reads supporting (aligned to +/- strand):  new base (10/0):  ref base (22/13):  total (32/13)
Fisher's exact test strand distribution p-value = 4.24e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTCGGTTGAACGCGAATGGTG  >  REL606/43092‑43159
                                 |                                  
tatTTAACGACAGCGCGTGCAAATTGCGGCGGggga                                  >  1:2231674/1‑36 (MQ=255)
tatTTAACGACAGCGCGTGCAAATTGCGGCGGgaga                                  <  1:7464751/36‑1 (MQ=255)
tatTTAACGACAGCGCGTGCAAATTGCGGCGGgaga                                  >  1:1554687/1‑36 (MQ=255)
 atTTAACGACAGCGCGTGCAAATTGCGGCGGGggag                                 >  1:7745014/1‑36 (MQ=255)
  tttAACGACAGCGCGTGCAAATTGCGGCGGGAGAGc                                <  1:3108512/36‑1 (MQ=255)
   ttAACGACAGCGCGTGCAAATTGCGGCGGGggaggt                               >  1:1293633/1‑34 (MQ=38)
   ttAACGACAGCGCGTGCAAATTGCGGCGGGGGAGCt                               >  1:1222267/1‑36 (MQ=255)
    tAACGACAGCGCGTGCAAATTGCGGCGGGAGAGCTg                              <  1:2339511/36‑1 (MQ=255)
    tAACGACAGCGCGTGCAAATTGCGGCGGGAGAGCTg                              <  1:405394/36‑1 (MQ=255)
     aaCGACAGCGCGTGCAAATTGCGGCGGGGGAGGTgg                             >  1:671058/1‑35 (MQ=37)
     aaCGACAGCGCGTGCAAATTGCGGCGGGGGAGGTgg                             >  1:2881510/1‑35 (MQ=37)
      aCGACAGCGCGCGCAAATTGCGGCGGGGGAGCTgag                            >  1:7292718/1‑34 (MQ=37)
       cGACAGCGCGTGCAAATTGCGGCGGGAGAGCTGCgg                           >  1:1131079/1‑36 (MQ=255)
       cGACAGCGCGTGCAAATTGCGGCGGGAGAGCTGCgg                           >  1:5343062/1‑36 (MQ=255)
        gACAGCGCGTGCAAATTGCGGCGGGGGAGCTGCggg                          >  1:2983786/1‑36 (MQ=255)
        gACAGCGCGTGCAAATTGCGGCGGGGGAGCTGCggg                          >  1:5820515/1‑36 (MQ=255)
        gACAGCGCGTGCAAATTGCGGCGGGGGAGCTGCggg                          >  1:1023867/1‑36 (MQ=255)
        gACAGCGCGTGCAAATTGCGGCGGGAGAGGTGCggg                          >  1:2904187/1‑36 (MQ=255)
        gACAGCGCGTGCAAATTGCGGCGGGAGAGCTGCggg                          <  1:2758375/36‑1 (MQ=255)
         aCAGCGCGTGCAAATTGCGGCGGGGGAGGTGCgggg                         >  1:4200001/1‑35 (MQ=37)
         aCAGCGCGTGCAAATTGCGGCGGGAGAGCTGAgggg                         >  1:4160406/1‑35 (MQ=255)
          cAGCGCGTGCAAATTGCGGCGGGGGAGCTGCggggc                        >  1:7973139/1‑34 (MQ=38)
          cAGCGCGTGCAAATTGCGGCGGGAGAGCTGCgggga                        >  1:7887207/1‑34 (MQ=255)
           aGCGCGTGCAAATTGCGGCGGGGGAGCTGCggggcg                       >  1:2401345/1‑33 (MQ=37)
              gcgTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTc                    <  1:6169597/36‑1 (MQ=255)
                gTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTCgg                  >  1:1348562/1‑36 (MQ=255)
                  gCAAATTGCGGCGGGAGAGCTGCGGGTAGTTCGGtt                <  1:3461791/36‑1 (MQ=255)
                  gCAAATTGCGGCGGGAGAGCTGCGGGTAGTTCGGtt                <  1:260091/36‑1 (MQ=255)
                   cAAATTGCGGCGGGAGAGCTGCGGGTAGTTCGGTTg               <  1:4168681/36‑1 (MQ=255)
                    aaaTTGCGGCGGGAGAGCTGCGGGTAGTTTGGTTGa              >  1:5353847/1‑36 (MQ=255)
                    aaaTTGCGGCGGGAGAGCTGCGGGTAGTTCGGTTGa              <  1:5498077/36‑1 (MQ=255)
                     aaTTGCGGCGGGAGAGGTGCGGGTTGTTTGGGTgta             >  1:169946/1‑34 (MQ=16)
                     aaTTGCGGCGGGAGAGCTGCGGGGAGTTCGGTTTaa             >  1:928164/1‑36 (MQ=37)
                      aTTGCGGCGGGAGAGCTGCGGGTAGTTCGGTTGAAc            <  1:7107837/36‑1 (MQ=255)
                       ttGCGGCGGGAGAGCTGCGGGGAGTTTGGTTTaaag           >  1:1223776/1‑34 (MQ=21)
                        tGCGGCGGGAGAGCTGCGGGTAGTTCGGTTGAAcgc          >  1:401051/1‑36 (MQ=255)
                         gcggcgGGAGAGCTGCGGGTAGTTCGGTTGAAcgcg         <  1:7035939/36‑1 (MQ=255)
                          cggcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGa        >  1:6099490/1‑36 (MQ=255)
                          cggcggGAGAGCTGCGGGGAGTTCGGTTGAACGCGa        >  1:5664213/1‑36 (MQ=255)
                           ggcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGaa       >  1:5479470/1‑36 (MQ=255)
                           ggcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGaa       <  1:4009081/36‑1 (MQ=255)
                           ggcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGaa       >  1:2513915/1‑36 (MQ=255)
                            gcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAAt      >  1:4955443/1‑36 (MQ=255)
                            gcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAAt      >  1:5391675/1‑36 (MQ=255)
                            gcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAAt      >  1:2250532/1‑36 (MQ=255)
                             cggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAATg     >  1:8146594/1‑36 (MQ=255)
                             cggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAATg     <  1:5339606/36‑1 (MQ=255)
                              gggAGAGCTGCGGGTAGTTCGGTTGAACGCGAATgg    >  1:2588378/1‑36 (MQ=255)
                              gggAGAGCTGCGGGTAGTTCGGTTGAACGCGAATgg    >  1:1186747/1‑36 (MQ=255)
                               ggAGAGCTGCGGGTAGTTCGGGTGAACGCGAATGgt   >  1:2905131/1‑36 (MQ=255)
                                gagaGCTGCGGGTAGTTCGGTTGAACGCGAATGgtg  <  1:4729108/36‑1 (MQ=255)
                                gagaGCTGCGGGTAGTTCGGTTGAACGCGAATGgtg  <  1:5683064/36‑1 (MQ=255)
                                gagaGCTGCGGGTAGTTCGGTTGAACGCGAATGgtg  <  1:4730919/36‑1 (MQ=255)
                                 |                                  
TATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTCGGTTGAACGCGAATGGTG  >  REL606/43092‑43159

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