Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,660,2300T→G100.0% 7.3 15T332P (ACC→CCC) xylAxylose isomerase
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (9/0):  ref base (4/2):  total (13/2)
Fisher's exact test strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GTACTTTGGCATCGAAGTTCAGACCACCGGTGGTGAAACCGCCTGCTTTGAGAATTTCATACATCACC  >  REL606/3660197‑3660264
                                 |                                  
gTACTTTGGCATCGAAGTTCAGCCCACCGgtggtg                                   <  1:7351476/35‑1 (MQ=255)
gTACTTTGGCATCGAAGTTCAGCCCACCGgtggtg                                   <  1:6816980/35‑1 (MQ=255)
   cTTTGGCATCGAAGTTCAGACCACCGGTGGGGaaa                                >  1:16914230/1‑35 (MQ=255)
   cTTTGGCATCGAAGTTCAGACCACCGGTGGGGaaa                                >  1:5058150/1‑35 (MQ=255)
   cTTGGGCATCGAAGTTCAGCCCACCGGTGGTGaaa                                <  1:1824826/35‑1 (MQ=38)
    tttGGCATCGAAGTTCAGCCCACCGGTGGTGAAAc                               <  1:12262739/35‑1 (MQ=255)
    tttGGCATCGAAGTTCAGACCACCGGTGGGGaaaa                               >  1:760180/1‑34 (MQ=255)
    tttGGCATCGAAGTTCAGACCACCGGTGGGGGAAc                               >  1:12202052/1‑35 (MQ=37)
    tttGGCATCGAAGTTCAGACCACCGGTGGGGAAAc                               >  1:9083929/1‑35 (MQ=255)
     ttGGCATCGAAGTTCAGACCACCGGTGGGGGAAcc                              >  1:12074276/1‑35 (MQ=37)
     ttGGCATCGAAGTTCAGACCACCGGTGGGGGAAcc                              >  1:7155803/1‑35 (MQ=37)
     ttGGCATCGAAGTTCAGACCACCGGTGGGGAAAca                              >  1:6007984/1‑34 (MQ=255)
       ggCATCGAAGTTCAGACCACCGGTGGGGAAACCGc                            >  1:18006598/1‑35 (MQ=255)
       ggCATCGAAGTTCAGACCACCGGTGGAGGAAACGc                            >  1:15147792/1‑35 (MQ=25)
        gCATCGAAGTTCAGACCACCGGTGGGGGAACCGcc                           >  1:5493418/1‑35 (MQ=37)
        gCATCGAAGTTCAGACCACCGGTGGGGAAAACGcc                           >  1:2126261/1‑35 (MQ=37)
        gCATCGAAGTTCAGACCACCGGTGGGGAAAACGcc                           >  1:14463657/1‑35 (MQ=37)
        gCATCGAAGTTCAGACCACCGGTGGGGAAAACGcc                           >  1:10448651/1‑35 (MQ=37)
         cATCGAAGTTCAGACCACCGGTGGGGGAACCGccc                          >  1:17954832/1‑34 (MQ=37)
         cATCGAAGTTCAGACCACCGGTGGGGAAACCGCCt                          >  1:17441631/1‑35 (MQ=255)
         cATCGAAGTTCAGACCACCGGTGGGGAAACCCccc                          >  1:16103988/1‑34 (MQ=37)
         cATCGAAGTTCAGACCACCGGTGGGGAAAACGCCt                          >  1:6301507/1‑35 (MQ=37)
          aTCGAAGTTCAGACCACCGGTGGGGGAACCGCCTg                         >  1:14793622/1‑35 (MQ=37)
           tCGAAGTTCAGACCACCGGTGGGGGAAACGCCTGc                        >  1:14123183/1‑35 (MQ=25)
           tCGAAGTTCAGACCACCGGTGGGGAAAACGCCTGc                        >  1:3977700/1‑35 (MQ=37)
            cGAAGTTCAGACCACCGGTGGGGAAACCGCCTGCt                       >  1:5864991/1‑35 (MQ=255)
            cGAAGTTCAGACCACCGGTGGGGAAACCGCCTGCt                       >  1:20075575/1‑35 (MQ=255)
            cGAAGTTCAGACCACCGGTGGGGAAACCGCCTGCt                       >  1:6806547/1‑35 (MQ=255)
             gAAGTTCAGACCACCGGTGGGGGAACCGCCTGCtt                      >  1:831635/1‑35 (MQ=37)
             gAAGTTCAGACCACCGGTGGGGGAACCGCCTGCtt                      >  1:11602797/1‑35 (MQ=37)
              aaGTTCAGACCACCGGTGGGGAAACCGCCTGCttt                     >  1:3922162/1‑35 (MQ=255)
               aGTTCAGACCACCTGTGGTGAAACCGCCTGCTTTg                    >  1:18018905/1‑35 (MQ=255)
               aGTTCAGACCACCGGTGGGGGAACCGCCTGCtttt                    >  1:9099438/1‑34 (MQ=37)
               aGTTCAGACCACCGGTGGGGGAAACCCCTGCtttt                    >  1:6731719/1‑34 (MQ=16)
                gTTCAGACCACCGGTGGGGAAACCCCCGCCtttag                   >  1:6484091/1‑33 (MQ=16)
                 ttCAGACCACCGGTGGGGGAACCGCCCGCTTTtag                  >  1:7873503/1‑35 (MQ=19)
                    aGACCACCGGTGGGGAAACCGCCTGCTTTGAGAAt               >  1:3678373/1‑35 (MQ=255)
                     gACCACCGGTGTTGAAACCGCCTGCTTTGAGAAtt              >  1:10980168/1‑35 (MQ=255)
                     gACCACCGGTGGGGAAACCGCCTGCTTTTAGAAtt              >  1:14039436/1‑35 (MQ=37)
                          ccGGTGGTGAAACCGCCTGCTTTGAGAATTTCATa         >  1:1206900/1‑35 (MQ=255)
                          ccGGTGGTGAAACCGCCTGCTTTGAGAATTTCATa         >  1:8431158/1‑35 (MQ=255)
                            ggtggtGAAACCGCCTGCTTTGAGAATTTCATACa       >  1:4295019/1‑35 (MQ=255)
                            ggtggaGAACCCGCCTGCTTTGAGAATTTCATACa       <  1:1775487/35‑1 (MQ=37)
                            gggggtaAACCCGCCTGCTTTGAGAATTTCATACa       <  1:9967918/32‑1 (MQ=25)
                             gtggtgAAACCGCCTGCTTTGAGAATTTCATACAt      >  1:9843137/1‑35 (MQ=255)
                              tggtgAAACCGCCTGCTTTGAGAATTTCATACATc     >  1:2138006/1‑35 (MQ=255)
                              tggtgAAACCGCCTGCTTTGAGAATTTCATACATc     >  1:1216119/1‑35 (MQ=255)
                              tggagAACCCGCCTGCTTTGAGAATTTCATACATc     <  1:18063526/35‑1 (MQ=37)
                               ggtgAACCCGCCTGCTTTGAGAATTTCATACATCa    <  1:20206518/35‑1 (MQ=255)
                               ggtgAAACCGCCTGCTTTGAGAATTTCATACATCa    >  1:18241310/1‑35 (MQ=255)
                               ggtgAAACCGCCTGCTTTGAGAATTTCATACATCa    >  1:16615970/1‑35 (MQ=255)
                               ggtgAAACCGCCTGCTTTGAGAATTTCATACATCa    >  1:1044075/1‑35 (MQ=255)
                                gtgAAACCGCCTGCTTTGAGAATTTCATACATCAc   >  1:19305241/1‑35 (MQ=255)
                                gtgAAACCGCCTGCTTTGAGAATTTCATACATCAc   >  1:19266666/1‑35 (MQ=255)
                                gtgAAACCGCCTGCTTTGAGAATTTCATACATCAc   >  1:8240841/1‑35 (MQ=255)
                                gtgAAACCGCCTGCTTTGAGAATTTCATACATCAc   >  1:15103237/1‑35 (MQ=255)
                                gtgAAACCGCCTGCTTTGAGAATTTCATACATCAc   >  1:12399001/1‑35 (MQ=255)
                                gtgAAACCGCATGCTTTGAGAATTTCATACATCAc   <  1:4939667/35‑1 (MQ=255)
                                 tgAAACCGCCTGCTTTGAGAATTTCATACATCAcc  >  1:581769/1‑35 (MQ=255)
                                 tgAAACCGCCTGCTTTGAGAATTTCATACATCAcc  >  1:14790025/1‑35 (MQ=255)
                                 |                                  
GTACTTTGGCATCGAAGTTCAGACCACCGGTGGTGAAACCGCCTGCTTTGAGAATTTCATACATCACC  >  REL606/3660197‑3660264

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