New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1 =0 (0.000)39 (0.870) 31/70 NT 100% intergenic (–/‑189) –/thrL –/thr operon leader peptide
?REL606 = 4629812 0 (0.000)intergenic (+24/–) lasT/– predicted rRNA methyltransferase/–

GACATATTGCCCGTTGCAGTCAGAATGAAAAGCT...................................  <  REL606/34‑1
..................................GAAAAATACTTACTAAGGCGTTTTTTATTTGGTGA  <  REL606/4629812‑4629778
                                                                     
GACATATTGCCCGTTGCAGTCAGAATGAAAAGCTGA                                   >  1:3245904/1‑36 (MQ=255)
GACATATTGCCCGTTGCAGTCAGAATGAAAAGCTGA                                   >  1:2128837/1‑36 (MQ=255)
 ACATATTGCCCGTTGCAGTCAGAATGAAAAGCTGAA                                  <  1:5525579/36‑1 (MQ=255)
 ACATATTGCCCGTTGCAGTCAGAATGAAAAGCTGAA                                  >  1:680225/1‑36 (MQ=255)
  CATATTGCCCGTTGCAGTCAGAATGAAAAGCTGAAA                                 <  1:444800/36‑1 (MQ=255)
  CATATTGCCCGTTGCAGTCAGAATGAAAAGCTGAAA                                 <  1:5659477/36‑1 (MQ=255)
   ATATTGCCCGTTGCAGTCAGAATGAAAAGCTGAAAA                                >  1:450956/1‑36 (MQ=255)
     ATTGCCCGTTGCAGTCAGAATGAAAAGCTGAAAAAT                              <  1:2887342/36‑1 (MQ=255)
       TGCCCGTTGCAGTCAGAATGAAAAGCTGAAAAATAC                            <  1:4277365/36‑1 (MQ=255)
        GCCCGTTGCAGTCAGAATGAAAAGCTGAAAAATACT                           <  1:4467497/36‑1 (MQ=40)
         CCCGTTGCAGTCAGAATGAAAAGCTGAAAAATACTT                          <  1:3109698/36‑1 (MQ=38)
         CCCGTTGCAGTCAGAATGAAAAGCTGAAAAATACTT                          <  1:1942765/36‑1 (MQ=38)
          CCGTTGCAGTCAGAATGAAAAGCTGAAAAATACTTA                         >  1:5930862/1‑36 (MQ=35)
          CCGTTGCAGTCAGAATGAAAAGCTGAAAAATACTTA                         <  1:618313/36‑1 (MQ=35)
            GTTGCAGTCAGAATGAAAAGCTGAAAAATACTTACT                       >  1:4174205/1‑36 (MQ=1)
            GTTGCAGTCAGAATGAAAAGCTGAAAAATACTTACT                       <  1:5854212/36‑1 (MQ=1)
              TGCAGTCAGAATGAAAAGCTGAAAAATACTTACTAA                     <  1:3254594/36‑1 (MQ=1)
               GCAGTCAGAATGAAAAGCTGAAAAATACTTACTAAG                    <  1:2487674/36‑1 (MQ=1)
               GCAGTCAGAATGAAAAGCTGAAAAATACTTACTAAG                    >  1:3088407/1‑36 (MQ=1)
                CAGTCAGAATGAAAAGCTGAAAAATACTTACTAAGG                   <  1:698088/36‑1 (MQ=1)
                 AGTCAGAATGAAAAGCTGAAAAATACTTACTAAGGC                  <  1:3671556/36‑1 (MQ=1)
                  GTCAGAATGAAAAGCTGAAAAATACTTACTAAGGCG                 >  1:2188994/1‑36 (MQ=1)
                   TCAGAATGAAAAGCTGAAAAATACTTACTAAGGCGT                <  1:2263167/36‑1 (MQ=1)
                   TCAGAATGAAAAGCTGAAAAATACTTACTAAGGCGT                <  1:4660931/36‑1 (MQ=1)
                      GAATGAAAAGCTGAAAAATACTTACTAAGGCGTTTT             <  1:4484142/36‑1 (MQ=1)
                      GAATGAAAAGCTGAAAAATACTTACTAAGGCGTTTT             <  1:1707841/36‑1 (MQ=1)
                       AATGAAAAGCTGAAAAATACTTACTAAGGCGTTTTT            <  1:5826044/36‑1 (MQ=1)
                        ATGAAAAGCTGAAAAATACTTACTAAGGCGTTTTTa           >  1:4300760/1‑35 (MQ=1)
                        ATGAAAAGCTGAAAAATACTTACTAAGGCGTTTTTT           >  1:5710874/1‑36 (MQ=1)
                         TGAAAAGCTGAAAAATACTTACTAAGGCGTTTTTTA          <  1:5394531/36‑1 (MQ=1)
                           AAAAGCTGAAAAATACTTACTAAGGCGTTTTTTATT        <  1:980736/36‑1 (MQ=1)
                           AAAAGCTGAAAAATACTTACTAAGGCGTTTTTTATT        <  1:5758557/36‑1 (MQ=1)
                            AAAGCTGAAAAATACTTACTAAGGCGTTTTTTATTT       <  1:575708/36‑1 (MQ=1)
                            AAAGCTGAAAAATACTTACTAAGGCGTTTTTTATTT       <  1:5260462/36‑1 (MQ=1)
                              AGCTGAAAAATACTTACTAAGGCGTTTTTTATTTGG     <  1:4078219/36‑1 (MQ=1)
                              AGCTGAAAAATACTTACTAAGGCGTTTTTTATTTGG     >  1:1183299/1‑36 (MQ=1)
                               GCTGAAAAATACTTACTAAGGCGTTTTTTATTTGGT    >  1:3847650/1‑36 (MQ=1)
                                CTGAAAAATACTTACTAAGGCGTTTTTTATTTGGTG   >  1:5942538/1‑36 (MQ=1)
                                 TGAAAAATACTTACTAAGGCGTTTTTTATTTGGTGA  <  1:2094290/36‑1 (MQ=1)
                                                                     
GACATATTGCCCGTTGCAGTCAGAATGAAAAGCT...................................  <  REL606/34‑1
..................................GAAAAATACTTACTAAGGCGTTTTTTATTTGGTGA  <  REL606/4629812‑4629778

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