Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,864,0660C→T12.8% 4.3 ‑7.1 0.061intergenic (‑469/‑216)ygdP/mutHdinucleoside polyphosphate hydrolase/DNA mismatch repair protein
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (9/0):  ref base (21/31):  total (30/31)
Fisher's exact test strand distribution p-value = 8.25e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATAACATGGTGT  >  REL606/2864031‑2864101
                                   |                                   
cAAGGCTGTAAATGGAAACAGTAGTGGCGGTTTTTc                                     <  1:605187/36‑1 (MQ=255)
cAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTc                                     <  1:5796273/36‑1 (MQ=255)
cAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTc                                     <  1:3500671/36‑1 (MQ=255)
cAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTc                                     <  1:4219496/36‑1 (MQ=255)
  aGGTTGTAAATGGAAACAGTAGTGTAGGTTTTTcac                                   <  1:4353628/36‑1 (MQ=38)
  aGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTtac                                   >  1:5847607/1‑36 (MQ=255)
  aGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTcac                                   <  1:1490564/36‑1 (MQ=255)
    gCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAg                                 <  1:3398548/36‑1 (MQ=255)
    gCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAg                                 <  1:3845054/36‑1 (MQ=255)
      tGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGtt                               <  1:549586/36‑1 (MQ=255)
       gTAAATGGAAACAGTAGTGGAGGTTTTTTACCGttt                              >  1:4758980/1‑35 (MQ=37)
       gTAAATGGAAACAGTAGTGGAGGTTTTTTACAGTTa                              >  1:3811422/1‑36 (MQ=255)
       gTAAATGGAAACAGTAGTGGAGGTTTTTTACAAttt                              >  1:3911040/1‑35 (MQ=37)
       gTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTa                              <  1:4194575/36‑1 (MQ=255)
        tAAATGGAAACCGTAGTGGAGGTTTTTCACAGTTAt                             <  1:2518081/36‑1 (MQ=255)
        tAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTAt                             <  1:2554100/36‑1 (MQ=255)
         aaaTGGAAACAGTAGTGGAGGTTTTTCACAGTTATc                            >  1:341081/1‑36 (MQ=255)
         aaaTGGAAACAGTAGTGGAGGTTTTTCACAGTTATc                            <  1:3546576/36‑1 (MQ=255)
         aaaTGGAAACAGTAGTGGAGGTTTTTCACAGTTATc                            <  1:605545/36‑1 (MQ=255)
         aaaTGGAAACAGTAGTGGAGGGTTTTCACAGTTATc                            >  1:274957/1‑36 (MQ=255)
         aaaTGGAAACAGGAGTGGAGGGTGTTTACAGTTATc                            >  1:75140/1‑36 (MQ=19)
          aaTGGAAACAGTAGTGGAGGTTTTTCACAGTTATcc                           <  1:5769673/36‑1 (MQ=255)
           aTGGAAACAGTAGTGGAGGTTTTTTACAGTTTTccc                          >  1:2790795/1‑36 (MQ=38)
           aTGGAAACAGTAGTGGAGGTTTTTCACAGTTATccc                          <  1:3488309/36‑1 (MQ=255)
           aTGGAAACAGTAGTGCAGGTTTTTCACATTTATccc                          <  1:2309977/36‑1 (MQ=37)
            tGGAAAAAGTAGTGGAGGATTTTTGCCGTTATGCCa                         >  1:1400957/1‑36 (MQ=9)
                aaCAGTAGTGGAGGTTTTTCACAGTTATCCCAGCtt                     <  1:1224433/36‑1 (MQ=255)
                 aCAGTCGTGGAGGTTTTTCACAGTTATCCCAGCttt                    <  1:2885867/36‑1 (MQ=255)
                 aCAGTAGTGGAGGTTTTTTACAGTTATCCCAGCttt                    >  1:458513/1‑36 (MQ=255)
                 aCAGTAGTGGAGGTTTTTCACAGTTATNNCAGCttt                    >  1:54296/1‑36 (MQ=255)
                 aCAGTAGTGGAGGGTTTTCACAGTTATTCCAACttt                    >  1:4408555/1‑36 (MQ=37)
                  cAGTAGTGGATGTTTTTCACAGTTATCCCAGCTTTc                   <  1:5770711/36‑1 (MQ=255)
                  cAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTc                   <  1:5861069/36‑1 (MQ=255)
                   agtagtGGAGGTTTTTCGCAGTTATCCCAGCTTTCt                  <  1:4875608/36‑1 (MQ=255)
                     tagtGGAGGTTTTTTAAAGTTATTCCAGGTTTCttt                >  1:2067689/1‑34 (MQ=16)
                     tagtGGAGGTTTTTCACAGTTATCCCAGCTTTCtgt                <  1:5708893/36‑1 (MQ=255)
                      agtGGAGGTTTTTCACAGTTATCCCAGCTTTCtgtg               <  1:5189661/36‑1 (MQ=255)
                      agtGGAGGTTTTTCACAGTTATCCCAGCTTTCtgtg               <  1:4274463/36‑1 (MQ=255)
                       gtGGAGGTTTTTTACAGTTATCCCAGCTTTTTGTgg              >  1:2750582/1‑36 (MQ=38)
                       gtGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTgg              >  1:2058973/1‑36 (MQ=255)
                       gtGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTgg              <  1:4739219/36‑1 (MQ=255)
                       gtGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTgg              >  1:890885/1‑36 (MQ=255)
                        tGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGa             >  1:1831605/1‑36 (MQ=255)
                        tGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGa             >  1:2205270/1‑36 (MQ=255)
                        tGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGa             <  1:754704/36‑1 (MQ=255)
                        tGGAGGTTTTTCACAGTTATCCCAGCCTTCTGTGGa             <  1:947404/36‑1 (MQ=255)
                         ggAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGAt            >  1:1872060/1‑36 (MQ=255)
                         ggAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGAt            >  1:4048851/1‑36 (MQ=255)
                         ggAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGAt            <  1:4537687/36‑1 (MQ=255)
                         ggAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGAt            <  1:2699644/36‑1 (MQ=255)
                          gAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATa           <  1:1590889/36‑1 (MQ=255)
                          gAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATa           >  1:3981989/1‑36 (MQ=255)
                           aGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATaa          <  1:1109124/36‑1 (MQ=255)
                           aGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATaa          <  1:740937/36‑1 (MQ=255)
                           aGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATaa          <  1:3512216/36‑1 (MQ=255)
                           aGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATaa          >  1:3639489/1‑36 (MQ=255)
                            ggTTTTTCACAGTTATCCCAGCTTTCTGTGGATAAc         >  1:1551128/1‑36 (MQ=255)
                            ggTTTTTCACAGTTATCCCAGCTTTCTGTGGATAAc         >  1:3541176/1‑36 (MQ=255)
                            ggTTTTTCACAGTTATCCCAGCTTTCTGTGGATAAc         >  1:5265422/1‑36 (MQ=255)
                             gTTTTTCACAGTTATCCCAGCTTTCTGTGGATAACa        >  1:5545923/1‑36 (MQ=255)
                             gTTTTTCACAGTTATCCCAGCTTTCTGTGGATAACa        >  1:5544078/1‑36 (MQ=255)
                             gTTTTTCACAGTTATCCCAGCTTTCTGTGGATAACa        >  1:840865/1‑36 (MQ=255)
                              tttttCACAGTTATCCCTGCTTTCTGTGGATAACAt       <  1:5944689/36‑1 (MQ=255)
                               ttttCACAGTTATCCCAGCTTTCTGTGGATAACATg      <  1:2221794/36‑1 (MQ=255)
                                tttCACAGTTATCCCAGCTTTCTGTGGATAACATgg     >  1:3341548/1‑36 (MQ=255)
                                tttCACAGTTATCCCAGCTTTCTGTGGATAACATgg     <  1:4261158/36‑1 (MQ=255)
                                  tCACAGTTATCCCAGCTTTCTGTGGATAACATGgtg   >  1:3649961/1‑36 (MQ=255)
                                  tCACAGTTATCCCAGCTTTCTGTGGATAACATGgtg   >  1:3387614/1‑36 (MQ=255)
                                   cacaGTTATCCCAGCTTTCTGTGGATAACATGgtgt  <  1:3749862/36‑1 (MQ=255)
                                   cacaGTTATCCCAGCTTTCTGTGGATAACATGgtgt  >  1:657607/1‑36 (MQ=255)
                                   |                                   
CAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATAACATGGTGT  >  REL606/2864031‑2864101

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