Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,739,309 T→C 15.8% K102E (AAG→GAG)  hycF ← formate hydrogenlyase complex iron‑sulfur protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,739,3090T→C15.8% 4.9 ‑3.8 ‑0.037304hycFformate hydrogenlyase complex iron‑sulfur protein
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (17/14):  total (17/20)
Fisher's exact test strand distribution p-value = 2.20e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.10e-01

GCGAAGCGGGACTGTTGCAGGAAGTCTTCTTTCTTCCACACCGCCAGTTCGTACTCTTGCGACAGTTTG  >  REL606/2739275‑2739343
                                  |                                  
gcgAAGCGGGACTGTTGCAGGAAGTCTTCTttcttc                                   >  1:4007773/1‑36 (MQ=255)
gcgAAGCGGGACTGTTGCAGGAAGTCTTCTttcttc                                   >  1:2858453/1‑36 (MQ=255)
gcgAAGCGGGACTGTTGCAGGAAGTCTTCTttcttc                                   <  1:7738190/36‑1 (MQ=255)
gcgAAGCGGGACTGTTGCAGGAAGTCTTCTttcttc                                   >  1:7755921/1‑36 (MQ=255)
gcgAAGCGGGACTGTTGCAGGAAGTCTTCTttcttc                                   <  1:1853309/36‑1 (MQ=255)
 cgAAGCGGGACTGTTGCAGGAAGTCTTCTTTCTTcc                                  <  1:977358/36‑1 (MQ=255)
  gAAGCGGGACTGTTGCAGGAAGTCTTCTTTCTTCca                                 <  1:4753621/36‑1 (MQ=255)
  gAAGCGGGACTGTTGCAGGAAGTCTTCTTTCTTCca                                 >  1:7834550/1‑36 (MQ=255)
   aaGCGGGACTGTTGCAGGAAGTCTTCTTTCTTCcac                                <  1:3539658/36‑1 (MQ=255)
   aaGCGGGACTGTTGCAGGAAGTCTTCTTTCTTCcac                                >  1:1581979/1‑36 (MQ=255)
   aaGCGGGACTGTTGCAGGAAGTCTTCTTTCTTCcac                                >  1:1850234/1‑36 (MQ=255)
    aGCGGGACTGTTGCAGGAAGTCTTCTTTCTTCcaca                               <  1:2260062/36‑1 (MQ=255)
     gCGGGACTGTTGCAGGAAGTCTTCTTTCTTCCacac                              <  1:157936/36‑1 (MQ=255)
      cGGGACTGTTGCAGGACGTCGTCTTTCTTCCACAcc                             >  1:2132504/1‑36 (MQ=37)
      cGGGACTGTTGCAGGAAGTCTTCTTTCTTCCACAcc                             <  1:2233962/36‑1 (MQ=255)
       gggACTGTTGCAGTAAGTCTTCTTTCTTCCACACCg                            >  1:3344873/1‑36 (MQ=255)
       gggACTGTTGCAGGAAGTCTTCTTTCTTCCACACCg                            <  1:4418870/36‑1 (MQ=255)
        ggACTGTTGCAGGAAGTCTTCTTTCTTCCACACCGc                           >  1:5061896/1‑36 (MQ=255)
        ggACTGTTGCAGGAAGTCTTCTTTCTTCCACACCGc                           <  1:3423650/36‑1 (MQ=255)
        ggACTGTTGCAGGAAGTCTTCTTTCTTCCACACCGc                           <  1:3671942/36‑1 (MQ=255)
         gACTGTTGCAGGAAGTCTTCTTTCTTCCACACCGcc                          >  1:1388353/1‑36 (MQ=255)
         gACTGTTGCAGGAAGTCTTCTTTCTTCCACACCGcc                          <  1:3385539/36‑1 (MQ=255)
           cTGTTGCAGGAAGTCTTCTTTCTTCCACACCGCCAg                        >  1:415746/1‑36 (MQ=255)
             gTTGCAGGAAGTCTTCTTTCTTCCACACCGCCAGtt                      >  1:4813385/1‑36 (MQ=255)
             gTTGCAGGAAGTCTTCTTTCTTCCACACCGCCAGtt                      >  1:3258839/1‑36 (MQ=255)
              ttGCAGGAAGTCTTCTTTCTTCCACACCGCCAGTTc                     >  1:7456112/1‑36 (MQ=255)
              ttGCAGGAAGTCTTCTTTCTTCCACACCGCCAGTTc                     >  1:1549579/1‑36 (MQ=255)
              ttGCAGAAAGTCTTCTTTCTTCCACACCGCCAGTTc                     <  1:1066859/36‑1 (MQ=255)
                gcgggAAGTCTTCTTTCTCCCACACCGCCAGTTCGt                   <  1:6363017/33‑1 (MQ=25)
                gCAGGAAGTCTTCTTTCTTCCACACCGCCAGTTCGt                   <  1:7335279/36‑1 (MQ=255)
                 cAGGAAGTCTTCTTTCTTCCACACCGCAAGTTCGTa                  <  1:777773/36‑1 (MQ=255)
                  gggAAGTCTTCTTTCTCCCCCACCGCCAGTTCGTAc                 <  1:3020199/35‑1 (MQ=37)
                    gAAGTCTTCTTTCTTCCACACCGCCAGTTCGTActc               >  1:2487071/1‑36 (MQ=255)
                    gAAGTCTTCTTTCTTCCACACCGCCAGTTCGTActc               >  1:4646032/1‑36 (MQ=255)
                    gAAGTCTTCTTTCTTCAACACCGCCAGTTCGTActc               >  1:2939086/1‑36 (MQ=255)
                     aaGTCTTCTTTCTCCCACACCGCCAGTTCGTActct              <  1:3788233/36‑1 (MQ=255)
                       gTCTTCTTTCTTCCACACCGCCAGTTCGTACTCTTg            >  1:605200/1‑36 (MQ=255)
                        ttttctttctCCCACACCGCCAGTTCGTACTCTTGc           <  1:6725467/34‑1 (MQ=37)
                        tcttctTTCTCCCACACCGCCAGTTCGTACTCTTGc           <  1:5139629/36‑1 (MQ=255)
                        tcttctTCCTCCCACACCGCCAGTTCGTACTCTTGc           <  1:3973919/36‑1 (MQ=37)
                         cttcttTCTTCCACACCGCCAGTTCGTACTCTTGCg          >  1:5687207/1‑36 (MQ=255)
                          ttctttctTCCACGCCGCCAGTTCGTACTCTTGCGa         <  1:277675/36‑1 (MQ=255)
                          tcctttcttCCACACCGCCAGTTCGTACTCTTGCGa         >  1:3832445/3‑36 (MQ=255)
                           tctttcttCCACACCCCCAGTTCGTACTCTTGCGAc        >  1:938034/1‑36 (MQ=255)
                            ctttcccCCACACCGCCAGTTCGTACTCTTGCGACa       <  1:4413839/36‑1 (MQ=37)
                             tttcttCCACACCGCCAGTTCGTACTCTTGCGACAg      <  1:1501607/36‑1 (MQ=255)
                              ttcttcCACACCGCCAGTTCGTACTCTTGCGACAGt     >  1:871804/1‑36 (MQ=255)
                              ttctccCACACCGCCAGTTCGTACTCTTGCGACAGt     <  1:809353/36‑1 (MQ=255)
                               tcttcCACACCGCCAGTTCGTACTCTTGCGACAGtt    >  1:3282423/1‑36 (MQ=255)
                               tcttcCACACCGCCAGTTCGTACTCTTGCGACAGtt    >  1:2470207/1‑36 (MQ=255)
                                 ttcCACACCGCCAGTTCGTACTCTTGCGACAGTTTg  <  1:7636769/36‑1 (MQ=255)
                                  |                                  
GCGAAGCGGGACTGTTGCAGGAAGTCTTCTTTCTTCCACACCGCCAGTTCGTACTCTTGCGACAGTTTG  >  REL606/2739275‑2739343

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG