Predicted mutation
evidence seq id position mutation annotation gene description
JC JC REL606 3,537,178 +GTACTGACCCCAAAA :: IS150 (–) :: +CAA coding (40/2199 nt) zntA → zinc, cobalt and lead efflux system

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 666133NA (NA)67 (1.110)
+TGGGGTCAGTAC
43/76 0.1 100% pseudogene (355/387 nt) gltL ATP‑binding protein of glutamate/aspartate transport system; b0652_1
?REL606 = 3537178 0 (0.000)coding (40/2199 nt) zntA zinc, cobalt and lead efflux system
* ? REL606 3537179 =0 (0.000)82 (1.100)
+TTG
52/94 0.1 100% coding (41/2199 nt) zntA zinc, cobalt and lead efflux system
?REL606 3651091 = NA (NA)noncoding (1/1443 nt) IS150 repeat region

GCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCA  >  REL606/666092‑666180
                                          |                                               
gCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCtt                                         <  1:997936/51‑1 (MQ=14)
gCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCtt                                         <  1:448683/51‑1 (MQ=14)
 ccTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttt                                        >  1:1434745/1‑51 (MQ=14)
  cTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCtttt                                       >  1:2881652/1‑51 (MQ=17)
   tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttttt                                      <  1:6615816/51‑1 (MQ=17)
   tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttttt                                      <  1:5099631/51‑1 (MQ=17)
   tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttttt                                      >  1:7179897/1‑51 (MQ=17)
   tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttttt                                      >  1:1934377/1‑51 (MQ=17)
    cGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTc                                     >  1:6792164/1‑51 (MQ=17)
      tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                                   >  1:4827349/1‑51 (MQ=18)
      tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                                   <  1:5168681/51‑1 (MQ=18)
      tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                                   <  1:2012718/51‑1 (MQ=18)
      tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                                   >  1:2384194/1‑51 (MQ=18)
      tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                                   >  1:333166/1‑51 (MQ=18)
      tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                                   >  1:6601617/1‑51 (MQ=18)
       gtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACg                                  >  1:84205/1‑51 (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:6864998/46‑1 (MQ=18)
         ttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAg                                <  1:4552814/51‑1 (MQ=21)
         ttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAg                                >  1:834322/1‑51 (MQ=21)
         ttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAg                                >  1:1961249/1‑51 (MQ=21)
          tAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGc                               >  1:4106993/1‑51 (MQ=21)
          tAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGc                               >  1:1603083/1‑51 (MQ=21)
            aCTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCg                             <  1:5959410/51‑1 (MQ=25)
              tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCgg                            <  1:5264045/50‑1 (MQ=25)
              tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                           >  1:3125213/1‑51 (MQ=25)
              tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                           <  1:3638213/51‑1 (MQ=25)
              tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                           <  1:4179080/51‑1 (MQ=25)
              tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                           <  1:2395099/51‑1 (MQ=25)
              tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                           <  1:6169288/51‑1 (MQ=25)
              tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                           >  1:4247414/1‑51 (MQ=25)
               gTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgcg                          >  1:886532/1‑51 (MQ=25)
               gTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgcg                          >  1:4146389/1‑51 (MQ=25)
               gTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgcg                          >  1:7000270/1‑51 (MQ=25)
                   aaCTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCttt                      <  1:1170625/51‑1 (MQ=38)
                    aCTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTa                     >  1:4812027/1‑51 (MQ=38)
                     cTTTTTGGGGTCAGTACATTTTCAGGGTTTTTTCACGAGCCGGCGCttt                      <  1:2709305/49‑1 (MQ=37)
                      tttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATc                   >  1:4367519/1‑51 (MQ=38)
                      tttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATc                   >  1:6456867/1‑51 (MQ=38)
                       ttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcg                  >  1:4822934/1‑51 (MQ=39)
                       ttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcg                  >  1:5429993/1‑51 (MQ=39)
                       ttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcg                  <  1:6296160/51‑1 (MQ=39)
                        tttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcgc                 >  1:881823/1‑51 (MQ=255)
                         ttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcgcg                >  1:6521583/1‑51 (MQ=255)
                          tGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGt               <  1:6921465/51‑1 (MQ=255)
                             ggTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTa            <  1:4372546/51‑1 (MQ=255)
                              gTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTaa           <  1:1118518/51‑1 (MQ=255)
                               tCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAg          >  1:1194321/1‑51 (MQ=255)
                               tCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAg          <  1:3471592/51‑1 (MQ=255)
                               tCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAg          <  1:7199754/51‑1 (MQ=255)
                                cAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGc         >  1:158141/1‑51 (MQ=255)
                                cAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGc         <  1:366655/51‑1 (MQ=255)
                                cAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGc         >  1:1214990/1‑51 (MQ=255)
                                  gTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCAc       <  1:1828050/51‑1 (MQ=255)
                                   tacattttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACt      >  1:5053943/6‑51 (MQ=255)
                                    acattttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACtt     <  1:6154015/47‑1 (MQ=255)
                                    acattttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACtt     >  1:879123/5‑51 (MQ=255)
                                     cattttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACttt    <  1:6553703/48‑1 (MQ=255)
                                     cattttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACttt    >  1:6191575/4‑51 (MQ=255)
                                     cattttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACttt    >  1:334635/4‑51 (MQ=255)
                                      attttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTc   <  1:136342/49‑1 (MQ=255)
                                      attttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTc   >  1:6512795/3‑51 (MQ=255)
                                       ttttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCa  <  1:6338214/50‑1 (MQ=255)
                                       ttttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCa  <  1:5904792/50‑1 (MQ=255)
                                       ttttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCa  >  1:593944/2‑51 (MQ=255)
                                          |                                               
GCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCA  >  REL606/666092‑666180

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.