Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,724,4230A→C 8.5% 8.7 ‑47.9 0.0263T173P (ACC→CCC) srlRDNA‑bindng transcriptional repressor
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Rejected: Polymorphic base substitution creates a homopolymer stretch.
Reads supporting (aligned to +/- strand):  new base (0/39):  ref base (171/52):  total (171/92)
Fisher's exact test strand distribution p-value = 1.60e-21
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 299 total aligned reads displayed.

CGAGCATTTCACCTTCGATAAATTGTTTATGGGCACCGACGGCATCGATCTCAATGCGGGCGTAACCAC  >  REL606/2724389‑2724457
                                  |                                  
cGAGCATTTCACCTTCGATAAATTGTTTATGGGCa                                    <  1:1828396/35‑1 (MQ=255)
cGAGCATTTCACCTTCGATAAATTGTTTATGGGCa                                    >  1:15758580/1‑35 (MQ=255)
cGAGCATTTCACCTTCGATAAATTGTTTATGGGCa                                    >  1:9153987/1‑35 (MQ=255)
cGAGCATTTCACCTTCGATAAATTGTTTATGGGCa                                    <  1:2577720/35‑1 (MQ=255)
cGAGCATTTCACCTTCGATAAATTGTTTATGGGCa                                    >  1:4184760/1‑35 (MQ=255)
cGAGCATTTCACCTTCGATAAATTGTTTATGGGCa                                    >  1:6472668/1‑35 (MQ=255)
cGAGCATTTCAACTTCGATAAATTGTTTATGGGCa                                    <  1:13936453/35‑1 (MQ=255)
 gggCATTTCCCCTTCGATAAATTGTTTATGGGCAc                                   <  1:23214/33‑1 (MQ=38)
 gAGCATTTCACCTTCGATAAATTGTTTATGGGCAc                                   >  1:18672201/1‑35 (MQ=255)
 gAGCATTTCACCTTCGATAAATTGTTTATGGGCAc                                   >  1:19922742/1‑35 (MQ=255)
 gAGCATTTCACCTTCGATAAATTGTTTATGGGCAc                                   >  1:16542815/1‑35 (MQ=255)
 gAGCATTTCACCTTCGATAAATTGTTTATGGGCAc                                   <  1:13202474/35‑1 (MQ=255)
 gAGCATTTCACCTTCGATAAATTGTTTATGGGCAc                                   >  1:10047716/1‑35 (MQ=255)
  aGCATTTCCCCTTCGATAAATTGTTTATGGGCAcc                                  <  1:909628/35‑1 (MQ=255)
  aGCATTTCCCCTTCGATAAATTGTTTATGGGCAcc                                  <  1:1268961/35‑1 (MQ=255)
  aGCATTTCACCTTCGATAAATTGTTTATGGGCAcc                                  <  1:6243484/35‑1 (MQ=255)
  aGCATTTCACCTTCGATAAATTGTTTATGGGCAcc                                  <  1:14291428/35‑1 (MQ=255)
   gCATTTCCCCTTCGATAAATTGTTTATGGGCACCg                                 <  1:20385350/35‑1 (MQ=255)
   gCATTTCCCCTTCGATAAATTGTTTATGGGCACCg                                 <  1:16088396/35‑1 (MQ=255)
   gCATTTCACCTTCGATAAATTGTTTATGGTCACCg                                 >  1:1669307/1‑35 (MQ=255)
   gCATTTCACCTTCGATAAATTGTTTATGGGCACCg                                 >  1:15959378/1‑35 (MQ=255)
   gCATTTCACCTTCGATAAATTGTTTATGGGCACCg                                 >  1:9340219/1‑35 (MQ=255)
   gCATTTCACCTTCGATAAATTGTTTATGGGCACCg                                 >  1:15659702/1‑35 (MQ=255)
   gCATTTCACCTTCGATAAATTGTTTATGGGCACCg                                 >  1:6764591/1‑35 (MQ=255)
   gCATTTCACCTTCGATAAATTGTTTATGGGCACCg                                 <  1:12472625/35‑1 (MQ=255)
    cATTTCCCCTTCGATAAATTGTTTATGGGCACCGa                                <  1:19145836/35‑1 (MQ=255)
    cATTTCACCTTCGATAAATTGTTTATGGGCACCGa                                <  1:1201948/35‑1 (MQ=255)
    cATTTCACCTTCGATAAATTGTTTATGGGCACCGa                                >  1:16560350/1‑35 (MQ=255)
    cATTTCACCTTCGATAAATTGTTTATGGGCACCGa                                >  1:7218006/1‑35 (MQ=255)
    cATTTCAACTTCGATAAATTGTTTATGGGCACCGa                                <  1:14834217/35‑1 (MQ=255)
     aTTTCCCCTTCGATAAATTGTTTATGGGCACCGAc                               <  1:13580489/35‑1 (MQ=255)
      tttCACCTTCGATAAATTGTTTATGGGCACCGCCg                              >  1:10292230/1‑35 (MQ=255)
       ttCACCTTCGATAAATTGTTTATGGGCACCGACgg                             >  1:6375954/1‑35 (MQ=255)
       ttCACCTTCGATAAATTGTTTATGGGCACCGACgg                             >  1:5551480/1‑35 (MQ=255)
         ccccTTCGATAAATTGTTTATGGGCACCGACGGCa                           <  1:17397909/33‑1 (MQ=255)
         cACCTTCGATAAATTGTTTATGGGCACCGACGGCa                           >  1:13118329/1‑35 (MQ=255)
           ccTTCGATAAATTGTTTATGGGCACCGACGGCATc                         >  1:8715081/1‑35 (MQ=255)
           ccTTCGATAAATTGTTTATGGGCACCGACGGCATc                         >  1:10202168/1‑35 (MQ=255)
            cTTCGATAAATTGTTTATGGGCACCGACGGCATCg                        >  1:10775840/1‑35 (MQ=255)
            cTTCGATAAATTGTTTATGGGCACCGACGGCATCg                        >  1:4697738/1‑35 (MQ=255)
               cGATAAATTGTTTATGGGCACCGCCCGCATTGaga                     <  1:9755981/35‑3 (MQ=21)
               cGATAAATTGTTTATGGGCACCGACGGCATCGAtc                     >  1:4953406/1‑35 (MQ=255)
               cGATAAATTGTTTATGGGCACCGACGGCATCGAtc                     >  1:7352697/1‑35 (MQ=255)
               cGAAAAATTGTTTATGGGCACCGACGGCATCGAtc                     >  1:16811313/1‑35 (MQ=255)
                gATAAATTTTTTGTGGGCCCCGACGGCATCGAtct                    <  1:20073535/35‑1 (MQ=25)
                gATAAATTGTTTATGGGCCCCGACGGCATCGAtct                    <  1:5996820/35‑1 (MQ=255)
                gATAAATTGTTTATGGGCACCGACGGCATCGAtct                    >  1:5101506/1‑35 (MQ=255)
                gATAAATTGTTTATGGGCACCGACGGCATCGAtct                    >  1:16329218/1‑35 (MQ=255)
                gATAAATTGTTTATGGGCACCGACGGCATCGAtct                    >  1:17311218/1‑35 (MQ=255)
                gATAAATTGTTTATGGGCACCGACGGCATCGAtct                    >  1:18492599/1‑35 (MQ=255)
                gATAAATTGTTTATGGGCACCGACGGCATCGAtct                    >  1:18151980/1‑35 (MQ=255)
                 aTAAATTGTTTGTGGGCCCCGACGGCATCGAtctc                   <  1:18977863/35‑1 (MQ=37)
                 aTAAATTGTTTATGGGCACCGACGGCATCGAtctc                   >  1:9026694/1‑35 (MQ=255)
                 aTAAATTGTTTATGGGCACCGACGGCATCGAtctc                   >  1:13793969/1‑35 (MQ=255)
                 aTAAATTGTTTATGGGCACCGACGGCATCGAtctc                   >  1:12741056/1‑35 (MQ=255)
                 aTAAATTGTTTATGGGCACCGACGGCATCGAtctc                   >  1:5890729/1‑35 (MQ=255)
                  tAAATTGTTTATGGGCACCGACGGCATCGATCTCa                  >  1:17699489/1‑35 (MQ=255)
                  tAAATTGTTTATGGCCCCCGACGGCATCGATCTCa                  <  1:11738498/35‑1 (MQ=37)
                   aaaTTGTTTATGGGCACCGACGGCATCGATCTCaa                 >  1:15531428/1‑35 (MQ=255)
                    aaTTGTTTATGGGCCCCGACGGCATCGATCTCAAt                <  1:1232660/35‑1 (MQ=255)
                    aaTTGTTTATGGGCCCCGACGGCATCGATCTAAAt                <  1:9401354/35‑1 (MQ=37)
                      ttttttATGGGCCCCGACGGCATCGATCTCAATGc              <  1:4109452/32‑1 (MQ=37)
                      ttGTTTATGGGCACCGACGGCATCGATCTCAATGc              >  1:3509572/1‑35 (MQ=255)
                       tGTTTATGGGCCCCGACGGCATCGATCTCAATGCg             <  1:12155031/35‑1 (MQ=255)
                       tGTTTATGGGCACCGACGGCATCGATCTCAATGCg             >  1:3749236/1‑35 (MQ=255)
                       tGTTTATGGGCACCGACGGCATCGATCTCAATGCg             >  1:3910147/1‑35 (MQ=255)
                       tGTTTATGGACCCCGACGGCATCGATCTCAATGCg             <  1:13718249/35‑1 (MQ=37)
                       tGTTTAGGGGCACCGACGGCATCGATCTCAATGCg             <  1:5276307/35‑1 (MQ=255)
                        ttttAGGGGCACCGACGGCATCGATCTCAATGCgg            <  1:8548709/34‑1 (MQ=255)
                        gTTTATGGGCACCGACGGCATCGATCTCAATGCgg            >  1:9522087/1‑35 (MQ=255)
                         tttATGTGCACCGACGGCATCGATCTCAATGCgtg           >  1:1437176/1‑33 (MQ=37)
                         tttATGGGCACCGACGGCATCGATCTCAATGCggg           >  1:7872341/1‑35 (MQ=255)
                         tttATGGGCACCGACGGCATCGATCTCAATGCggg           >  1:7577906/1‑35 (MQ=255)
                         tttATGGGCACCGACGGCATCGATCTCAATGCggg           >  1:20177643/1‑35 (MQ=255)
                         tttATGGGCACCGACGGCATCGATCTCAATGCggg           <  1:2882302/35‑1 (MQ=255)
                          ttATGGGCACCGACGGCATCGATCTCAATGCGGGc          >  1:17517451/1‑35 (MQ=255)
                          ttAGGGGCCCCGACGGCATCGATCTCAATGCGGGc          <  1:18835196/35‑1 (MQ=37)
                           tATGGGCACCGACGGCATCGATCTCAATGCGGGCg         >  1:16229046/1‑35 (MQ=255)
                           tATGGGCACCGACGGCATCGATCTCAATGCGGGCg         >  1:3275102/1‑35 (MQ=255)
                           tATGGGCACCGACGGCATCGATCTCAATGCGGGCg         >  1:18064025/1‑35 (MQ=255)
                           tATGGGCACCGACGGCATCGATCTCAATGCGGGCg         >  1:17892015/1‑35 (MQ=255)
                            gtGGGCCCCGACGGCATCGATCTCAATGCGGGCGt        <  1:4383321/34‑1 (MQ=255)
                            attggCACCGACGGCATCGATCTCAATGCGGGCGt        >  1:2123233/4‑35 (MQ=255)
                            aTGGGCCCCGACGGCATCGATCTCAATGCGGGCGt        <  1:17249833/35‑1 (MQ=255)
                            aTGGGCACCGACTGCATCGATCTCAATGCGGGCGt        >  1:19339393/1‑35 (MQ=255)
                            aTGGGCACCGACGGCATCGATCTCAATGCGGGCGt        >  1:14523774/1‑35 (MQ=255)
                             tGGGCACCGACGGCATCGATCTCAATGCGGGCgca       >  1:2289124/1‑33 (MQ=255)
                             tGGGCACCGACGGCATCGATCTCAATGCGGGCGTa       >  1:10243953/1‑35 (MQ=255)
                             tGGGCACCGACGGCATCGATCTCAATGCGGGCGTa       >  1:11627559/1‑35 (MQ=255)
                             tGGGCACCGACGGCATCGATCTCAATGCGGGCGTa       >  1:11163847/1‑35 (MQ=255)
                             tGGGCACCGACGGCATCGATCTCAATGCGGGCGTa       >  1:15285990/1‑35 (MQ=255)
                             ggggCCCCGACGGCATCGATCTCAATGCGGGCGTa       <  1:6650286/34‑1 (MQ=255)
                              gggCACCGACGGCATCGATCTCAATGCGGGCGTaa      >  1:13007088/1‑35 (MQ=255)
                               ggCACCGACGGCATCGATCTCAATGCGGGCGTAac     >  1:801130/1‑35 (MQ=255)
                               ggCACCGACGGCATCGATCTCAATGCGGGCGTAac     <  1:3133704/35‑1 (MQ=255)
                                 cACCGACGGCATCGATCTCAATGCGGGCGTAacca   >  1:8231965/1‑35 (MQ=255)
                                 cACCGACGGCATCGATCTCAATGCGGGCGTAacca   >  1:4560641/1‑35 (MQ=255)
                                 cACCGACGGCATCGATCTCAATGCGGGCGTAacca   >  1:3620290/1‑35 (MQ=255)
                                  aCCGACGGCATCGATCTCAATGCGGGCGTAaccac  >  1:19657384/1‑35 (MQ=255)
                                  aCCGACGGCATCGATCTCAATGCGGGCGTAaccac  >  1:17365023/1‑35 (MQ=255)
                                  |                                  
CGAGCATTTCACCTTCGATAAATTGTTTATGGGCACCGACGGCATCGATCTCAATGCGGGCGTAACCAC  >  REL606/2724389‑2724457

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.