Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,619,834 IS150 (–) +3 bp 100% intergenic (‑71/+139) hokD ← / ← ECB_01533 small toxic polypeptide/conserved hypothetical protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 664688 =NA (NA)89 (0.990) 56/76 NT 98.8% noncoding (1/1443 nt) IS150 repeat region
?REL606 1619834 = 1 (0.010)intergenic (‑71/+141) hokD/ECB_01533 small toxic polypeptide/conserved hypothetical protein
* ? REL606 = 666130NA (NA)70 (0.760) 46/78 NT 98.5% noncoding (1443/1443 nt) IS150 repeat region
?REL606 = 1619836 1 (0.010)intergenic (‑73/+139) hokD/ECB_01533 small toxic polypeptide/conserved hypothetical protein

GACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCAC  >  REL606/666083‑666175
                                               |                                              
gACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAttg                                            <  1:5342077/51‑2 (MQ=18)
         gCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCtt                                    <  1:997936/51‑1 (MQ=14)
         gCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCtt                                    <  1:448683/51‑1 (MQ=14)
          ccTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttt                                   >  1:1434745/1‑51 (MQ=14)
           cTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCtttt                                  >  1:9280813/1‑51 (MQ=17)
           cTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCtttt                                  >  1:2881652/1‑51 (MQ=17)
            tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttttt                                 <  1:8275728/51‑1 (MQ=17)
            tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttttt                                 <  1:5099631/51‑1 (MQ=17)
            tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttttt                                 >  1:1934377/1‑51 (MQ=17)
            tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttttt                                 <  1:6615816/51‑1 (MQ=17)
            tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttttt                                 >  1:7179897/1‑51 (MQ=17)
             cGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTc                                >  1:6792164/1‑51 (MQ=17)
              gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCa                               <  1:7610303/51‑1 (MQ=18)
               tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                              <  1:7342744/51‑1 (MQ=18)
               tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                              <  1:10024178/51‑1 (MQ=18)
               tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                              >  1:333166/1‑51 (MQ=18)
               tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                              <  1:8990722/51‑1 (MQ=18)
               tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                              >  1:6601617/1‑51 (MQ=18)
               tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                              <  1:2012718/51‑1 (MQ=18)
               tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                              >  1:2384194/1‑51 (MQ=18)
               tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                              <  1:5168681/51‑1 (MQ=18)
               tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                              <  1:10822822/51‑1 (MQ=18)
               tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                              >  1:4827349/1‑51 (MQ=18)
                gtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACg                             >  1:84205/1‑51 (MQ=18)
                gtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACg                             <  1:7360990/51‑1 (MQ=18)
                 tttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACg                             >  1:7558621/1‑50 (MQ=18)
                 tttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGa                            <  1:1023650/51‑1 (MQ=18)
                 tttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGa                            <  1:3217572/51‑1 (MQ=18)
                 tttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGa                            >  1:4479159/1‑51 (MQ=18)
                  ttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTc                                <  1:9848575/46‑1 (MQ=18)
                  ttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTc                                <  1:6864998/46‑1 (MQ=18)
                  ttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAg                           >  1:9123523/1‑51 (MQ=21)
                  ttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAg                           >  1:834322/1‑51 (MQ=21)
                  ttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAg                           >  1:1961249/1‑51 (MQ=21)
                  ttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAg                           <  1:4552814/51‑1 (MQ=21)
                  ttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAg                           >  1:7783829/1‑51 (MQ=21)
                   tAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGc                          >  1:4106993/1‑51 (MQ=21)
                   tAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGc                          >  1:8273448/1‑51 (MQ=21)
                   tAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGc                          >  1:1603083/1‑51 (MQ=21)
                    aaCTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGcc                         >  1:9622666/1‑51 (MQ=21)
                    aaCTGTCCAACTTTTTGGGGTCACTACATTTTCAGGGCTTTTTCACGAGcc                         <  1:7996587/51‑1 (MQ=21)
                     aCTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCg                        <  1:5959410/51‑1 (MQ=25)
                     aCTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCg                        <  1:7403951/51‑1 (MQ=25)
                      cTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCgg                       <  1:10493974/51‑1 (MQ=25)
                       tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCgg                       <  1:5264045/50‑1 (MQ=25)
                       tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                      <  1:3638213/51‑1 (MQ=25)
                       tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                      <  1:2395099/51‑1 (MQ=25)
                       tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                      <  1:6169288/51‑1 (MQ=25)
                       tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                      <  1:4179080/51‑1 (MQ=25)
                       tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                      >  1:10788621/1‑51 (MQ=25)
                       tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                      >  1:4247414/1‑51 (MQ=25)
                       tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                      >  1:3125213/1‑51 (MQ=25)
                        gTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgcg                     >  1:7000270/1‑51 (MQ=25)
                        gTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgcg                     >  1:4146389/1‑51 (MQ=25)
                        gTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgcg                     >  1:886532/1‑51 (MQ=25)
                         tCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGcc                         <  1:8390772/46‑1 (MQ=25)
                         tCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgcgc                    >  1:7793163/1‑51 (MQ=37)
                          ccAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCt                   >  1:8752331/1‑51 (MQ=37)
                           cAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCtt                  <  1:7382956/51‑1 (MQ=37)
                            aaCTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCttt                 <  1:1170625/51‑1 (MQ=38)
                             aCTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTa                >  1:4812027/1‑51 (MQ=38)
                              cTTTTTGGGGTCAGTACATTTTCAGGGTTTTTTCACGAGCCGGCGCttt                 <  1:2709305/49‑1 (MQ=37)
                               tttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATc              <  1:9823213/51‑1 (MQ=38)
                               tttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATc              >  1:8673015/1‑51 (MQ=38)
                               tttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATc              >  1:6456867/1‑51 (MQ=38)
                               tttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATc              >  1:4367519/1‑51 (MQ=38)
                                ttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcg             >  1:5429993/1‑51 (MQ=39)
                                ttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcg             <  1:6296160/51‑1 (MQ=39)
                                ttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcg             >  1:4822934/1‑51 (MQ=39)
                                 tttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcgc            <  1:7808204/51‑1 (MQ=255)
                                 tttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcgc            >  1:881823/1‑51 (MQ=255)
                                  ttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcgcg           >  1:7919436/1‑51 (MQ=255)
                                  ttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcgcg           >  1:6521583/1‑51 (MQ=255)
                                   tGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGt          <  1:6921465/51‑1 (MQ=255)
                                      ggTCAGTACATTTTCAGGGTTTTTTCACGAGCCGGCGCTTTATCGCGTTTa       <  1:10238063/51‑1 (MQ=255)
                                      ggTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTa       >  1:9745683/1‑51 (MQ=255)
                                      ggTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTa       <  1:8469558/51‑1 (MQ=255)
                                      ggTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTa       <  1:4372546/51‑1 (MQ=255)
                                      ggTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTa       <  1:10186258/51‑1 (MQ=255)
                                       gTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTaa      <  1:1118518/51‑1 (MQ=255)
                                        tCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAg     <  1:3471592/51‑1 (MQ=255)
                                        tCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAg     <  1:7199754/51‑1 (MQ=255)
                                        tCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAg     >  1:1194321/1‑51 (MQ=255)
                                        tCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAg     >  1:7889880/1‑51 (MQ=255)
                                         cAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGc    >  1:158141/1‑51 (MQ=255)
                                         cAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGc    >  1:1214990/1‑51 (MQ=255)
                                         cAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGc    <  1:366655/51‑1 (MQ=255)
                                          aGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCa   >  1:10343001/1‑51 (MQ=255)
                                           gTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCAc  <  1:1828050/51‑1 (MQ=255)
                                               |                                              
GACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCAC  >  REL606/666083‑666175

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