Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,621,349 C→A 5.8% intergenic (+50/+35) trpR → / ← yjjX Trp operon repressor/NTPase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,621,3490C→A 5.8% 4.7 ‑0.3 ‑0.3129intergenic (+50/+35)trpR/yjjXTrp operon repressor/NTPase
Reads supporting (aligned to +/- strand):  new base (3/5):  ref base (57/64):  total (60/69)
Fisher's exact test strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.67e-01

Only 100 of 135 total aligned reads displayed.

AAGACGCGGCAGCGTCGCATCAGGCATCGTGCACCGAATGCCGGATGCGGCGTGAACGCCTTATCCGTC  >  REL606/4621315‑4621383
                                  |                                  
aaGACGCGGCAGCGTCGCATCAGGCATCGTGCAcc                                    >  1:14309323/1‑35 (MQ=38)
 aGACGCGGCAGCGTCGCATCAGGCATCGTGCACCg                                   <  1:13450511/35‑1 (MQ=40)
 aGACGCGGCAGCGTCGCATCAGGCATCGTGCACCg                                   <  1:14349443/35‑1 (MQ=40)
 aGACGCGGCAGCGTCGCATCAGGCATCGTGCACCg                                   <  1:18434816/35‑1 (MQ=40)
    cgcgGCAGCGTCGCATCAGGCATCGTGCACCGAAt                                <  1:11322415/35‑1 (MQ=255)
    cgcgGCAGCGTCGCATCAGGCATCGTGCACAGAAt                                <  1:14774207/35‑1 (MQ=255)
      cgGCAGCGGCGCATCAGGCATCGTGCACCGAATGc                              >  1:16902292/1‑35 (MQ=255)
       ggCAGCGTCGCATCAGGCATCGTGCACCGAGTGcc                             <  1:19171331/35‑1 (MQ=255)
       ggCAGCGTCGCATCAGGCATCGTGCACCGAATGcc                             >  1:7336409/1‑35 (MQ=255)
       ggCAGCGTCGCATCAGGCATCGTGCACCGAATGcc                             <  1:11977620/35‑1 (MQ=255)
       ggCAGCGTCGCATCAGGCATCGTGCACCGAATGcc                             <  1:12239136/35‑1 (MQ=255)
        gCAGCGTCGCATCAGGCATCGTGCACCGAATGcct                            >  1:16040442/1‑34 (MQ=255)
        gCAGCGTCGCATCAGGCATCGTGCACCGAATGCCg                            >  1:17234061/1‑35 (MQ=255)
        gCAGCGTCGCATCAGGCATCGTGCACCGAATGCCg                            >  1:7843335/1‑35 (MQ=255)
         cAGCGTCGCATCAGGCATCGTGCACCGAATGCCgg                           >  1:16611338/1‑35 (MQ=40)
         cAGCGTCGCATCAGGCATCGTGCACCGAATGCCgg                           <  1:11100704/35‑1 (MQ=40)
         cAGCGTCGCATCAGGCATCGTGCACCGAATGCCgg                           <  1:1090875/35‑1 (MQ=40)
          aGCGTCGCATCAGGCATCGTGCACCGAATGCCGGa                          >  1:2856124/1‑35 (MQ=38)
          aGCGTCGCATCAGGCATCGTGCACCGAATGCCGGa                          >  1:8987437/1‑35 (MQ=38)
           gCGTCGCATCAGGCATCGTGCACCGAATGCCGGAt                         >  1:20173754/1‑35 (MQ=38)
           gCGTCGCATCAGGCATCGTGCACCGAATGCCGGAt                         >  1:17546447/1‑35 (MQ=38)
           gCGTCGCATCAGGCATCGTGCACCGAATGCCGGAt                         <  1:16610985/35‑1 (MQ=38)
           gCGTCGCATCAGGCATCGTGCACCGAATGCCGGAt                         <  1:4971849/35‑1 (MQ=38)
            cGTCGCATCAGGCATCGTGCACCGAATGCCGGATg                        >  1:7433888/1‑35 (MQ=38)
            cGTCGCATCAGGCATCGTGCACCGAATGCCGGATg                        <  1:7203996/35‑1 (MQ=38)
            cGTCGCATCAGGCATCGTGCACCGAATGCCGGATg                        >  1:19845519/1‑35 (MQ=38)
             gTCGCATCGGGCATCGTGTACCGAATGCCGGATgc                       <  1:9151616/35‑1 (MQ=14)
             gTCGCATCAGGCATCGTGCACCGAATGCCGGATgc                       >  1:14899749/1‑35 (MQ=38)
             gTCGCATCAGGCATCGTGCACCGAATGCCGGATgc                       <  1:17564424/35‑1 (MQ=38)
             gTCGCATCAGGCATCGTGCACCGAATGCCGGATgc                       >  1:4434978/1‑35 (MQ=38)
             gTCGCATCAGGCATCGTGCACCGAATGCCGGATgc                       <  1:437018/35‑1 (MQ=38)
             gTCGCATCAGGCATCGTGCACCGAATGCCGGATgc                       <  1:2485931/35‑1 (MQ=38)
              tcccATCAGGCATCGTGCACCGAATGCCGGATgcg                      <  1:20353723/32‑1 (MQ=255)
              tagCATCAGGCATCGTGCACCGAATGCCGGATgcg                      <  1:12014294/33‑1 (MQ=255)
              tCGCATCAGGCATCGTGCACCGAATGCCGGATgcg                      <  1:18800519/35‑1 (MQ=38)
              tCGCATCAGGCATCGTGCACCGAATGCCGGATgcg                      >  1:1454876/1‑35 (MQ=38)
              tCGCATCAGGCATCGTGCACCGAATGCCGGATgcg                      >  1:4514888/1‑35 (MQ=38)
              tCGCATCAGGCATCGTGCACCGAATGCCGGATgcg                      <  1:12579413/35‑1 (MQ=38)
               cGCATCAGGCATCGTGCACCGAATGCCGGATgcgg                     >  1:5073236/1‑35 (MQ=38)
               cGCATCAGGCATCGTGCACCGAATGCCGGATgcgg                     <  1:288078/35‑1 (MQ=38)
               cGCATCAGGCATCGTGCACCGAATGCCGGATgcgg                     >  1:7573529/1‑35 (MQ=38)
               cGCATCAGGCATCGTGCACCGAATGCCGGATgcgg                     >  1:10080971/1‑35 (MQ=38)
               cGCATCAGGCATCGTGCACCGAATGCCGGATgcgc                     >  1:15148035/1‑34 (MQ=38)
                gCATCAGGCATTGTGCACAGAACGGCGGGGgcggg                    >  1:19002163/1‑34 (MQ=17)
                gCATCAGGCATCGTGCACCGAATGCCGGATgcggc                    <  1:19247952/35‑1 (MQ=38)
                 cATCAGGCATCGTGCACCGAATGCCGGATgcggcg                   >  1:13734076/1‑35 (MQ=38)
                 cATCAGGCATCGTGCACCGAATGCCGGATgcggcg                   >  1:360047/1‑35 (MQ=38)
                 cATCAGGCATCGTGCACCGAATGCCGGATgcggcg                   <  1:17657130/35‑1 (MQ=38)
                 cATCAGGCATCGTGCACCGAATGCCGGATgcggcg                   >  1:12479968/1‑35 (MQ=38)
                 cATCAGGCATCGTGCACCGAATGCCGGATgcggcg                   <  1:5002437/35‑1 (MQ=38)
                  aTCAGGCATCGTGCACCGAATGCCGGATGTGGCGt                  <  1:6984531/35‑1 (MQ=255)
                  aTCAGGCATCGTGCACCGAATGCCGGATGCGGCGt                  <  1:13398/35‑1 (MQ=38)
                   tCAGGCATCGTGCACCGAATGCCGGATGCGGCGTg                 >  1:4012187/1‑35 (MQ=38)
                   tCAGGCATCGTGCACCGAATGCCGGATGCGGCGTg                 >  1:2826233/1‑35 (MQ=38)
                   tCAGGCATCGTGCACCGAATGCCGGATGCGGCGTg                 >  1:14036900/1‑35 (MQ=38)
                   tCAGGCATCGTGCACCGAATGCCGGATGCGGCGTg                 <  1:11346616/35‑1 (MQ=38)
                   tCAGGCATCGTGCACAGAATGCCGGATGCGGCGTg                 <  1:2158021/35‑1 (MQ=255)
                    cAGGCATCGTGCACCGAATGCCGGATGCGGCGTGa                >  1:12078681/1‑35 (MQ=38)
                    cAGGCATCGTGCACCGAATGCCGGATGCGGCGTGa                <  1:14224993/35‑1 (MQ=38)
                    cAGGCATCGTGCACCGAATGCCGGATGCGGCGTGa                <  1:6981712/35‑1 (MQ=38)
                     aGGCATCGTGCACCGAATGCCGGATGCGGCGTGaa               <  1:11296337/35‑1 (MQ=38)
                     aGGCATCGTGCACCGAATGCCGGATGCGGCGTGaa               >  1:291775/1‑35 (MQ=38)
                     aGGCATCGTGCACCGAATGCCGGATGCGGCGTGaa               <  1:7976640/35‑1 (MQ=38)
                      ggCATCGTGCACCGAATGCCGGATGCGGCGTGAAc              <  1:10476558/35‑1 (MQ=38)
                      ggCATCGTGCACCGAATGCCGGATGCGGCGTGAAc              >  1:4068562/1‑35 (MQ=38)
                      ggCATCGTGCACCGAATGCCGGATGCGGCGTGAAc              <  1:4737277/35‑1 (MQ=38)
                       gCATCGTGCACCGAATGCCGGATGCGGCGTGAACg             <  1:6203471/35‑1 (MQ=38)
                       gCATCGTGCACCGAATGCCGGATGCGGCGTGAACg             <  1:12844720/35‑1 (MQ=38)
                       gCATCGTGCACCGAATGCCGGATGCGGCGTGAACg             <  1:14737174/35‑1 (MQ=38)
                        cATCGTGCACCGAATGCCGTATGCGGCGTGAACGc            >  1:10959033/1‑35 (MQ=255)
                        cATCGTGCACCGAATGCCGGATGCGGGGGGaaagg            >  1:1008115/1‑32 (MQ=14)
                        cATCGTGCACCGAATGCCGGATGCGGCGTGAACGc            <  1:12010655/35‑1 (MQ=38)
                        cATCGTGCACCGAATGCCGGATGCGGCGTGAACGc            <  1:15174901/35‑1 (MQ=38)
                        cATCGTGCACCGAATGCCGGATGCGGCGTGAACGc            <  1:7530173/35‑1 (MQ=38)
                         aTCGTGCACCGAATGCCGGATGCGGGGTGAACGcc           >  1:14512971/1‑35 (MQ=255)
                          tCGTGCACCGAATGCCGGATGCGGCGTGAACGCCt          >  1:15304508/1‑35 (MQ=38)
                          tCGTGCACCGAATGCCGGATGCGGCGTGAACGCCt          <  1:9178408/35‑1 (MQ=38)
                          tCGTGCACCGAATGCCGGATGCGGCGTGAACGCCt          >  1:18831084/1‑35 (MQ=38)
                          tCGTGCACCGAATACCGGATGCGGCGTGAACGCCt          <  1:19472937/35‑1 (MQ=255)
                           cGTGCACCGAATGCCGGATGCGGCGTGAACGCCtt         >  1:7263301/1‑35 (MQ=38)
                           cGTGCACCGAATGCCGGATGCGGCGTGAACGCCtt         <  1:7992646/35‑1 (MQ=38)
                           cGTGCACCGAATGCCGGATGCGGCGTGAACGCCtt         >  1:14502105/1‑35 (MQ=38)
                            gTGCACCGAATGCCGGATGCGGCGTGAACGCCTTa        <  1:10933221/35‑1 (MQ=38)
                             tGCACCGAATGCCGGATGCGGCGTGAACGCCTTAt       <  1:9775381/35‑1 (MQ=38)
                             tGCACCGAATGCCGGATGCGGCGTGAACGCCTTAt       >  1:6281174/1‑35 (MQ=38)
                              ggaCCGAATGCCGGATGCGGCGTGAACGCCTTATc      <  1:1699759/33‑1 (MQ=255)
                              gCACCGAATGCCGGATGCGGCGTGAACGCCTTATc      <  1:5928076/35‑1 (MQ=38)
                              gCACCGAATGCCGGATGCGGCGTGAACGCCTTATc      <  1:5515440/35‑1 (MQ=38)
                              gCACCGAATGCCGGATGCGGCGTGAACGCCTTATc      <  1:2984078/35‑1 (MQ=38)
                              gCACCGAATGCCGGATGCGGCGTGAACGCCTTATc      >  1:16802147/1‑35 (MQ=38)
                              gCACAGAATGCCGGATGCGGCGTGAACGCCTTATc      >  1:15905455/1‑35 (MQ=255)
                               cACCGAATGCCGGATGCGGCGTGAACGCCTTATcc     <  1:11505308/35‑1 (MQ=38)
                               cACCGAATGCCGGATGCGGCGTGAACGCCTTATcc     <  1:10930340/35‑1 (MQ=38)
                               cACCGAATGCCGGATGCGGCGTGAACGCCTTATcc     <  1:13915361/35‑1 (MQ=38)
                               cACCGAATGCCGGATGCGGCGTGAACGCCTTATcc     <  1:9556542/35‑1 (MQ=38)
                               cACAGAATGCCGGATGCGGCGTGAACGCCTTATcc     <  1:583079/35‑1 (MQ=255)
                                aCCGAATGCCGGATGCGGCGTGAACGCCTTATCCg    <  1:7016371/35‑1 (MQ=38)
                                aCCGAATGCCGGATGCGGCGTGAACGCCTTATCCg    >  1:6130048/1‑35 (MQ=38)
                                aCCGAATGCCGGATGCGGCGTGAACGCCTTATCCg    >  1:4834435/1‑35 (MQ=38)
                                  cGAATGCCGGATGCGGCGTGAACGCCTTATCCGTc  >  1:13286470/1‑35 (MQ=40)
                                  |                                  
AAGACGCGGCAGCGTCGCATCAGGCATCGTGCACCGAATGCCGGATGCGGCGTGAACGCCTTATCCGTC  >  REL606/4621315‑4621383

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