Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,691,510 G→A 36.0% intergenic (+45/‑61) ydgR → / → gst putative tripeptide transporter permease/glutathionine S‑transferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,691,5100G→A36.0% 133.5 ‑0.3 ‑0.2150intergenic (+45/‑61)ydgR/gstputative tripeptide transporter permease/glutathionine S‑transferase
Reads supporting (aligned to +/- strand):  new base (22/32):  ref base (43/53):  total (65/85)
Fisher's exact test strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.43e-01

Only 100 of 166 total aligned reads displayed.

AACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCACTACCATACTT  >  REL606/1691476‑1691544
                                  |                                  
acctctTTAAAAAGCCGCTAACTTTTCGTTAGCgg                                    <  1:10618943/33‑1 (MQ=25)
aaCTCTTTTACAAGCCGCTAACTTTTCGTTAGCgg                                    >  1:19562068/1‑35 (MQ=255)
aaCTCTTTTACAAGCCGCTAACTTTTCGTTAGCgg                                    >  1:5820812/1‑35 (MQ=255)
 cccgtttctctggcaataatcgcggttaattCGGc                                   <  1:937810‑M1/4‑1 (MQ=255)
 cccgtttctctggcaataatcgcggttaattCGGc                                   >  1:3806001‑M1/32‑35 (MQ=255)
 aCTCTTTTACAAGCCGCTAACTTTTCGCTAGCGGc                                   >  1:18143526/1‑35 (MQ=255)
 aCTCTTTAAAAAGCCGCTAACTTTTCGTTAGCGGc                                   <  1:13731814/35‑1 (MQ=37)
  ctctTTTACAAGCCGCTAACTTTTCGTTAGCGGCt                                  >  1:11510685/1‑35 (MQ=255)
  ccgtttctctggcaataatcgcggttaattCGGCt                                  >  1:16204962‑M1/31‑35 (MQ=255)
  ccgtttctctggcaataatcgcggttaattCGGCt                                  >  1:15141085‑M1/31‑35 (MQ=255)
  ccgtttctctggcaataatcgcggttaattCGGCt                                  >  1:855929‑M1/31‑35 (MQ=255)
   tctTTTACAAGCCGCTAACTTTTCGTTAGCGGCtt                                 >  1:15516921/1‑35 (MQ=255)
   tctTTAAAAAGCCGCTAACTTTTCGTTAGCGACtt                                 <  1:14493550/35‑1 (MQ=25)
   cgtttctctggcaataatcgcggttaattTGGCtt                                 <  1:13597224‑M1/6‑1 (MQ=255)
   cgtttctctggcaataatcgcggttaattCGGCtt                                 <  1:19339706‑M1/6‑1 (MQ=255)
   cgtttctctggcaataatcgcggttaattCGGCtt                                 <  1:2924524‑M1/6‑1 (MQ=255)
   cgtttctctggcaataatcgcggttaattCGGCtt                                 <  1:10048728‑M1/6‑1 (MQ=255)
   cgtttctctggcaataatcgcggttaattCGGCtt                                 >  1:13019330‑M1/30‑35 (MQ=255)
    gtttctctggcaataatcgcgggtaattCGGCttt                                <  1:15061296‑M1/7‑1 (MQ=255)
    ctTTTACAAGCCGCTAACTTTTCGTTAGCGGCttt                                >  1:9213871/1‑35 (MQ=255)
    ctTTTACAAGCCGCTAACTTTTCGTTAGCGGCttt                                >  1:1287510/1‑35 (MQ=255)
     ttttACAAGCCGCTAACTTTTCGTTAGCGGCtttt                               >  1:5083523/1‑35 (MQ=255)
     ttttACAAGCCGCTAACTTTTCGTTAGCGGCtttt                               >  1:13567555/1‑35 (MQ=255)
     tttctctggcaataatcgcggttaattCGGCtttt                               <  1:4052707‑M1/8‑1 (MQ=255)
      ttctctggcaataatcgcggttaattCGGCttttt                              <  1:6579571‑M1/9‑1 (MQ=255)
       ttACAAGCCGCTAACTTTTCGTTAGCGGCtttttt                             >  1:16295465/1‑35 (MQ=255)
       tctctggcaataatcgcggttaattCGGCtttttt                             <  1:5205646‑M1/10‑1 (MQ=255)
       tctctggcaataatcgcggttaattCGGCtttttt                             <  1:10465239‑M1/10‑1 (MQ=255)
        tACAAGCCGCTAACTTTTCGTTAGCGACttttttt                            >  1:8377365/1‑35 (MQ=255)
        ctctggcaataatcgcggttaattCGGCttttttt                            <  1:2182189‑M1/11‑1 (MQ=255)
        aaaaaGCCGCTAACTTTTCGTTAGCGACttttttt                            <  1:5981617/32‑1 (MQ=37)
        aaCAAGCCGCTAACTTTTCGTTAGCGACttttttt                            <  1:13748388/34‑1 (MQ=255)
         aCAAGCCGCTAACTTTTCGTTAGCGGCtttttttt                           >  1:1561176/1‑35 (MQ=255)
         aCAAGCAGCTAACTTTTCGTTAGCGACtttttttt                           <  1:5983265/35‑1 (MQ=37)
          ctggcaataatcgcggttaattCGgctttttttct                          <  1:464778‑M1/13‑3 (MQ=255)
          ctggcaataatcgcggttaattCGGcttttttttc                          <  1:2201791‑M1/13‑2 (MQ=255)
          cAAGCCGCTAACTTTTCGTTAGCGGCttttttttt                          >  1:16991555/1‑35 (MQ=255)
          cAAGCCGCTAACTTTTCGTTAGCGGCttttttttt                          >  1:14550017/1‑35 (MQ=255)
          cAAGCCGCTAACTTTTCGTTAGCGACttttttttt                          >  1:13503693/1‑35 (MQ=255)
           aaGCCGCTAACTTTTCGTTAGCGGCtttttttttt                         >  1:2430966/1‑34 (MQ=255)
           aaGCCGCTAACTTTTCGTTAGCGGCtttttttttt                         >  1:3236574/1‑34 (MQ=255)
           aaGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTg                         <  1:6738152/35‑1 (MQ=255)
           aaGCCGCTAACTTTTAGTTAGCGGCTTTTTTTTTg                         <  1:10061142/35‑1 (MQ=255)
            ggcaataatcgcggttaattCGgcttttttttttt                        <  1:8636261‑M1/15‑3 (MQ=255)
            aGCCGCTAACTTTTCGTTAGCGACTTTTTTTTTGt                        <  1:20327418/35‑1 (MQ=255)
            aGCCGCAACTTTTCGTTAGCGACTTTTTTTTTGtt                       >  1:5463120/1‑35 (MQ=37)
             gCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTTtt                       >  1:16340801/1‑35 (MQ=255)
             gCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGtt                       <  1:17617428/35‑1 (MQ=255)
             gCCGCTAACTTTTCGTTAGCGACTTTTTTTTTGtt                       >  1:11104609/1‑35 (MQ=255)
             gCCGCTAACTTTTCGTTAGCGACTTTTTTTTTGtt                       <  1:11818772/35‑1 (MQ=255)
              ccGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTc                      <  1:15199448/35‑1 (MQ=255)
              ccGCTAACTTTTCGTTAGCGACTTTTTTTTTGTTc                      >  1:2796091/1‑35 (MQ=255)
              caataatcgcggttaattCGGCTTTTTTTTTTTTc                      <  1:400980‑M1/17‑1 (MQ=255)
               cGCTAACTTTTCGTTAGCGGCTTTTTTTTTTtttt                     >  1:20178019/1‑33 (MQ=38)
               cGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCa                     <  1:15605052/35‑1 (MQ=255)
               cGCTAACTTTTCGTTAGCGACTTTTTTTTTGTTCa                     >  1:8057876/1‑35 (MQ=255)
               cGCTAACTTTTCGTTAGCGACTTTTTTTTTGTTCa                     <  1:19714470/35‑1 (MQ=255)
               aataatcgcggttaattCGGCTTTTTTTTTTTtcc                     <  1:16092123‑M1/18‑2 (MQ=255)
                gCTAACTTTTCGTTAGCGACTTTTTTTTTGTTCAg                    <  1:9639550/35‑1 (MQ=255)
                 cTAACTTTTCGTTAGCGGCTTTTTTTTTTTTCAGc                   >  1:1739134/1‑35 (MQ=255)
                 cTAACTTTTCGTTAGCGGCTTTTTTTTTGTTGAGc                   <  1:17111001/35‑1 (MQ=255)
                 cTAACTTTTCGTTAGCGAC‑‑TTTTTTTGTTCAGCaa                 <  1:3278988/35‑1 (MQ=21)
                  tAACTTTTCGTTAGCGGCTTTTTTTTTTTTTAGCa                  >  1:13185316/1‑35 (MQ=37)
                  tAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCa                  <  1:1029311/35‑1 (MQ=255)
                  tAACTTTTCGTTAGCGACTTTTTTTTTGTTCAGCa                  <  1:18885967/35‑1 (MQ=255)
                  tAACTTTTCGTTAGCGACTTTTTTTTTGTTCAGCa                  <  1:924871/35‑1 (MQ=255)
                   atcTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCaa                 <  1:11565504/33‑1 (MQ=255)
                   aaCTTTTCGTTAGCTGCTTTTTTTTTGTTCAGCaa                 >  1:7011866/1‑35 (MQ=255)
                   aaCTTTTCGTTAGCGACTTTTTTTTTGTTCATCaa                 >  1:5562415/1‑35 (MQ=38)
                    aCTTTTCGTTAGCGGCTTTTTTTTTTTTAGCaaac               >  1:19142394/1‑33 (MQ=21)
                    aCTTTTCGTTAGCGACTTTTTTTTTGTTCAGCAAc                <  1:13383075/35‑1 (MQ=255)
                     cTTTTCGTTAGCGACTTTTTTTTTGTTCAGCAACg               >  1:1086708/1‑35 (MQ=255)
                      ttttCGTTAGCGGCTTTTTTTTTGTTCAGCAACGc              <  1:18093844/35‑1 (MQ=255)
                      ttttCGTTAGCGGCTTTTTTTTTGTTCAGCAACGc              <  1:4024713/35‑1 (MQ=255)
                      ttttCGTTAGCGGCTTTTTTTTTGTTCAGCAACGc              <  1:3079814/35‑1 (MQ=255)
                       tttCGTTAGCGGCTTTTTTTTTTTTTTGCAACGCa             >  1:12938480/1‑35 (MQ=25)
                       tttCGTTAGCGACTTTTTTTTTGTTCAGCAACGCa             <  1:8306201/35‑1 (MQ=255)
                        ttCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCAc            <  1:4448331/35‑1 (MQ=255)
                        ttCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCAc            <  1:6321255/35‑1 (MQ=255)
                        ttCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCAc            <  1:14534151/35‑1 (MQ=255)
                         tCGTTAGCGACTTTTTTTTTGTTCAGCAACGCACt           >  1:14326701/1‑35 (MQ=255)
                         tCGTTAGCGACTTTTTTTTTGTTCAGCAACGCACt           <  1:16185490/35‑1 (MQ=255)
                         tCGTTAGCGACTTTTTTTTTGTTCAGCAACGCACt           <  1:16854952/35‑1 (MQ=255)
                         tCGTTAGCGACTTTTTTTTTGTTCAGCAACGAACt           >  1:11370712/1‑35 (MQ=37)
                          cGTTAGCGGCTTTTTTTTTGTTCAGCAACGCACTa          <  1:14748810/35‑1 (MQ=255)
                          cGTTAGCGACTTTTTTTTTGTTCAGCAACGCACTa          <  1:20572063/35‑1 (MQ=255)
                          cGTTAGCGACTTTTTTTTTGTTCAGCAACGCACTa          <  1:7145175/35‑1 (MQ=255)
                           gTTAGCGGCTTTTTTTTGTTCAGCAACGCACTAcc        <  1:759106/35‑1 (MQ=37)
                             tAGCGACTTTTTTTTTGTTCAGCAACGCACTACCa       >  1:16387864/1‑35 (MQ=255)
                             tAGCGACTTTTTTTTTGTTCAGCAACGCACTACCa       <  1:18896293/35‑1 (MQ=255)
                              ngCGGTNTTTTTTTTGTTCAGCAACGCACTACCAt      <  1:5204116/34‑1 (MQ=37)
                              aGCGGCTTTTTTTTTGTTCAGCAACGCACTACCAt      <  1:19669074/35‑1 (MQ=255)
                              aGCGACTTTTTTTTTGTTCAGCAACGCACTACGAt      <  1:7936319/35‑1 (MQ=25)
                              aGCGACTTTTTTTTTGTTCAGCAACGCACTACCAt      <  1:12558159/35‑1 (MQ=255)
                               gctGCTTTTTTTTTGTTCAGCAACGCACTACCATa     <  1:5399452‑M1/32‑1 (MQ=255)
                               gCGGCTTTTTTTTTTTTCAGCAACGCACTACCATa     >  1:20135222/1‑35 (MQ=255)
                                CGACTTTTTTTTTGTTCAGCAACGCACTACCATAc    >  1:18682827‑M1/1‑35 (MQ=255)
                                 GGCTTTTTTTTTGTTCAGCAACGCACTACCATACt   >  1:18188522‑M1/1‑35 (MQ=255)
                                 GGCTTTTTTTTTGTTCAGCAACGCACTACCATACt   <  1:17002942‑M1/35‑1 (MQ=255)
                                 GGCTTTTTTTTTGTTCAGCAACGCACTACCATACt   <  1:14628923‑M1/35‑1 (MQ=255)
                                  |                                  
AACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCACTACCATACTT  >  REL606/1691476‑1691544

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.