Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,724,1560G→C 3.4% 2.4 ‑14.5 0.0241A84P (GCT→CCT) 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.
Reads supporting (aligned to +/- strand):  new base (0/9):  ref base (189/42):  total (190/51)
Fisher's exact test strand distribution p-value = 4.86e-07
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

CATCAACACCCACAAGAAAGAGCTGATTGCAGAAGCTGCCGTTAGTTTTATTCATGATGGCGACTCAAT  >  REL606/2724122‑2724190
                                  |                                  
catcaccacccACAAGAAAGAGCTGATTGCAGAAg                                    <  1:10436221‑M2/24‑1 (MQ=255)
catcaACACCCACAAGAAAGAGCTGATTGCAGAAg                                    >  1:15570704/1‑35 (MQ=255)
 accaccacccACAAGAAAGAGCTGATTGCAGAAGc                                   <  1:4316565‑M2/25‑1 (MQ=255)
 aacaccacccACAAGAAAGAGCTGATTGCAGAAGc                                   <  1:11242965‑M2/25‑1 (MQ=255)
  tcaACACCCACAAGAAAGAGCTGATTGCAGAAGCt                                  >  1:3717568/1‑35 (MQ=255)
  tcaACACCCACAAGAAAGAGCTGATTGCAGAAGCt                                  >  1:122798/1‑35 (MQ=255)
  tcaACACCCACAAGAAAGAGCTGATTGCAGAAGCt                                  >  1:19533621/1‑35 (MQ=255)
  cccccacccACAAGAAAGAGCTGATTGCAGAAGCt                                  <  1:15432459‑M2/26‑1 (MQ=255)
  ccaccgcccACAAGAAAGAGCTGATTGCAGAAGCt                                  <  1:16485849‑M2/26‑1 (MQ=255)
  ccaccacccACAAGAAAGAGCTGATTGCAGAAGCt                                  <  1:4913352‑M2/26‑1 (MQ=255)
  ccaccacccACAAGAAAGAGCTGATTGCAGAAGCt                                  <  1:12977909‑M2/26‑1 (MQ=255)
  acaccacccACAAGAAAGAGCTGATTGCAGAAGCt                                  <  1:19724864‑M2/26‑1 (MQ=255)
   ccccacccACAAGAAAGAGCTGATTGCAGAAGCTg                                 <  1:7554545‑M2/27‑1 (MQ=255)
   caccacccACAAGAAAGAGCTGATTGCAGAAGCTg                                 <  1:7390228‑M2/27‑1 (MQ=255)
   caACACCCACAAGAAAGAGCTGATTGCAGAAGCTg                                 >  1:14743605/1‑35 (MQ=255)
   caACACCCACAAGAAAGAGCTGATTGCAGAAGCTg                                 >  1:17690919/1‑35 (MQ=255)
   caACACCCACAAGAAAGAGCTGATTGCAGAAGCTg                                 >  1:7743436/1‑35 (MQ=255)
   caACACCCACAAGAAAGAGCTGATTGCAGAAGCTg                                 >  1:1993311/1‑35 (MQ=255)
    cccaaccACAAGAAAGAGCTGATTGCAGAAGCTGc                                <  1:3407889‑M2/28‑1 (MQ=255)
    accacccACAAGAAAGAGCTGATTGCAGAAGCTGc                                <  1:3594885‑M2/28‑1 (MQ=255)
    accacccACAAGAAAGAGCTGATTGCAGAAGCTGc                                <  1:7156724‑M2/28‑1 (MQ=255)
    accacccACAAGAAAGAGCTGATTGCAGAAGCTGc                                <  1:8769321‑M2/28‑1 (MQ=255)
    accacccACAAGAAAGAGCTGATTGCAGAAGCTGc                                <  1:1148497‑M2/28‑1 (MQ=255)
    accacccACAAGAAAGAGCTGATTGCAGAAGCTGc                                <  1:18071998‑M2/28‑1 (MQ=255)
    accacccACAAGAAAGAGCTGATTGCAGAAGATGc                                <  1:16804091‑M2/28‑1 (MQ=255)
    aaCACCCACAAGAAAGAGCTGATTGCAGAAGCTGc                                >  1:1268672/1‑35 (MQ=255)
    aaCACCCACAAGAAAGAGCTGATTGCAGAAGCTGc                                >  1:9792171/1‑35 (MQ=255)
    aaCACCCACAAGAAAGAGCTGATTGCAGAAGCTGc                                >  1:10641320/1‑35 (MQ=255)
    aaCACCCACAAGAAAGAGCTGATTGCAGAAGCTGc                                >  1:3214797/1‑35 (MQ=255)
    aaCACCCACAAGAAAGAGCTGATTGCAGAAGCTGc                                >  1:10114650/1‑35 (MQ=255)
     ccccccACAAGAAAGAGCTGATTGCAGAAGCTGcc                               <  1:6362212‑M2/29‑1 (MQ=255)
     ccacccACAAGAAAGAGCTGATTGCAGAAGCTGcc                               <  1:15895606‑M2/29‑1 (MQ=255)
     ccacccACAAGAAAGAGCTGATTGCAGAAGCTGcc                               <  1:8252350‑M2/29‑1 (MQ=255)
     acacCCACAAGAAAGAGCTGATTGCAGAAGCTGcc                               >  1:20576068/1‑35 (MQ=255)
     acacCCACAAGAAAGAGCTGATTGCAGAAGCTGcc                               >  1:20059891/1‑35 (MQ=255)
      cacccACAAGAAAGAGCTGATTGCAGAAGCTGCCg                              >  1:16153222‑M2/6‑35 (MQ=255)
      cacccACAAGAAAGAGCTGATTGCAGAAGCTGCCg                              >  1:4678642‑M2/6‑35 (MQ=255)
       acccACAAGAAAGAGCTGATTGCAGAATCTGCCGt                             >  1:4581886‑M2/5‑35 (MQ=255)
       acccACAAGAAAGAGCTGATTGCAGAAGCTGCCGt                             >  1:18415485‑M2/5‑35 (MQ=255)
         ccACAAGAAAGAGCTGATTGCAGAAGCTGCCGTTa                           >  1:6633581‑M2/3‑35 (MQ=255)
         ccACAAGAAAGAGCTGATTGCAGAAGCTGCCGTTa                           >  1:9917473‑M2/3‑35 (MQ=255)
         ccACAAGAAAGAGCTGATTGCAGAAGCTGCAGTTa                           <  1:19820487‑M2/33‑1 (MQ=255)
          cACAAGAAAGAGCTGATTGCAGAAGCTGCCGTTAg                          >  1:3776169‑M2/2‑35 (MQ=255)
            caAGAAAGAGCTGATTGCAGAAGCTGCAGTTAGtt                        <  1:2114827/35‑1 (MQ=255)
             aagaaagaGCTGATTGCAGAAGCTGCCGTTAGttt                       >  1:5338063/1‑35 (MQ=255)
             aagaaagaGCTGATTGCAGAAGCTGCCGTTAGttt                       >  1:852527/1‑35 (MQ=255)
             aagaaagaGCTGATTGCAGAAGCTGCCGTTAGttt                       >  1:9133353/1‑35 (MQ=255)
              agaaagaGCTGATTGCAGAAGCTGCCGTTAGtttt                      >  1:6263986/1‑35 (MQ=255)
              agaaagaGCTGATTGCAGAAGCTGCCGTTAGtttt                      >  1:16201619/1‑35 (MQ=255)
               gaaagaGCTGATTGCAGAAGCTGCCGTTCGCTTTa                     >  1:16628520/1‑35 (MQ=37)
               gaaagaGCTGATTGCAGAAGCTGCCGTTAGTTTTa                     >  1:19979551/1‑35 (MQ=255)
               gaaagaGCTGATTGCAGAAGCTGCCGTTAGTTTTa                     >  1:19336120/1‑35 (MQ=255)
               gaaagaGCTGATTGCAGAAGCTGCCGTTAGTTTTa                     >  1:12344291/1‑35 (MQ=255)
               gaaagaGCTGATTGCAGAAGCTGCCGTTAGTTTTa                     >  1:16411484/1‑35 (MQ=255)
                aaagaGCTGATTGCAGAAGCTGCCGTTAGTTTTAt                    >  1:15802077/1‑35 (MQ=255)
                aaagaGCTGATTGCAGAAGCTGCCGTTAGTTTTAt                    >  1:13144714/1‑35 (MQ=255)
                 aagaGCTGATTGCAGAAGCTGCCGGTAGTTTTAtt                   >  1:18813260/1‑35 (MQ=255)
                  agagCTGATTGCAGAAGCTGCCGTTAGTTTTTTTc                  >  1:9612886/1‑35 (MQ=255)
                  agagCTGATTGCAGAAGCTGCCGTTAGTTTTATTc                  >  1:1160611/1‑35 (MQ=255)
                  agagCTGATTGCAGAAGCTGCCGTTAGTTTTATTc                  >  1:14898161/1‑35 (MQ=255)
                   gagCTGATTGCAGAAGCTGCCGTTAGTTTTATTCa                 >  1:90200/1‑35 (MQ=255)
                   gagCTGATTGCAGAAGCTGCCGTTAGTTTTATTCa                 >  1:8543943/1‑35 (MQ=255)
                    agCTGATTGCAGACCCTGCCGTTAGTTTTATTCat                <  1:18238180/35‑1 (MQ=37)
                    agCTGATTGCAGAAGCTGCCGTTAGTTTTATTCat                >  1:12031524/1‑35 (MQ=255)
                    agCTGATTGCAGAAGCTGCCGTTAGTTTTATTCat                >  1:10085723/1‑35 (MQ=255)
                     gCTGATTTCAGAAGCTGCCGTTAGTTTTATTCatg               >  1:598237/1‑35 (MQ=255)
                     gCTGATTGCAGACCCTGCCGTTAGTTTTATTCatg               <  1:19003478/35‑1 (MQ=37)
                      cTGATTGCAGAAGCTGCCGTTAGTTTTATTCatga              >  1:20288598/1‑35 (MQ=255)
                      cTGATTGCAGAAGCTGCCGTTAGTTTTATTCatga              >  1:2221757/1‑35 (MQ=255)
                       tGATTGCAGAAGCTGCCGTTAGTTTTATTCatgat             >  1:9782864/1‑35 (MQ=255)
                       tGATTGCAGAAGCTGCCGTTAGTTTTATTCatgat             >  1:19608350/1‑35 (MQ=255)
                       tGATTGCAGAACCTGCCGTTAGTTTTATTCatgat             <  1:8031432/35‑1 (MQ=255)
                        gATTGCAGACCCTGCCGTTAGTTTTATTCatgatg            <  1:731997/35‑1 (MQ=37)
                        gATTGCAGAAGCTGCCGTTAGTTTTATTCatgatg            >  1:5653454/1‑35 (MQ=255)
                        gATTGCAGAAGCTGCCGTTAGTTTTATTCatgatg            >  1:700608/1‑35 (MQ=255)
                        gATTGCAGAACCTGCCGTTAGTTTTATTCatgatg            <  1:6675460/35‑1 (MQ=255)
                         aTTGCAGAAGCTGCCGTTAGTTTTATTCATGATgg           >  1:13258844/1‑35 (MQ=255)
                          ttGCAGAAGCTGCCGTTAGTTTTATTCATGATggg          >  1:17085235/1‑34 (MQ=255)
                          ttGCAGAAGCTGCCGTTAGTTTTATTCATGATGGc          <  1:11102233/35‑1 (MQ=255)
                          ttGCAGAACCTGCCGTTAGTTTTATTCATGATGGc          <  1:15184220/35‑1 (MQ=255)
                           tGCAGAAGCTGCCGTTAGTTTTATTCATGATGGCg         >  1:17833052/1‑35 (MQ=255)
                           tGCAGAAGCTGCCGTTAGTTTTATTCATGATGGCg         >  1:254482/1‑35 (MQ=255)
                           tGCAGAAGCTGCCGTTAGGTTTATTCATGATGGCg         >  1:8213768/1‑35 (MQ=255)
                            gCAGACCCTGCCGTTAGTTTTATTCATGATGGCGa        <  1:14482664/35‑1 (MQ=37)
                            gCAGAAGCTGCCGTTAGTTTTATTCATGATGGCGa        >  1:3618076/1‑35 (MQ=255)
                            gCAGAAGCTGCCGTTAGTTTTATTCATGATGGCGa        >  1:20425119/1‑35 (MQ=255)
                            gCAGAAGCTGCCGTTAGTTTTATTCATGATGGCGa        >  1:2674981/1‑35 (MQ=255)
                             cAGAACCTGCCGTTAGTTTTATTCATGATGGCGAc       <  1:4982537/35‑1 (MQ=255)
                             cAGAACCTGCCGTTAGTTTTATTCATGATGGCGAc       <  1:17540248/35‑1 (MQ=255)
                              aGAAGCTGCCGTTAGTTTTATTCATGATGGCGACt      >  1:2843002/1‑35 (MQ=255)
                              aGAAGCTGCCGTTAGTTTTATTCATGATGGCGACt      >  1:17427381/1‑35 (MQ=255)
                              aGAACCTGCCGTTAGTTTTATTCATGATGGCGACt      <  1:11770844/35‑1 (MQ=255)
                                aaGCTGCCGTTAGTTTTATTCATGATGGCGACTCa    >  1:13985156/1‑35 (MQ=255)
                                aaGCTGCCGTTAGTTTTATTCATGATGGCGACTCa    >  1:16887181/1‑35 (MQ=255)
                                 aGCTGCCGTTAGTTTTATTCATGATGGCGACTCaa   >  1:3084946/1‑35 (MQ=255)
                                 aGCTGCCGTTAGTTTTATTCATGATGGCGACTCaa   >  1:2313065/1‑35 (MQ=255)
                                  gCTGCCGTTAGTTTTATTCATGATGGCGCCTCaat  >  1:13585171/1‑35 (MQ=255)
                                  gCTGCCGTTAGTTTTATTCATGATGGCGACTCaat  >  1:4122528/1‑35 (MQ=255)
                                  gCTGCCGTTAGTTTTATTCATGATGGCGACTCaat  >  1:14834275/1‑35 (MQ=255)
                                  gCTGCCGTTAGTTTTATTCATGATGGCGACTCaat  <  1:15507151/35‑1 (MQ=255)
                                  |                                  
CATCAACACCCACAAGAAAGAGCTGATTGCAGAAGCTGCCGTTAGTTTTATTCATGATGGCGACTCAAT  >  REL606/2724122‑2724190

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ 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.