Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 896,463 IS150 (+) +3 bp 100% coding (159‑161/432 nt) ECB_00838 → putative phage tail protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 666130NA (NA)75 (0.790) 47/80 NT 100% noncoding (1443/1443 nt) IS150 repeat region
?REL606 896463 = 0 (0.000)coding (159/432 nt) ECB_00838 putative phage tail protein
* ? REL606 = 8964650 (0.000)87 (0.970) 49/76 NT 100% coding (161/432 nt) ECB_00838 putative phage tail protein
?REL606 = 2775877 NA (NA)noncoding (1/1443 nt) IS150 repeat region

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.