Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,581,9660C→A15.1% 7.3 ‑18.5 ‑0.0981075tsrmethyl‑accepting chemotaxis protein I, serine sensor receptor
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (0/17):  ref base (58/21):  total (58/40)
Fisher's exact test strand distribution p-value = 9.24e-09
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.99e-01

Only 100 of 105 total aligned reads displayed.

TGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGC  >  REL606/4581931‑4582001
                                   |                                   
tggtgCAGACTATGCGCGATATCTCCACCAGTTCGc                                     <  1:3547626/36‑1 (MQ=255)
tggtgCAGACTATGCGCGATATCTCCACCAGTTCGc                                     <  1:2304724/36‑1 (MQ=255)
tggtgCAGACTATGCGCGATATCTCCACCAGTTCGc                                     >  1:4754933/1‑36 (MQ=255)
tggtgCAGACTATGCGCGATATCTCCACCAGTTCGc                                     <  1:5360582/36‑1 (MQ=255)
 ggtgCAGCCTATGCGCGATATCTCCACCAGTTCGCa                                    <  1:6567599/36‑1 (MQ=255)
 ggtgCAGACTATGCGCGATATCTCCACCAGTTCGCa                                    >  1:2869533/1‑36 (MQ=255)
 ggtgCAGACTATGCGCGATATCTCCACCAGTTCGCa                                    <  1:753278/36‑1 (MQ=255)
 ggtgCAGACTATGCGCGATATCTCCACCAGTTCGCa                                    >  1:1854586/1‑36 (MQ=255)
  gtgCAGACTATGCGCGATATCTCCACCAGTTCGCAg                                   <  1:6669507/36‑1 (MQ=255)
    gCAGACTATGCGCGATATCTCCACCAGTTCGCAGaa                                 <  1:2676837/36‑1 (MQ=255)
    gCAGACTATGCGCGATATCTCCACCAGTTCGCAGaa                                 <  1:6416714/36‑1 (MQ=255)
     cAGACTATGCGCGATATCTCCACCAGTTCGCAGaaa                                >  1:7348775/1‑36 (MQ=255)
     cAGACTATGCGCGATATCTCCACCAGTTCGCAGaaa                                >  1:3770941/1‑36 (MQ=255)
      aGACTATGCGCGATATCTCCACCAGTTCGCCGaaaa                               >  1:4448389/1‑36 (MQ=255)
      aGACTATGCGCGATATCTCCACCAGTTCGCAGaaaa                               >  1:2584826/1‑36 (MQ=255)
      aGACTATGCGCGATATCTCCACCAGTTCGCAGaaaa                               >  1:6773215/1‑36 (MQ=255)
      aGACTATGCGCGATATCTCCACCAGTTCGCAGaaaa                               >  1:2567921/1‑36 (MQ=255)
       gACTATGCGCGATATCTCCACCAGTTCGCAGAAAAt                              <  1:1114440/36‑1 (MQ=255)
       gACTATGCGCGATATCTCCACCAGTTCGCAGAAAAt                              >  1:2558409/1‑36 (MQ=255)
       gACTATGCGCGATATCTCCACCAGTTCGCAGAAAAt                              >  1:1575420/1‑36 (MQ=255)
        aCTATGCGCGATATCTCCCCCAGTTCTCAGGAAAtt                             >  1:208604/1‑35 (MQ=21)
        aCTATGCGCGATATCTCCACCAGTTCGCAGAAAATc                             >  1:1400103/1‑36 (MQ=255)
        aCTATGCGCGATATCTCCACCAGTTCGCAGAAAATc                             >  1:3955203/1‑36 (MQ=255)
        aCTATGCGCGATATCTCCACCAGTTCGCAGAAAATc                             >  1:696085/1‑36 (MQ=255)
        aCTATGCGCGATATCTCCACCAGTTCGCAGAAAATc                             <  1:4944100/36‑1 (MQ=255)
        aCTATGCGCGATATCTCCACCAGTTCGCAGAAAATc                             <  1:1926430/36‑1 (MQ=255)
         cTATGCGCGATATCTCCACCAGTTCGCAGAAAATCg                            <  1:4085916/36‑1 (MQ=255)
         cTATGCGCGATATCTCCACCAGTTCGCAGAAAATCg                            <  1:2257930/36‑1 (MQ=255)
         cTATGCGCGATATCTCCACCAGTTCGCAGAAAATCg                            >  1:179939/1‑36 (MQ=255)
          tATGCGCGATATCTCCACCAGTTCGCAGAAAATCGc                           >  1:809853/1‑36 (MQ=255)
          tATGCGCGATATCTCCACCAGTTCGCAGAAAATCGc                           <  1:2863594/36‑1 (MQ=255)
          tATGCGCGATATCTCCACCAGTTCGCAGAAAATCGc                           >  1:2724328/1‑36 (MQ=255)
          tATGCGCGATATCTCCACCAGTTCGCAGAAAATCGc                           <  1:2048751/36‑1 (MQ=255)
          tATGCGCGATATCTCCACCAGTTCGCAGAAAATCGc                           >  1:1667905/1‑36 (MQ=255)
           aTGCGCGATATCTCCACCAGTTCGCAGAAAATCGcc                          >  1:2850915/1‑36 (MQ=255)
           aTGCGCGATATCTCCACCAGTTCGCAGAAAATCGcc                          >  1:2135152/1‑36 (MQ=255)
           aTGCGCGATATCTCCACCAGTTCGCAGAAAATCGcc                          <  1:5591742/36‑1 (MQ=255)
           aTGCGCGATATCTCCACCAGTTCGCAGAAAATCGcc                          >  1:2920754/1‑36 (MQ=255)
           aTGCGCGATATCTCCACCAGTTCGCAGAAAATCGcc                          >  1:7136672/1‑36 (MQ=255)
            tGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCg                         >  1:4978383/1‑36 (MQ=255)
            tGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCg                         <  1:6344815/36‑1 (MQ=255)
            tGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCg                         >  1:6971286/1‑36 (MQ=255)
              cgcgATATCTCCACCAGTTCGCAGAAAATCGCCGat                       >  1:615700/1‑36 (MQ=255)
              cgcgATATCTCCACCAGTTCGCAGAAAATCGCCGat                       >  1:7210631/1‑36 (MQ=255)
              cgcgATACCTCCACAGGTTCAAAAAAAATCGCCGat                       <  1:630635/36‑1 (MQ=9)
               gcgATATCTCCACCAGTTCGCAGAAAATCGCCGata                      >  1:5834726/1‑36 (MQ=255)
               gcgATATCTCCACCAGTTCGCAGAAAATCGCCGata                      >  1:357860/1‑36 (MQ=255)
               gcgATATCTCCACCAGTTCGCAGAAAATCGCCGata                      >  1:6852089/1‑36 (MQ=255)
                cgATATCTCCACCAGTTCGCAGAAAATCGCCGatat                     >  1:4502035/1‑36 (MQ=255)
                cgATATCTCCACCAGTTCGCAGAAAATCGCCGatat                     >  1:6669122/1‑36 (MQ=255)
                cgATATCTCCACCAGTTCGCAGAAAATCGCCGatat                     >  1:1159858/1‑36 (MQ=255)
                cgATATCTCCACCAGTTCGAAGAAAATCGCCGatat                     <  1:1325440/36‑1 (MQ=255)
                 gATATCTCCACCAGTTCGCAGAAAATCGCCGAtatt                    >  1:1413513/1‑36 (MQ=255)
                  atatCTCCAGCAGTTCGCCGAAAATCGCCGAtatta                   <  1:1436027/36‑1 (MQ=37)
                  atatCTCCACCAGTTCGCAGAAAATCGCCGAtatta                   >  1:3138305/1‑36 (MQ=255)
                  atatCTCCACCAGTTCGCAGAAAATCGCCGAtatta                   >  1:5500435/1‑36 (MQ=255)
                  atatCTCCACCAGTTCGAAAAAAATCGCCGAtatta                   <  1:2169216/36‑1 (MQ=25)
                   tatCTCCACCAGTTCGCAGAAAATCGCCGAtattat                  >  1:6644872/1‑36 (MQ=255)
                   tatCTCCAACAGTTCGTAGAAACTCGCCGAtattat                  <  1:2565870/36‑1 (MQ=37)
                   tatCTCAACCAGTTCGAATAAAATCGCCGAtattat                  <  1:3728065/36‑1 (MQ=25)
                    atCTCCACCAGTTCGCAGAAAATCGCCGATATTATc                 >  1:6291820/1‑36 (MQ=255)
                    atCTCCACCAGTTCGCAGAAAATCGCCGATATTATc                 >  1:157744/1‑36 (MQ=255)
                    atCTCCACCAGTTCGCAGAAAATCGCCGATATTATc                 >  1:5646476/1‑36 (MQ=255)
                    atCTCCACCAGTTCGCAGAAAATCGCCGATATTATc                 >  1:6370931/1‑36 (MQ=255)
                    atCTCCACCAGTTCGAAGAAAATCGCCGATATTATc                 <  1:1787547/36‑1 (MQ=255)
                     cccccaccaGTTCGCAAAAAATCGCCGATATTATCa                <  1:2323871/33‑1 (MQ=18)
                      ctcCACCAGTTCGCAGAAAATCGCCGATATTATCAg               >  1:6743332/1‑36 (MQ=255)
                      ctcCACAAGTTCGAAAAAAGTCGCCGATATTATCAg               <  1:2040142/36‑1 (MQ=19)
                       tcCACTATTGCGGAGAAAATCGCCGATATTATCAGc              <  1:3548870/36‑1 (MQ=16)
                       tcCACCAGTTCGAAGAAAATCGCCGATATTATCAGc              <  1:4957547/36‑1 (MQ=255)
                       cccaccaGTTCGAAAAAAATCGCCGATATTATCAGc              <  1:2582140/35‑1 (MQ=18)
                        ccaccaGTTCGCAGAAAATCGCCGATATTATCAGCg             >  1:2903772/1‑36 (MQ=255)
                        ccacaaGTTCGCAGAAAATCGCCGATATTATCAGCg             <  1:5869065/36‑1 (MQ=255)
                         caccaTTTCGAAGAAAATCGCCGATATTATCAGCGt            <  1:2418641/36‑1 (MQ=17)
                          accaGTTCGCAGAAAATCGCCGATATTATCAGCGTa           >  1:4442088/1‑36 (MQ=255)
                           ccaGTTCGCAGAAAATCGCCGATATTATCAGCGTaa          >  1:127629/1‑36 (MQ=255)
                           ccaGTTCGCAGAAAATCGCCGATATTATCAGCGTaa          <  1:465562/36‑1 (MQ=255)
                           ccaGTTCGAAGAAAATCGCCGATATTATCAGCGTaa          <  1:6517126/36‑1 (MQ=255)
                           ccaGTTCGAAGAAAATCGCCGATATTATCAGCGTaa          <  1:5851482/36‑1 (MQ=255)
                            caGTTCGAAGAAAATCGCCGATATTATCAGCGTAAt         <  1:3567326/36‑1 (MQ=255)
                            caGTTCGAAGAAAATCGCCGATATTATCAGCGTAAt         <  1:4685834/36‑1 (MQ=255)
                             aGTTCGCAGAAAATCGCCGATATTATCAGCGTAAtt        >  1:7352738/1‑36 (MQ=255)
                              gTTCGCAGAAAATCGCCGATATTATCAGCGTAATTg       >  1:2593260/1‑36 (MQ=255)
                              gTTCGCAGAAAATCGCCGATATTATCAGCGTAATTg       >  1:2363466/1‑36 (MQ=255)
                              gTTCGCAGAAAATCGCCGATATTATCAGCGTAATTg       >  1:7009596/1‑36 (MQ=255)
                              gTTCGCAGAAAATCGCCGATATTATCAGCGTAATTg       >  1:490978/1‑36 (MQ=255)
                              gTTCGCAGAAAATCGCCGATATTATCAGCGTAATTg       >  1:5193786/1‑36 (MQ=255)
                              gTTCGCAGAAAATCGCCGATATTATCAGCGTAATTg       >  1:6079202/1‑36 (MQ=255)
                               ttCGCAGAANATCGCCGATATTATCAGCGTAATTGa      <  1:4566097/36‑1 (MQ=255)
                               ttCGCAGAAAATCGCCGATATTATCAGCGTAATTGa      >  1:1801523/1‑36 (MQ=255)
                               ttCGCAGAAAATCGCCGATATGATAAGAGTAAttaa      >  1:6250802/1‑34 (MQ=21)
                               ttCGAAGAAAATCGCCGATATTATCAGCGTAATTGa      <  1:7120395/36‑1 (MQ=255)
                               ttCGAAGAAAATCGCCGATATTATCAGCGTAATTGa      <  1:6417144/36‑1 (MQ=255)
                                tCGCAGAAAATCGCCGATATTATCAGCGTAATTGAc     >  1:6249949/1‑36 (MQ=255)
                                 cGCAGAAAATCGCCGATATTATCAGCGTAATTGACg    <  1:5791229/36‑1 (MQ=255)
                                 cGCAGAAAATCGCCGATATTATCAGCGTAATTGACg    >  1:7497356/1‑36 (MQ=255)
                                  gCAGAAAATCGCCGATATTATCAGCGTAATTGACgg   >  1:2871131/1‑36 (MQ=255)
                                  gCAGAAAATCGACGATATTATCAGCGTAATTGACgg   >  1:6091886/1‑36 (MQ=255)
                                   cAGAAAATCGCCGATATTATCAGCGTAATTGACGGc  >  1:685378/1‑36 (MQ=255)
                                   cAGAAAATCGCCGATATTATCAGCGTAATTGACGGc  >  1:6189693/1‑36 (MQ=255)
                                   |                                   
TGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGC  >  REL606/4581931‑4582001

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