New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 417727 =31 (0.650)6 (0.130) 5/54 NT 19.7% coding (682/891 nt) cyoE protoheme IX farnesyltransferase
?REL606 1333598 = 18 (0.290)pseudogene (67/129 nt) yciX predicted protein; b4522(b4574)

GCCGCGGCGACCACCAGATATTTATATCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAA  >  REL606/417693‑417762
                                  |                                   
gccgcGGCGACCACCAGATATTTATATCCAGCGTaa                                    >  1:618876/1‑36 (MQ=255)
gccgcGGCGACCACCAGATATTTATATCCAGCGTaa                                    <  1:4289897/36‑1 (MQ=255)
cgacaaccacccagaaaagtgcagaaccaagaaaaa                                    >  1:4859928‑M1/35‑36 (MQ=255)
 gacaaccacccagaaaagtgcagaaccaagaaaAAc                                   <  1:7013072‑M1/3‑1 (MQ=255)
 gacaaccacccagaaaagtgcagaaccaagaaaAAc                                   <  1:4703759‑M1/3‑1 (MQ=255)
 gacaaccacccagaaaagtgcagaaccaagaaaAAc                                   <  1:4622036‑M1/3‑1 (MQ=255)
 ccgcGGCGACCACCAGAACCTTATATCCAGCGTAAc                                   >  1:4559277/1‑36 (MQ=37)
  cgcgGCGACCACCAGATATTTATATCCAGCGTAAcc                                  <  1:1182143/36‑1 (MQ=255)
   gcggcgACCACCAGATATTTATATCCAGCGTAACCg                                 >  1:7677981/1‑36 (MQ=255)
     ggcgACCACCAGATATTTATATCCAGCGTAACCGcc                               >  1:7854846/1‑36 (MQ=255)
     ggcgACCACCAGATATTTATATCCAGCGTAACCGcc                               >  1:7106469/1‑36 (MQ=255)
      gcgACCACCAGATATTTATATCCAGCGTAACCGCCa                              >  1:2561812/1‑36 (MQ=255)
         accaccAGATATTTATATCCAGCGTAACCGCCAaga                           >  1:2706311/1‑36 (MQ=255)
           caccaGATATTTATATCCAGCGTAACCGCCAagaga                         >  1:5597676/1‑36 (MQ=255)
           caccaGATATTTATATCCAGCGTAACCGCCAagaga                         >  1:6628672/1‑36 (MQ=255)
             ccaGATATTTATATCCAGCGTAACCGCCAagagaga                       >  1:2496114/1‑36 (MQ=255)
              caGATATTTATATCCAGCGTAACCGCCAagagagag                      >  1:5235450/1‑36 (MQ=255)
              caGATATTTATATCCAGCGTAACCGCCAagagagag                      >  1:2053921/1‑36 (MQ=255)
              aaaatatTAATACCAAGCGAAACCGCCAagagagag                      <  1:3360113/33‑1 (MQ=12)
               aGATATTTATATCCAGCGTAACCGCCAAGAGAGAGc                     >  1:3458162/1‑36 (MQ=255)
               aGATATTTATATCCAGCGTAACCGCCAAGAGAGAGc                     >  1:3978246/1‑36 (MQ=255)
                gATATTTATATCCAGCGTAACCGCCAAGAGAGAGCa                    >  1:6139917/1‑36 (MQ=255)
                gATATTTATATCCAACGTATCCCCCAATATCGAGCa                    >  1:547071/1‑36 (MQ=255)
                  tatTTATATCCAGCGTAACCGCCAAGAGAGAGCATc                  >  1:5827523/1‑36 (MQ=255)
                   atTTATATCCAGCGTAACCGCCAAGAGAGAGCATCa                 >  1:7400120/1‑36 (MQ=255)
                   atTTATATCCAGCGTAACCGCCAAGAGAGAGCATCa                 >  1:8162514/1‑36 (MQ=255)
                    tttATATCCAGCGTAACCGCCAAGAGAGAGCATCAg                >  1:7865170/1‑36 (MQ=255)
                    tttATATCCAGCGTAACCGCCAAGAGAGAGCATCAg                >  1:5411143/1‑36 (MQ=255)
                     ttATATCCAGCGTAACCGCCAAGAGAGAGCATCAGc               >  1:2194868/1‑36 (MQ=255)
                      tataTCCAGCGTAACCGCCAAGAGAGAGCATCAGCg              >  1:8036925/1‑36 (MQ=255)
                        tatCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTg            >  1:2478136/1‑36 (MQ=255)
                        aaccaagaaaAACCGCCAAGAGAGAGCATCAGCGTg            >  1:2070191‑M1/11‑36 (MQ=255)
                         atCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTgg           >  1:7336873/1‑36 (MQ=255)
                          tCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGc          >  1:3913360/1‑36 (MQ=255)
                          cccagaaaAACCGCCAAGAGAGAGCATCAGCGTGGc          >  1:2116308‑M1/9‑36 (MQ=255)
                           ccAGCGTAACCGCCACGAGAGAGCATCAGCGTGGCa         >  1:2678777/1‑36 (MQ=255)
                           ccAGAGAAACCGCCAAGAGAGAGCATCAGCGTGGCa         <  1:3103202/36‑1 (MQ=37)
                           caaGCGAAACCGCCAAGAGAGAGCATCAGCGTGGCa         <  1:2057384/34‑1 (MQ=37)
                            cAGCGAAACCGCCAAGAGAGAGCATCAGCGTGGCaa        <  1:2679131/36‑1 (MQ=255)
                             agaaaAACCGCCAAGAGAGAGCATCAGCGTGGCAAc       >  1:3143853‑M1/6‑36 (MQ=255)
                             aGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAAc       >  1:4407536/1‑36 (MQ=255)
                             aGCGAAACCGCCAAGAGAGAGCATCAGCGTGGCAAc       <  1:5659639/36‑1 (MQ=255)
                              gCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACg      >  1:8015183/1‑36 (MQ=255)
                              aaaaAACCGCCAAGAGAGAGCATCAGCGTGGCAACg      >  1:3700259‑M1/5‑36 (MQ=255)
                               cGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACgg     >  1:6860720/1‑36 (MQ=255)
                               aaaAACCGCCAAGAGAGAGCATCAGCGTGGCAACgg     >  1:2850123‑M1/4‑36 (MQ=255)
                                gTAACCGCCAAGAGAGAGCATCAGCGTGGCAAAGGc    >  1:3645636/1‑36 (MQ=255)
                                 tAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCa   >  1:5066076/1‑36 (MQ=255)
                                 aAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCa   >  1:1364483‑M1/2‑36 (MQ=255)
                                  AACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCaa  >  1:7381852‑M1/1‑36 (MQ=255)
                                  |                                   
GCCGCGGCGACCACCAGATATTTATATCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAA  >  REL606/417693‑417762

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