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)32 (1.310)
+TGGGGTCAGTAC
25/254 NT 97.2% pseudogene (355/387 nt) gltL ATP‑binding protein of glutamate/aspartate transport system; b0652_1
?REL606 = 3537178 1 (0.040)coding (40/2199 nt) zntA zinc, cobalt and lead efflux system
* ? REL606 3537179 =1 (0.040)25 (0.960)
+TTG
15/272 NT 96.2% coding (41/2199 nt) zntA zinc, cobalt and lead efflux system
?REL606 3651091 = NA (NA)noncoding (1/1443 nt) IS150 repeat region

GGTCACTTCCTGATCGTCCGCTCGCTGTATCTCTGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCGACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAACGCTCTCTGAAAACGTCTCCGGCTCCCGCTATAGCTGGAAAGTCAGCGGCATGGACTGCGCCGCCTGTG  >  REL606/3537030‑3537325
                                                                                                                                                     |                                                                                                                                                  
accaaaacgttccgagaagttcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaat                                                                                                                                                    >  1:651146‑M1/150‑150 (MQ=255)
accaaaacgtgccgagatgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaatggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaat                                                                                                                                                    >  1:102356‑M1/150‑150 (MQ=255)
accaaaacgtgccgagatgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaat                                                                                                                                                    >  2:548036‑M1/150‑150 (MQ=255)
accaaaacgtgccgagatgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaat                                                                                                                                                    >  2:159490‑M1/150‑150 (MQ=255)
accaaaacgtgccgagatgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaat                                                                                                                                                    >  2:103275‑M1/150‑150 (MQ=255)
     aacgtgccgagatgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTtgctg                                                                                                                                               >  1:109848‑M1/145‑150 (MQ=255)
              agatgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAAc                                                                                                                                      >  1:339517‑M1/136‑150 (MQ=255)
                      ctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAAcc                                                                                                                              >  2:128868‑M1/128‑150 (MQ=255)
                      ctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAAcc                                                                                                                              >  2:113846‑M1/128‑150 (MQ=255)
                                  tcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAAc                                                                                                                  >  2:39647‑M1/116‑150 (MQ=255)
                                  caggttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAAc                                                                                                                  >  1:133173‑M1/116‑150 (MQ=255)
                                     gttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAACCCGCTAACCACGGTACAGAACGcc                                                                                                               <  1:54009‑M1/38‑1 (MQ=255)
                                     gttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAgaacctg                                                                                                               <  2:133557‑M1/38‑4 (MQ=255)
                                     gttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGcc                                                                                                               <  2:566740‑M1/38‑1 (MQ=255)
                                     gttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagcgtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTACAAACCACTAACCACGGTACAGAACGcc                                                                                                               <  1:599961‑M1/38‑1 (MQ=255)
                                       tgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCaa                                                                                                             >  2:154088‑M1/111‑150 (MQ=255)
                                             gtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTg                                                                                                       <  2:455458‑M1/46‑1 (MQ=255)
                                                                      ttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCCAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGc                                                                              >  2:423842‑M1/80‑150 (MQ=255)
                                                                                ggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAACGc                                                                    <  1:344519‑M1/81‑1 (MQ=255)
                                                                                          tgtataagagacaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAACGCTCTCTGaaaa                                                          >  1:546855‑M1/60‑150 (MQ=255)
                                                                                            ttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAACGCTCTCTGAAAACg                                                        >  2:321185‑M1/58‑150 (MQ=255)
                                                                                                   attaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAACGCTCTCTGAAAACGTCTCCgg                                                 <  2:184563‑M1/100‑1 (MQ=255)
                                                                                                   attaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAACGCTCTCTGAAAACGTCTCCgg                                                 <  1:507839‑M1/100‑1 (MQ=255)
                                                                                                   attaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAACGCTCTCTGAAAACGTCTCCgg                                                 >  1:435364‑M1/51‑150 (MQ=255)
                                                                                                       gaTGCCATGTCGACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAACGCTCTCTGAAAACGTCTCCGGCTcc                                             >  2:115989/1‑150 (MQ=255)
                                                                                                        gggctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAAcgctctctgaaaactgtctcttataca                                            <  2:548440‑M1/105‑14 (MQ=255)
                                                                                                           ctattccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAACGCTCTCTGAAAACGTCTCCGGCTCCCGCt                                         >  1:454167‑M1/43‑150 (MQ=255)
                                                                                                                ccatttcatcgtccaacaaaatgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAACGCTCTCTGAAAACGTCTCCGGCTCCCGCTATAGc                                    <  2:483450‑M1/113‑1 (MQ=255)
                                                                                                                                     tgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAACGCTCTCTGAAAACGTCTCCGGCTCCCGCTATAGCTGGAAAGTCAGCGGCATGGAc               >  2:75240‑M1/17‑150 (MQ=255)
                                                                                                                                     tgggtgcagtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAACGCTCTCTGAAAACGTCTCCGGCTCCCGCTATAGCTGGAAAGTCAGCGGCATGGAc               >  2:207131‑M1/17‑150 (MQ=255)
                                                                                                                                     tgggtgaagtacacaaTTGCTGCGTTCAAACCGTTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAACGCTCTCTGAAAACGTCTCCGGCTCCCGCTATAGCTGGAAAGTCAGCGGCATGGAc               >  1:565092‑M1/17‑150 (MQ=255)
                                                                                                                                             gtacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAACGCTCTCTGAAAACGTCTCCGGCTCCCGCTATAGCTGGAAAGTCAGCGGCATGGACTGcgccgc       <  2:57455‑M1/142‑1 (MQ=255)
                                                                                                                                              tacacaaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAACGCTCTCTGAAAACGTCTCCGGCTCCCGCTATAGCTGGAAAGTCAGCGGCATGGACTGCgccgcc      >  2:653375‑M1/8‑150 (MQ=255)
                                                                                                                                                  caaTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAACGCTCTCTGAAAACGTCTCCGGCTCCCGCTATAGCTGGAAAGTCAGCGGCATGGACTGCGCCGCCtgtg  <  1:64127‑M1/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGTCACTTCCTGATCGTCCGCTCGCTGTATCTCTGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCGACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCTCGCCAACGCTCTCTGAAAACGTCTCCGGCTCCCGCTATAGCTGGAAAGTCAGCGGCATGGACTGCGCCGCCTGTG  >  REL606/3537030‑3537325

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