New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 3995530 =46 (0.870)10 (0.190) 4/60 NT 17.7% coding (58/768 nt) ECB_03724 putative transcriptional regulator
?REL606 = 4141016 47 (0.760)coding (89/705 nt) yijC DNA‑binding transcriptional repressor

GTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAT  >  REL606/4140982‑4141051
                                  |                                   
gTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACc                                    <  1:4787648‑M2/36‑2 (MQ=255)
gTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACc                                    >  1:1767827‑M2/1‑35 (MQ=255)
gTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACc                                    >  1:7341501‑M2/1‑35 (MQ=255)
 tGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACcc                                   <  1:656620‑M2/36‑3 (MQ=255)
 tGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACcc                                   <  1:6852849‑M2/36‑3 (MQ=255)
 tGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACcc                                   <  1:2923902‑M2/36‑3 (MQ=255)
    tGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgcc                                <  1:5571892‑M2/36‑6 (MQ=255)
    tGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgcc                                >  1:7433392‑M2/1‑31 (MQ=255)
        cGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTc                            >  1:819558/1‑36 (MQ=255)
        agTAAGAGCAAAACAAAAAGAAAAAACCCGCCGTTc                            <  1:7841874/35‑1 (MQ=25)
         gTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCg                           >  1:3447303/1‑36 (MQ=255)
           aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCt                         >  1:1930025/1‑36 (MQ=255)
           aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCt                         >  1:8334920/1‑36 (MQ=255)
            agagCGAAAAAAAAAGAAAAACCCCGCCGTTCGCtg                        <  1:6090784/36‑1 (MQ=25)
             gagCGCAACAAAAAGAAAAAACCCGCCGTTCGCtgg                       >  1:4956375/1‑36 (MQ=255)
              agCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggt                      >  1:6638661/1‑36 (MQ=255)
              agCGAAAAAAAAAAAAAAAACCCGCCGTTCGCtggt                      <  1:5434703/36‑1 (MQ=18)
               gcgcAACAAAAAGAAAAAACCCGCCGTTCGCtggtg                     >  1:621075/1‑36 (MQ=255)
                cgcAACAAAAAGAAAAAACCCGCCGTTCGCtggtgg                    >  1:809768/1‑36 (MQ=255)
                cgcAACAAAAAGAAAAAACCCGCCGTTCGCtggtgg                    >  1:2747246/1‑36 (MQ=255)
                 gaaacaaAAAGAAAAACCCCGCCGTTCGCTGGTGGa                   <  1:873195/34‑1 (MQ=38)
                  caacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGaa                  >  1:8032097/1‑36 (MQ=255)
                   aaaaaaaaGAAAAAACCCGCCGTTCGCTGGTGGAAg                 <  1:1817160/33‑1 (MQ=255)
                    acaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGc                >  1:2777313/1‑36 (MQ=255)
                    acaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGc                >  1:4757174/1‑36 (MQ=255)
                    acaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGc                >  1:3377767/1‑36 (MQ=255)
                    aaaaaaaGAAAAAACCCGCCGTTCGCTGGTGGAAGc                <  1:849823/34‑1 (MQ=255)
                    aaaaaaaAAAAAAACCCGCCGTTCGCTGGTGGAAGc                <  1:8007084/34‑1 (MQ=37)
                    aaaaaaaAAAAAAACCCGCCGTTCGCTGGTGGAAGc                <  1:5901898/34‑1 (MQ=37)
                     caaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGcc               >  1:5843321/1‑36 (MQ=255)
                     caaAAAGAAAAAACCCGCCGCTCGCTGGTGGAAGcc               >  1:8301022/1‑36 (MQ=255)
                     caaAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGcc               <  1:5913146/36‑1 (MQ=255)
                     aaaaaaGAAAAACCCCGCCGTTCGCTGGTGGAAGcc               <  1:3153254/35‑1 (MQ=255)
                     aaaaaaAAAAAAACCCGCCGTTCGCTGGTGGAAGcc               <  1:6589166/35‑1 (MQ=255)
                      aaaaaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCg              >  1:4411566/1‑36 (MQ=255)
                      aaaaaAAAAAAACccgccagtccacgtaataacgtt              >  1:6205237‑M2/1‑13 (MQ=255)
                      aaaaaAAAAAAACccgccagtccacgtaataacgtt              >  1:3203353‑M2/1‑13 (MQ=255)
                      aaaaaAAAAAAACCCGCCGTTCGCTGGTGGAAGCCg              <  1:2574093/36‑1 (MQ=255)
                       gaaaaAAAAAACccgccgaaaactgtaccgagaaca             >  1:1713578‑M2/2‑12 (MQ=255)
                       aaaaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGc             >  1:4498743/1‑36 (MQ=255)
                       aaaaGAAAAAACCCGACGTTCGCTGGTGGAAGCCGc             >  1:7363172/1‑36 (MQ=255)
                       aaaaAAAAAAACccgccgtaccagattaaggtattc             >  1:2278046‑M2/1‑12 (MQ=255)
                       aaaaAAAAAAACccgccgaaaactgtaccgagaaca             >  1:2798616‑M2/1‑12 (MQ=255)
                       aaaaAAAAAAACccgccagtccacgtaataacgttt             >  1:3633073‑M2/1‑12 (MQ=255)
                       aaaaAAAAAAACccgccagtccacgtaataacgttt             >  1:1408672‑M2/1‑12 (MQ=255)
                        aaaGAAAAAACCCGCCGTGCGCTGGTGGAAGCCGCa            >  1:2614074/1‑36 (MQ=255)
                        aaaAAAAAAACccgccgaaaactgtaccgagaacaa            >  1:4701667‑M2/1‑11 (MQ=255)
                         aaaaaaaacCCCGCCGTTCGCTGGTGGAAGCCGCAt           <  1:8430524/33‑1 (MQ=37)
                         aaaaaaaaaCCCGCCGTTCGCTGGTGGAAGCCGCAt           <  1:5356793/33‑1 (MQ=255)
                         aaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt           >  1:5964286/1‑36 (MQ=255)
                         aaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt           >  1:4585948/1‑36 (MQ=255)
                          ggaAAGAACccgcctctgcgctggctcatgcggtac          >  1:2980193‑M2/2‑9 (MQ=255)
                          ggaAAAAAAccgcctctgcgctggctcatgcggtac          >  1:1686106‑M2/2‑9 (MQ=255)
                          aGAAAAAACccgcctctgcgctggctcatgcggtac          >  1:3421879‑M2/1‑9 (MQ=255)
                           gAAAAAACccgcctctgcgctggctcatgcggtaca         >  1:2043722‑M2/1‑8 (MQ=255)
                           gAAAAAACCCGCCGTTCGCTGGGGGAAGCCGCAttt         >  1:3239008/1‑36 (MQ=255)
                           aaaaaaccccgccagtccacgtaataacgtttgtgc         >  1:6322589‑M2/2‑8 (MQ=255)
                           aaaaaacCCCGCCGTTCGCTGGTGGAAGCCGCAttt         <  1:2018795/35‑1 (MQ=255)
                           aaaaaaacccgccgaaaactgtaccgagaacaacca         >  1:6338915‑M2/2‑8 (MQ=255)
                           aaaaaaacccgccagtccacgtaataacgtttgtgc         >  1:31253‑M2/2‑8 (MQ=255)
                            aaaaaaGccgcctctgcgctggctcatgcggtacaa        >  1:7890023‑M2/1‑7 (MQ=255)
                            aaaaaaCccgccgtaccagattaaggtattccagcc        >  1:437479‑M2/1‑7 (MQ=255)
                            aaaaaaCCCGCCGTTCGCTGGTGGAAGCCGCATTTa        >  1:2606491/1‑36 (MQ=255)
                            aaaaaaCCCGCCGTTCGCTGGTGGAAGCCGCATTTa        >  1:440891/1‑36 (MQ=255)
                             aaacaCccgccagtccacgtaataacgtttgtgcgc       >  1:3196237‑M2/1‑6 (MQ=255)
                             aaaacCccgccgtaccagattaaggtattccagccg       >  1:7202282‑M2/1‑6 (MQ=255)
                             aaaacCCCGCCGTTCGCTGGTGGAAGCCGCATTTAg       <  1:7558666/36‑1 (MQ=255)
                             aaaaaCCCGCCGTTCGCTGGTGGAAGCCGCATTTAg       >  1:2434732/1‑36 (MQ=255)
                             aaaaaCCCGCCGTTCGCTGGTGGAAGCCGCATTTAg       <  1:4806931/36‑1 (MQ=255)
                             aaaaaCCCGCCGTGCGCTGGTGCAAGCAGCATTTAg       >  1:2091679/1‑36 (MQ=25)
                             aaaaaAccgcctctgcgctggctcatgcggtacaaa       >  1:3998721‑M2/1‑6 (MQ=255)
                             aaaaaAccgccagtccacgtaataacgtttgtgcgc       >  1:7673818‑M2/1‑6 (MQ=255)
                              aaaaCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGc      <  1:185560/36‑1 (MQ=255)
                              aaaaCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGc      >  1:2525575/1‑36 (MQ=255)
                              aaaaAccgccgaaaactgtaccgagaacaaccaggt      >  1:4500260‑M2/1‑5 (MQ=255)
                              aaaaAccgccgaaaactgtaccgagaacaaccaggt      >  1:6997068‑M2/1‑5 (MQ=255)
                               aaccccGCCGTTCGCTGGTGGAAGCCGCATTTAGcc     <  1:2121173/33‑1 (MQ=255)
                               aaaCccgcctctgcgctggctcatgcggtacaaagc     >  1:3501407‑M2/1‑4 (MQ=255)
                               aaaCccgccgaaaactgtaccgagaacaaccaggta     >  1:5485810‑M2/1‑4 (MQ=255)
                                accccgccgtaccagattaaggtattccagccgcct    >  1:4074373‑M2/3‑3 (MQ=255)
                                aaCccgcctctgcgctggctcatgcggtacaaagcg    >  1:2674592‑M2/1‑3 (MQ=255)
                                 aCCCGCCGTTCGCTGGTGGAAGCCGTATTTAGCCaa   >  1:6218303/1‑36 (MQ=255)
                                 aCCCGCCGTTCGCTGGTGGAAGCCACATTTAGCCaa   >  1:5012111/1‑36 (MQ=255)
                                  cccgccgtaccagattaaggtattccagccgccttt  >  1:2604152‑M2/1‑1 (MQ=255)
                                  cccgccgtaccagattaaggtattccagccgccttt  >  1:1237391‑M2/1‑1 (MQ=255)
                                  cccgccagtccacgtaataacgtttgtgcgcctgcg  >  1:3603617‑M2/1‑1 (MQ=255)
                                  cccgccagtccacgtaataacgtttgtgcgcctgcg  >  1:5758017‑M2/1‑1 (MQ=255)
                                  cccGCCGTTCGCTGGTGGAACCCGCATTTAGCCAAt  <  1:3209081/36‑1 (MQ=255)
                                  |                                   
GTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAT  >  REL606/4140982‑4141051

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG