New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 885072 =27 (0.510)9 (0.170) 5/60 NT 25.7% coding (35/858 nt) ECB_00824 protein similar to retron EC67 DNA adenine methylase
?REL606 3944634 = 25 (0.400)intergenic (+23/‑124) proM/aslB tRNA‑Pro/predicted regulator of arylsulfatase activity

TATGAGCACCATCCTGAAATGGGCGGGTAATAAAACCGCCATTATGTCCGAACTGAAAAAATACCTTCCTG  >  REL606/885037‑885107
                                   |                                   
tATGAGCACCATCCTGAAATGGGCGGGTAATAAAAc                                     >  1:470106/1‑36 (MQ=255)
gggaaagtgacgaaatgcacagaaaacaaaaaaccc                                     <  1:439797‑M1/1‑1 (MQ=255)
 ggaaagtgacgaaatgcacagaaaacaaaaaacccc                                    >  1:1337399‑M1/35‑36 (MQ=255)
 ggaaagtgacgaaatgcacagaaaacaaaaaacccc                                    >  1:1509820‑M1/35‑36 (MQ=255)
 ggaaagtgacgaaatgcacagaaaacaaaaaacccc                                    >  1:6203280‑M1/35‑36 (MQ=255)
 ggaaagtgacgaaatgcacagaaaacaaaaaacccc                                    >  1:6144908‑M1/35‑36 (MQ=255)
 ggaaagtgacgaaatgcacagaaaacaaaaaacccc                                    >  1:4390377‑M1/35‑36 (MQ=255)
  tGAGCACCATCCTGAAATGGGCGGGTAATAAAACCg                                   <  1:1475381/36‑1 (MQ=255)
  tGAGCACCATCCTGAAATGGGCGGGTAATAAAACCg                                   >  1:2082691/1‑36 (MQ=255)
  gaaagtgacgaaatgcacagaaaacaaaaaaccCCg                                   >  1:7702630‑M1/34‑36 (MQ=255)
  gaaagtgaagaaatgcacagaaaacaaaaaaccCCg                                   >  1:6882411‑M1/34‑36 (MQ=255)
   gAGCACCATCCTGAAATGGGCGGGTAATAAAACCGc                                  >  1:5516644/1‑36 (MQ=255)
   gAGCACCATCCTGAAATGGGCGGGTAATAAAACCGc                                  >  1:2208095/1‑36 (MQ=255)
   gAGCACCATCCTGAAATGGGCGGGTAATAAAACCGc                                  <  1:29680/36‑1 (MQ=255)
   aaagtgacgaaatgcacagaaaacaaaaaaccCCGc                                  >  1:1088977‑M1/33‑36 (MQ=255)
      cACCATCCTGAAATGGGCGGGTAATAAAACCGCCAt                               >  1:7146196/1‑36 (MQ=255)
        ccATCCTGAAATGGGCGGGTAATAAAACCGCCATTa                             >  1:310378/1‑36 (MQ=255)
           tCCTGAAATGGGCGGGTAATAAAACCGCCATTATGt                          >  1:453186/1‑36 (MQ=255)
           tCCTGAAATGGGCGGGTAATAAAACCGCCATTATGt                          >  1:2152037/1‑36 (MQ=255)
            ccagAAAGGGGCGAAAAAAAAACCCGCCATTATGTc                         <  1:1938157/33‑1 (MQ=0)
             cTGAAATGGGCGGGTAATAAAACCGCCATTATGTcc                        >  1:5943847/1‑36 (MQ=255)
                aaaTGGGCGGGTAATAAAACCGCCATTTGGTCCgta                     >  1:6204865/1‑34 (MQ=37)
                 aaTGGGCGGGTAATAAAACCGCCATTATGTCCGAAc                    >  1:493/1‑36 (MQ=255)
                  aTGGGCGGGTAATAAAACCGCCATTCNNNCCGAACt                   >  1:2789407/1‑36 (MQ=37)
                  aTGGGCGGGTAATAAAACCGCCATTATGTCCGAACt                   >  1:4842355/1‑36 (MQ=255)
                    ggaaaaaaaaaaaccCCGCCATTATGTCCGAACTGa                 >  1:3772690‑M1/16‑36 (MQ=255)
                     ggCGGGTAATAAAACCGCCATTATGTCCGAACTGaa                >  1:3223445/1‑36 (MQ=255)
                     ggCGGGTAATAAAACCGCCATTATGTCCGAACTGaa                >  1:2821833/1‑36 (MQ=255)
                      gCGGGTAATAAAACCGCCATTATGTCCGAACTGaaa               >  1:2281646/1‑36 (MQ=255)
                       cGGGTAATAAAACCGCCATTATGTCCGAACTGaaaa              >  1:6115806/1‑36 (MQ=255)
                        gggTAATAAAACCGCCATTATGTCCGAACTGaaaaa             >  1:2682927/1‑36 (MQ=255)
                        aaaaaaaaaccCCGCCATTATGTCCGAACTGaaaaa             >  1:5484827‑M1/12‑36 (MQ=255)
                          gTAATAAAACCGCCATTATGTCCGAACTGAAAAAAt           >  1:3555398/1‑36 (MQ=255)
                           aaaaaacaCCGCCATTATGTCCGAACTGAAAAAATa          >  1:2067614‑M1/9‑36 (MQ=255)
                            aataaAACCGCCATTATGTCCGAACTGAAAAAATAc         >  1:4670815/1‑36 (MQ=255)
                            aataaAACCGCCATTATGTCCGAACTGAAAAAATAc         >  1:6548365/1‑36 (MQ=255)
                            aataaAACCGCCATGATGGCCGAACTGAAAAAATAc         >  1:7079913/1‑36 (MQ=37)
                            aaaaaacCCGCCATTATGTCCGAACTGAAAAAATAc         >  1:4566477‑M1/8‑36 (MQ=255)
                            aaaaaacCCGCCATTATGTCCGAACTGAAAAAATAc         >  1:2597505‑M1/8‑36 (MQ=255)
                             ataaAACCGCCATTATGTCCGAACTGAAAAAATAcc        >  1:414145/1‑36 (MQ=255)
                              taaAACCGCCATTATGTCCGAACTGAAAAAATACCt       >  1:4549977/1‑36 (MQ=255)
                               aatcCCGCCATTATGTCCGAACTGAAAAAATACCtt      >  1:6032381‑M1/5‑36 (MQ=255)
                               aaccCCGCCATTATGTCCGAACTGAAAAAATACCtt      >  1:624798‑M1/5‑36 (MQ=255)
                               aaccCCGCCATTATGTCCGAACTGAAAAAATACCtt      >  1:7409053‑M1/5‑36 (MQ=255)
                                aaaCCGCCATTATGTCCGAACTGAAAAAATACCTTc     >  1:666135/1‑36 (MQ=255)
                                  cCCGCCATTATGTCCGAACTGAAAAAATACCTTCct   >  1:7832785‑M1/2‑36 (MQ=255)
                                   CCGCCATTATGTCCGAACTGAAAAAATACCTTCctg  >  1:6282527‑M1/1‑36 (MQ=255)
                                   |                                   
TATGAGCACCATCCTGAAATGGGCGGGTAATAAAACCGCCATTATGTCCGAACTGAAAAAATACCTTCCTG  >  REL606/885037‑885107

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG