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)89 (1.570)
+TGGGGTCAGTAC
66/176 0.0 99.0% pseudogene (355/387 nt) gltL ATP‑binding protein of glutamate/aspartate transport system; b0652_1
?REL606 = 3537178 1 (0.020)coding (40/2199 nt) zntA zinc, cobalt and lead efflux system
* ? REL606 3537179 =1 (0.020)67 (1.070)
+TTG
52/194 0.2 98.6% coding (41/2199 nt) zntA zinc, cobalt and lead efflux system
?REL606 3651091 = NA (NA)noncoding (1/1443 nt) IS150 repeat region

TGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACAGCTCG  >  REL606/666043‑666230
                                                                                           |                                                                                                 
tGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttt                                                                                          >  2:775918/1‑101 (MQ=255)
tGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttt                                                                                          >  2:1150422/1‑101 (MQ=255)
tGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttt                                                                                          <  1:872176/101‑1 (MQ=255)
 gAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCtttt                                                                                         <  2:1321153/101‑1 (MQ=255)
  aaaTTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttttt                                                                                        <  2:54763/101‑1 (MQ=255)
  aaaTTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttttt                                                                                        <  1:1486536/101‑1 (MQ=255)
  aaaTTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttttt                                                                                        <  2:45971/101‑1 (MQ=255)
    aTTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCa                                                                                      >  1:534251/1‑101 (MQ=255)
        aaaGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAg                                                                                  <  1:1021004/101‑1 (MQ=255)
        aaaGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAg                                                                                  >  2:473770/1‑101 (MQ=255)
          aGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGcc                                                                                >  1:386256/1‑101 (MQ=255)
          aGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGcc                                                                                <  1:514019/101‑1 (MQ=255)
                gACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCt                                                                          >  2:754383/1‑101 (MQ=255)
                     cAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATc                                                                     >  2:495947/1‑101 (MQ=255)
                     cAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATc                                                                     >  2:798676/1‑101 (MQ=255)
                      aaTTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcg                                                                    >  1:1227492/1‑101 (MQ=255)
                      aaTTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcg                                                                    >  1:1235210/1‑101 (MQ=255)
                       aTTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcgc                                                                   >  1:206869/1‑101 (MQ=255)
                         tGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGt                                                                 >  1:798510/1‑101 (MQ=255)
                          gAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGtt                                                                >  1:1229718/1‑101 (MQ=255)
                           aaTATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGttt                                                               <  1:150360/101‑1 (MQ=255)
                             tatCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTaa                                                             >  1:1126712/1‑101 (MQ=255)
                                 ggAATCAGTCCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCAc                                                         <  1:495952/101‑1 (MQ=255)
                                  gAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACt                                                        <  1:183604/101‑1 (MQ=255)
                                   aaTCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACtt                                                       >  1:1560540/1‑101 (MQ=255)
                                    aTCAGACCTATATGGCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACttt                                                      <  2:1143716/101‑1 (MQ=255)
                                    aTCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACttt                                                      >  2:68144/1‑101 (MQ=255)
                                    aTCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACttt                                                      >  2:532109/1‑101 (MQ=255)
                                       aGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCAc                                                   >  1:1219144/1‑101 (MQ=255)
                                        gACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCAcc                                                  <  1:1246411/101‑1 (MQ=255)
                                         aCCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCt                                                 >  2:1273793/1‑101 (MQ=255)
                                          ccTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTg                                                >  1:164424/1‑101 (MQ=255)
                                           cTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTgc                                               >  2:532788/1‑101 (MQ=255)
                                              tatGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGcc                                            <  1:775923/101‑1 (MQ=255)
                                               atGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCa                                           >  2:367352/1‑101 (MQ=255)
                                                tGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAg                                          <  2:248679/101‑1 (MQ=255)
                                                  ccTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAggg                                        >  1:1658054/1‑101 (MQ=255)
                                                   cTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGt                                       >  2:910633/1‑101 (MQ=255)
                                                    tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGATc                                      <  2:1071544/101‑1 (MQ=255)
                                                      gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAg                                    >  2:944486/1‑101 (MQ=255)
                                                      gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAg                                    >  1:595948/1‑101 (MQ=255)
                                                       tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAgg                                   <  2:265284/101‑1 (MQ=255)
                                                       tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAgg                                   >  2:1210726/1‑101 (MQ=255)
                                                       tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAgg                                   >  1:106661/1‑101 (MQ=255)
                                                        gtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGt                                  >  1:146578/1‑101 (MQ=255)
                                                         tttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGTCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGtt                                 <  2:206923/101‑1 (MQ=255)
                                                            aaCTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTc                              >  1:736755/1‑101 (MQ=255)
                                                             aCTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCg                             <  1:798681/101‑1 (MQ=255)
                                                               tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGAt                           >  2:523231/1‑101 (MQ=255)
                                                               tGTCCAACTTTTTGGGGGCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGAt                           >  2:1428502/1‑101 (MQ=255)
                                                                 tCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATaa                         >  2:1503475/1‑101 (MQ=255)
                                                                 tCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATaa                         <  2:456252/101‑1 (MQ=255)
                                                                 tCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATaa                         <  2:1210062/101‑1 (MQ=255)
                                                                  ccAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAAt                        <  2:274111/101‑1 (MQ=255)
                                                                  ccAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAAt                        >  1:1257312/1‑101 (MQ=255)
                                                                    aaCTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCg                      <  1:331095/101‑1 (MQ=255)
                                                                        ttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAg                  <  2:1265962/101‑1 (MQ=255)
                                                                        ttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAg                  <  1:532793/101‑1 (MQ=255)
                                                                             gggTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAg             <  1:532114/101‑1 (MQ=255)
                                                                              ggTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAgg            >  1:418779/1‑101 (MQ=255)
                                                                               gTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAggg           >  1:1144591/1‑101 (MQ=255)
                                                                                tCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGa          >  1:1543606/1‑101 (MQ=255)
                                                                                tCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGa          >  2:898719/1‑101 (MQ=255)
                                                                                  aGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAAc        >  1:347107/1‑101 (MQ=255)
                                                                                  aGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAAc        >  2:1637490/1‑101 (MQ=255)
                                                                                  aGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAAc        <  2:1126707/101‑1 (MQ=255)
                                                                                   gttcattttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACa       <  2:1444095/95‑1 (MQ=255)
                                                                                   gTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACa       <  1:1579741/101‑1 (MQ=255)
                                                                                   gTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACa       <  1:569397/101‑1 (MQ=255)
                                                                                    tacattttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACAg      <  2:164424/96‑1 (MQ=255)
                                                                                       attttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACAGCTc   >  2:1250192/3‑101 (MQ=255)
                                                                                       attttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACAGCTc   >  2:935674/3‑101 (MQ=255)
                                                                                       attttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAACCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACAGCTc   >  2:260394/3‑101 (MQ=255)
                                                                                        ttttCAGGTTTTGTTCAGGAGCCGGCGATTTACCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGTGAACAGCTCg  <  2:1131293/100‑1 (MQ=255)
                                                                                           |                                                                                                 
TGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACAGCTCG  >  REL606/666043‑666230

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