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)

ATCCTGAATTAAAACGTAAAGCCTGGCTGGCGGTTTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTCG  >  REL606/1333563‑1333631
                                   |                                 
tgccgttgccacgctgatgctctctcttggcggttt                                   >  1:1364483‑M2/36‑36 (MQ=255)
aTCCTGAATTAAAACGTAAAGCCTGGCTGGCGGttt                                   <  1:1760490/36‑1 (MQ=255)
 tCCTGAATTAAGACGTAAAGCCTGGCTGGCGGtttt                                  <  1:5471872/36‑1 (MQ=255)
  ccgttgccacgctgatgctctctcttggcggttttt                                 >  1:2850123‑M2/34‑36 (MQ=255)
  ccTGAATTAAAACGTAAAGCCTGGCTGGCGGttttt                                 <  1:3396901/36‑1 (MQ=255)
   cgttgccacgctgatgctctctcttggcggtttttt                                >  1:3700259‑M2/33‑35 (MQ=255)
   cTGAATTAAAACGTAAAGCCTGGCTGGCGGttgttt                                >  1:3530103/1‑35 (MQ=255)
   cTGAATTAAAACGTAAAGCCTGGCTGGCGGTTTTTc                                <  1:615555/36‑1 (MQ=255)
    gttgccacgctgatgctctctcttggcggttTTTCt                               >  1:3143853‑M2/32‑36 (MQ=255)
     gAATTAAAACGTAAAGCCTGGCTGGCGGTTTTTTtt                              >  1:8346573/1‑36 (MQ=255)
     gAATTAAAACGTAAAGCCTGGCTGGCGGGTTTTTtt                              >  1:2623050/1‑36 (MQ=37)
       gccacgctgatgctctctcttggcggttTTTCTGgg                            >  1:2116308‑M2/29‑36 (MQ=255)
       aTTAAAACGTAAAGCCTGGCTGGCGGTTTTTCtttg                            >  1:338296/1‑34 (MQ=255)
       aTTAAAACGTAAAGCCTGGCTGGCGGTTTTTCTTgg                            <  1:1974142/36‑1 (MQ=255)
       aTTAAAACGTAAAGCCTGGCTGGCGGTTTTTCTTgg                            <  1:4703686/36‑1 (MQ=255)
       aTTAAAACGTAAAGCCTGGCTGGCGGGTTTTTtttg                            >  1:2438515/1‑34 (MQ=25)
         tAAAACGTAAAGCCTGGCTGGCGGTTTTTCTTGGtt                          <  1:661092/36‑1 (MQ=255)
         tAAAACGTAAAGCCTGGCTGGCGGTTTTTCTTGGtt                          <  1:7418670/36‑1 (MQ=255)
         cacgctgatgctctctcttggcggttTTTCTTGGtt                          >  1:2070191‑M2/27‑36 (MQ=255)
            aaCGTAAAGCCTGGCTGGCGGTTTTTCTTGGTTCTg                       <  1:2422911/36‑1 (MQ=255)
             aCGTAAAGCCTGGCTGGCGGTTTTTCTTGGTTCTGc                      <  1:3500995/36‑1 (MQ=255)
                tacaGCCTGGCTGGCGGTTTTTCTTGGTTCTGCACt                   <  1:4758529/33‑1 (MQ=255)
                  aaGCCTGGCTGGCGGTTTTTTTTTGTTTTTCACttt                 >  1:3035067/1‑36 (MQ=16)
                   aGCCTGGCTGGCGGTTTTTTTTGGTTTTTCATtttt                >  1:2868779/1‑36 (MQ=16)
                   aGCCTGGCTGGCGGTTTTTCTTGGTTCTGCACtttt                <  1:3900372/36‑1 (MQ=255)
                   aGCCTGGCTGGCGGTTTTTCTTGGTTCTGCACtttt                <  1:5529101/36‑1 (MQ=255)
                    gCCTGGCTGGCGGTTTTTTTTTTTTTTTCACttttt               >  1:5807925/1‑35 (MQ=9)
                    gCCTGGCTGGCGGTTTTTCTTTGTTCTTCACTTTTc               >  1:6069063/1‑36 (MQ=37)
                    gCCTGGCTGGCGGTTTTTCTTTGTTCTGCACttttt               >  1:3495707/1‑35 (MQ=255)
                    gCCTGGCTGGCGGTTTTTCTTGGTTCTGCACTTTTc               <  1:7178966/36‑1 (MQ=255)
                    gCCTGGCTGGCGGTTTTTCTTGGTTCTGCACTTTTc               <  1:2204397/36‑1 (MQ=255)
                     ccTGGCTGGCGGTTTTTCTTGGTTCTGCACTTTTCt              <  1:5188285/36‑1 (MQ=255)
                      ctggctggCGGTTTTTCTTTGTTCTGCACctttttt             >  1:2762833/1‑33 (MQ=25)
                        ggctggcGGTTTTTCTTGGTTTTGCACTTTTCTggg           >  1:2899412/1‑36 (MQ=255)
                            ggcGGTTTTTCTTGGGTTTGCACTTTTCTTggtggt       >  1:284962/1‑36 (MQ=25)
                                gttTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTc   <  1:7013072‑M2/33‑1 (MQ=255)
                                gttTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTc   <  1:4703759‑M2/33‑1 (MQ=255)
                                gttTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTc   <  1:4622036‑M2/33‑1 (MQ=255)
                                 ttTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTCg  >  1:4859928‑M2/3‑36 (MQ=255)
                                   |                                 
ATCCTGAATTAAAACGTAAAGCCTGGCTGGCGGTTTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTCG  >  REL606/1333563‑1333631

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 19 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG