Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,341,8140C→T 7.5% 0.6 ‑3.7 ‑0.060216yjdCpredicted transcriptional regulator
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (4/1):  ref base (14/41):  total (18/42)
Fisher's exact test strand distribution p-value = 2.51e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.88e-01

TGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACAT  >  REL606/4341779‑4341848
                                   |                                  
tGATAACGTGCCCGGCGGTTTTGTTCGGCGGTTTGc                                    <  1:4634125/36‑1 (MQ=19)
  aTAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGt                                  >  1:682120/1‑36 (MQ=255)
   taccGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGtt                                 <  1:4667256/33‑1 (MQ=255)
   tAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGtt                                 <  1:2633079/36‑1 (MQ=255)
   tAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGtt                                 <  1:549455/36‑1 (MQ=255)
   tAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGtt                                 <  1:5913712/36‑1 (MQ=255)
   tAACGTGCCAGCAGCTTTTGTTCGGCGGCTTGCGtt                                 >  1:151375/1‑36 (MQ=255)
   tAACATGCCAGCAGCTTTTGTTCGGCGGTTTGCGtt                                 <  1:4442835/36‑1 (MQ=255)
    aaCGTGCCCGTAGCTTTTGTTCGGCGGTTTGCGttt                                <  1:4895806/36‑1 (MQ=37)
    aaCGTGCCCGCAGCTTTTGTTCGGCGGTTTGCGttt                                <  1:4302609/36‑1 (MQ=255)
    aaCGTGCCAGCAGCTTTTGTTCGGCGGTTTTCGttt                                >  1:4658710/1‑36 (MQ=255)
     aCGTGCCAGCATCTTTTGTTCGGCGGTTTGCGTTTc                               <  1:4606125/36‑1 (MQ=255)
     aCGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTc                               <  1:5148732/36‑1 (MQ=255)
     aCGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTc                               <  1:4561536/36‑1 (MQ=255)
     aCGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTc                               <  1:4052953/36‑1 (MQ=255)
      cGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCg                              <  1:2478966/36‑1 (MQ=255)
       gTGCCAGCAGCTTTTGTTCGGCGGTTTTTGttttta                             >  1:4818712/1‑33 (MQ=25)
       gTGCCAGCAGCTTTTGTTCGGCGGTTTGTGTTTCGt                             <  1:1392938/36‑1 (MQ=255)
       gTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGt                             <  1:4308284/36‑1 (MQ=255)
       gTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGt                             <  1:3679075/36‑1 (MQ=255)
        tGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTc                            <  1:2623906/36‑1 (MQ=255)
        tGCCAGCAGCTTTTGTTCGGCGGGTTTTGTTTTGtt                            >  1:3134941/1‑35 (MQ=16)
         gCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTcc                           <  1:2143777/36‑1 (MQ=255)
         gCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTcc                           <  1:53814/36‑1 (MQ=255)
         gCCAGCAGCTTTTGTTCGGCGGGTTGTGTTTCTTcc                           >  1:1041406/1‑36 (MQ=18)
          ccAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCa                          <  1:2865715/36‑1 (MQ=255)
             gcagcTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCa                       <  1:4387380/36‑1 (MQ=255)
             gcagcTTTTGTTCGGCGGTTTGCGTTTCGTCCAACa                       <  1:1656155/36‑1 (MQ=255)
               cgcTTTTGTTCGGCGTTTTGCGTTTCGTCCAGCATc                     <  1:3799972/35‑1 (MQ=255)
               agcTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATc                     <  1:4160204/36‑1 (MQ=255)
                gcTTTTGTTCTGCGGTTTGCGTTTCGTCCAGCATCa                    <  1:3660208/36‑1 (MQ=255)
                gcTTTTGTTCGGCGGTTTTCGTTTTGTCCCGCATCa                    >  1:5868251/1‑36 (MQ=25)
                  ttttGTTCGGCGGTTTTTGTTTTTTCCAGCATCatc                  >  1:3206483/1‑34 (MQ=16)
                  ttttGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGc                  <  1:1552348/36‑1 (MQ=255)
                  ttttGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGc                  <  1:4351586/36‑1 (MQ=255)
                   tttGTTCGGCGGTTTTCGTTTTGTTCAGCACCAGCt                 >  1:5304495/1‑36 (MQ=16)
                   tttGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCt                 <  1:1811878/36‑1 (MQ=255)
                   tttGTTCGGCGGGTTTCGTTTTTTCCCGCATCAGCt                 >  1:881583/1‑36 (MQ=16)
                   tttGTCCGGCGGTTTGCGTTTCGTCCAGCATCAGcc                 <  1:3660933/36‑2 (MQ=255)
                     tGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTga               <  1:5404175/36‑1 (MQ=255)
                      gTTCGGCGGTTTTCGTTTCGTCCAGCATCCGCTgtc              >  1:4350505/1‑34 (MQ=25)
                      gTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTgac              <  1:107020/36‑1 (MQ=255)
                         cggcggTTTGCGTTTCGTCCATCATCAGCTgccgac           <  1:3785952/36‑1 (MQ=37)
                         cggcggTTTGCGTTTCGTCCAGCATCAGCTgacgac           <  1:54306/36‑1 (MQ=255)
                         cggcggTTTGCGTTTCGTCCAGCATCAGCTgacgac           <  1:1010708/36‑1 (MQ=255)
                         cggcggTTTGCGTTTCGTCCAGCATCAGCTgacgac           <  1:2370646/36‑1 (MQ=255)
                         cggcggTTTGCGTTTCGTCCAGCATCAGCTgacgac           <  1:3347159/36‑1 (MQ=255)
                         cggcggTTTGCGTTTCGTCCAGCATCAGCTgacgac           <  1:4065148/36‑1 (MQ=255)
                         cggcggTTTGCGTTTCGTCCAGCATCAGCTgacgac           <  1:4405815/36‑1 (MQ=255)
                         cggcggTTTGCGTTTCGTCCAGCATCAGCTgacgac           <  1:467471/36‑1 (MQ=255)
                          ggcggTTTTCGTTTCGTGCCCCATTCGCTGacgacg          >  1:5914584/1‑36 (MQ=9)
                           gcggTTTGCGTTTCGTCCAGCATCAGCTGACGACGc         <  1:1341333/36‑1 (MQ=255)
                            cggTTTGCGTTTCGTCCAGCATCAGCTGACGACGcc        <  1:514416/36‑1 (MQ=255)
                             ggTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCa       >  1:1620007/1‑36 (MQ=255)
                              gTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAg      >  1:3694666/1‑36 (MQ=255)
                               tttGCGTTTCGTCCAGCATCAGCTGACGACGCCAGa     <  1:4957975/36‑1 (MQ=255)
                               tttGCGTTTCGTCCAGCATCAGCTGACGACGCCAGa     >  1:2877744/1‑36 (MQ=255)
                               tttGCGTTTCGTCCAGCATCAGCTGACGACGCCAGa     >  1:1584561/1‑36 (MQ=255)
                               tttGCGTTTCGTCCAGAATCAGCTGACGACGCCAGa     >  1:5831168/1‑36 (MQ=255)
                                ttGCGTTTCGTCCAGCATCAGCTGACGACGCCAGAc    >  1:3001845/1‑36 (MQ=255)
                                  gCGTTTCGTCCAGCATCAGCTGACGACGCCAGACat  <  1:1884181/36‑1 (MQ=255)
                                   |                                  
TGATAACGTGCCAGCAGCTTTTGTTCGGCGGTTTGCGTTTCGTCCAGCATCAGCTGACGACGCCAGACAT  >  REL606/4341779‑4341848

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 16 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG