Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 3,990,387 IS150 (+) +3 bp 100% coding (1250‑1252/1491 nt) ECB_03719 ← conserved hypothetical protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 3652533NA (NA)124 (1.380) 59/76 NT 100% noncoding (1443/1443 nt) IS150 repeat region
?REL606 3990387 = 0 (0.000)coding (1252/1491 nt) ECB_03719 conserved hypothetical protein
* ? REL606 3893554 =NA (NA)97 (1.050) 60/78 NT 100% noncoding (1/1443 nt) IS150 repeat region
?REL606 = 3990389 0 (0.000)coding (1250/1491 nt) ECB_03719 conserved hypothetical protein

Only 100 of 283 total aligned reads displayed.

CCGGATATCCTCCTTTCTGCATGAAATAGCCCACCACGCCTGCAATCAGCATAATGCCAACATCAAAAAGGCTATTGTTCAACGCATAAGTCCCGACCACA  >  REL606/3990339‑3990439
                                                  |                                                  
gacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGCa                                                    >  1:4652918‑M2/49‑51 (MQ=255)
gacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGCa                                                    <  1:9794313‑M2/3‑1 (MQ=255)
gacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGCa                                                    >  1:10687928‑M2/49‑51 (MQ=255)
 cGGATATCCTCCTTTCTGCATGAAATAGCCCACCACGCCTGCAATCAGCAt                                                   <  1:2078416‑M2/51‑2 (MQ=255)
  cctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGCATa                                                  <  1:10498271‑M2/5‑1 (MQ=255)
  cctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGCATa                                                  <  1:1955380‑M2/5‑1 (MQ=255)
  cctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGCATa                                                  <  1:10407325‑M2/5‑1 (MQ=255)
   gATATCCTCCTTTCTGCATGAAATAGCCCACCACGCCTGCAATCAGCAtgt                                                 <  1:7257911‑M2/51‑4 (MQ=255)
    tatatgcctcgtgtttaactgtccaactttttggggtcagtacaGCATAAt                                                <  1:4376567‑M2/7‑1 (MQ=255)
    tatatgcctcgtgtttaactgtccaactttttggggtcagtacaGCATAAt                                                >  1:2621534‑M2/45‑51 (MQ=255)
     atatgcctcgtgtttaactgtccaactttttggggtcagtacaGCATAATg                                               <  1:9454791‑M2/8‑1 (MQ=255)
     atatgcctcgtgtttaactgtccaactttttggggtcagtacaGCATAATg                                               <  1:4905633‑M2/8‑1 (MQ=255)
       tcctccTTTCTGCATGAAATAGCCCACCACGCCTGCAATCAGCAtgtactg                                             >  1:5868087‑M2/1‑44 (MQ=255)
       tcctccTTTCTGCATGAAATAGCCCACCACGCCTGCAATCAGCAtgtactg                                             <  1:4784922‑M2/51‑8 (MQ=255)
        tgcctcgtgtttaactgtccaactttttggggtcagtacaGCATAATGCCa                                            <  1:9296270‑M2/11‑1 (MQ=255)
        tgcctcgtgtttaactgtccaactttttggggtcagtacaGCATAATGCCa                                            <  1:3367358‑M2/11‑1 (MQ=255)
           ctcgtgtttaactgtccaactttttggggtcagtacaGCATAATGCCAACa                                         >  1:8100884‑M2/38‑51 (MQ=255)
           ctcgtgtttaactgtccaactttttggggtcagtacaGCATAATGCCAACa                                         <  1:5704136‑M2/14‑1 (MQ=255)
           ctcgtgtttaactgtccaactttttggggtcagtacaGCATAATGCCAACa                                         <  1:6186414‑M2/14‑1 (MQ=255)
            tcgtgtttaactgtccaactttttggggtcagtacaGCATAATTCCAACAt                                        >  1:1658123‑M2/37‑51 (MQ=255)
            tcgtgtttaactgtccaactttttggggtcagtacaGCATAATGCCAACAt                                        >  1:7523424‑M2/37‑51 (MQ=255)
            ctTTCTGCATGAAATAGCCCACCACGCCTGCAATCAGCAtgtactgcaccc                                        <  1:4618310‑M2/51‑13 (MQ=255)
            ctTTCTGCATGAAATAGCCCACCACGCCTGCAATCAGCAtgtactgcaccc                                        <  1:7728659‑M2/51‑13 (MQ=255)
            ctTTCTGCATGAAATAGCCCACCACGCCTGCAATCAGCAtgtactgcaccc                                        >  1:3083727‑M2/1‑39 (MQ=255)
             tttCTGCATGAAATAGCCCACCACGCCTGCAATCAGCAtgtactgcaccca                                       >  1:10452833‑M2/1‑38 (MQ=255)
             tttCTGCATGAAATAGCCCACCACGCCTGCAATCAGCAtgtactgcaccca                                       >  1:5574906‑M2/1‑38 (MQ=255)
             tttCTGCATGAAATAGCCCACCACGCCTGCAATCAGCAtgtactgcaccca                                       <  1:2810329‑M2/51‑14 (MQ=255)
               tgtttaactgtccaactttttggggtcagtacaGCATAATGCCAACATCaa                                     <  1:10058576‑M2/18‑1 (MQ=255)
               tgtttaactgtccaactttttggggtcagtacaGCATAATGCCAACATCaa                                     >  1:8311465‑M2/34‑51 (MQ=255)
               tgtttaactgtccaactttttggggtcagtacaGCATAATGCCAACATCaa                                     >  1:10207341‑M2/34‑51 (MQ=255)
               tgtttaactgtccaactttttggggtcagtacaGCATAATGCCAACATCaa                                     >  1:3506074‑M2/34‑51 (MQ=255)
                 tttaactgtccaactttttggggtcagtacaGCATAATGCCAACATCaaaa                                   >  1:9238124‑M2/32‑51 (MQ=255)
                 tttaactgtccaactttttggggtcagtacaGCATAATGCCAACATCaaaa                                   >  1:2937364‑M2/32‑51 (MQ=255)
                 tGCATGAAATAGCCCACCACGCCTGCAATCAGCAtgtactgcacccatttt                                   <  1:5286775‑M2/51‑18 (MQ=255)
                 tGCATGAAATAGCCCACCACGCCTGCAATCAGCAtgtactgcacccatttt                                   <  1:3648997‑M2/51‑18 (MQ=255)
                  ttaactgtccaactttttggggtcagtacaGCATAATGCCAACATc                                       <  1:6477674‑M2/16‑1 (MQ=255)
                  ttaactgtccaactttttggggtcagtacaGCATAATGCCAACATc                                       <  1:5093132‑M2/16‑1 (MQ=255)
                  ttaactgtccaactttttggggtcagtacaGCATAATGCCAACATCaaaaa                                  >  1:3798885‑M2/31‑51 (MQ=255)
                   taactgtccaactttttggggtcagtacaGCATAATGCCAACATCAAAAAg                                 >  1:4665843‑M2/30‑51 (MQ=255)
                   cATGAAATAGCCCACCACGCCTGCAATCAGCAtgtactgcacccattttgt                                 >  1:3318481‑M2/1‑32 (MQ=255)
                    aactgtccaactttttggggtcagtacaGCATAATGCCAACATCAAAAAgg                                <  1:2371765‑M2/23‑1 (MQ=255)
                    aactgtccaactttttggggtcagtacaGCATAATGCCAACATCAAAAAgg                                <  1:10427048‑M2/23‑1 (MQ=255)
                     tGAAATAGCCCACCACGCCTGCAATCAGCAtgtactgcacccattttgttg                               <  1:3206234‑M2/51‑22 (MQ=255)
                      gAAATAGCCCACCACGCCTGCAATCAGCAtgtactgcacccattttgttgg                              >  1:7839584‑M2/1‑29 (MQ=255)
                      ctgtccaactttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCt                              >  1:7156114‑M2/27‑51 (MQ=255)
                      ctgtccaactttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCt                              >  1:5688408‑M2/27‑51 (MQ=255)
                      ctgtccaactttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCt                              >  1:7279987‑M2/27‑51 (MQ=255)
                       aaaTAGCCCACCACGCCTGCAATCAGCAtgtactgcacccattttgttgga                             <  1:8244302‑M2/51‑24 (MQ=255)
                        aaTAGCCCACCACGCCTGCAATCAGCAtgtactgcacccattttgttggac                            >  1:8955824‑M2/1‑27 (MQ=255)
                         tccaactttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTAtt                           <  1:4479395‑M2/28‑1 (MQ=255)
                         tccaactttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTAtt                           <  1:9429017‑M2/28‑1 (MQ=255)
                         tccaactttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTAtt                           >  1:2577430‑M2/24‑51 (MQ=255)
                         tccaactttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTAtt                           <  1:6011933‑M2/28‑1 (MQ=255)
                         tccaactttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTAtt                           >  1:5332622‑M2/24‑51 (MQ=255)
                         aTAGCCCACCACGCCTGCAATCAGCAtgtactgcacccattttgttggacg                           >  1:3261930‑M2/1‑26 (MQ=255)
                         aTAGCCCACCACGCCTGCAATCAGCAtgtactgcacccattttgttggacg                           >  1:1420102‑M2/1‑26 (MQ=255)
                          ccaactttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTg                          <  1:2310346‑M2/29‑1 (MQ=255)
                           aGCCCACCACGCCTGCAATCAGCAtgtactgcacccattttgttggacgat                         <  1:4262077‑M2/51‑28 (MQ=255)
                            gCCCACCACGCCTGCAATCAGCAtgtactgcacccattttgttggacgatg                        >  1:815618‑M2/1‑23 (MQ=255)
                            aactttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTGtt                        >  1:3979676‑M2/21‑51 (MQ=255)
                            aactttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTGtt                        >  1:3389217‑M2/21‑51 (MQ=255)
                             cccTCCACGCCTGCAATCAGCAtgtactgcacccattttgttggacgatga                       >  1:1130327‑M2/1‑22 (MQ=255)
                              ctttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTGTTNa                      <  1:1567544‑M2/33‑1 (MQ=255)
                              ctttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTGTTCa                      <  1:8012789‑M2/33‑1 (MQ=255)
                              ccaccaCGCCTGCAATCAGCAtgtactgcacccattttgttggacgatgaa                      >  1:3623385‑M2/1‑21 (MQ=255)
                               tttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTGTTCaa                     >  1:4413552‑M2/18‑51 (MQ=255)
                               tttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTGTTCaa                     <  1:2679845‑M2/34‑1 (MQ=255)
                                ttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTGTTCAAc                    <  1:9928304‑M2/35‑1 (MQ=255)
                                ttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTGTTCAAc                    >  1:1781966‑M2/17‑51 (MQ=255)
                                ttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTGTTCAAc                    <  1:7688580‑M2/35‑1 (MQ=255)
                                ttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTGTTCAAc                    <  1:9402906‑M2/35‑1 (MQ=255)
                                ttttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTGTTCAAc                    <  1:8607066‑M2/35‑1 (MQ=255)
                                accacGCCTGCAATCAGCAtgtactgcacccattttgttggacgatgaaat                    <  1:8538188‑M2/51‑33 (MQ=255)
                                 tttggggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTGTTCAACg                   <  1:10853249‑M2/36‑1 (MQ=255)
                                 ccacGCCTGCAATCAGCAtgtactgcacccattttgttggacgatgaaatg                   <  1:3902555‑M2/51‑34 (MQ=255)
                                  cacGCCTGCAATCAGCAtgtactgcacccattttgttggacgatgaaatgg                  >  1:6792771‑M2/1‑17 (MQ=255)
                                  cacGCCTGCAATCAGCAtgtactgcacccattttgttggacgatgaaatgg                  <  1:6715836‑M2/51‑35 (MQ=255)
                                   tggggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTGTTCAACGCa                 <  1:1248927‑M2/38‑1 (MQ=255)
                                   tggggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTGTTCAACGCa                 >  1:5010036‑M2/14‑51 (MQ=255)
                                    ggggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTGTTCAACGCAt                >  1:10329636‑M2/13‑51 (MQ=255)
                                     gCCTGCAATCAGCAtgtactgcacccattttgttggacgatgaaatggaat               >  1:6412519‑M2/1‑14 (MQ=255)
                                     gCCTGCAATCAGCAtgtactgcacccattttgttggacgatgaaatggaat               >  1:6899006‑M2/1‑14 (MQ=255)
                                     gCCTGCAATCAGCAtgtactgcacccattttgttggacgatgaaatggaat               <  1:4188369‑M2/51‑38 (MQ=255)
                                     gCCTGCAATCAGCAtgtactgcacccattttgttggacgatgaaatggaat               <  1:2118948‑M2/51‑38 (MQ=255)
                                     gCCTGCAATCAGCAtgtactgcacccattttgttggacgatgaaatggaat               <  1:84296‑M2/51‑38 (MQ=255)
                                     gCCTGCAATCAGCAtgtactgcacccattttgttggacgatgaaatggaat               >  1:8806137‑M2/1‑14 (MQ=255)
                                      ggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTGTTCAACGCATaa              >  1:2009809‑M2/11‑51 (MQ=255)
                                      ggtcagtacaGCATAATGCCAACATCAAAAAGGCTATTGTTCAACGCATaa              <  1:5741926‑M2/41‑1 (MQ=255)
                                          cAATCAGCAtgtactgcacccattttgttggacgatgaaatggaatagccc          <  1:9749259‑M2/51‑43 (MQ=255)
                                          cAATCAGCAtgtactgcacccattttgttggacgatgaaatggaatagccc          <  1:886978‑M2/51‑43 (MQ=255)
                                          agtacaGCATAATGCCAACATCAAAAAGGCTATTGTTCAACGCATAAGTcc          <  1:7056688‑M2/45‑1 (MQ=255)
                                          agtacaGCATAATGCCAACATCAAAAAGGCTATTGTTCAACGCATAAGTcc          <  1:5193868‑M2/45‑1 (MQ=255)
                                            tacaGCATAATGCCAACATCAAAAAGGCTATTGTTCAACGCATAAGTCCCg        <  1:9629546‑M2/47‑1 (MQ=255)
                                             acaGCATAATGCCAACATCAAAAAGGCTATTGTTCAACGCATAAGTCCCGa       <  1:387899‑M2/48‑1 (MQ=255)
                                              cAGCAtgtactgcacccattttgttggacgatgaaatggaatagcccctaa      <  1:684779‑M2/51‑47 (MQ=255)
                                              cAGCAtgtactgcacccattttgttggacgatgaaatggaatagcccctaa      >  1:10570552‑M2/1‑5 (MQ=255)
                                               aGCAtgtactgcacccattttgttggacgatgaaatggaatagcccctaat     >  1:5494800‑M2/1‑4 (MQ=255)
                                                gCATAATGCCAACATCAAAAAGGCTATTGTTCAACGCATAAGTCCCGACca    >  1:69904/1‑51 (MQ=255)
                                                gCATAATGCCAACATCAAAAAGGCTATTGTTCAACGCATAAGTCCCGACca    >  1:8781200/1‑51 (MQ=255)
                                                 cAtgtactgcacccattttgttggacgatgaaatggaatagcccctaatat   <  1:1221628‑M2/51‑50 (MQ=255)
                                                  |                                                  
CCGGATATCCTCCTTTCTGCATGAAATAGCCCACCACGCCTGCAATCAGCATAATGCCAACATCAAAAAGGCTATTGTTCAACGCATAAGTCCCGACCACA  >  REL606/3990339‑3990439

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