Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,741,525 +G 6.7% coding (748/924 nt) hycD ← hydrogenase 3, membrane subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,741,5251.→G 6.7% 4.3 ‑2.7 ‑0.530coding (748/924 nt)hycDhydrogenase 3, membrane subunit
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (6/22):  total (8/22)
Fisher's exact test strand distribution p-value = 6.44e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.96e-01

GATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGC.GGTG.AA.GGTTTCCATTTGTCCCCACGGAATAAAC  >  REL606/2741490‑2741559
                                    |                                    
gATGGCAATCCCCAGCGCCAGCAGCAGTCCACctgc                                       <  1:6158505/36‑1 (MQ=37)
gATGGCAATCACCAGCGCCATCAGCAGTCCACcggc                                       <  1:4720916/36‑1 (MQ=255)
gATGGCAATCACCAGCGCCAGCAGCAGTCCACcggc                                       >  1:1782527/1‑36 (MQ=255)
gATGGCAATCACCAGCGCCAGCAGCAGTCCACcggc                                       >  1:6473583/1‑36 (MQ=255)
 aTGGCAATCACCAGCGCCAGCAGCAGTCCACcggc.g                                     <  1:5549546/36‑1 (MQ=255)
 aTGGCAATCACCAGCGCCAGCAGCAGTCCACcggc.g                                     >  1:3332077/1‑36 (MQ=255)
 aTGGCAATCACCAGCGCCAGCAGCAGTACACcggc.g                                     >  1:208010/1‑36 (MQ=255)
  tGGCAATCACCAGCGCCAGCAGCAGTCCACcggc.gg                                    <  1:1099898/36‑1 (MQ=255)
   ggCAATCACCAGCGCCAGCAGCAGTCCGCCGGC.GGt                                   <  1:2719749/36‑1 (MQ=255)
     cAATCACCAGCGCCAGCAGCAGTCCACCGGC.GGTgg                                 >  1:4670904/1‑35 (MQ=255)
      aaccACCAGCGCCAGCAGCAGTCCACCGGC.GGTG.aa                               <  1:4984671/33‑1 (MQ=255)
      aaTCACCAGCGCCAGCAGCAGTCCACCGGC.GGTG.aa                               <  1:6810360/36‑1 (MQ=255)
          aCCAGCGCCAGCAGCAGTCCACCGGC.GGTG.AA.GGtt                          <  1:5283984/36‑1 (MQ=255)
          aCCAGCGCCAGCAGCAGTCCACCGGC.GGGG.GA.gggg                          >  1:3981269/1‑34 (MQ=21)
           ccAGCGCCAGCAGCAGTCCACCGGC.GGTG.AA.GGttt                         <  1:3885371/36‑1 (MQ=255)
           ccAGCGCCAGCAGCAGTCCACCGGC.GGTG.AA.GGttt                         <  1:5105565/36‑1 (MQ=255)
            cAGCGCCAGCAGCAGTCCACCGGC.GGTG.AA.GGTTTc                        <  1:3834657/36‑1 (MQ=255)
            cAGCGCCAGCAGCAGTCCACCGGC.GGGGAAA.GGttt                         >  1:265337/1‑36 (MQ=21)
             aGCGCCAGCAGCAGTCCACCGGC.GGTG.AA.GGTTTcc                       <  1:7473419/36‑1 (MQ=255)
             aGCGCCAGCAGCAGTCCACCGGC.GGTG.AA.GGTTTcc                       <  1:6328310/36‑1 (MQ=255)
               cgcCAGCAGCAGTCCACCGGC.GGTG.AA.GGTTTCCAt                     <  1:2049948/36‑1 (MQ=255)
                gcCAGCACCAGTCCTCCGGC.GGTG.AA.GGTTTCCAtt                    <  1:294222/36‑1 (MQ=37)
                 ccAGCAGCAGTCCACCGGC.GGTG.AA.GGTTTCCAttt                   <  1:4127954/36‑1 (MQ=255)
                 ccAGCAGCAGTCCACCGGC.GGTG.AA.GGTTTCCAttt                   <  1:5338504/36‑1 (MQ=255)
                  cagcagcagTCCACCGGC.GGTG.AA.GGTTTCCATTTg                  <  1:3687061/36‑1 (MQ=255)
                  cagcagcagTCCACCGGC.GGTG.AA.GGTTTCCATTTg                  <  1:6035491/36‑1 (MQ=255)
                   agcagcagTCCACCGGCGGGGG.AG.GGTTTCCTtttt                  >  1:2918116/1‑35 (MQ=12)
                    gcagcagTCCACCGGCGGGGG.AA.GGTTTCCTttttt                 >  1:4457472/1‑34 (MQ=16)
                     cagcagTCCACCGGC.GGTG.AA.GGTTTCCATTTGTcc               <  1:1138812/36‑1 (MQ=255)
                       gcagTCCACCGGC.GGTG.AA.GGTTTCCATTTGTcccc             <  1:2594359/36‑1 (MQ=255)
                       gcagCCCCCCGGC.GGTG.AA.GGTTTCCATTTGTcccc             <  1:5932109/36‑1 (MQ=37)
                        cagTCCACCGGC.GGTG.AA.GGTTTCCATTTGTCCCCa            <  1:2257278/36‑1 (MQ=255)
                         agTCCACCGGC.GGTG.AA.GGTTTCCATTTGTCCCCAc           <  1:2509464/36‑1 (MQ=255)
                            ccACCGGC.GGTG.AAGGGTTTCCTTTTTTCCCCCCgc         >  1:3245360/1‑35 (MQ=12)
                            ccACCGGC.GGTG.AA.GGTTTCCATTTGTCCCCACGGa        <  1:7419778/36‑1 (MQ=255)
                              aCCGGC.GGTG.AA.GGTTTCCATTTGTCCCCACGGAAt      <  1:3956253/36‑1 (MQ=255)
                               ccGGC.GGTG.AA.GGTTTCCATTTGTCCCCACGGAATa     >  1:3617013/1‑36 (MQ=255)
                               ccGGC.GGTG.AA.GGTTTCCATTTGTCCCCACGGAATa     >  1:7202245/1‑36 (MQ=255)
                                cggc.ggTG.AA.GGTTTCCATTTGTCCCCACGGAATaa    <  1:6275984/36‑1 (MQ=255)
                                  gc.ggTG.AA.GGTTTCCATTTGTCCCCACGGAATAAac  >  1:2941707/1‑36 (MQ=255)
                                    |                                    
GATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGC.GGTG.AA.GGTTTCCATTTGTCCCCACGGAATAAAC  >  REL606/2741490‑2741559

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