Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,730,1880G→A10.2% 3.1 ‑10.6 0.0972098hypFcarbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand):  new base (1/11):  ref base (62/23):  total (63/34)
Fisher's exact test strand distribution p-value = 2.61e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCCCGCCGCTAAATACCAGCGTAGTGATAC  >  REL606/2730153‑2730223
                                   |                                   
aTAATGCGCCAGACGTGCACGCAGCAAACGGTTATg                                     >  1:2694348/1‑36 (MQ=255)
 tAATGCGCCAGACGTGCACGCAGCAAACGGTTATGa                                    <  1:7205247/36‑1 (MQ=255)
 tAATGCGCCAGACGTGCACGCAGCAAACGGTTATGa                                    >  1:3243117/1‑36 (MQ=255)
 tAATGCGCCAGACGTGCACGCAGCAAACGGTTATGa                                    >  1:645827/1‑36 (MQ=255)
  aaTGCGCCAGACGTGCACGCAGCAAACGGTTATGaa                                   <  1:7624395/36‑1 (MQ=255)
  aaTGCGCCAGACGTGCACGCAGCAAACGGTTATGaa                                   <  1:5232107/36‑1 (MQ=255)
  aaTGCGCCAGACGTGCACGAAGCAAACGGTTATGaa                                   <  1:5565360/36‑1 (MQ=255)
   aTGCGCCAGACGTGCACGGAGCAAACGGTTATGAAt                                  >  1:5440429/1‑36 (MQ=255)
   aTGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt                                  <  1:3083495/36‑1 (MQ=255)
   aTGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt                                  >  1:7986335/1‑36 (MQ=255)
   aTGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt                                  >  1:6046320/1‑36 (MQ=255)
   aTGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt                                  >  1:6148075/1‑36 (MQ=255)
   aTGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt                                  >  1:2356353/1‑36 (MQ=255)
   aTGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt                                  >  1:2264760/1‑36 (MQ=255)
   aTGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt                                  <  1:4103923/36‑1 (MQ=255)
   aTGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt                                  <  1:4406359/36‑1 (MQ=255)
   aTGCGCCAGACGTGCACGCAGCAAACGGGTATGAAt                                  >  1:6396434/1‑36 (MQ=255)
    tGCGCCAGACGTGCACGCAGCAAACGGTTATGAATa                                 >  1:3168689/1‑36 (MQ=255)
     gcgcCAGACGTGCACGCAGCAAACGGTTATGAATaa                                <  1:5454741/36‑1 (MQ=255)
     gcgcCAGACGTGCACGCAGCAAACGGTTATGAATaa                                >  1:33621/1‑36 (MQ=255)
       gcCAGACGTGCACGCAGCAAACGGTTATGAATAAcc                              >  1:5825479/1‑36 (MQ=255)
        ccAGACGTGCACGCAGCAAACGGTTATGAATAAccc                             <  1:3195962/36‑1 (MQ=255)
        ccAGACGTGCACGCAGCAAACGGTTATGAATAAccc                             <  1:5810501/36‑1 (MQ=255)
        ccAGACGTGCACGCAGCAAACGGTTATGAATAAccc                             >  1:3446615/1‑36 (MQ=255)
        ccAGACGTGCACGCAGCAAACGGTTATGAATAAccc                             >  1:5836494/1‑36 (MQ=255)
         cAGACGTGCACGCAGCAAACGGTTATGAATAAcccc                            >  1:3091385/1‑36 (MQ=255)
         cAGACGTGCACGCAGCAAACGGTTATGAATAAcccc                            <  1:2396399/36‑1 (MQ=255)
         cAGACGTGCACGCAGCAAACGGTTATGAATAAcccc                            <  1:2368406/36‑1 (MQ=255)
         cAGACGTGCACGCAGCAAACGGTTATGAATAAcccc                            <  1:6933777/36‑1 (MQ=255)
         cAGACGTGCACGCAGCAAACGGTTATGAATAAcccc                            <  1:1203135/36‑1 (MQ=255)
         cAGACGTGCACGCAGCAAACGGTTATGAATAAcccc                            >  1:5168893/1‑36 (MQ=255)
         cAGACGTGCACGAAGCAAACGGTTATGAATAAcccc                            <  1:5833441/36‑1 (MQ=255)
          aGACGTGCACGCAGCAAACGGTTATGAATAACCccg                           <  1:7129924/36‑1 (MQ=255)
          aGACGTGCACGCAGCAAACGGTTATGAATAACCccg                           >  1:6986060/1‑36 (MQ=255)
          aGACGTGCACGCAGCAAACGGTTATGAATAACCccg                           <  1:6003955/36‑1 (MQ=255)
          aGACGTGCACGCAGCAAACGGTTATGAATAACCccg                           >  1:8080330/1‑36 (MQ=255)
          aGACGTGCACGCAGCAAACGGTTATGAATAACCccg                           <  1:2245692/36‑1 (MQ=255)
          aGACGTGCACGCAGCAAACGGTTATGAATAACCccg                           <  1:6620244/36‑1 (MQ=255)
           gACGTGCACGCAGCAAACGGTTATGAATAACCccgc                          >  1:2256144/1‑36 (MQ=255)
           gACGTGCACGCAGCAAACGGTTATGAATAACCccgc                          <  1:2132675/36‑1 (MQ=255)
            aCGTGCACGCAGCAAACGGTTATGAATAACCccgcc                         <  1:7280489/36‑1 (MQ=255)
            aCGTGCACGCAGCAAACGGGTATGAATAACCccgcc                         >  1:6391047/1‑36 (MQ=255)
            aCGTGCACGAAGCAAACGGTTATGAATAACCccgcc                         <  1:7602496/36‑1 (MQ=255)
             cGTGCACGCAGCAAACGGTTATGAATAACCccgccg                        >  1:2761372/1‑36 (MQ=255)
             cGTGCACGCAGAAAACGGTTATGAATAACCccgccg                        <  1:1401010/36‑1 (MQ=255)
              gTGCACGCAGCAAACGGTTATGAATAACCCcgccgc                       >  1:8106939/1‑36 (MQ=255)
               tGCACGCAGCAAACGGTTATGAATAACCCCGCCGCt                      >  1:4805776/1‑36 (MQ=255)
                gCACGCAGCAAACGGTTATGAATAACCCCGCCGCTa                     >  1:7262257/1‑36 (MQ=255)
                gCACGCAGCAAACGGTTATAAATAACCCCGCCGCTa                     >  1:8184735/1‑36 (MQ=255)
                 cACGCAGCAAACGGTTATGAATAACCCCGCCGCTaa                    >  1:6544376/1‑36 (MQ=255)
                 cACGCAGCAAACGGTTATGAATAACCCCGCCGCTaa                    >  1:7300164/1‑36 (MQ=255)
                  aCGCAGCAAACGGTTATGAATAACCCCGCCGCTaaa                   >  1:7161363/1‑36 (MQ=255)
                    gcagcaAACGGTTATGAATAACCCCGCCGCTaaaca                 >  1:3877326/1‑34 (MQ=255)
                    gcagcaAACGGTTATGAATAACCCCGCCGCTAAATa                 >  1:8363715/1‑36 (MQ=255)
                     cagcaAACGGTTATGAATAACCCCGCCGCTAAATAc                >  1:1813910/1‑36 (MQ=255)
                     cagcaAACGGTTATGAATAACCCCGCCGCTAAATAc                >  1:2662436/1‑36 (MQ=255)
                     cagcaAACGGTTATGAATAACCCCGCCGCTAAATAc                >  1:4114079/1‑36 (MQ=255)
                      agcaAACGGTTATGAATAACCCCGCCGCTAAATAcc               >  1:8153725/1‑36 (MQ=255)
                      agcaAACGGTTATGAATAACCCCGCCGCTAAATAcc               >  1:2870555/1‑36 (MQ=255)
                      agcaAACGGTTATGAATAACCCCGCCGCTAAATAcc               >  1:3900174/1‑36 (MQ=255)
                      agcaAACGGTTATGAATAACCCCGCCGCTAAATAcc               >  1:4067528/1‑36 (MQ=255)
                       gcaAACGGTTATGAATAACCCCGCCGCTAAATACCa              >  1:4485348/1‑36 (MQ=255)
                       gcaAACGGTTATGAATAACCCCGCCGCTAAATACCa              >  1:2406566/1‑36 (MQ=255)
                       gcaAACGGTTATGAATAACCCCGCCGCTAAATACCa              >  1:6155646/1‑36 (MQ=255)
                       gcaAACGGTTATGAATAACCCCGCCGCTAAATACCa              >  1:5430528/1‑36 (MQ=255)
                       gaaaaCGGTTATAAAAAACCCCGCCGCTAAATACCa              <  1:1720970/34‑1 (MQ=25)
                       aaaaaGGGTTATAAAAAACCCCGCCGCTAAATACCa              <  1:1878477/34‑1 (MQ=21)
                        caAACGGTTATGAATAACCCCGCCGCTAAATACCAg             >  1:2931677/1‑36 (MQ=255)
                        caAACGGTTATGAATAACCCCGCCGCTAAATACCAg             >  1:5909481/1‑36 (MQ=255)
                        caAACGGTTATGAATAACCCCGCCGCTAAATACCAg             >  1:5107312/1‑36 (MQ=255)
                        aaaaCGTTTATAAATAACCCCGCCGCTAAATACCAg             <  1:2150293/35‑1 (MQ=37)
                        aaaaCGTTTATAAAAAACCCCGCCGCTAAATACCAg             <  1:5840656/35‑1 (MQ=21)
                         aaaCGGTTATGAATAACCCCGCCGCTAAATACCAGc            >  1:2781100/1‑36 (MQ=255)
                         aaaCGGTTATGAATAACCCCGCCGCTAAATACCAGc            >  1:4831770/1‑36 (MQ=255)
                         aaaCGGTTATAAATAACCCCGCCGCTAAATACCAGc            <  1:1190200/36‑1 (MQ=255)
                          aaCGGTTATGAATAACCCCGCCGCTAAATACCAGCg           >  1:2517254/1‑36 (MQ=255)
                          aaCGGTTATGAATAACCCCGCCGCTAAATACCAGCg           >  1:7903269/1‑36 (MQ=255)
                           aCGGTTATGAATAACCCCGCCGCTAAATACCAGCGt          >  1:1743421/1‑36 (MQ=255)
                           aCGGTTATGAATAACCCCGCCGCTAAATACCAGCGt          >  1:7815693/1‑36 (MQ=255)
                            cgtttATAAAAAACCCCGCCGCTAAATACCAGCGTa         <  1:7807256/33‑1 (MQ=21)
                            cGGTTATAAAAAACCCCGCCGCTAAATACCAGCGTa         <  1:2236492/36‑1 (MQ=37)
                             ggTTATGAATAACCCCGCCGCTAAATACCAGCGTAg        >  1:8327851/1‑36 (MQ=255)
                             ggTTATGAATAACCCCGCCGCTAAATACCAGCGTAg        >  1:6508528/1‑36 (MQ=255)
                             ggTTATGAATAACCCCGCCGCTAAATACCAGCGTAg        >  1:7860099/1‑36 (MQ=255)
                             ggTTATGAATAACCCCGCCGCTAAATACCAGCATAg        >  1:128444/1‑36 (MQ=255)
                             ggTTAAAAATAACCCCGCCGCTAAATACCAGCGTAg        <  1:3955380/36‑1 (MQ=37)
                              gTTATGAATAACCCCGCCGCTAAATACCAGCGTAGt       >  1:5888046/1‑36 (MQ=255)
                               ttATAAATAACCCCGCCGCTAAATACCAGCGTAGTg      <  1:4985828/36‑1 (MQ=255)
                               ttATAAAAAACCCCGCCGCTAAATACCAGCGTAGTg      <  1:2752921/36‑1 (MQ=37)
                                tATGAATAACCCCGCCGCTAAATACCAGCGTAGTGa     >  1:7896713/1‑36 (MQ=255)
                                tATGAATAACCCCGCCGCTAAATACCAGCGTAGTGa     >  1:5480410/1‑36 (MQ=255)
                                tATGAATAACCCCGCCGCTAAATACCAGCGTAGTGa     >  1:6509416/1‑36 (MQ=255)
                                tATGAATAACCCCGCCGCTAAATACCAGCGTAGTGa     >  1:6531673/1‑36 (MQ=255)
                                tATAAATAACCCCGCCGCTAAATACCAGCGTAGTGa     <  1:5239886/36‑1 (MQ=255)
                                 aTGAATAACCCCGCCGCTAAATACCAGCGTAGTGAt    >  1:7500001/1‑36 (MQ=255)
                                 aTGAATAACCCCGCCGCTAAATACCAGCGTAGTGAt    >  1:5366843/1‑36 (MQ=255)
                                 aTGAATAACCCCGCCGCTAAATACCAGCGTAGTGAt    >  1:6472944/1‑36 (MQ=255)
                                  tGAATAACCCCGCCGCTAAATACCAGCGTAGTGATa   >  1:5666058/1‑36 (MQ=255)
                                   gAATAACCCCGCCGCTAAATACCAGCGTAGTGATAc  <  1:3420138/36‑1 (MQ=255)
                                   gAATAACCCCGCCGCTAAATACCAGCGTAGTGATAc  >  1:3881958/1‑36 (MQ=255)
                                   |                                   
ATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCCCGCCGCTAAATACCAGCGTAGTGATAC  >  REL606/2730153‑2730223

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG