Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,847,2521.→A 1.0% 2.2 ‑1.5 ‑0.0209coding (3371/3543 nt)recBexonuclease V (RecBCD complex), beta subunit
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (95/112):  total (97/112)
Fisher's exact test strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.39e-01

Only 100 of 237 total aligned reads displayed.

AAAGTGACGATCATAGTCGTAATCAGCAATGCGATGGCGCAGATAACGATGCAGCGCCAGGGTATAAAG  >  REL606/2847218‑2847286
                                   |                                  
aaaGTGACGATCATAGTCGTAATCAGCAATTCGAt                                     >  1:10003850/1‑35 (MQ=255)
aaaGTGACGATCATAGTCGTAATCAGCAATGCGAt                                     >  1:667230/1‑35 (MQ=255)
aaaGTGACGATCATAGTCGTAATCAGCAATGCGAt                                     <  1:6039144/35‑1 (MQ=255)
aaaGTGACGATCATAGTCGTAATCAGCAATGCGAt                                     >  1:19909370/1‑35 (MQ=255)
aaaGTGACGATCATAGTCGTAATCAGCAATGCGAt                                     <  1:19732975/35‑1 (MQ=255)
 aaGTGACGATCATAGTCGTAATCAGCAATGCGATg                                   <  1:2242049/35‑1 (MQ=255)
  aGTGACGATCATAGTCGTAATCAGCAATGCGATgg                                  <  1:10504658/35‑1 (MQ=255)
  aGTGACGATCATAGTCGTAATCAGCAATGCGATgg                                  >  1:17827915/1‑35 (MQ=255)
   gTGACGATCATAGTCGTAATCAGCAATGCGATGgc                                 >  1:15481398/1‑35 (MQ=255)
   gTGACGATCATAGTCGTAATCAGCAATGCGATGgc                                 <  1:13038652/35‑1 (MQ=255)
    tGACGATCATAGTCGTAATCAGCAATGCGATGgcg                                >  1:9719378/1‑35 (MQ=255)
    tGACGATCATAGTCGTAATCAGCAATGCGATGgcg                                >  1:4188421/1‑35 (MQ=255)
    tGACGATCATAGTCGTAATCAGCAATGCGATGgcg                                >  1:2699291/1‑35 (MQ=255)
    tGACGATCATAGTCGTAATCAGCAATGCGATGgcg                                <  1:20410707/35‑1 (MQ=255)
    tGACGATCATAGTCGTAATCAGCAATGCGATGgcg                                >  1:14796892/1‑35 (MQ=255)
     gACGATCATAGTCGTAATCAGCAATGCGATGgcgc                               >  1:9459188/1‑35 (MQ=255)
     gACGATCATAGTCGTAATCAGCAATGCGATGgcgc                               <  1:8407569/35‑1 (MQ=255)
     gACGATCATAGTCGTAATCAGCAATGCGATGgcgc                               <  1:686298/35‑1 (MQ=255)
     gACGATCATAGTCGTAATCAGCAATGCGATGgcgc                               >  1:6357415/1‑35 (MQ=255)
      actaTCATAGTCGTAATCAGCAATGCGATGGCGCa                              >  1:8514323/4‑35 (MQ=255)
      aCGATCATAGTCGTAATCAGCAATGCGATGGCGCa                              <  1:166202/35‑1 (MQ=255)
      aCGATCATAGTCGTAATCAGCAATGCGATGGCGCa                              >  1:15191230/1‑35 (MQ=255)
       cGATCATAGTCGTAATCAGCAATGCGATGGCGCAg                             <  1:17462844/35‑1 (MQ=255)
       cGATCATAGTCGTAATCAGCAATGCGATGGCGCAg                             >  1:4379689/1‑35 (MQ=255)
       cGATCATAGTCGTAATCAGCAATGCGATGGCGCAg                             >  1:12147678/1‑35 (MQ=255)
       cGATCATAGTCGTAATCAGCAATGCGATGGCGCAg                             >  1:14723747/1‑35 (MQ=255)
        gATCATAGTCGTAATCAGCAATGCGATGGCGCAGa                            <  1:16440451/35‑1 (MQ=255)
        gATCATAGTCGTAATCAGCAATGCGATGGCGCAGa                            >  1:10686192/1‑35 (MQ=255)
         aTCATAGTCGTAATCAGCAATGCGATGGCGCAGAt                           <  1:10212342/35‑1 (MQ=255)
          tCATAGTCGTAATCAGCAATGCGATGGCGCAGATa                          >  1:3105104/1‑35 (MQ=255)
          tCATAGTCGTAATCAGCAATGCGATGGCGCAGATa                          <  1:1627146/35‑1 (MQ=255)
           cATATTCGTAATCAGCAATGCGATGGCGCAGATaa                         >  1:15726712/1‑35 (MQ=255)
           cATAGTCGTAATCAGCAATGCGATGGCGCAGATaa                         <  1:14496646/35‑1 (MQ=255)
           cATAGTCGTAATCAGCAATGCGATGGCGCAGATaa                         <  1:20332772/35‑1 (MQ=255)
           cATAGTCGTAATCAGCAATGCGATGGCGCAGATaa                         <  1:11743291/35‑1 (MQ=255)
            aTAGTCGTAATCAGCAATGCGATGGCGCAGATAAn                        <  1:18340812/35‑1 (MQ=255)
            aTAGTCGTAATCAGCAATGCGATGGCGCAGATAAc                        >  1:5771924/1‑35 (MQ=255)
            aTAGTCGTAATCAGCAATGCGATGGCGCAGATAAc                        <  1:14975856/35‑1 (MQ=255)
             tAGTCGTAATCAGCAATGCGATGGCGCAGATAACg                       <  1:7358116/35‑1 (MQ=255)
             tAGTCGTAATCAGCAATGCGATGGCGCAGATAACg                       >  1:13585393/1‑35 (MQ=255)
             tAGTCGTAATCAGCAATGCGATGGCGCAGATAACg                       <  1:2500755/35‑1 (MQ=255)
             tAGTCGTAATCAGCAATGCGATGGCGCAGATAACg                       <  1:13161667/35‑1 (MQ=255)
              gcacGTAATCAGCAAGCGATATACGCAGATAACGa                      >  1:19741412/4‑35 (MQ=1)
              aGTCGTAATCAGCAATGCGATGGCGCAGATAACGa                      <  1:4882409/35‑1 (MQ=255)
              aGTCGTAATCAGCAATGCGATGGCGCAGATAACGa                      <  1:7715741/35‑1 (MQ=255)
              aGTCGTAATCAGCAATGCGATGGCGCAGATAACGa                      >  1:12648565/1‑35 (MQ=255)
               gTCGTAATCAGCAATGCGATGGCGCAGATAACGAt                     >  1:10100531/1‑35 (MQ=255)
               gTCGTAATCAGCAATGCGATGGCGCAGATAACGAt                     <  1:19017272/35‑1 (MQ=255)
               gTCGTAATCAGCAATGCGATGGCGCAGATAACGAt                     <  1:654420/35‑1 (MQ=255)
               cacGTAATCAGCAAGCGATATACGCAGATAACGAt                     >  1:11665264/3‑35 (MQ=1)
                tCGTAATCAGCAATGCGATGGCGCAGATAACGATg                    <  1:18612895/35‑1 (MQ=255)
                tCGTAATCAGCAATGCGATGGCGCAGATAACGATg                    <  1:8657266/35‑1 (MQ=255)
                tCGTAATCAGCAATGCGATGGCGCAGATAACGATg                    <  1:11913821/35‑1 (MQ=255)
                 cGTAATCAGCAATGCGATGGCGCAGATAACGATGc                   >  1:921250/1‑35 (MQ=255)
                 cGTAATCAGCAATGCGATGGCGCAGATAACGATGc                   >  1:18843682/1‑35 (MQ=255)
                 cGTAATCAGCAATGCGATGGCGCAGATAACGATGc                   >  1:8885685/1‑35 (MQ=255)
                 cGTAATCAGCAATGCGATGGCGCAGATAACGATGc                   <  1:12076036/35‑1 (MQ=255)
                  gTAATCAGCAATGCGATGGCGCAGATAACGATGCa                  <  1:16845259/35‑1 (MQ=255)
                  gTAATCAGCAATGCGATGGCGCAGATAACGATGCa                  <  1:3384676/35‑1 (MQ=255)
                  gTAATCAGCAATGCGATGGCGCAGATAACGATGCa                  >  1:3853917/1‑35 (MQ=255)
                  gTAATCAGCAATGCGATGGCGCAGATAACGATGCa                  <  1:16927060/35‑1 (MQ=255)
                  gTAATCAGCAATGCGATGGCGCAGATAACGATGCa                  <  1:4937063/35‑1 (MQ=255)
                  gTAATCAGCAATGCGATGGCGCAGATAACGATGCa                  <  1:12259927/35‑1 (MQ=255)
                  gTAATCAGCAATGCGACGGCGCAGATAACGATGCa                  >  1:6986239/1‑35 (MQ=255)
                   tAATCAGCAATGCGATGGCGCAGATAACGATGCAg                 >  1:10844825/1‑35 (MQ=255)
                   tAATCAGCAATGCGATGGCGCAGATAACGATGCAg                 <  1:16062913/35‑1 (MQ=255)
                   tAATCAGCAATGCGATGGCGCAGATAACGATGCAg                 <  1:18269266/35‑1 (MQ=255)
                    aaTCAGCAATGCGATGGCGCAGATAACGATGCacc                <  1:7119669/35‑3 (MQ=255)
                    aaTCAGCAATGCGATGGCGCAGATAACGATGCAgc                <  1:20531193/35‑1 (MQ=255)
                    aaTCAGCAATGCGATGGCGCAGATAACGATGCAgc                <  1:7788743/35‑1 (MQ=255)
                     aTCAGCAATGCGATGGCGCAGATAACGATGCAgcg               >  1:4023235/1‑35 (MQ=255)
                     aTCAGCAATGCGATGGCGCAGATAACGATGCAgcg               >  1:12625513/1‑35 (MQ=255)
                      tCAGCAATGCGATGGCGCAGATAACGATGCAgcgc              >  1:2155801/1‑35 (MQ=255)
                      tCAGCAATGCGATGGCGCAGATAACGATGCAgcgc              <  1:12885743/35‑1 (MQ=255)
                       cAGCAATGCGATGGCGCAGATAACGATGCAGCGcc             >  1:1911670/1‑35 (MQ=255)
                       cAGCAATGCGATGGCGCAGATAACGATGCAGCGcc             <  1:8791289/35‑1 (MQ=255)
                       cAGCAATGCGATGGCGCAGATAACGATGCAGCGcc             >  1:15115457/1‑35 (MQ=255)
                       cAGCAATGCGATGGCGCAGATAACGATGCAGCGcc             <  1:14370769/35‑1 (MQ=255)
                        aGCAATGCGATGGCGCAGATAACGATGCAGCGCCa            >  1:13997909/1‑35 (MQ=255)
                        aGCAATGCGATGGAGCAGATAACGATGCAGCGcct            <  1:5662884/35‑2 (MQ=255)
                         gCAATGCGATGGCGCAGATAACGATGCAGCGCCAg           <  1:8981897/35‑1 (MQ=255)
                         gCAATGCGATGGCGCAGATAACGATGCAGCGCCAg           <  1:14135305/35‑1 (MQ=255)
                         gCAATGCGATGGCGCAGATAACGACGCAGCGCCAg           <  1:9587978/35‑1 (MQ=255)
                         gCAATGAGATGGCGCAGATAACGATGCAGCGCCAg           <  1:17175493/35‑1 (MQ=255)
                          cAATGCGATGGCGCAGATAACGATGCAGCGCCAgg          <  1:8056761/35‑1 (MQ=255)
                          cAATGCGATGGCGCAGATAACGATGCAGCGCCAgg          >  1:18708/1‑35 (MQ=255)
                          cAATGCGATGGCGCAGATAACGATGCAGCGCCAgg          <  1:17708488/35‑1 (MQ=255)
                          cAATGCGATGGCGCAGATAACGATGCAGCGCCAgg          <  1:20208967/35‑1 (MQ=255)
                          cAATGCGATGGCGCAGATAACGATGCAGCGCCAgg          <  1:1623596/35‑1 (MQ=255)
                           aaTGCGATGGCGCAGATAACGATGCAGCGCCAggg         >  1:14064001/1‑35 (MQ=255)
                           aaTGCGATGGCGCAGATAACGATGCAGCGCCAggg         <  1:11507775/35‑1 (MQ=255)
                             tGCGATGGCGCAGATAACGATGCAGCGCCAGGGta       >  1:11110575/1‑35 (MQ=255)
                             ggCGATGGCGCAGATAACGATGCAGCGCCAGGGta       <  1:19563508/34‑1 (MQ=255)
                              gCGATGGCGCAGATAACGATGCAGCGCCAGGGtat      <  1:1568964/35‑1 (MQ=255)
                               cGATGGCGCAGATAACGATGCAGCGCCAGGGcata     >  1:3927051/1‑35 (MQ=255)
                                 aTGGCGCAGATAACGATGCAGCGCCAGGGTATaaa   <  1:4513023/35‑1 (MQ=255)
                                 aTGGCGCAGATAACGATGCAGCGCCAGGGTATaaa   >  1:17953210/1‑35 (MQ=255)
                                  tGGCGCAGATAACGATGCAGCGCCAGTGTATAAAg  >  1:4903539/1‑35 (MQ=255)
                                  tGGCGCAGATAACGATGCAGCGCCAGGGTATAAAg  >  1:7059055/1‑35 (MQ=255)
                                  tGGCGCAGATAACGATGCAGCGCCAGGGTATAAAg  >  1:2865473/1‑35 (MQ=255)
                                   |                                  
AAAGTGACGATCATAGTCGTAATCAGCAATGCGATGGCGCAGATAACGATGCAGCGCCAGGGTATAAAG  >  REL606/2847218‑2847286

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