New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 230935240 (0.770)5 (0.100) 3/60 NT 13.7% coding (795/1203 nt) yfaY competence damage‑inducible protein A
?REL606 3944634 = 23 (0.380)intergenic (+23/‑124) proM/aslB tRNA‑Pro/predicted regulator of arylsulfatase activity

CTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTCCCT  >  REL606/3944599‑3944669
                                   |                                   
ctctcGCCGACCAATTTTGAACCCCGCTTCGGCggg                                     >  1:1500704/1‑36 (MQ=255)
ctctcGCCGACCAATTTTGAACCCCGCTTCGGCggg                                     <  1:206216/36‑1 (MQ=255)
ctctcGCCGACCAATTTTGAACCCCGCTTCGGCggg                                     <  1:2609635/36‑1 (MQ=255)
 gaattgctcgctcaacgtcaggctgaactggcgggt                                    >  1:5178745‑M2/35‑35 (MQ=255)
  ctcGCCGACCAATTTTGAACCCCGCTTCGGCGGGGt                                   <  1:6063098/36‑1 (MQ=255)
   tcGCCGACCAATTTTGAACCCCGCTTCGGCGGGGtt                                  >  1:7137647/1‑36 (MQ=255)
   tcGCCGACCAATTTTGAACCCCGCTTCGGCGGGGtt                                  <  1:2377332/36‑1 (MQ=255)
   tcGCCGACCAATTTTGAACCCCGCTTCGGCGGGGtt                                  >  1:3683381/1‑36 (MQ=255)
    cGCCGACCAATTTTGAACCCCGCTTCGGCGGGGttt                                 >  1:5385111/1‑36 (MQ=255)
      gctcgctcaacgtcaggctgaactggcggTGttttt                               >  1:423406‑M2/30‑36 (MQ=255)
       ctcgctcaacgtcaggctgaactggcggGGtttttt                              >  1:2032517‑M2/29‑36 (MQ=255)
       ctcgctcaacgtcaggctgaactggcggGCtttttt                              >  1:6613612‑M2/29‑36 (MQ=255)
       cGACCAATTTTGAACCCCGCTTCGGCGGGGtttttt                              >  1:198866/1‑36 (MQ=255)
            tcaacgtcaggctgaactggcggGGTTTTTTTtttt                         >  1:3053623‑M2/24‑36 (MQ=255)
            aaTTTTGAACCCCGCTTCGGCGGGGTTTTTTGtttt                         <  1:6552640/36‑1 (MQ=255)
            aaTTTTGAACCCCGCTTCGGCGGGGTTTTTTGtttt                         >  1:527621/1‑36 (MQ=255)
            aaTTTTGAACCCCGCTTCGGCGGGGTTTTTTGtttt                         >  1:2947712/1‑36 (MQ=255)
             aTTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTc                        <  1:4441040/36‑1 (MQ=255)
             aTTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTc                        <  1:451560/36‑1 (MQ=255)
                ttGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCtgt                     <  1:3839218/36‑1 (MQ=255)
                ttGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCtgt                     <  1:2791341/36‑1 (MQ=255)
                 tGAACCCCGCTTCGGCGGGGGTTTTTGTTTTCtgtg                    >  1:6252440/1‑36 (MQ=255)
                  gAACCCCGCTGCGGCGGGGGTTTTTTTGTTCTGTGc                   >  1:7646440/1‑36 (MQ=9)
                      cccGCTTCGGCGGGGTTTTTTGTTTTCTGTGCAttt               >  1:342072/1‑36 (MQ=255)
                      cccGCTTCGGCGGGGTTTTTTGTTTTCTGTGCAttt               >  1:284088/1‑36 (MQ=255)
                      cccGCTTCGGCGGGGTTTTTTGTTTTCTGTGCAttt               >  1:5715029/1‑36 (MQ=255)
                       ccGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTc              <  1:5824988/36‑1 (MQ=255)
                        cGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCg             >  1:6272430/1‑36 (MQ=255)
                        cGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCg             >  1:2275914/1‑36 (MQ=255)
                          cTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTc           <  1:4148155/36‑1 (MQ=255)
                          cTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTc           <  1:5023388/36‑1 (MQ=255)
                                cggGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTc     >  1:5712339‑M2/4‑36 (MQ=255)
                                cggGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTc     >  1:4387668‑M2/4‑36 (MQ=255)
                                 ggGGTTTTTTGTTTTCTGTGCATTTCGTCCCTTTcc    <  1:7500632‑M2/34‑1 (MQ=255)
                                 ggGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTcc    >  1:6394775‑M2/3‑36 (MQ=255)
                                  gGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTccc   >  1:905893‑M2/2‑36 (MQ=255)
                                  gGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTccc   >  1:595873‑M2/2‑36 (MQ=255)
                                  gGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTccc   <  1:2957668‑M2/35‑1 (MQ=255)
                                   ggTTTTTTGTTTTCTGTGCATTTCGTCACTTTCCCt  >  1:3592842/1‑36 (MQ=255)
                                   ggTTTTTTGTTTTCTGTGCATTTCGTCACTTTCCCt  >  1:111475/1‑36 (MQ=255)
                                   |                                   
CTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTCCCT  >  REL606/3944599‑3944669

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