New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 3367495 =25 (0.470)4 (0.080) 4/60 NT 12.5% intergenic (‑81/+26) yhdN/rplQ hypothetical protein/50S ribosomal protein L17
?REL606 3725714 = 31 (0.500)coding (68/432 nt) yibN predicted rhodanese‑related sulfurtransferase

GTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCA  >  REL606/3367461‑3367530
                                  |                                   
gTGGGGATACTACGGGTATAAAAAAACCCGCCgggg                                    <  1:1823165/36‑1 (MQ=255)
ccatactgagtatcgcctggatcgcgttactggtgg                                    <  1:7462574‑M1/2‑1 (MQ=255)
ccatactgagtatcgcctggatcgcgttactggtgg                                    >  1:1040412‑M1/35‑36 (MQ=255)
ccatactgagtatcgcctggatcgcgttactggtgg                                    <  1:2324755‑M1/2‑1 (MQ=255)
 tGGTGATACTACGGGTATAAAAAAACCCGCCGGGGc                                   <  1:2367754/36‑1 (MQ=255)
 tGGGGATACTACGGGTATAAAAAAACCCGCCGGGGc                                   <  1:5215868/36‑1 (MQ=255)
 tGGGGATACTACGGGTATAAAAAAACCCGCCGGGGc                                   <  1:4691228/36‑1 (MQ=255)
 catactgagtatcgcctggatcgcgttactggtGGc                                   <  1:4635656‑M1/3‑1 (MQ=255)
 catactgagtatcgcctggatcgcgttactggtGGc                                   <  1:3596840‑M1/3‑1 (MQ=255)
  ggggATACTACGGGTATAAAAAAACCCGCCGGGGCg                                  >  1:4867671/1‑36 (MQ=255)
  ggggATACTACGGGTATAAAAAAACCCGCCGGGGCg                                  >  1:2112534/1‑36 (MQ=255)
  atactgagtatcgcctggatcgcgttactggtGGCg                                  <  1:6481352‑M1/4‑1 (MQ=255)
  atactgagtatcgcctggatcgcgttactggtGGCg                                  <  1:5524367‑M1/4‑1 (MQ=255)
  atactgagtatcgcctggatcgcgttactggtGGCg                                  <  1:424769‑M1/4‑1 (MQ=255)
  atactgagtatcgcctggatcgcgttactggtGGCg                                  <  1:3373305‑M1/4‑1 (MQ=255)
  atactgagtatcgcctggatcgcgttactggtGGCg                                  >  1:1037231‑M1/33‑36 (MQ=255)
     ctgagtatcgcctggatcgcgttactggtGGCGGGt                               <  1:3255767‑M1/7‑1 (MQ=255)
      aTACTACGGGTATAAAAAAACCCGCCGGGGCGGGtt                              >  1:4359202/1‑36 (MQ=255)
         ctacGGGTATAAAAAAACCCGCCGGGGCGGGttttt                           <  1:1184654/36‑1 (MQ=255)
          tatcgcctggatcgcgttactggtGGCGGGtctttt                          <  1:6864921‑M1/12‑1 (MQ=255)
           atcgcctggatcgcgttactggtGGCGGGttttttt                         <  1:4581517‑M1/13‑1 (MQ=255)
           acGGGTATAAAAAAACCCGCCGGGGCGGGttttttt                         >  1:2609531/1‑36 (MQ=255)
             gggTAAAAAAAAACCCGCCGGGGCGGGTTTTTTTAc                       <  1:354488/36‑1 (MQ=255)
             cgcctggatcgcgttactggtGGCGGGTTTTTTTAc                       <  1:8470611‑M1/15‑1 (MQ=255)
              ggTATAAAAAACCCCGCCGGGGCGGGTTTTTTTACg                      <  1:2252943/36‑1 (MQ=255)
              ggTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACg                      <  1:4198205/36‑1 (MQ=255)
               gTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGt                     <  1:832421/36‑1 (MQ=255)
                tataAAAAAACCCGCCGGGGCGGGTTTTTTTACGtt                    >  1:2059014/1‑36 (MQ=255)
                tataAAAAAACCCGCCGGGGCGGGTTTTTTTACGtt                    >  1:8448454/1‑36 (MQ=255)
                  taAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGc                  >  1:3845039/1‑36 (MQ=255)
                  taAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGc                  <  1:5718283/36‑1 (MQ=255)
                    aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCtt                >  1:8234416/1‑36 (MQ=255)
                    aaaaaaCCCGCCGGGGCGGGTTTTTTTACGGTGCtt                >  1:2183664/1‑36 (MQ=255)
                     aacaaCCCGCCGGGGCGGGTTTTTTTACGTTGCGTc               <  1:8353500/33‑1 (MQ=25)
                      tctaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCa              <  1:4531767/33‑1 (MQ=255)
                       aaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAg             >  1:2461380/1‑36 (MQ=255)
                         aCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGAt           >  1:2963621/1‑36 (MQ=255)
                         aCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGAt           >  1:4503191/1‑36 (MQ=255)
                           ccgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTa         >  1:7043713/1‑36 (MQ=255)
                           ccgccgGGGCGGGGTTTTTTACGCTGATTCAGATTa         >  1:7894392/1‑36 (MQ=18)
                            cgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTAc        >  1:7706474/1‑36 (MQ=255)
                            cgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTAc        <  1:1613797/36‑1 (MQ=255)
                              ccgGGGCGGGTTTTTTTACGTTTCTTCAGATTActc      >  1:6442380/1‑36 (MQ=255)
                                 gggCGGGTTTTTTTACGTTGCTTCAGATTACTCTgc   <  1:2366990/36‑1 (MQ=255)
                                 gggCGGGTTATTTTACGTTGCTTCAGATTACTCTgc   <  1:3560874/36‑1 (MQ=255)
                                  GGCGGGTTTTTTTACGTTTCTTCAGATTACTCTgca  <  1:1938599‑M1/36‑1 (MQ=255)
                                  GGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgca  <  1:680412‑M1/36‑1 (MQ=255)
                                  GGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgca  <  1:1519910‑M1/36‑1 (MQ=255)
                                  |                                   
GTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCA  >  REL606/3367461‑3367530

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG