Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 43,125 A→G 17.5% 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→G17.5% 5.1 ‑2.5 ‑0.145pseudogene (184/420 nt)caiBcrotonobetainyl‑CoA:carnitineCoA‑transferase; b0038_2
Reads supporting (aligned to +/- strand):  new base (8/0):  ref base (23/13):  total (32/13)
Fisher's exact test strand distribution p-value = 8.20e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.77e-01

TATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTCGGTTGAACGCGAATGGTGT  >  REL606/43092‑43160
                                 |                                   
tatTTAACGACAGCGCGTGCAAATTGCGGCGGgaga                                   >  1:7840020/1‑36 (MQ=255)
tatTTAACGACAGCGCGTGCAAATTGCGGCGGgaga                                   <  1:5442039/36‑1 (MQ=255)
tatTTAACGACAGCGCGTGCAAATTGCGGCGGgaga                                   <  1:6978909/36‑1 (MQ=255)
 atTTAACGACAGCGCGTGCAAATTGCGGCGGGggag                                  >  1:3337430/1‑36 (MQ=255)
  tttAACGACAGCGCGTGCAAATTGCGGCGGGGGAGc                                 >  1:4029417/1‑36 (MQ=255)
   ttAACGACAGCGCGTGCAAATTGCGGCGGGAGAGCt                                >  1:6834904/1‑36 (MQ=255)
    tAACGACAGCGCGTGCAAATTGCGGCGGGGGAGCTg                               >  1:5222489/1‑36 (MQ=255)
     aaCGACAGCGCGTGCAAATTGCGGCGGGAGAGCTGc                              <  1:4938438/36‑1 (MQ=255)
      aCGACAGCGCGTGCAAATTGCGGCGGGGGGGCTGCg                             >  1:21011/1‑36 (MQ=38)
      aCGACAGCGCGTGCAAATTGCGGCGGGGGAGCTGCg                             >  1:5638474/1‑36 (MQ=255)
      aCGACAGCGCGTGCAAATTGCGGCGGGAGAGGTGCg                             >  1:6070470/1‑36 (MQ=255)
      aCGACAGCGCGTGCAAATTGCGGCGGGAGAGCTGct                             <  1:4975068/36‑2 (MQ=255)
       cGACAGCGCGTGCAAATTGCGGCGGGTGAGCCGCgg                            >  1:1813915/1‑36 (MQ=38)
       cGACAGCGCGTGCAAATTGCGGCGGGGGAGGCGCgg                            >  1:3485333/1‑36 (MQ=21)
        gACAGCGCGTGCAAATTGCGGCGGTGGAGCCTCggg                           >  1:2659819/1‑36 (MQ=19)
        gACAGCGCGTGCAAATTGCGGCGGGGGAGGTGCggg                           >  1:2732905/1‑36 (MQ=37)
        gACAGCGCGTGCAAATTGCGGCGGGGGAGGTGCggg                           >  1:3290011/1‑36 (MQ=37)
        gACAGCGCGTGCAAATTGCGGCGGGAGAGGTGAggg                           >  1:5432402/1‑36 (MQ=37)
        gACAGCGCGTGCAAATTGCGGCGGGAGAGCTGCggg                           >  1:3374722/1‑36 (MQ=255)
          cAGCGCGTGCAAATTGCGGCGGGAGAGCTGCggggg                         >  1:3871197/1‑34 (MQ=255)
           aGCGCGTGCAAATTGCGGCGGGAGAGCTGCGGGTAg                        >  1:6231787/1‑36 (MQ=255)
             cgcgTGCAAATTGCGGCGGGAGAGCTGCGGGTAGtt                      >  1:577408/1‑36 (MQ=255)
              gcgTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTc                     <  1:4344905/36‑1 (MQ=255)
              gcgTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTc                     <  1:2814275/36‑1 (MQ=255)
              gcgTGCAAATTGCGGCGGGAGAGCTGCGGGGGGTTc                     >  1:164532/1‑36 (MQ=37)
                gTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTTgg                   >  1:6337435/1‑36 (MQ=255)
                  tcAAATTGCGGCGGGAGAGCTGCGGGTAGTTCGGtt                 <  1:3821515/35‑1 (MQ=255)
                  gCAAATTGCGGCGGGAGAGCTGGGGGTAGTTTgggt                 >  1:3862627/1‑34 (MQ=25)
                  gCAAATTGCGGCGGGAGAGCTGGCGGGAGGTCGGtt                 >  1:5264637/1‑36 (MQ=16)
                  gCAAATTGCGGCGGGAGAGCTGCGGGGAGTTCGGtt                 >  1:3366802/1‑36 (MQ=255)
                    aaaTTGCGGCGGGAGAGCTGCGGATAGTTCGGTTGa               <  1:204931/36‑1 (MQ=255)
                      aTTGCGGCGGGAGAGCTGCGGGTAGTTTGGTTGaag             >  1:6229322/1‑35 (MQ=255)
                      aTTGCGGCGGGACAGCTGCGGGTGGTTCGGGTGatc             >  1:1134957/1‑34 (MQ=25)
                       ttGCGGCGGGAGAGCTGCCGGTAGGTCGGGTGAAcg            >  1:7383682/1‑36 (MQ=25)
                         gcggcgGGAGAGCTGCGGGTAGGTTGGTTGAAcgcg          >  1:5444135/1‑36 (MQ=37)
                          cggcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGa         >  1:7776457/1‑36 (MQ=255)
                          cggcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGa         <  1:3623427/36‑1 (MQ=255)
                          cggcggGAGAGCTGCAGGTAGTTCGGTTGAACGCGa         <  1:3151600/36‑1 (MQ=255)
                           ggggggAGTGCTGCGGGTAGTTCGGTTGAACGCGaa        <  1:256805/33‑1 (MQ=38)
                           ggcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGaa        >  1:2959747/1‑36 (MQ=255)
                           ggcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGaa        >  1:2487032/1‑36 (MQ=255)
                            gcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAAt       <  1:5155652/36‑1 (MQ=255)
                            gcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAAt       <  1:2314663/36‑1 (MQ=255)
                            gcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAAt       >  1:1079708/1‑36 (MQ=255)
                             cggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAATg      <  1:5956666/36‑1 (MQ=255)
                              gggAGAGCTGCGGGTAGTTCGGTTGAACGCGAATgg     >  1:5262359/1‑36 (MQ=255)
                              gggAGAGCTGCGGGTAGTTCGGTTGAACGCGAATgg     >  1:4607032/1‑36 (MQ=255)
                               ggAGAGCTGCGGGTAGTTCGGTTTAACGCGAATGgt    >  1:1232837/1‑36 (MQ=255)
                               ggAGAGCTGCGGGTAGTTCGGTTGAACGCGAATGgt    <  1:7872151/36‑1 (MQ=255)
                                 agagCTGCGGGTAGTTCGGTTGAACGCGAATGgtgt  <  1:2298587/36‑1 (MQ=255)
                                 agagCTGCGGGTAGTTCGGTTGAACGCGAATGgtgt  >  1:813008/1‑36 (MQ=255)
                                 |                                   
TATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTCGGTTGAACGCGAATGGTGT  >  REL606/43092‑43160

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