Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,750,8160A→C 4.4% 0.7 ‑30.1 0.0232Y540S (TAC→TCC) fhlADNA‑binding transcriptional activator
Rejected: E-value exceeds prediction threshold.
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/16):  ref base (190/26):  total (190/42)
Fisher's exact test strand distribution p-value = 8.40e-14
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

GTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCA  >  REL606/2750782‑2750850
                                  |                                  
gTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTa                                    <  1:8784092/35‑1 (MQ=255)
gTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTa                                    >  1:19860446/1‑35 (MQ=255)
gTCGCCGACCGTGAGTTCCGTAGAGATCTCTATTa                                    <  1:2195168/35‑1 (MQ=255)
 tCGCCGACCGTGAGTTCCGTAGCGATCTCTATTAc                                   >  1:3655817/1‑35 (MQ=255)
 tCGCCGACCGTGAGTTCCGTAGCGATCTCTATTAc                                   <  1:4930739/35‑1 (MQ=255)
 tCGCCGACCGTGAGTTCCGTAGCGATCTCTATTAc                                   <  1:6563476/35‑1 (MQ=255)
 tCGCCGACCGTGAGTTCCGTAGCGATCTCTATTAc                                   >  1:19770100/1‑35 (MQ=255)
 tCGCCGACCGTGAGTTCCGTAGCGATCTCTATTAc                                   <  1:18921557/35‑1 (MQ=255)
 tCGCCGACCGTGAGTTCCGTAGCGATCTCTATTAc                                   <  1:1008013/35‑1 (MQ=255)
 tCGCCGACCGTGAGTTCCGTAGCGATCTCTATTAc                                   <  1:17772135/35‑1 (MQ=255)
  cGCCGACCGTTAGTTCCGTAGCGATCTCTATTAcc                                  >  1:3146076/1‑35 (MQ=255)
  cGCCGACCGTGAGTTCCGTAGCGATCTCTATTAcc                                  >  1:3381714/1‑35 (MQ=255)
  cGCCGACCGTGAGTTCCGTAGCGATCTCTATTAcc                                  >  1:5150833/1‑35 (MQ=255)
  cGCCGACCGTGAGTTCCGTAGCGATCTCTATTAcc                                  >  1:13937704/1‑35 (MQ=255)
  cGCCGACCGTGAGTTCCGTAGCGATCTCTATTAcc                                  <  1:13829340/35‑1 (MQ=255)
   gCCGACCGTGAGTTCCGTAGCGATCTCTATTACCg                                 <  1:7468148/35‑1 (MQ=255)
   gCCGACCGTGAGTTCCGTAGCGATCTCTATTACCg                                 >  1:754058/1‑35 (MQ=255)
   gCCGACCGTGAGTTCCGTAGCGATCTCTATTACCg                                 <  1:6194970/35‑1 (MQ=255)
    ccGCCCGTGAGTTCCGTAGCGATCTCTATTACCGc                                <  1:19105114/35‑1 (MQ=255)
    ccGACCGTGAGTTCCGTAGCGATCTCTATTACCGc                                >  1:4795532/1‑35 (MQ=255)
    ccGACCGTGAGTTCCGTAGCGATCTCTATTACCGc                                >  1:6017864/1‑35 (MQ=255)
     cGACCGTGAGTTCCGTAGCGATCTCTATTACCGcc                               >  1:3093451/1‑35 (MQ=255)
     cGACCGTGAGTTCCGTAGCGATCTCTATTACCGcc                               >  1:11777490/1‑35 (MQ=255)
     cGACCGTGAGTTCCGTAGCGATCTCTATTACCGcc                               >  1:7848754/1‑35 (MQ=255)
      gACCGTGAGTTCCGTAGCGATCTCTATTACCGCCt                              >  1:4711865/1‑35 (MQ=255)
      gACCGTGAGTTCCGTAGCGATCTCTATTACCGCCt                              >  1:12212220/1‑35 (MQ=255)
      gACCGTGAGTTCCGTAGCGATCTCTATTACCGCCt                              >  1:8003024/1‑35 (MQ=255)
       aCCGTGAGTTCCGTAGCGATCTCTATTACCGCCTg                             >  1:390874/1‑35 (MQ=255)
       aCCGTGAGTTCCGTAGCGATCTCTATTACCGCCTg                             >  1:5433256/1‑35 (MQ=255)
       aCCGTGAGTTCCGTAGCGATCTCTATTACCGCCTg                             >  1:20421212/1‑35 (MQ=255)
        ccGTGAGTTCCGTAGCGATCTCTATTACCGCCTGa                            >  1:6654030/1‑35 (MQ=255)
        ccGTGAGTTCCGTAGCGATCTCTATTACCGCCTGa                            >  1:9554751/1‑35 (MQ=255)
        ccGTGAGTTCCGTAGCGATCTCTATTACCGCCTGa                            >  1:10902437/1‑35 (MQ=255)
        ccGTGAGTTCCGTAGCGATCTCTATTACCGCCTGa                            >  1:16911773/1‑35 (MQ=255)
         cGTGAGTTCCGTATCGATCTCTATTACCGCCTGaa                           >  1:14652934/1‑35 (MQ=255)
         cGTGAGTTCCGTAGCGATCTCTATTACCGCCTGaa                           >  1:18403051/1‑35 (MQ=255)
          gTGAGTTCCGTAGCGATCTCTATTACCGCCTGAAc                          >  1:558058/1‑35 (MQ=255)
          gTGAGTTCCGTAGCGATCTCTATTACCGCCTGAAc                          >  1:3312353/1‑35 (MQ=255)
          gTGAGTTCCGTAGCGATCTCTATTACCGCCTGAAc                          >  1:8529125/1‑35 (MQ=255)
           tGAGTTCCGTAGCGATCTCTATTACCGCCTGAACg                         >  1:11944487/1‑35 (MQ=255)
           tGAGTTCCGTAGCGATCTCTATTACCGCCTGAACg                         >  1:8166572/1‑35 (MQ=255)
            gAGTTCCGTAGCGATCTCTATTACCGCCTGAACGt                        >  1:17854657/1‑35 (MQ=255)
            gAGTTCCGTAGCGATCTCTATTACCGCATGAACGt                        <  1:225565/35‑1 (MQ=255)
             aGTTCCGTAGCGATCTCTATTACCGCCTGAACGTa                       >  1:18529035/1‑35 (MQ=255)
             aGTTCCGTAGCGATCTCTATTACCGCCTGAACGTa                       >  1:12417063/1‑35 (MQ=255)
             aGTTCCGTAGCGATCTCTATTACCGCCTGAACGTa                       >  1:1603349/1‑35 (MQ=255)
             aGTTCCGTAGCGATCTCTATTACCGCCTGAACGTa                       >  1:7390132/1‑35 (MQ=255)
             aGTTCCGTAGCGATCTCTATTACCGCCTGAACGTa                       >  1:2984608/1‑35 (MQ=255)
               ttCCGTAGCGATCTCTATTACCGCCTGAACGTAtt                     >  1:2763191/1‑35 (MQ=255)
                 ccGTAGCGATCTCTATTACCGCCTGAACGTATTcc                   >  1:9320687/1‑35 (MQ=255)
                 ccGTAGCGATCTCTATTACCGCCTGAACGTATTcc                   >  1:17577953/1‑35 (MQ=255)
                  cGTAGCGATCTCTATTACCGCCTGAACGTATTccc                  >  1:19496014/1‑35 (MQ=255)
                  cGTAGCGATCTCTATTACCGCCTGAACGTATTccc                  >  1:16186466/1‑35 (MQ=255)
                   gTAGCGATCTCTATTACCGCCTGAACGTATTCCCg                 >  1:11058281/1‑35 (MQ=255)
                   gTAGCGATCTCTATTACCGCCTGAACGTATTCCCg                 >  1:5853831/1‑35 (MQ=255)
                    tAGCGATCTCTATTACCGCCTGAACGTATTCCCGa                >  1:10693319/1‑35 (MQ=255)
                    tAGCGATCTCTATTACCGCCTGAACGTATTCCCGa                >  1:1305870/1‑35 (MQ=255)
                    tAGCGATCTCTATTACCGCCTGAACGTATTCCCGa                >  1:6144740/1‑35 (MQ=255)
                     aGCGTTCTTTTTTCCCGCCTGAACGTATTCCCGAt               <  1:4255953/35‑1 (MQ=16)
                     aGCGATCTCTATTACCGCCTGAACGTATTCCCGAt               >  1:17228024/1‑35 (MQ=255)
                     aGCGATCTCTATTACCGCCTGAACGTATTCCCGAt               >  1:12037561/1‑35 (MQ=255)
                     aGCGATCTCTATTACCGCCTGAACGTATTCCCGAt               >  1:4320537/1‑35 (MQ=255)
                      gCGATCTCTATTCCAGCCTGAACGTATTCCCGAtt              <  1:18266408/35‑1 (MQ=37)
                      gCGATCTCTATTACCGCCTGAACGTATTCCCGAtt              >  1:6839034/1‑35 (MQ=255)
                      gCGATCTCTATTACCGCCTGAACGTATTCCCGAtt              >  1:14350798/1‑35 (MQ=255)
                       cGATCTCTATTCCCGCCTGAACGTATTCCCGATTc             <  1:7099606/35‑1 (MQ=255)
                       cGATCTCTATTACCGCCTGAACGTATTCCCGATTc             >  1:11343185/1‑35 (MQ=255)
                       cGATCTCTATTACCGCCTGAACGTATTCCCGATTc             >  1:2368706/1‑35 (MQ=255)
                        gATTTCTTTTCCCGCCTGAACGTATTCCCGATTCa            <  1:19974452/35‑1 (MQ=21)
                        gATCTCTATTACCGCCTGAACGTATTCCCGATTCa            >  1:6397641/1‑35 (MQ=255)
                        gATCTCTATTACCGCCTGAACGTATTCCCGATTCa            >  1:4423550/1‑35 (MQ=255)
                        gATCTCTATTACCGCCTGAACGTATTCCCGATTCa            >  1:18072988/1‑35 (MQ=255)
                         tcttctTTTCCCGCCTGAACGTATTCCCGATTCAc           <  1:12967297/32‑1 (MQ=21)
                         aTCTCTATTACCGCCTGAACGTATTCCCGATTCAc           >  1:6920583/1‑35 (MQ=255)
                         aTCTCTATTACCGCCTGAACGTATTCCCGATTCAc           >  1:10002627/1‑35 (MQ=255)
                          tctcTATTACCGCCTGAACGTATTCCCGATTCAcc          >  1:10595106/1‑35 (MQ=255)
                           ctttATTCCCGCCTGAACGTATTCCCGATTCACCt         <  1:9115727/32‑1 (MQ=25)
                           ctttATTCCCGCCTGAACGTATTCCCGATTCACCt         <  1:11151828/32‑1 (MQ=25)
                           ctctATTCCCGCCTGAACGTATTCCCGATTCACCt         <  1:19413268/35‑1 (MQ=255)
                           ctctATTACCGCCTGAACGTATTCCCGATTCACCt         >  1:8268196/1‑35 (MQ=255)
                           ctctATTACCGCCTGAACGTATTCCCGATTCACCt         >  1:8677001/1‑35 (MQ=255)
                           ctctATTACCGCCTGAACGTATTCCCGATTCACCt         >  1:2127141/1‑35 (MQ=255)
                            tctATTACCGCCTGAACGTATTCCCGATTCACCTg        >  1:17975079/1‑35 (MQ=255)
                            cctATTACCGCCTGAACGTATTCCCGATTCACCTg        >  1:16465849/2‑35 (MQ=255)
                             ctATTACCGCCTGAACGTATTCCCGATTCACCTgc       >  1:5942607/1‑35 (MQ=255)
                             ctATTACCGCCTGAACGTATTCCCGATTCACCTgc       >  1:12485529/1‑35 (MQ=255)
                             ctATTACCGCCTGAACGTATTCCCGATTCACCTgc       >  1:20283211/1‑35 (MQ=255)
                              ttttACCGCCTGAACGTATTCCCGATTCACGTgcc      <  1:14963496/33‑1 (MQ=37)
                              tATTACCGCCTGAACGTATTCCCGATTCACCTgcc      >  1:13444336/1‑35 (MQ=255)
                              tATTACCGCCTGAACGTATTCCCGACTCATCTgcc      >  1:16109855/1‑35 (MQ=25)
                               aTTACCGCCTGAACGTATTCCCGATTCACCTgccg     >  1:3752167/1‑35 (MQ=255)
                                ttACCGCCTGAACGTATTCCCGATTCACCTgccgc    >  1:15857354/1‑35 (MQ=255)
                                 tACCGCCTGAACGTATTCCCGATTCACCTgccgcc   >  1:17397111/1‑35 (MQ=255)
                                 tACCGCCTGAACGTATTCCCGATTCACCTgccgcc   >  1:13703072/1‑35 (MQ=255)
                                  aCCGCCTGAACGTATTCCCGATTCACCTGCCGCCa  >  1:17379054/1‑35 (MQ=255)
                                  aCCGCCTGAACGTATTCCCGATTCACCTGCCGCCa  >  1:7274107/1‑35 (MQ=255)
                                  aCCGCCTGAACGTATTCCCGATTCACCTGCCGCCa  >  1:2706620/1‑35 (MQ=255)
                                  aCCGCCTGAACGTATTCCCGATTCACCTGCCGCCa  >  1:15370106/1‑35 (MQ=255)
                                  aCCGCCTGAACGTATTCCCGATTCACCTGCCGCCa  >  1:11206063/1‑35 (MQ=255)
                                  aCCGCCTGAACGTATTCCCGATTCACCTGCCGCCa  >  1:5317774/1‑35 (MQ=255)
                                  |                                  
GTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCA  >  REL606/2750782‑2750850

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