Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,766,9581.→C 1.5% 5.8 0.0 ‑0.2194coding (394/711 nt)ispD2‑C‑methyl‑D‑erythritol 4‑phosphate cytidylyltransferase
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Rejected: Polymorphic indel expands or contracts a homopolymer stretch.
Reads supporting (aligned to +/- strand):  new base (2/1):  ref base (106/85):  total (108/86)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.63e-01

Only 100 of 239 total aligned reads displayed.

CGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCCAAC  >  REL606/2766924‑2766992
                                   |                                  
cGTTTCATAGTATCGCGCACTGGTGCGGCGAGGAt                                     <  1:16286370/35‑1 (MQ=255)
 gTTTCATAGTATCGCGCACTGGTGCGGCGAGGATc                                   <  1:18062561/35‑1 (MQ=255)
 gTTTCATAGTATCGCGCACTGGTGCGGCGAGGATc                                   <  1:1579840/35‑1 (MQ=255)
 gTTTCATAGTATAGCGCACTGGTGCGGCGAGGATc                                   <  1:17687320/35‑1 (MQ=255)
  tttCATAGTATCGCGCACTGGTGCGGCGAGGATcc                                  <  1:6370505/35‑1 (MQ=255)
     cATATTATCGCGCACTGGTGCGGCGAGGATccccc                               <  1:7781124/35‑1 (MQ=255)
     cATAGTATCGCGCACTGGTGCGGCGAGGAAccccc                               >  1:3004864/1‑35 (MQ=255)
      aTAGTATCGCGCACTGGTGCGGCGAGGATcccccc                              <  1:19218204/35‑1 (MQ=255)
      aTAGTATCGCGCACTGGTGCGGCGAGGATcccccc                              <  1:18132207/35‑1 (MQ=255)
      aTAGTATCGCGCACTGGTGCGGCGAGGATcccccc                              <  1:15415644/35‑1 (MQ=255)
      aTAGTATCGCGCACTGGTGCGGCGAGGATcccccc                              <  1:11987926/35‑1 (MQ=255)
       tAGTATCGCGCACTGGTGCGGCGAGGATccccccc                             <  1:19174086/35‑1 (MQ=255)
       tAGTATCGCGCACTGGTGCGGCGAGGAAccccccc                             >  1:2385676/1‑35 (MQ=255)
       tAGGATCGCGCACTGGTGCGGCGAGGATccccccc                             <  1:316813/35‑1 (MQ=255)
        aGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCg                            >  1:20356626/1‑35 (MQ=255)
        aGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCg                            <  1:15180614/35‑1 (MQ=255)
        aGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCg                            >  1:7625682/1‑35 (MQ=255)
        aGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCg                            >  1:113724/1‑35 (MQ=255)
        aGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCg                            >  1:10585397/1‑35 (MQ=255)
         gTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGt                           >  1:7211183/1‑35 (MQ=255)
          tATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTg                          <  1:16472673/35‑1 (MQ=255)
           aTCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGc                         >  1:3773285/1‑35 (MQ=255)
           aTCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGc                         >  1:2745208/1‑35 (MQ=255)
           aTCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGc                         <  1:9628448/35‑1 (MQ=255)
            tCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCg                        <  1:17099643/35‑1 (MQ=255)
            tCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCg                        <  1:17269693/35‑1 (MQ=255)
            tCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCg                        <  1:6804407/35‑1 (MQ=255)
             cgcgCACTGGTGCGGCGAGGATCCCCCCCGTGCgg                       >  1:8919144/1‑35 (MQ=255)
              gcgcACTGGTGCGGCGAGGATCCCCCCCGTGCGGc                      <  1:12425086/35‑1 (MQ=255)
              gcgcACTGGTGCGGCGAGGATCCCCCCCGTGCGGc                      <  1:10571102/35‑1 (MQ=255)
               cgcACTGGTGCGGCGAGGATCCCCCCCGTGCGGCt                     <  1:5207385/35‑1 (MQ=255)
                gcACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTg                    <  1:7671825/35‑1 (MQ=255)
                gcACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTg                    <  1:15940727/35‑1 (MQ=255)
                gcACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTg                    >  1:10651071/1‑35 (MQ=255)
                gcACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTg                    >  1:1861345/1‑35 (MQ=255)
                 cACTGGTGCGGCGAGGATCCCCCCCGTGCTGCTgg                   >  1:6721158/1‑35 (MQ=255)
                 cACTGGTGCGGCGAGGATCCCCCCCGTGCGGCtcg                   >  1:12309459/1‑33 (MQ=255)
                 cACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTgg                   >  1:6229571/1‑35 (MQ=255)
                 cACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTgg                   >  1:8696497/1‑35 (MQ=255)
                 cACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTgg                   <  1:9979467/35‑1 (MQ=255)
                 cACTGGTGCGGCGAGGATCCCCCCCGCGCGGCTgg                   >  1:20281290/1‑35 (MQ=255)
                 cACTGGTGCGGCGACGATCCCCCCCGTGCGGCTgg                   <  1:4786713/35‑1 (MQ=255)
                  aCTGGTGCGGCGAGGATCCCCCCCGTGCGGCTggg                  >  1:16739332/1‑34 (MQ=255)
                  aCTGGTGCGGCGAGGATCCCCCCCGTGCGGCTggg                  >  1:8095031/1‑34 (MQ=255)
                  aCTGGTGCGGCGAGGATCCCCCCCGTGCGGCTggg                  >  1:8228553/1‑34 (MQ=255)
                  aCTGGTGCGGCGAGGATCCCCCCCGTGCGGCTggg                  >  1:9527488/1‑34 (MQ=255)
                  aCTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGt                  >  1:1985857/1‑35 (MQ=255)
                   cTGGTGCGGCGAGGATCCCCCCCGTGCGGCTgggt                 >  1:4324675/1‑33 (MQ=255)
                   cTGGTGCGGCGAGGATCCCCCCCGTGCGGCTgggt                 >  1:1397442/1‑33 (MQ=255)
                   cTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGtt                 <  1:10103735/35‑1 (MQ=255)
                   cTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGtt                 <  1:15520766/35‑1 (MQ=255)
                    tGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGttt                <  1:15908648/35‑1 (MQ=255)
                    tGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGgtt                >  1:4879704/1‑35 (MQ=255)
                     ggTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTc               <  1:13013395/35‑1 (MQ=255)
                     ggTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTc               <  1:4083037/35‑1 (MQ=255)
                     ggTGCGGCGAGGATCCCCCCCGTGCGGCTGGGTTc               >  1:1324979/1‑35 (MQ=255)
                      gTGCGGCGAGGATCCCCCCCGTGCGGCTGGgtttc              >  1:2098377/1‑33 (MQ=38)
                       tGCGGCGAGGATCCCCCCCGTGCGGCTGGgtttgg             >  1:1183903/1‑32 (MQ=38)
                       tGCGGCGAGGATCCCCCCCGTGCGGCTGGgtttgg             >  1:5370680/1‑32 (MQ=38)
                       tGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGc             <  1:5122185/35‑1 (MQ=255)
                       tGCGGCGAGGATCCCCCCCGTGCGGCTGGGTTCGc             >  1:8549284/1‑35 (MQ=255)
                        gcggcgAGGATCCCCCCCGTGCGGCTGGTTTCGCt            <  1:188814/35‑1 (MQ=255)
                        gcggcgAGGATCCCCCCCGTGCGGCTGGTTTCGCt            <  1:9903321/35‑1 (MQ=255)
                         cggcgAGGATCCCCCCCCGTGCGGCTGGGTTCggt            >  1:336843/1‑33 (MQ=25)
                         cggcgAGGATCCCCCCCGTGCGGCTGGTTTCGCTc           <  1:6096241/35‑1 (MQ=255)
                         cggcgAGGATCCCCCCCGTGCGGCTGGTTTCGCTc           <  1:11593708/35‑1 (MQ=255)
                         cggcgAGGATCCCCCCCGTGCGGCTGGTTTCGCTc           <  1:4606845/35‑1 (MQ=255)
                         cggcgAGGATCCCCCCCGTGCGGCTGGTTTCGCTc           <  1:14795552/35‑1 (MQ=255)
                         cggcgAGGATCCCCCCCGTGCGGCTGGGTTTGCTc           >  1:19769286/1‑35 (MQ=37)
                         cggcgAGGATCCCCCCCGTGCGGCTGGGTTCGCTc           >  1:8433079/1‑35 (MQ=255)
                         cggcgAGGATCCCCCCCGTGCGGCTGGGTTCGCTc           >  1:9093764/1‑35 (MQ=255)
                         cggcgAGGATCCCCCCCGTGCGGCTGGGTTCGCTc           >  1:7269014/1‑35 (MQ=255)
                          ggcgAGGATCCCCCCCGTGCGGCTGGTTTCGCTCa          >  1:15488984/1‑35 (MQ=255)
                          ggcgAGGATCCCCCCCGTGCGGCTGGGTTCGGTCa          >  1:7037499/1‑35 (MQ=38)
                          ggcgAGGATCCCCCCCGTGCGGCTGGGTTCGCTCa          >  1:388602/1‑35 (MQ=255)
                           gcgAGGATCCCCCCCGTGCGGCTGGTTTCGCTCaa         <  1:4452564/35‑1 (MQ=255)
                           gcgAGGATCCCCCCCGTGCGGCTGGTTTCGCTCaa         >  1:10742866/1‑35 (MQ=255)
                           gcgAGGATCCCCCCCGTGCGGCTGGGTTCGCTCaa         >  1:5624010/1‑35 (MQ=255)
                            cgAGGATCCCCCCCGTGCGGCTGGTTTCGCTCCAc        >  1:18791314/1‑35 (MQ=255)
                            cgAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAAc        >  1:18300255/1‑35 (MQ=255)
                            cgAGGATCCCCCCCGTGCGGCTGGGTTTGCTCAAc        >  1:13551128/1‑35 (MQ=37)
                            cgAGGATCCCCCCAGTGCGGCTGGTTTCGCTCAAc        <  1:1695423/35‑1 (MQ=255)
                             gAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACg       <  1:11073628/35‑1 (MQ=255)
                             gAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACg       <  1:17520929/35‑1 (MQ=255)
                             gAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACg       <  1:1744264/35‑1 (MQ=255)
                             gAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACg       <  1:17200707/35‑1 (MQ=255)
                             gAGGATCCCCCCCGTGCGGCTGGGTTTGCTCAACg       >  1:15585026/1‑35 (MQ=37)
                              aGGATCCCCCCCGTGCGGCTGGGTTCGCTCAACGc      >  1:14727754/1‑35 (MQ=255)
                              aGGATCCCCCCCGTGCGGCTGGGTTAGCGCAACgg      >  1:12202032/1‑34 (MQ=25)
                                gATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCca    <  1:10358331/35‑1 (MQ=255)
                                gATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCca    <  1:19618225/35‑1 (MQ=255)
                                gATCCCCCCCGTGCGGCTGGTTTCGCTAAACGCca    <  1:18435675/35‑1 (MQ=255)
                                gATCCCCCCCGTGCGGCTGGGTTTGCTCCACGCca    >  1:14899546/1‑35 (MQ=25)
                                 aTCCCCCCCGTGCTGCTGGTTTCGCTCAACGCcaa   >  1:16206181/1‑35 (MQ=255)
                                 aTCCCCCCCGTGCGGCTGGTTTCGCTCAACGCcaa   <  1:13429087/35‑1 (MQ=255)
                                 aTCCCCCCCGTGCGGCTGGTTTCGCTCAACGCcaa   >  1:9711125/1‑35 (MQ=255)
                                  tCCCCCCCGTGCGGCTGGTTTCGCTCAACGCcaac  <  1:13634383/35‑1 (MQ=255)
                                  tCCCCCCCGTGCGGCTGGTTTCGCTCAACGCcaac  <  1:14144363/35‑1 (MQ=255)
                                  tCCCCCCCGTGCGGCTGGTTTCGCTCAACGCcaac  <  1:10898452/35‑1 (MQ=255)
                                  tCCCCCCCGTGCGGCTGGTTTCGCTCAACGCcaac  <  1:10406808/35‑1 (MQ=255)
                                   |                                  
CGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCCAAC  >  REL606/2766924‑2766992

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.