Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,893,550 Δ1 bp 29.5% intergenic (+5/‑51) kup → / → insJ‑5 potassium transporter/IS150 hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,893,5500A→.29.5% 145.4 ‑3.0 0.0129intergenic (+5/‑51)kup/insJ‑5potassium transporter/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  new base (26/12):  ref base (44/47):  total (70/59)
Fisher's exact test strand distribution p-value = 5.20e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

ATCGAACTGGGTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTTGGACGATGAAA  >  REL606/3893516‑3893584
                                    |                                    
aTCGAACTGGGTACTCAGGTCGAAATCTAACGCC‑‑a                                      <  1:3524980/35‑1 (MQ=255)
  cGAACTGGGTACTCAGGTCGAAATCTAACGCC‑‑A‑‑Ga                                  >  1:8567957/1‑35 (MQ=255)
  cGAACTGGGTACTCAGGTCGAAATCTAACGCC‑‑A‑‑Ga                                  >  1:13199/1‑35 (MQ=255)
   gAACTGGGTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GAc                                 >  1:8466115/1‑35 (MQ=255)
   gAACTGGGTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GAc                                 <  1:15125974/35‑1 (MQ=255)
   gAACTGGGTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GAc                                 <  1:13000793/35‑1 (MQ=255)
    aaCTGGGTACTCAGGTCGAAATCTAACGCC‑‑A‑‑ggac                                <  1:15948582/35‑4 (MQ=255)
    aaCTGGGTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACt                                <  1:2567460/35‑1 (MQ=255)
    aaCTGGGTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACt                                <  1:12547314/35‑1 (MQ=255)
    aaCTGGGTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACt                                >  1:13942699/1‑35 (MQ=255)
     aCTGGGTACTCAGTTCGAAATCTAACGCC‑‑A‑‑GACTg                               >  1:5880668/1‑35 (MQ=255)
     aCTGGGTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTg                               <  1:18522639/35‑1 (MQ=255)
     aCTGGGTACTCAGGTCGAAATCTAACGCC‑‑‑‑GACTGt                              >  1:3397272/1‑35 (MQ=37)
     aCTGGGTACTCAGGTCGAAATCTAACGCC‑‑‑‑GACTGt                              >  1:6001015/1‑35 (MQ=37)
     aCTGGGTACTCAGGTCGAAATCTAACGCC‑‑‑‑GACTGt                              >  1:10403113/1‑35 (MQ=37)
      cTGGGTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGt                              <  1:8037827/35‑1 (MQ=255)
      cTGGGTACTCAGGTCGAAATCTAACGCC‑‑‑‑GACTGTa                             >  1:387057/1‑35 (MQ=37)
      cTGGGTACTCAGGTCGAAATCTAACGCC‑‑‑‑GACTGTa                             >  1:43534/1‑35 (MQ=37)
      cTGGGTACTCAGGTCGAAATCTAACGCC‑‑‑‑GACTGTa                             >  1:8872603/1‑35 (MQ=37)
      cTGGGTACTCAGGTCGAAATCTAACGCC‑‑‑‑GACTGTa                             <  1:16408655/35‑1 (MQ=37)
       tGGGTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTa                             <  1:2436759/35‑1 (MQ=255)
       tGGGTACTCAGGTCGAAATCTAACGCC‑‑‑‑GACTGTAc                            >  1:12404993/1‑35 (MQ=38)
       tGGGTACTCAGGTCGAAATCTAACGCC‑‑‑‑GACTGTAc                            >  1:407990/1‑35 (MQ=38)
       tGGGTACTCAGGTCGAAATCTAACGCC‑‑‑‑GACTGTAc                            <  1:1877561/35‑1 (MQ=38)
       tGGGTACTCAGGTCGAAATCTAACGCC‑‑‑‑GACTGTAc                            >  1:7838440/1‑35 (MQ=38)
       tGGGTACTCAGGTCGAAATCTAACGCC‑‑‑‑GACTGTAc                            >  1:7079601/1‑35 (MQ=38)
        gggTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTAc                            <  1:12590829/35‑1 (MQ=255)
        gggTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTAc                            <  1:3814987/35‑1 (MQ=255)
         ggTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACt                           <  1:14041386/35‑1 (MQ=255)
         ggTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACt                           <  1:4809957/35‑1 (MQ=255)
         ggTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACt                           <  1:17311730/35‑1 (MQ=255)
         ggTACTCAGGTCGAAATCTAACGCC‑‑‑‑GACTGTACTg                          >  1:11720369/1‑35 (MQ=38)
          gTACTCATGTCGAAATCTAACGCC‑‑‑‑GACTGTACTGc                         >  1:10712625/1‑35 (MQ=37)
          gTACTCAGGTCGAAATCTAACGCC‑‑AGGACTGTACt                           >  1:1687942/1‑35 (MQ=37)
          gTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACTg                          <  1:7954058/35‑1 (MQ=255)
          gTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACTg                          <  1:18982549/35‑1 (MQ=255)
          gTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACTg                          <  1:4901076/35‑1 (MQ=255)
           tACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACTGc                         >  1:13813666/1‑35 (MQ=255)
           tACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACTGc                         >  1:1105679/1‑35 (MQ=255)
            aCTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACTGCa                        >  1:6268134/1‑35 (MQ=255)
            aCTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACTGCa                        <  1:17349458/35‑1 (MQ=255)
            aCTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACTGCa                        <  1:8589558/35‑1 (MQ=255)
            aCTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACTGCa                        >  1:15360774/1‑35 (MQ=255)
            aCTCAGGTCGAAATCTAACGCC‑‑‑‑GACTGTACTGCAc                       <  1:2493460/35‑1 (MQ=38)
             cTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACTGCAc                       >  1:12056090/1‑35 (MQ=255)
             cTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACTGCAc                       <  1:13449550/35‑1 (MQ=255)
             cTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACTGCAc                       <  1:10172757/35‑1 (MQ=255)
             cTCAGGTCGAAATCTAACGCC‑‑‑‑GACTGTACTGCAcc                      <  1:5129361/35‑1 (MQ=38)
              tCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACTGCAcc                      <  1:12221382/35‑1 (MQ=255)
              tCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACTGCAcc                      <  1:13422155/35‑1 (MQ=255)
              tCAGGTCGAAATCTAACGCC‑‑‑‑GACTGTACTGCAccc                     >  1:9276157/1‑35 (MQ=38)
              tCAGGTCGAAATCTAACGCC‑‑‑‑GACTGTACTGCAccc                     <  1:20604119/35‑1 (MQ=38)
                aGGTCGAAATCTAACGCC‑‑AGGACTGTACTGCActc                     >  1:5969988/1‑33 (MQ=25)
                 ggTCGAAATCTAACGCC‑‑AGGGACTGTACTGCAccc                     <  1:2156930/35‑1 (MQ=25)
                 ggTCGAAATCTAACGCC‑‑A‑‑GACTGTACTGCACCCAt                   >  1:4080697/1‑35 (MQ=255)
                 ggTCGAAATCTAACGCC‑‑‑‑GACTGTACTGAACCCAtt                  <  1:2324917/35‑1 (MQ=37)
                  gTCGAAATCTAACGCC‑‑A‑‑GACTGTACTGCACCCAtt                  <  1:3058667/35‑1 (MQ=255)
                  gTCGAAATCTAACGCC‑‑A‑‑GACTGTACTGCACCCAtt                  >  1:11448860/1‑35 (MQ=255)
                  gTCGAAATCTAACGCC‑‑A‑‑GACTGTACTGCACCCAtt                  >  1:18337572/1‑35 (MQ=255)
                   tCGAAATCTAACGCC‑‑AGGACTGTACTGCACCCAtt                  <  1:11300379/35‑1 (MQ=25)
                   tCGAAATCTAACGCC‑‑‑‑GACTGTACTGCACCCAtttt                >  1:7868169/1‑35 (MQ=25)
                   tCGAAATCTAACGCC‑‑‑‑GACTGTACTGCACCCAtttt                >  1:12577838/1‑35 (MQ=25)
                    cGAAATCTAACGCC‑‑A‑‑GACTGTACTGCACCCAtttt                <  1:12837366/35‑1 (MQ=255)
                    cGAAATCTAACGCC‑‑A‑‑GACTGTACTGCACCCAtttt                >  1:19607082/1‑35 (MQ=255)
                    cGAAATCTAACGCC‑‑A‑‑GACTGTACTGCACCCAtttt                >  1:6247719/1‑35 (MQ=255)
                    cGAAATCTAACGCC‑‑‑‑GACTGTACTGCACCCATTttg               <  1:11696013/35‑1 (MQ=25)
                    cGAAATCTAACGC‑‑A‑‑GACTGTACTGCACCCATTttg               <  1:6078691/35‑1 (MQ=25)
                     gAAATCTAACGCCAGA‑‑GACTGTACTGCACCCAttt                 <  1:8936458/35‑1 (MQ=18)
                     gAAATCTAACGCCAGA‑‑GACTGTACTGCACCCAttt                 <  1:15502501/35‑1 (MQ=18)
                     gAAATCTAACGCC‑‑A‑‑GACTGTACTGCACCCATTttg               <  1:12475066/35‑1 (MQ=255)
                     gAAATCTAACGCC‑‑A‑‑GACTGTACTGCACCCATTttg               >  1:6636512/1‑35 (MQ=255)
                     gAAATCTAACGCC‑‑A‑‑GACTGTACTGCACCCATTttg               >  1:7463737/1‑35 (MQ=255)
                     gAAATCTAACGCC‑‑‑‑GACTGTACTGCACCCATTttgt              <  1:12325836/35‑1 (MQ=21)
                     aaaaTCTAACGCC‑‑‑‑GACTGTACTGCACCCATTttgt              <  1:13333545/34‑1 (MQ=21)
                      aaaTCTAACGCC‑‑A‑‑GACTGTACTGCACCCATTttgt              >  1:3448004/1‑35 (MQ=255)
                       aaTCTAACGCC‑‑A‑‑GACTGTACTGCACCCATTttgtt             <  1:934360/35‑1 (MQ=255)
                       aaTCTAACGCC‑‑A‑‑GACTGTACTGCACCCATTttgtt             >  1:17124981/1‑35 (MQ=255)
                        aTCTAACGCC‑‑‑‑GACTGTACTGCACCCATTTTGTTgg           >  1:18215698/1‑35 (MQ=18)
                         tCTAACGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTTgg           >  1:10230123/1‑35 (MQ=255)
                         tCTAACGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTTgg           >  1:8016303/1‑35 (MQ=255)
                         tCTAACGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTTgg           >  1:241903/1‑35 (MQ=255)
                         tCTAACGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTTgg           >  1:19676540/1‑35 (MQ=255)
                          cTAACGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTTGGa          <  1:9375308/35‑1 (MQ=255)
                          cTAACGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTTGGa          <  1:12577486/35‑1 (MQ=255)
                          cTAACGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTTGGa          <  1:4759200/35‑1 (MQ=255)
                          cTAACGCC‑‑‑‑GACTGTACTGCACCCATTTTGTTGGAc         <  1:5101731/35‑1 (MQ=17)
                          cTAACGCC‑‑‑‑GACTGTACTGCACCCATTTTGTTGGAc         <  1:3646778/35‑1 (MQ=17)
                           tAACGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTGGGAc         <  1:10046849/35‑1 (MQ=255)
                            aaCGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTTGGACg        >  1:11718126/1‑35 (MQ=255)
                            aaCGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTTGGACg        >  1:15434413/1‑35 (MQ=255)
                            aaCGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTTGGACg        >  1:15502888/1‑35 (MQ=255)
                            aaCGCC‑‑‑‑GACTGTACTGCACCCATTTTGTTGGACGa       >  1:7885979/1‑35 (MQ=14)
                             aCGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTTGGACGa       >  1:678807/1‑35 (MQ=255)
                             aCGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTTGGACGa       >  1:7307590/1‑35 (MQ=255)
                              cGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTTGGACGAt      >  1:19050878/1‑35 (MQ=255)
                              cGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTTGGACGAt      <  1:15207157/35‑1 (MQ=255)
                              cGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTTGGACGAt      <  1:7521486/35‑1 (MQ=255)
                                cc‑‑A‑‑GACTGTACTGCACCCATTTTGTTGGACGATGa    >  1:11225130/1‑35 (MQ=255)
                                 c‑‑A‑‑GACTGTACTGCACCCATTTTGTTGGACGATGaa   >  1:3007875/1‑35 (MQ=40)
                                    a‑‑GACTGTACTGCACCCATTTTGTTGGACGATGaaa  >  1:514955/1‑35 (MQ=38)
                                    |                                    
ATCGAACTGGGTACTCAGGTCGAAATCTAACGCC‑‑A‑‑GACTGTACTGCACCCATTTTGTTGGACGATGAAA  >  REL606/3893516‑3893584

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.