Predicted mutation
evidence seq id position mutation annotation gene description
JC JC REL606 1,619,870 IS150 (+) +3 bp intergenic (‑107/+103) hokD ← / ← ECB_01533 small toxic polypeptide/conserved hypothetical protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 664688 =NA (NA)102 (1.830) 61/146 0.0 94.2% noncoding (1/1443 nt) IS150 repeat region
?REL606 = 1619872 6 (0.110)intergenic (‑109/+103) hokD/ECB_01533 small toxic polypeptide/conserved hypothetical protein
* ? REL606 = 666130NA (NA)74 (1.290) 43/150 0.3 92.2% noncoding (1443/1443 nt) IS150 repeat region
?REL606 1619870 = 6 (0.110)intergenic (‑107/+105) hokD/ECB_01533 small toxic polypeptide/conserved hypothetical protein

CTGCTGACCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAA  >  REL606/664616‑664761
                                                                        |                                                                         
ctgctgACCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTa                                                                        <  2:1370041/76‑1 (MQ=255)
 tgctgACCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTAc                                                                       <  1:1322698/76‑1 (MQ=255)
 tgctgACCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTAc                                                                       >  1:777364/1‑76 (MQ=255)
  gctgACCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACt                                                                      >  1:637573/1‑76 (MQ=255)
  gctgACCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACt                                                                      >  2:617821/1‑76 (MQ=255)
    tgACCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGc                                                                    >  2:1709623/1‑76 (MQ=255)
    tgACCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGc                                                                    <  2:464741/76‑1 (MQ=255)
      aCCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCAc                                                                  >  2:787029/1‑76 (MQ=255)
       ccgccgGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCAcc                                                                 >  2:965623/1‑76 (MQ=255)
        cgccgGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCAccc                                                                >  2:1694938/1‑76 (MQ=255)
        cgccgGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCAccc                                                                <  1:776807/76‑1 (MQ=255)
          ccgGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCAt                                                              <  2:1472124/76‑1 (MQ=255)
             gAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCAtttt                                                           <  2:1008206/76‑1 (MQ=255)
              agaGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTttg                                                          <  2:1069991/74‑1 (MQ=255)
              aaaGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTttg                                                          >  1:1713815/1‑76 (MQ=255)
              aaaGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTttg                                                          >  1:1602469/1‑76 (MQ=255)
                       ccGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGAt                                                 <  1:681092/76‑1 (MQ=255)
                       ccGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGAt                                                 <  1:1679772/76‑1 (MQ=255)
                       ccGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGAt                                                 <  2:1188331/76‑1 (MQ=255)
                        cGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATg                                                <  1:55872/76‑1 (MQ=255)
                         ggAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGa                                               <  2:688064/76‑1 (MQ=255)
                           aaaCTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGaaa                                             >  2:1309018/1‑76 (MQ=255)
                                tATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGaa                                        <  2:308727/76‑1 (MQ=255)
                                tATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGaa                                        >  1:1697036/1‑76 (MQ=255)
                                  ttGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCGCCCATTTTGTTGGACGATGAAATGGAATa                                      <  1:716500/76‑1 (MQ=255)
                                  ttGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa                                      <  2:1320467/76‑1 (MQ=255)
                                   tGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAg                                     <  2:185922/76‑1 (MQ=255)
                                   tGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAg                                     <  2:1734028/76‑1 (MQ=255)
                                   tGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAg                                     <  1:898412/76‑1 (MQ=255)
                                    ggCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGc                                    >  1:840699/1‑76 (MQ=255)
                                        tGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCt                                <  2:1289515/76‑1 (MQ=255)
                                        tGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCt                                <  2:1188418/76‑1 (MQ=255)
                                        tGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCt                                >  1:913774/1‑76 (MQ=255)
                                        tGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCt                                <  1:878025/76‑1 (MQ=255)
                                        tGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCt                                <  1:561718/76‑1 (MQ=255)
                                        tGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCt                                >  2:630597/1‑76 (MQ=255)
                                           gCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAat                             <  1:1211399/76‑1 (MQ=255)
                                            cAGAATGCACCACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAata                            <  2:233378/76‑1 (MQ=18)
                                            cAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGGCGATGAAATGGGATAGACACTAata                            >  2:1137209/1‑76 (MQ=21)
                                             aGAAAGCCCGACACGCTCCATCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAatat                           >  1:294340/1‑76 (MQ=255)
                                               aaaGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGt                         >  2:9783/1‑76 (MQ=255)
                                                aaGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGtn                        <  1:246977/76‑2 (MQ=255)
                                                    ccGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAg                    >  1:1263968/1‑76 (MQ=255)
                                                      gACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGcc                  >  1:697135/1‑76 (MQ=255)
                                                      gACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGcc                  <  1:352494/76‑1 (MQ=255)
                                                       acacGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCa                 <  2:23154/76‑1 (MQ=255)
                                                       acacGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCa                 <  1:1821663/76‑1 (MQ=255)
                                                        cacGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaa                <  2:1361776/76‑1 (MQ=255)
                                                        cacGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaa                <  1:90922/76‑1 (MQ=255)
                                                        cacGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaa                <  1:665150/76‑1 (MQ=255)
                                                            cTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATa            >  2:354102/1‑76 (MQ=255)
                                                             tCCAGCAGTTTTGTACTGCACCCATTTTGTTTGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAc           <  1:1026259/76‑1 (MQ=255)
                                                             tCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAc           <  1:1026260/76‑1 (MQ=255)
                                                              ccAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAcc          >  1:1470792/1‑76 (MQ=255)
                                                                agcagTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCt        >  2:1463321/1‑76 (MQ=40)
                                                                   agTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtttt     <  1:419855/76‑1 (MQ=35)
                                                                   agTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtttc     <  1:613054/76‑2 (MQ=25)
                                                                     ttttGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa   <  2:743/76‑1 (MQ=32)
                                                                     ttttGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa   >  2:903860/1‑76 (MQ=32)
                                                                     ttttGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa   >  2:1621167/1‑76 (MQ=32)
                                                                      tttGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGaa  <  1:787029/76‑1 (MQ=31)
                                                                      tttGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGaa  >  2:1030372/1‑76 (MQ=31)
                                                                        |                                                                         
CTGCTGACCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAA  >  REL606/664616‑664761

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