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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL60642,8190A→G12.5% 2.1 ‑3.7 0.070pseudogene (158/885 nt)caiBcrotonobetainyl‑CoA:carnitineCoA‑transferase; b0038_1
Reads supporting (aligned to +/- strand):  new base (10/0):  ref base (37/22):  total (48/22)
Fisher's exact test strand distribution p-value = 2.44e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GGTCAGGCCAGAAAAGTAATCGGCGGTATACGGGAAGGCAGGCATTGGCTGGTCAACATCACCGTTCTGA  >  REL606/42785‑42854
                                  |                                   
ggTCAGGCCAGAAAAGTAATCGGCGGTATACGGGaa                                    <  1:5169483/36‑1 (MQ=255)
ggTCAGGCCAGAAAAGTAATCGGCGGGATACGGGaa                                    >  1:3617135/1‑36 (MQ=255)
 gTCAGGCCAGAAAAGTAATCGGCGGGATACGGGGag                                   >  1:6957672/1‑36 (MQ=37)
 gTCAGGCCAGAAAAGTAATCGGCGGGATAAGGGGag                                   >  1:3566850/1‑36 (MQ=25)
  tCAGGCCAGAAAAGTAATCGGCGGTATTTGGGAtgg                                  >  1:7557966/1‑36 (MQ=25)
   cAGGCCAGAAAAGTAATCGGCGGTATACGGGGaggg                                 >  1:1110264/1‑35 (MQ=255)
   cAGGCCAGAAAAGTAATCGGCGGTATACGGGAaggg                                 >  1:5523732/1‑35 (MQ=255)
   cAGGCCAGAAAAGTAATCGGCGGTATACGGGAaggg                                 >  1:7309080/1‑35 (MQ=255)
    aGGCCAGAAAAGTAATCGGCGGTATACGGGGaggga                                >  1:3240494/1‑34 (MQ=37)
    aGGCCAGAAAAGTAATCGGCGGTATACGGGGaggcg                                >  1:6356733/1‑35 (MQ=255)
    aGGCCAGAAAAGAAATCGGCGGTGTAAGGGAaggga                                >  1:5056441/1‑34 (MQ=21)
     ggCCAGAAAAGTAATCGGCGGTATACGGGAagggag                               >  1:7834122/1‑36 (MQ=255)
      gCCAGAAAAGTAATCGGCGGTATTCCGGGgggcagg                              >  1:5013563/1‑36 (MQ=19)
      gCCAGAAAAGTAATCGGCGGTATTAGGGAagggagg                              >  1:4749539/1‑36 (MQ=21)
      gCCAGAAAAGTAATCGGCGGTATACGGGGaggcagg                              >  1:6717492/1‑36 (MQ=255)
      gCCAGAAAAGTAATCGGCGGTATACGGGAaggcagg                              <  1:7165625/36‑1 (MQ=255)
      gCCAGAAAAGTAATCGGCGGCATACGGGAagggagg                              >  1:622700/1‑36 (MQ=37)
       ccAGAAAAGTAATCGGCGGTATACGGGAaggcaggc                             <  1:1849885/36‑1 (MQ=255)
       ccAGAAAAGTAATCGGCGGTATAAGGGAagggaggg                             >  1:3443746/1‑35 (MQ=37)
        cAGAAAAGTAATCGGCGGTATACGGGGagggaggga                            >  1:3595079/1‑34 (MQ=37)
        cAGAAAAGTAATCGGCGGTATACGGGAagggaggga                            >  1:5092546/1‑34 (MQ=38)
        cAGAAAAGTAATCGGCGGTATAAGGGGagggagggg                            >  1:2986510/1‑34 (MQ=25)
          gAAAAGTAATCGGCGGTGTACGGGTCggctggcaag                          >  1:396123/1‑34 (MQ=16)
          gAAAAGTAATCGGCGGTATCCGGGAAGGCAGGCAtt                          >  1:6726877/1‑36 (MQ=255)
          gAAAAGTAATCGGCGGTATACGGGAAGGCAGGCAtt                          <  1:5710150/36‑1 (MQ=255)
          gAAAAGTAATCGGCGGTATACGGGAAGGCAGGCAtt                          <  1:4463128/36‑1 (MQ=255)
           aaaaGTAATCGGCGGTATACGGGAAGGCAGGCATTg                         <  1:1339706/36‑1 (MQ=255)
           aaaaGTAATCGGCGGTATACGGGAAGGCAGGCATTg                         <  1:2833144/36‑1 (MQ=255)
            aaaGTAATCGGCGGTGTATGGGGAGGGAGGGCTTgg                        >  1:4953082/1‑36 (MQ=255)
            aaaGTAATCGGCGGTATACGGGAAGGGAGGCGTTgg                        >  1:1283224/1‑36 (MQ=38)
            aaaGTAATCGGCGGTATAAGGGAAGGCAGGCATTgg                        >  1:1979955/1‑36 (MQ=255)
             aaGTAATCGGCGGTATACGGGAAGGCAGGCATTGGc                       <  1:3267947/36‑1 (MQ=255)
             aaGTAATCGGCGGTATACGGGAAGGCAGGCATTGGc                       <  1:4978870/36‑1 (MQ=255)
             aaGTAATCGGCGGTATACGGGAAGGCAGGCATTGGc                       <  1:4691002/36‑1 (MQ=255)
               gTAATCGGCGGTATACGGGAAGGGAGGCATTGGGTg                     >  1:531825/1‑36 (MQ=38)
               gTAATCGGCGGTATACGGGAAGGCAGGCATTGGGTg                     >  1:2701976/1‑36 (MQ=255)
                tAATCGGCGGTATACGGGAAGGCAGGAATTGGGTgg                    >  1:1086800/1‑36 (MQ=38)
                 gaTCGGCGGTATAGGGGAAGGCAGGCATTGGCTGGt                   <  1:1364007/35‑1 (MQ=255)
                 aaTCGGCGGTATACGGGAAGGCAGGCATTGGCTggg                   >  1:304829/1‑35 (MQ=255)
                  aTCGGCGTTATACGGGAAGGCAGGCATTGGCTGGTc                  <  1:2016068/36‑1 (MQ=255)
                   tCGGCGGTATACGGGAAGGCAGGCATTGGCTGGTCa                 <  1:5419428/36‑1 (MQ=255)
                    cggaggTATACGGGAAGGCAGGCATTGGCTGGTCaa                >  1:2825987/1‑36 (MQ=255)
                     ggcggTATACGGGAAGGCAGGCATTGGCTGGTGAAc               >  1:457478/1‑36 (MQ=255)
                     ggcggTATACGGGAAGGCAGGCATTGGCTGGTCAAc               <  1:4358042/36‑1 (MQ=255)
                     ggcggTATACGGGAAGGCAGGCATTGGCTGGTCAAc               <  1:6559340/36‑1 (MQ=255)
                      gcggTATACGGGAAGGCAGGCATTGGCTGGTCAACa              <  1:866498/36‑1 (MQ=255)
                      gcggTATACGGGAAGGCAGGCATTGGCTGGTCAACa              >  1:7042235/1‑36 (MQ=255)
                       cggTATACGGTAAGGCAGGCATTGGCTGGTCAACAt             <  1:5951850/36‑1 (MQ=255)
                        gggataCGGGAAGGCAGGCATTGGCTGGGGAACATc            >  1:4550854/4‑36 (MQ=25)
                        ggTATACGGGAAGGCAGGCATTGGCTGGTCAACATc            >  1:833771/1‑36 (MQ=255)
                        ggTATACGGGAAGGCAGGCATTGGCTGGTCAACATc            >  1:3980314/1‑36 (MQ=255)
                        ggTATACGGGAAGGCAGGCATTGGCTGGTCAACATc            <  1:4972483/36‑1 (MQ=255)
                        ggTATACGGGAAGGCAGGCATTGGCTGGTCAAAAtt            >  1:4303708/1‑35 (MQ=255)
                         gTATACGGGAAGGCAGGCATTGGCTGGTCAACATCa           >  1:4061460/1‑36 (MQ=255)
                         gTATACGGGAAGGCAGGCATTGGCTGGTCAACATCa           >  1:6302271/1‑36 (MQ=255)
                         gTATACGGGAAGGCAGGCATTGGCTGGTCAACATCa           >  1:6546386/1‑36 (MQ=255)
                         gTATACGGGAAGGCAGGCATTGGCTGGGCAACATCa           >  1:1603639/1‑36 (MQ=255)
                          tataCGGGAAGGCAGGCATTGGCTGGTCAACATCAc          >  1:2785202/1‑36 (MQ=255)
                           ataCGGGAAGGCAGGCATTGGCTGGTCAACATCAcc         <  1:7600811/36‑1 (MQ=255)
                            taCGGGAAGGCAGGCATTGGCTGGTCAACATCACCg        >  1:1804488/1‑36 (MQ=255)
                            taCGGGAAGGCAGGCATTGGCTGGGCAACATCACCg        >  1:269000/1‑36 (MQ=255)
                             aCGGGAAGGCAGGCATTGGCTGGTCAACATCACCGt       >  1:6366540/1‑36 (MQ=255)
                             aCGGGAAGGCAGGCATTGGCTGGTCAACATCACCGt       <  1:752384/36‑1 (MQ=255)
                             aCGGGAAGGCAGGCATTGGCTGGTCAACATCAACGt       >  1:6531864/1‑36 (MQ=255)
                             aCGGGAAGGCAGGCATTGGCTGGTCAACAGCACCGt       >  1:1213297/1‑36 (MQ=255)
                             aCGGGAAGGCAGGCATTGGCTGGTCAAAATCACCGt       >  1:5288364/1‑36 (MQ=255)
                              cGGGAAGGCAGGCATTGGCTGGTCAACATCACCGtt      >  1:2661801/1‑36 (MQ=255)
                              cGGGAAGGCAGGCATTGGCTGGTCAACATCACCGtt      <  1:5454397/36‑1 (MQ=255)
                               gggAAGGCAGGCATTGGCTGGTCAACATCTACGttt     >  1:1776714/1‑35 (MQ=37)
                               gggAAGGCAGGCATTGGCTGGTCAACATCACCGTTc     <  1:3653533/36‑1 (MQ=255)
                                ggAAGGCAGGCATTGGCTGGTCACCATCACCGTTCt    <  1:4448531/36‑1 (MQ=255)
                                ggAAGGCAGGCATTGGCTGGTCAACATCACCGTTCt    >  1:6628354/1‑36 (MQ=255)
                                  aaGGCAGGCATTGGCTGGTCAACATCACCGTTCTGa  <  1:4329447/36‑1 (MQ=255)
                                  aaGGCAGGCATTGGCTGGTCAACATCACCGTTCTGa  >  1:2227080/1‑36 (MQ=255)
                                  aaGGCAGGCATTGGCTGGTCAACATCACCGTTCTGa  <  1:7685716/36‑1 (MQ=255)
                                  aaGGCAGGCATTGGCTGGGCAACATCACCCTTCTGa  >  1:2304336/1‑36 (MQ=37)
                                  |                                   
GGTCAGGCCAGAAAAGTAATCGGCGGTATACGGGAAGGCAGGCATTGGCTGGTCAACATCACCGTTCTGA  >  REL606/42785‑42854

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