Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 3,537,178 +GTACTGACCCCAAAA :: IS150 (–) :: +CAA 100% 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)65 (0.990)
+TGGGGTCAGTAC
37/56 NT 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)105 (1.200)
+TTG
51/74 NT 100% coding (41/2199 nt) zntA zinc, cobalt and lead efflux system
?REL606 3651091 = NA (NA)noncoding (1/1443 nt) IS150 repeat region

Only 100 of 103 total aligned reads displayed.

GCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCT  >  REL606/666092‑666183
                                          |                                                  
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:1934377/1‑51 (MQ=17)
   tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttttt                                         <  1:8275728/51‑1 (MQ=17)
   tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttttt                                         <  1:5099631/51‑1 (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:4827349/1‑51 (MQ=18)
      tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                                      <  1:10024178/51‑1 (MQ=18)
      tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                                      <  1:10822822/51‑1 (MQ=18)
      tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                                      >  1:6601617/1‑51 (MQ=18)
      tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                                      <  1:5168681/51‑1 (MQ=18)
      tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                                      >  1:333166/1‑51 (MQ=18)
      tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                                      <  1:8990722/51‑1 (MQ=18)
      tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                                      >  1:2384194/1‑51 (MQ=18)
      tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCAc                                      <  1:2012718/51‑1 (MQ=18)
       gtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACg                                     <  1:7360990/51‑1 (MQ=18)
       gtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACg                                     >  1:84205/1‑51 (MQ=18)
        tttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACg                                     >  1:7558621/1‑50 (MQ=18)
        tttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGa                                    >  1:4479159/1‑51 (MQ=18)
        tttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGa                                    <  1:1023650/51‑1 (MQ=18)
        tttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGa                                    <  1:3217572/51‑1 (MQ=18)
         ttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTc                                        <  1:9848575/46‑1 (MQ=18)
         ttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTc                                        <  1:6864998/46‑1 (MQ=18)
         ttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAg                                   >  1:834322/1‑51 (MQ=21)
         ttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAg                                   >  1:7783829/1‑51 (MQ=21)
         ttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAg                                   >  1:9123523/1‑51 (MQ=21)
         ttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAg                                   >  1:1961249/1‑51 (MQ=21)
         ttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAg                                   <  1:4552814/51‑1 (MQ=21)
          tAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGc                                  >  1:1603083/1‑51 (MQ=21)
          tAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGc                                  >  1:4106993/1‑51 (MQ=21)
          tAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGc                                  >  1:8273448/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:10788621/1‑51 (MQ=25)
              tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                              <  1:6169288/51‑1 (MQ=25)
              tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                              >  1:4247414/1‑51 (MQ=25)
              tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                              <  1:2395099/51‑1 (MQ=25)
              tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                              >  1:3125213/1‑51 (MQ=25)
              tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgc                              <  1:3638213/51‑1 (MQ=25)
               gTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgcg                             >  1:7000270/1‑51 (MQ=25)
               gTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgcg                             >  1:886532/1‑51 (MQ=25)
               gTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGgcg                             >  1:4146389/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:8673015/1‑51 (MQ=38)
                      tttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATc                      <  1:9823213/51‑1 (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:4822934/1‑51 (MQ=39)
                       ttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcg                     <  1:6296160/51‑1 (MQ=39)
                        tttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcgc                    >  1:881823/1‑51 (MQ=255)
                        tttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcgc                    <  1:7808204/51‑1 (MQ=255)
                         ttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcgcg                   >  1:6521583/1‑51 (MQ=255)
                         ttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcgcg                   >  1:7919436/1‑51 (MQ=255)
                          tGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGt                  <  1:6921465/51‑1 (MQ=255)
                             ggTCAGTACATTTTCAGGGTTTTTTCACGAGCCGGCGCTTTATCGCGTTTa               <  1:10238063/51‑1 (MQ=255)
                             ggTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTa               <  1:10186258/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)
                              gTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTaa              <  1:1118518/51‑1 (MQ=255)
                               tCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAg             >  1:1194321/1‑51 (MQ=255)
                               tCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAg             >  1:7889880/1‑51 (MQ=255)
                               tCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAg             <  1:3471592/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)
                                 aGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCa           >  1:10343001/1‑51 (MQ=255)
                                  gTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCAc          <  1:1828050/51‑1 (MQ=255)
                                   tacattttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACt         >  1:5053943/6‑51 (MQ=255)
                                    acattttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACtt        >  1:879123/5‑51 (MQ=255)
                                    acattttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACtt        <  1:6154015/47‑1 (MQ=255)
                                     cattttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACttt       <  1:10805764/48‑1 (MQ=255)
                                     cattttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACttt       >  1:334635/4‑51 (MQ=255)
                                     cattttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACttt       >  1:6191575/4‑51 (MQ=255)
                                     cattttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACttt       <  1:6553703/48‑1 (MQ=255)
                                      attttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTc      <  1:136342/49‑1 (MQ=255)
                                      attttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTc      >  1:6512795/3‑51 (MQ=255)
                                       ttttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCa     <  1:5904792/50‑1 (MQ=255)
                                       ttttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCa     <  1:6338214/50‑1 (MQ=255)
                                       ttttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCa     >  1:593944/2‑51 (MQ=255)
                                          tCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCt  >  1:7812297/1‑51 (MQ=255)
                                          tCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCt  <  1:7793324/51‑1 (MQ=255)
                                          tCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCt  >  1:8200497/1‑51 (MQ=255)
                                          |                                                  
GCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCT  >  REL606/666092‑666183

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.