Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 300,681 +C 6.7% coding (702/1473 nt) betB ← betaine aldehyde dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606300,6811.→C 6.7% 6.4 ‑2.7 ‑0.045coding (702/1473 nt)betBbetaine aldehyde dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  new base (0/3):  ref base (26/16):  total (26/19)
Fisher's exact test strand distribution p-value = 6.83e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.92e-01

AGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTG.CCGCTGGCGACACCGCCGGTAAATGACACTTTGG  >  REL606/300646‑300715
                                    |                                  
aGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTg                                     >  1:5131623/1‑36 (MQ=255)
aGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTg                                     >  1:4627374/1‑36 (MQ=255)
aGGGAAGAGGCCGCCGAGCTAGCCATCACTTTTTTg                                     >  1:4898874/1‑36 (MQ=255)
 gggAAGAGGCCGCCGAGTTAGCCATCACTTTTTTG.c                                   >  1:5485663/1‑36 (MQ=255)
  ggAAGAGGCCGCCGAGTTAGCCATCACTTTTTTG.cc                                  <  1:3039040/36‑1 (MQ=255)
    aaGAGGCCGCCGAGTTAGCCATCACTTTTTTG.CCGc                                >  1:7417344/1‑36 (MQ=255)
     agagGCCGCCGAGTTAGCCATCACTTTTTTG.CCGCt                               <  1:4622599/36‑1 (MQ=255)
     agagGCCGCCGAGTTAGCCATCACTTTTTTG.CCGCt                               <  1:7029933/36‑1 (MQ=255)
      gagGCCGCCGAGTTAGCAATCACTTTTTTG.CCGCTg                              >  1:3418951/1‑36 (MQ=255)
       agtccgccgAGTTAGCCATCACTTTTTTG.CCGCTgg                             <  1:3457973/33‑1 (MQ=255)
        ggCCGCCGAGTTAGCCATCACTTTTTTG.CCGCTGGc                            <  1:1611910/36‑1 (MQ=255)
        ggCCGCCGAGTTAGCCATCACTTTTTTG.CCGCTGGc                            >  1:6975739/1‑36 (MQ=255)
         gccgccGAGTTAGCCATCACTTTTTTG.CCGCTGGCg                           <  1:5271135/36‑1 (MQ=255)
         gccgccGAGTTAGCCATCACTTTTTTG.CCGCTGGCg                           >  1:6896588/1‑36 (MQ=255)
         gccgccGAGTTAGCCATCACTTTTTTG.CCGCTGGCg                           >  1:3336813/1‑36 (MQ=255)
          ccgccgAGTTAGCCATAACTTTTTTG.CAGCTGGCGa                          >  1:2053025/1‑36 (MQ=37)
           cgccgAGTTAGCCATCACTTTTTTG.CCGCTGGCGac                         <  1:6768894/36‑1 (MQ=255)
            gccgAGTTAGCCATCACTTTTTTG.CCGCTGGCGaca                        <  1:6929967/36‑1 (MQ=255)
            gccgAGTTAGCCATCACTTTTTTG.CCGCTGGCGaca                        >  1:5193686/1‑36 (MQ=255)
             ccgAGTTAGCCATCACTTTTTTG.CCGCTGGCGacac                       <  1:2480416/36‑1 (MQ=255)
             ccgAGTTAGCCATCACTTTCTTG.CCGCTGGCGacac                       >  1:645843/1‑36 (MQ=255)
              cgAGTTAGCCATCACTTTTTTG.CCGCTGGCGACAcc                      <  1:4985649/36‑1 (MQ=255)
               gAGTTAGCCATCACTTTTTTG.CCGCTGGCGACAccg                     >  1:6550686/1‑36 (MQ=255)
                aGTTAGCCATCACTTTTTTG.CCGCTGGCGACAccgc                    >  1:7433508/1‑36 (MQ=255)
                 gTTAGCCATCACTTTTTTG.CCGCTGGCGACAccgcc                   <  1:7359768/36‑1 (MQ=255)
                  ttAGCCATCACTTTTTTG.CCGCTGGCGACAccgccg                  >  1:6517889/1‑36 (MQ=255)
                  ttAGCCATCACTTTTTTG.CCGCTGGCGACAccgccg                  >  1:5508475/1‑36 (MQ=255)
                  ttAGCCATCACTTTTTTG.CCGCTGGCGACAccgccg                  >  1:871359/1‑36 (MQ=255)
                  ttAGCCATCACTTTTTTG.CCGCTGGCGACAccgccg                  <  1:5019443/36‑1 (MQ=255)
                    aGCCATCACTTTTTTG.CCGCTGGCGACACCGCCGGt                >  1:5551218/1‑36 (MQ=255)
                     gCCATCACTTTTTTG.CCGCTGGCGACACCGCCGGTa               >  1:783832/1‑36 (MQ=255)
                         tCACTTTTTTG.CCGCTGGCTACACCCCGGGTAAATg           >  1:6085398/1‑36 (MQ=25)
                         tCACTTTTTTG.CCGCTGGCGACACCGCCGGTAAATg           >  1:6836855/1‑36 (MQ=255)
                          ctttttttcC.CCGCGGGCAAGCCCGCCGGTAAATGa          <  1:5105772/33‑1 (MQ=2)
                          cACTTTTTTG.CCGCTGGCCACACCGCCGGTATATGa          >  1:3696402/1‑36 (MQ=37)
                           cattttttGCCCGCGGGCGACACCGCCGGTAAATGa          <  1:1159618/34‑1 (MQ=25)
                           aCTTTTTTG.CCGCTGGCGACCCCGCCGGTAAATGac         >  1:3399890/1‑36 (MQ=255)
                            cTTTTTGC.CCGCGGGCGACACCGCCGGTAAATGaca        <  1:7619647/36‑1 (MQ=37)
                             ttttttGCCCGCTGGCGACACCGCCGGTAAATGaca        <  1:5914831/36‑1 (MQ=37)
                             tttttgC.CCGCTGGCACCACCGCCGGTAAATGacac       <  1:7333778/36‑1 (MQ=16)
                               ttttGCCCGCGGGCAACCCCGCCGGTAAATGACACt      <  1:446447/36‑1 (MQ=16)
                               ttttG.CCGCTGGCGACACCGCCGGTAAATGACACtt     >  1:4833376/1‑36 (MQ=255)
                               ttttG.CCGCTGGCGACACCGCCGGTAAATGACACtt     >  1:1478264/1‑36 (MQ=255)
                               ttttG.CCGCTGGCGACACCGCCGGTAAATGACACtt     >  1:131743/1‑36 (MQ=255)
                                tttG.CCGCTGGCGACACCGCCGGTAAATTACACttt    >  1:3764801/1‑36 (MQ=255)
                                tttC.CCGCTGGAAACACCGCCGGTAAATGACACttt    <  1:4125830/36‑1 (MQ=25)
                                tgtG.CCGCGGGCAACACCGCCGGTAAATGACACttt    <  1:584181/34‑1 (MQ=25)
                                  tG.CCGCTGGCGACACCGCCGGTAAATGACACTTTgg  >  1:4515222/1‑36 (MQ=255)
                                  tG.CCGCTGGCGACACCGCCGGTAAATGACACTTTgg  >  1:6927993/1‑36 (MQ=255)
                                  tG.CCGCTGGCGACACCGCCGGTAAATGACACTTTgg  >  1:4064036/1‑36 (MQ=255)
                                    |                                  
AGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTG.CCGCTGGCGACACCGCCGGTAAATGACACTTTGG  >  REL606/300646‑300715

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 15 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG