Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,721,8310C→T 8.1% 25.8 ‑42.0 ‑0.0233S194F (TCT→TTT) srlEglucitol/sorbitol‑specific enzyme IIB component of PTS
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 (19/0):  ref base (15/198):  total (34/199)
Fisher's exact test strand distribution p-value = 5.47e-19
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.97e-01

Only 100 of 248 total aligned reads displayed.

TTCGTCTCGGCGCTCATTGGCATCATTATGGCTTCTGGCCTTGGTGACTGGATTGCCCACGGCCTTGCT  >  REL606/2721797‑2721865
                                  |                                  
 tCGTCTCGGCGCTCATTGGCATCATTATGGCTTCt                                   <  1:16522782/35‑1 (MQ=255)
 tCGTCTCGGCGCTCATTGGCATCATTATGGCTTCt                                   <  1:1585676/35‑1 (MQ=255)
 tCGTCTCGGCGCTCATTGGCATCATTATGGCTTCt                                   <  1:4155697/35‑1 (MQ=255)
 tCGTCTCGGCGCTCATTGGCATCATTATGGCTTCt                                   <  1:5776457/35‑1 (MQ=255)
  cGTCTCGGCGCTCATTGGCATCATTATGGCTTCTg                                  <  1:17697498/35‑1 (MQ=255)
  cGTCTCGGCGCTCATTGGCATCATTATGGCTTCTg                                  <  1:3617979/35‑1 (MQ=255)
  cGTCTCGGCGCTCATTGGCATCATTATGGCTTCTg                                  <  1:4740266/35‑1 (MQ=255)
   gTCTCGGCGCTCATTGGCATCATTATGGCTTCTgg                                 <  1:16940916/35‑1 (MQ=255)
   gTCTCGGCGCTCATTGGCATCATTATGGCTTCTgg                                 <  1:7999051/35‑1 (MQ=255)
   gTCTCGGCGCTCATTGGCATCATTATGGCTTCTgg                                 <  1:9730548/35‑1 (MQ=255)
    tctcGGCGCTCATTGGCATCATTATGGCTTCTGGc                                <  1:850770/35‑1 (MQ=255)
    tctcGGCGCTCATTGGCATCATTATGGCTTCTGGc                                <  1:10253179/35‑1 (MQ=255)
     ctcGGCGCTCATTGGCATCATTATGGCTTCTGGcc                               <  1:15684439/35‑1 (MQ=255)
     ctcGGCGCTCATTGGCATCATTATGGCTTCTGGcc                               <  1:20373034/35‑1 (MQ=255)
     ctcGGCGCTCATTAGCATCATTATGGCTTTTggtc                               >  1:7290535/1‑33 (MQ=25)
      tcGGCGCTCATTGGCATCATTATGGCTTTTTGCCt                              >  1:11097531/1‑35 (MQ=25)
      tcGGCGCTCATTGGCATCATTATGGCTTCTGGCCt                              <  1:16115750/35‑1 (MQ=255)
      tcGGCGCTCATTGGCATCATTATGGCTTCTGGCCt                              <  1:8426705/35‑1 (MQ=255)
      tcGGCGCTCATTGGCATCATTATGGCTTCTGGCCt                              <  1:952214/35‑1 (MQ=255)
      tcGGCGCTCATTGGCATCATTATGGCTTCTGGCCt                              <  1:19787271/35‑1 (MQ=255)
      tcGGAGCTCATTGGCATCATTATGGCTTCTGGCCt                              <  1:3880046/35‑1 (MQ=255)
       cGGCGCTCATTGGCATCATTATGGCTTCTGGCCtt                             <  1:10379926/35‑1 (MQ=255)
       cGGCGCTCATTGGCATCATTATGGCTTCTGGCCtt                             <  1:6880924/35‑1 (MQ=255)
       cGGCGCTCATTGGCATCATTATGGCTTCTGGCCtt                             <  1:7617623/35‑1 (MQ=255)
       cGGCGCTCATTGGCATCATTATGGCTTCTGGCCtt                             <  1:19249570/35‑1 (MQ=255)
        ggCGCTCATTGGCATCATTATGGCTTCTGGCCTTg                            <  1:6058828/35‑1 (MQ=255)
        ggCGCTCATTGGCATCATTATGGCTTCTGGCCTTg                            <  1:11973403/35‑1 (MQ=255)
        ggCGCCCATTGGCATCATTATGGCTTCTGGCCTTg                            <  1:2127065/35‑1 (MQ=255)
        ggCACTCATTGGCATCATTATGGCTTCTGGCCTTg                            <  1:1650613/35‑1 (MQ=255)
         gcgcTCATTGGCATCATTATGGCTTTTTTCCTTgt                           >  1:8264990/1‑34 (MQ=21)
         gcgcTCATTGGCATCATTATGGCTTCTGGCCTTgg                           <  1:7511504/35‑1 (MQ=255)
          cgcTCATTTGCATCATTATGACTTCTGGCCTTGGt                          >  1:18697506/1‑35 (MQ=37)
          cgcTCATTGGCATCATTATGGCTTTTTGTCTTGGt                          >  1:9156774/1‑35 (MQ=25)
          cgcTCATTGGCATCATTATGGCTTTTGGCTTTTGt                          >  1:17498398/1‑35 (MQ=25)
          cgcTCATTGGCATCATTATGGCTTCTGGCCTTGGt                          <  1:9238110/35‑1 (MQ=255)
          cgcTCATTGGCATCATTATGGCTTCTGGCCTTGGt                          <  1:11265568/35‑1 (MQ=255)
          cgcTCATTGGCATCATTATGGCTTCTGGCCTTGGt                          <  1:18895271/35‑1 (MQ=255)
           gcTCATTGGCATCATTATGGCTTTTTTGCTTGGTg                         >  1:5365318/1‑35 (MQ=16)
           gcTCATTGGCATCATTATGGCTTTTTGTCTTGGTg                         >  1:16794915/1‑35 (MQ=25)
           gcTCATTGGCATCATTATGGCTTCTGGCCTTGGTg                         <  1:7018638/35‑1 (MQ=255)
           gcTCATTGGCATCATTATGGCTTCTGGCATTGGTg                         <  1:11015207/35‑1 (MQ=255)
            cTCATTGGCATCATTATGGCTTTTTTGCTTGGTgg                        >  1:8984643/1‑34 (MQ=12)
            cTCATTGGCATCATTATGGCTTTTTTCCTTTGTgg                        >  1:18418136/1‑34 (MQ=12)
            cTCATTGGCATCATTATGGCTTTTTGTCTTGGTgg                        >  1:6218229/1‑34 (MQ=21)
            cTCATTGGCATCATTATGGCTTTTTGTCTTGGTgg                        >  1:19639015/1‑34 (MQ=21)
            cTCATTGGCATCATTATGGCTTCTGGCCTTGGTGa                        <  1:10145221/35‑1 (MQ=255)
             tCATTGGCATCATTATGGCTTTTTTCCTTTGTggc                       >  1:10023283/1‑33 (MQ=12)
             tCATTGGCATCATTATGGCTTCTGGCCTTGGTGAc                       <  1:7722391/35‑1 (MQ=255)
             tCATTGGCATCATTATGGCTTCTGGCCTTGGTGAc                       <  1:5864297/35‑1 (MQ=255)
              cATTGGCATCATTATGGCTTCTGGCCTTGGTGACt                      <  1:13852819/35‑1 (MQ=255)
              cATTGGCATCATTATGGCTTCTGGCCTTGGTGACt                      <  1:529766/35‑1 (MQ=255)
               aTTGGCATCATTATGGCTTCTGGCCTTGGTGACTg                     <  1:10654723/35‑1 (MQ=255)
               aTTGGCATCATTATGGCTTCTGGCCTTGGTGACTg                     <  1:6660529/35‑1 (MQ=255)
               aTTGGCATCATTATGGCTTCTGGCCTTGGTGACTg                     <  1:12908146/35‑1 (MQ=255)
               aTTGGCATCATTATGGCTTCTGGCCTTGGTGACTg                     <  1:7675790/35‑1 (MQ=255)
                ttGGCATCATTATGGCTTGTGGCCTTGGTGACTgg                    <  1:1337674/35‑1 (MQ=255)
                ttGGCATCATTATGGCTTCTGGCCTTGGTGACTgg                    <  1:6181735/35‑1 (MQ=255)
                ttGGCATCATTATGGCTTCTGGCCTTGGTGACTgg                    <  1:18870857/35‑1 (MQ=255)
                 tGGCATCATTATGGCTTCTGGCCTTGGTGACTGGa                   <  1:8856489/35‑1 (MQ=255)
                 tGGCATCATTATGGCTTCTGGCCTTGGTGACTGGa                   <  1:19602222/35‑1 (MQ=255)
                 tGGCATCATTATGGCTTCTGGCCTTGGTGACTGGa                   <  1:12421375/35‑1 (MQ=255)
                   gCATCATTATGGCTTCTGGCCTTGGTGACTGGAtt                 <  1:860316/35‑1 (MQ=255)
                   gCATCATTATGGCTTCTGGCCTTGGTGACTGGAtt                 <  1:17280627/35‑1 (MQ=255)
                    catcatTATGGCTTCTGGCCTTGGTGACTGGATTg                <  1:20195101/35‑1 (MQ=255)
                    catcatTATGGCTTCTGGCCTTGGTGACTGGATTg                <  1:14435560/35‑1 (MQ=255)
                     atcatTATGGCTTCTGGCCTTGGTGACTGGATTGc               <  1:11211089/35‑1 (MQ=255)
                      tcatTATGGCTTCTGGCCTTGGTGACTGGATTGcc              <  1:3244670/35‑1 (MQ=255)
                      tcatTATGGCTTCTGGCCTTGGTGACTGGATTGcc              <  1:13476664/35‑1 (MQ=255)
                       catTATGGCTTCTGGCCTTGGTGACTGGATTGccc             <  1:12361639/35‑1 (MQ=255)
                       catTATGGCTTCTGGCCTTGGTGACTGGATTGccc             <  1:2279130/35‑1 (MQ=255)
                       catTATGGCTTCTGGCCTTGGTGACTGGATTGccc             <  1:14815886/35‑1 (MQ=255)
                         ttATGGCTTCTGGCCTTGGTGACTGGATTGCCCAc           <  1:18482040/35‑1 (MQ=255)
                         ttATGGCTTCTGGCCTTGGTGACTGGATTGCCCAc           <  1:20538499/35‑1 (MQ=255)
                         ttATGGCTTCTGGCCTTGGTGACTGGATTGACCAc           <  1:3726412/35‑1 (MQ=255)
                          tATGGCTTCTGGCCTTGGTGACTGGATTGCCCACg          <  1:12493144/35‑1 (MQ=255)
                          tATGGCTTCTGGCCTTGGTGACTGGATTGCCCACg          <  1:17538316/35‑1 (MQ=255)
                          tATGGCTTCTGGCCTTGGTGACTGGATTGCCCACg          <  1:7833842/35‑1 (MQ=255)
                           aTGGCTTCTGGCCTTGGTGACTGGATTGCCCACgg         <  1:17859702/35‑1 (MQ=255)
                           aTGGCTTCTGGCCTTGGTGACTGGATTGCCCACgg         <  1:2478219/35‑1 (MQ=255)
                            tGGCTTCTGGCCTTGGTGACTGGATTGCCCACGGc        <  1:18074797/35‑1 (MQ=255)
                            tGGCTTCTGGCCTTGGTGACTGGATTGCCCACGGc        <  1:11427101/35‑1 (MQ=255)
                            tGGCTTCTGGCCTTGGTGACTGGATTGCCCACGGc        <  1:3048668/35‑1 (MQ=255)
                              gCTTCTGGCCTTGGTGACTGGATTGCCCACGGCCt      >  1:17800035/1‑35 (MQ=255)
                              gCTTCTGGCCTTGGTGACTGGATTGCCCACGGCCt      <  1:11715720/35‑1 (MQ=255)
                              gCTTCTGGCCTTGGTGACTGGATTGCCCACGGCCt      <  1:4459313/35‑1 (MQ=255)
                              gCTTCTGGCCTTGGTGACTGGATTGCCCACGGCCt      <  1:5876735/35‑1 (MQ=255)
                              gCTTCTGGCCTTGGTGACTGGATTGCCCACGGCCt      >  1:13983359/1‑35 (MQ=255)
                              gCTTCTGGCCTTGGTGACTGGATTGCCAACGGCCt      <  1:11932991/35‑1 (MQ=255)
                               cTTCTGGCCTTGGTGACTGGATTGCCCACGGCCtt     <  1:3487921/35‑1 (MQ=255)
                                ttCTGGCCTTGGTGACTGGATTGCCCACGGCCTTg    >  1:16200483/1‑35 (MQ=255)
                                 tCTGGCCTTGGTGACTGGATTGCCCACGGCCTTGc   <  1:15114823/35‑1 (MQ=255)
                                 tCTGGCCTTGGTGACTGGATTGCCCACGGCCTTGc   >  1:1544898/1‑35 (MQ=255)
                                 tCTGGCCTTGGTGACTGGATTGCCCACGGCCTTGc   <  1:5646340/35‑1 (MQ=255)
                                 tCTGGCCTTGGTGACTGGATTGCCCACGGCATTGc   <  1:11581785/35‑1 (MQ=255)
                                  cTGGCCTTTGTGACTGGATTGCCCACGGCCTTGCt  >  1:1421663/1‑35 (MQ=255)
                                  cTGGCCTTGGTGACTGGATTGCCCACGGCCTTGCt  <  1:7591368/35‑1 (MQ=255)
                                  cTGGCCTTGGTGACTGGATTGCCCACGGCCTTGCt  >  1:12070199/1‑35 (MQ=255)
                                  cTGGCCTTGGTGACTGGATTGCCCACGGCCTTGCt  <  1:734532/35‑1 (MQ=255)
                                  cTGGCCTTGGTGACTGGATTGCCCACGGCCTTGCt  <  1:17341292/35‑1 (MQ=255)
                                  cTGGCCTTGGTGACGGGATTGCCCACGGCCTTGCt  <  1:2817052/35‑1 (MQ=255)
                                  |                                  
TTCGTCTCGGCGCTCATTGGCATCATTATGGCTTCTGGCCTTGGTGACTGGATTGCCCACGGCCTTGCT  >  REL606/2721797‑2721865

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.