Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,732,6700C→A11.7% 6.1 ‑17.2 0.0141296hydNformate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (1/21):  ref base (79/40):  total (80/61)
Fisher's exact test strand distribution p-value = 3.29e-08
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

CCGGACGTACCACCACTTCCATCGCGCCATACGGGCAAGCCACAACGCAGGTTTTGCAACCAATGCAACGT  >  REL606/2732635‑2732705
                                   |                                   
ccGGACGTACCACCACTTCCATCGCGCCATACGGGc                                     >  1:6534986/1‑36 (MQ=255)
 cGGACGTTCCACCACTTCCATCGCGCCATACGGGCa                                    >  1:2423160/1‑36 (MQ=255)
 cGGACGTACCACCACTTCCATCGCGCCATACGGGCa                                    <  1:5752578/36‑1 (MQ=255)
  ggACGTACCACCACTTCCATCGCGCCATACGGGCaa                                   >  1:7022823/1‑36 (MQ=255)
  ggACGTACCACCACTTCCATCGCGCCATACGGGCaa                                   >  1:46695/1‑36 (MQ=255)
   gaagTACCACCACTTCCATCGCGCCATACGGGCAAg                                  >  1:5567185/4‑36 (MQ=255)
   gACGTACCACCACTTCCATCGCGCCATACGGGCAAg                                  >  1:5810179/1‑36 (MQ=255)
    aCGTACCACCACTTCCATCGCGCCATACGGGCAAGc                                 <  1:6062307/36‑1 (MQ=255)
     cGTACCACCACTTCCATCGCGCCATACGGGCAAGcc                                >  1:1042728/1‑36 (MQ=255)
     cGTACCACCACTTCCATCGCGCCATACGGGCAAGcc                                >  1:3855963/1‑36 (MQ=255)
     cGTACCACCACTTCCATCGCGCCATACGGGCAAGcc                                >  1:238288/1‑36 (MQ=255)
      gTACCACCACTTCCATCGCGCCATACGGGCAAGCca                               <  1:4782091/36‑1 (MQ=255)
       tACCACCACTTCCATCGCGCCATACGGGCAAGCcac                              <  1:8436353/36‑1 (MQ=255)
       tACCACCACTTCCATCGCGCCATACGGGCAAGCcac                              >  1:5835284/1‑36 (MQ=255)
       tACCACCACTTCCATCGCGCCATACGGGCAAGCcac                              <  1:3067381/36‑1 (MQ=255)
       tACCACCACTTCCATCGCGCCATACGGGCAAGCcac                              >  1:2882855/1‑36 (MQ=255)
       tACCACCACTTCCATCGCGCCATACGGGCAAGCcac                              >  1:1033782/1‑36 (MQ=255)
        cccaccaCTCCCATCGCGCAAAACGGGCAAGCcaca                             <  1:1484482/35‑1 (MQ=18)
        accaccACTTCCATCGCGCCATACGGGCAAGCcaca                             >  1:4845364/1‑36 (MQ=255)
         ccaccaCTTCCATCGCGCCATACGGGCAAGCCACaa                            >  1:2828817/1‑36 (MQ=255)
          caccacTTCCATCGCGCCATACGGGCAAGCCACAAc                           >  1:6194607/1‑36 (MQ=255)
           accacTTCCATCGCGCCATACGGGCAAGCCACAACg                          >  1:6036870/1‑36 (MQ=255)
           accacTTCCATCGCGCCATACGGGCAAGCCACAACg                          >  1:2216094/1‑36 (MQ=255)
              ccTTCCACCGCGCAATACGTTCAAGCCACAACGCAg                       <  1:2937732/35‑1 (MQ=16)
              acTTCCATCGCGCCATACGGGCAAGCCACAACGCAg                       >  1:5117246/1‑36 (MQ=255)
              acTTCCATCGCGCCATACGGGCAAGCCACAACGCAg                       >  1:4690963/1‑36 (MQ=255)
               cTTCCATCGCGCCATACGGGCAAGCCACAACGCAgg                      >  1:5934852/1‑36 (MQ=255)
               cTTCCATCGCGCCATACGGGCAAGCCACAACGCAgg                      <  1:3925665/36‑1 (MQ=255)
               cTTCCATCGCGCAATACGGGCAAGCCACAACGCAgg                      <  1:4899721/36‑1 (MQ=255)
                ttCCATCGCGCCATAGGGGCAAGCCACAACGCAGGt                     <  1:6726823/36‑1 (MQ=255)
                ttCCATCGCGCAATCCGGGCAAGCCACAACGCAGGt                     <  1:673005/36‑1 (MQ=37)
                ttCCACCGCGCCATACGGGCAAGCCACAACGCAGGt                     <  1:2528178/36‑1 (MQ=255)
                tcccATCGCGCAATACGGGAAAGCCACAACGCAGGt                     <  1:2446559/34‑1 (MQ=25)
                 tCCATCGCGCCATACGGGCAAGCCACAACGCAGGtt                    >  1:4602187/1‑36 (MQ=255)
                 tCCATCGCGCCATACGGGCAAGCCACAACGCAGGtt                    >  1:4180646/1‑36 (MQ=255)
                 tCCATCGCGCCATACGGGCAAGCCACAACGCAGGtt                    >  1:1058331/1‑36 (MQ=255)
                  ccATCGCGCCATACGGGCAAGCCACAACGCAGGttt                   >  1:5753149/1‑36 (MQ=255)
                  ccATCGCGCCATACGGGCAAGCCACAACGCAGGttt                   <  1:2615518/36‑1 (MQ=255)
                  ccACCGCGCAATACGGGCAAGCCACAACGCAGGttt                   <  1:6496570/36‑1 (MQ=37)
                   cATCGCGCCATACGGGCAAGCCACAACGCCGGtttt                  >  1:5502185/1‑36 (MQ=255)
                   cATCGCGCCATACGGGCAAGCCACAACGCAGGtttt                  >  1:7956824/1‑36 (MQ=255)
                   cATCGCGCCATACGGGCAAGCCACAACGCAGGtttt                  >  1:4716332/1‑36 (MQ=255)
                   cATCGCGCCATACGGGCAAGCCACAACGCAGGtttt                  >  1:3948763/1‑36 (MQ=255)
                    ctCGCGCAATAGGGGCAAGCCACAACGCAGGTTTTg                 <  1:645221/35‑1 (MQ=37)
                    agcgccCCAAACGGGAAAGCCACAACGCAGGTTTTg                 <  1:3822959/34‑1 (MQ=21)
                    aTCGCGCCATACGGGCCCGCCACCACGCACGCTTTg                 >  1:3649387/1‑36 (MQ=16)
                    aTCGCGCCATACGGGCAAGCCACAACGCAGGTTTTg                 >  1:3453707/1‑36 (MQ=255)
                    aTCGCGACATACGGGCAAGCCACAACGCAGGTTTTg                 >  1:7646886/1‑36 (MQ=255)
                     tCGCGCCATACGGGCAAGCCACAACGCAGGTTTTGc                >  1:1894010/1‑36 (MQ=255)
                     tCGCGCCATACGGGCAAGCCACAACGCAGGTTTTGc                >  1:8041815/1‑36 (MQ=255)
                     tCGCGCCATACGGGCAAGCCACAACGCAGGTTTTGc                >  1:5233391/1‑36 (MQ=255)
                     tCGCGCCATAAGGGAAAGCCACAACGCAGGTTTTGc                >  1:5417998/1‑36 (MQ=37)
                     ccgcaCCATCCGGGCAAGCCACAACGCAGGTTTTGc                <  1:538536/35‑1 (MQ=37)
                       gcgcCATACGGGCAAGCCACAACGCAGGTTTTGCaa              >  1:6164620/1‑36 (MQ=255)
                       gcgcCATACGGGCAAGCCACAACGCAGGTTTTGCaa              >  1:1093165/1‑36 (MQ=255)
                       gcgcCATACGGGCAAGCCACAACGCAGGTTTTGCaa              >  1:2729688/1‑36 (MQ=255)
                       gcgcCATACGGGCAAGCCACAACGCAGGCTTTGCaa              >  1:8110286/1‑36 (MQ=255)
                       gcgcAATACGGGAAAGCCACAACGCAGGTTTTGCaa              <  1:1467753/36‑1 (MQ=37)
                        cgcCATACGGGCAAGCCTCAACGCAGGTTTTGCAAc             <  1:2134108/36‑1 (MQ=255)
                        cgcCATACGGGCAAGCCACAACGCAGGTTTTGCAAc             >  1:1129673/1‑36 (MQ=255)
                        cgcCACGCGGGAAAACCACAACGCAGGTTTTGCAAc             <  1:541712/36‑1 (MQ=16)
                         gcaaTACGGGCAAGCCACAACGCAGGTTTTGCAAcc            <  1:7330819/33‑1 (MQ=255)
                         gcCATACGGGCAAGCCACAACGCAGGTTTTGCAAcc            >  1:6317285/1‑36 (MQ=255)
                         gcCATACGGGCAAGCCACAACGCAGGTTTTGCAAcc            >  1:8018088/1‑36 (MQ=255)
                         gcCATACGGGCAAGCCACAACGCAGGTTTTGCAAcc            >  1:3748152/1‑36 (MQ=255)
                         gcCATACGGAAAAGCCACAACGCAGGTTTTGCAAcc            <  1:3574326/36‑1 (MQ=37)
                          ccATACGGGCAAGCCACAACGCAGGTTTTGCAACCa           >  1:5320589/1‑36 (MQ=255)
                          ccATACGGGCAAGCCACAACGCAGGTTTTGCAACCa           >  1:1554157/1‑36 (MQ=255)
                          ccATACGGGCAAGCCACAACGCAGGTTTTGCAACCa           >  1:1823725/1‑36 (MQ=255)
                          ccATACGGGAAAGCCACAACGCAGGTTTTGCAACCa           <  1:4993816/36‑1 (MQ=255)
                          caaTACGGGAAAGCCACAACGCAGGTTTTGCAACCa           <  1:6959251/34‑1 (MQ=37)
                          caaTACGGGAAAGCCACAACGCAGGTTTTGCAACCa           <  1:6655142/34‑1 (MQ=37)
                           cATACGGGCAAGCCACAACGCAGGTTTTGCAACCaa          >  1:7214556/1‑36 (MQ=255)
                           cATACGGGCAAGCCACAACGCAGGTTTTGCAACCaa          >  1:5239787/1‑36 (MQ=255)
                           cATACGGGCAAGCCACAACGCAGGTTTTGCAACCaa          >  1:2245643/1‑36 (MQ=255)
                           cATACGGGAAAGCCACAACGCAGGTTTTGCAACCaa          <  1:212673/36‑1 (MQ=255)
                           aaTACGGGCAAGCCACAACGCAGGTTTTGCAACCaa          <  1:6692763/35‑1 (MQ=255)
                            aaaCGGACAAGCCACAACGCAGGTTTTGCAACCAAt         <  1:4793806/34‑1 (MQ=37)
                            aTACGGGCAAGCCACAACGCAGGTTTTGCAACCAAt         >  1:5348250/1‑36 (MQ=255)
                            aTACGGGCAAGCCACAACGCAGGTTTTGCAACCAAt         >  1:3299505/1‑36 (MQ=255)
                            aTACGGGAAAGCCACAACGCAGGTTTTGCAACCAAt         <  1:3550132/36‑1 (MQ=255)
                            aTACGAGAAAGCCACAACGCAGGTTTTGCAACCAAt         <  1:1243082/36‑1 (MQ=37)
                             ttcGGGAAAGCCACAACGCAGGTTTTGCAACCAATg        <  1:323408/34‑1 (MQ=37)
                             tACGGGCAAGCCACAACGCAGGTTTTGCAACCAATg        >  1:5392392/1‑36 (MQ=255)
                             tACGGGCAAGCCACAACGCAGGTTTTGCAACCAATg        >  1:4517096/1‑36 (MQ=255)
                             aaCGGGCAAGCCACAACGCAGGTTTTGCAACCAATg        <  1:7469643/35‑1 (MQ=255)
                              aCGGGCAAGCCACAACGCAGGTTTTGCAACCAATGc       <  1:2300081/36‑1 (MQ=255)
                              aCGGGCAAGCCACAACGCAGGTTTTGCAACCAATGc       <  1:5488663/36‑1 (MQ=255)
                              aCGGGCAAGCCACAACGAAGGTTTTGCAACCAATGc       >  1:2018252/1‑36 (MQ=255)
                              aCGGGAAAGCCACAACGCAGGTTTTGCAACCAATGc       <  1:6341777/36‑1 (MQ=255)
                               cGGGCAAGCCACAACGCAGGTTTTGCAACCAATGCa      >  1:2479626/1‑36 (MQ=255)
                               cGGGAAAGCCACAACGCAGGTTTTGCAACCAATGCa      <  1:2846359/36‑1 (MQ=255)
                                gggCAAGCCACAACGCAGGTTTTGCAACCAATGCaa     >  1:1092755/1‑36 (MQ=255)
                                gggCAAGCCACAACGCAGGTTTTGCAACCAATGCaa     >  1:2380845/1‑36 (MQ=255)
                                gggCAAGCCACAACGCAGGTTTTGCAACCAATGCaa     >  1:5001108/1‑36 (MQ=255)
                                gggCAAGCCACAACGAAGGTTTTGCAACCAATGCaa     >  1:3247866/1‑36 (MQ=255)
                                 ggCAAGCCACAACGCAGGTTTTGCAACCAATGCAAc    >  1:5548321/1‑36 (MQ=255)
                                 ggCAAGCCACAACGCAGGTTTTGCAACCAATGCAAc    >  1:1865953/1‑36 (MQ=255)
                                 ggCAAGCCACAACGCAGGTTTTGCAACCAATGCAAc    <  1:3223223/36‑1 (MQ=255)
                                   cAAGCCACAACGCAGGTTTTGCAACCAATGCAACGt  >  1:5497886/1‑36 (MQ=255)
                                   |                                   
CCGGACGTACCACCACTTCCATCGCGCCATACGGGCAAGCCACAACGCAGGTTTTGCAACCAATGCAACGT  >  REL606/2732635‑2732705

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG