New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 588495 =NA (NA)40 (0.940) 30/66 NT 100% noncoding (1/1977 nt) IS150 repeat region
?REL606 1270660 = 0 (0.000)intergenic (‑167/+238) ldrC/chaA toxic polypeptide, small/calcium/sodium:proton antiporter

ATTTCATCGTCCAACAAAATGGGTGCAGTACA...................................  <  REL606/588526‑588495
..............................CACCGGATAATTCCGGCTTGGTGTGGATACTACTTCT  >  REL606/1270660‑1270694
                                                                   
ATTTCATCGTCCAACAAAATGGGGGCAGTACACCGG                                 >  1:4580688/1‑36 (MQ=0)
 TTTCATCGTCCAACAAAATGGGTGCAGTACACCGGA                                <  1:4613429/36‑1 (MQ=1)
  TTCATCGTCCAACAAAATGGGTGCAGTACACCGGAT                               <  1:5070007/36‑1 (MQ=1)
  TTCATCGTCCAACAAAATGGGTGCAGTACACCGGAT                               <  1:5594563/36‑1 (MQ=1)
     ATCGTCCAACAAAATGGGTGCAGTACACCGGATAAT                            >  1:1515675/1‑36 (MQ=1)
     ATCGTCCAACAAAATGGGTGCAGTACACCGGATAAT                            <  1:2170195/36‑1 (MQ=1)
       CGTCCAACAAAATGGGTGCAGTACACCGGATAATTC                          >  1:5873072/1‑36 (MQ=1)
        GTCCAACANNATGGGTGCAGTACACCGGATAATTCC                         <  1:433213/36‑1 (MQ=1)
        GTCCAACAAAATGGGTGCAGTACACCGGATAATTCC                         <  1:5536741/36‑1 (MQ=1)
        GTCCAACAAAATGGGGGCAGTACACCGGATAATTCC                         >  1:5623394/1‑36 (MQ=0)
         TCCAACAAAATGGGTGCAGTACACCGGATAATTCCG                        >  1:5605697/1‑36 (MQ=1)
         TCCAACAAAATGGGTGCAGTACACCGGATAATTCCG                        >  1:4315735/1‑36 (MQ=1)
         TCCAACAAAATGGGTGCAGTACACCGGATAATTCCG                        <  1:3888645/36‑1 (MQ=1)
          CCAACAAAATGGGTGCAGTACACCGGATAATTCCGG                       >  1:4497642/1‑36 (MQ=1)
          CCAACAAAATGGGTGCAGTACACCGGATAATTCCGG                       <  1:826857/36‑1 (MQ=1)
           CAACAAAATGGGTGCAGTACACCGGATAATTCCGGC                      <  1:128683/36‑1 (MQ=1)
             ACAAAATGGGTGCAGTACACCGGATAATTCCGGCTT                    <  1:5178675/36‑1 (MQ=1)
             ACAAAATGGGTGCAGTACACCGGATAATTCCGGCTT                    >  1:4768942/1‑36 (MQ=1)
             ACAAAATGGGTGCAGTACACCGGATAATTCCGGCTT                    <  1:4584963/36‑1 (MQ=1)
             ACAAAATGGGTGCAGTACACCGGATAATTCCGGCTT                    >  1:5751756/1‑36 (MQ=1)
             ACAAAATGGGTGCAGTACACCGGATAATTCCGGCTT                    <  1:4230018/36‑1 (MQ=1)
              CAAAATGGGTGCAGTACACCGGATAATTCCGGCTTG                   >  1:2356863/1‑36 (MQ=1)
              CAAAATGGGTGCAGTACACCGGATAATTCCGGCTTG                   >  1:5585103/1‑36 (MQ=1)
              CAAAATGGGTGCAGTACACCGGATAATTCCGGCTTG                   <  1:3536269/36‑1 (MQ=1)
                AAATGGGTTCAGTACACCGGATAATTCCGGCTTGGT                 <  1:1198653/36‑1 (MQ=0)
                AAATGGGTGCAGTACACCGGATAATTCCGGCTTGGT                 >  1:1446677/1‑36 (MQ=1)
                 AATGGGTGCAGTACACCGGATAATTCCGGCTTGGTG                <  1:3333517/36‑1 (MQ=1)
                  ATGGGTGCAGTACACCGGATAATTCCGGCTTGGTGT               >  1:5231706/1‑36 (MQ=1)
                   TGGGTGCAGTACACCGGATAATTCCGGCTTGGTGTG              <  1:525096/36‑1 (MQ=1)
                   TGGGTGCAGTACACCGGATAATTCCGGCTTGGTGTG              >  1:3775868/1‑36 (MQ=1)
                    GGGTGCAGTACACCGGATAATTCCGGCTTGGTGTGG             >  1:1496496/1‑36 (MQ=1)
                     GGTGCAGTACACCGGATAATTCCGGCTTGGTGTGGA            <  1:5405459/36‑1 (MQ=1)
                     GGTGCAGTACACCGGATAATTCCGGCTTGGTGTGGA            <  1:2093463/36‑1 (MQ=1)
                     GGTGCAGTACACCGGATAATTCCGGCTTGGTGTGGA            <  1:1523881/36‑1 (MQ=1)
                     GGTGCAGTACACCGGATAATTCCGGCTTGGTGTGGA            >  1:1189543/1‑36 (MQ=1)
                       TGCAGTACCCCGGATAATTCCGGCTTGGTGTGGATA          <  1:869381/36‑1 (MQ=0)
                          AGTACACCGGATAATTCCGGCTTGGTGTGGATACTA       >  1:4358576/1‑36 (MQ=1)
                           GTACACCGGATAATTCCGGCTTGGTGTGGATACTAC      <  1:4341235/36‑1 (MQ=1)
                           GTACACCGGATAATTCCGGCTTGGTGTGGATACTAC      <  1:2750286/36‑1 (MQ=1)
                             ACACCGGATAATTCCGGCTTGGTGTGGATACTACTT    >  1:4756720/1‑36 (MQ=1)
                                                                   
ATTTCATCGTCCAACAAAATGGGTGCAGTACA...................................  <  REL606/588526‑588495
..............................CACCGGATAATTCCGGCTTGGTGTGGATACTACTTCT  >  REL606/1270660‑1270694

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