New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1299694 =34 (0.740)11 (0.240) 5/52 NT 29.3% coding (919/1632 nt) oppA oligopeptide transporter subunit
?REL606 = 2574502 19 (0.310)intergenic (+344/+474) sseA/sseB 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate

ACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGC  >  REL606/1299659‑1299728
                                   |                                  
aCCTGTGCACTTACTATTACGAAATTAACAACCAGa                                    >  1:4158768/1‑36 (MQ=255)
 aaacatgtaaattcactttcctttaaaaaacaaaaa                                   <  1:7656191‑M1/2‑1 (MQ=255)
  cTGTGCACTTACTATTACGAAATTAACAACCAGaaa                                  <  1:6781377/36‑1 (MQ=255)
  cTGTGCACTTACTATTACGAAATTAACAACCAGaaa                                  <  1:6597322/36‑1 (MQ=255)
  cTGTGCACTTACTATTACGAAATTAACAACCAGaaa                                  >  1:4119067/1‑36 (MQ=255)
  aacatgtaaattcactttcctttaaaaaacaaaaaa                                  <  1:7657773‑M1/3‑1 (MQ=255)
  aacatgtaaattcactttcctttaaaaaacaaaaaa                                  >  1:2573126‑M1/34‑36 (MQ=255)
  aacatgtaaattcactttcctttaaaaaacaaaaaa                                  <  1:3205902‑M1/3‑1 (MQ=255)
   tgtgCACTTACTATTACGAAATTAACAACCAGAAAc                                 <  1:3146587/36‑1 (MQ=255)
   tgtgCACTTACTATTACGAAATTAACAAACAGAAAc                                 >  1:6965856/1‑36 (MQ=255)
   acatgtaaattcactttcctttaaaaaacaaaaaac                                 <  1:4406065‑M1/4‑1 (MQ=255)
    gtgCACTTACTATTACGAAATTAACAACCAGAAAcc                                >  1:7502927/1‑36 (MQ=255)
     atgtaaattcactttcctttaaaaaacaaaAAACCg                               <  1:5464243‑M1/6‑1 (MQ=255)
      tgtaaattcactttcctttaaaaaacaaaAAACCGc                              <  1:1279457‑M1/7‑1 (MQ=255)
       gtaaattcactttcctttaaaaaacaaaAAACCGcc                             <  1:4479403‑M1/8‑1 (MQ=255)
       cACTTACTATTACGAAATTAACAACCAGAAACCGcc                             <  1:3911650/36‑1 (MQ=255)
        taaattcactttcctttaaaaaacaaaAAACCGCCa                            <  1:3641118‑M1/9‑1 (MQ=255)
         cTTACTATTACGAAATTAACAACCAGAAACCGCCAt                           >  1:2932721/1‑36 (MQ=255)
          ttACTATTACGAAATTAACAACCAGAAACCGCCAtt                          >  1:5484541/1‑36 (MQ=255)
          ttACTATTACGAAATTAACAACCAGAAACCGCCAtt                          >  1:2400633/1‑36 (MQ=255)
              tATTACGAAATTAACAACCAGAAACCGCCATTCAAc                      >  1:7787863/1‑36 (MQ=255)
                ttACGAAATTAACAACCAGAAACCGCCATTCAACGa                    >  1:1632217/1‑36 (MQ=255)
                 tACGAAATTAACAACCAGAAACCGCCATTCAACGAt                   >  1:1980022/1‑36 (MQ=255)
                 tACGAAATTAACAACCAGAAACCGCCATTCAACGAt                   >  1:4932727/1‑36 (MQ=255)
                    aaaaTAAAAAACAAAAAACCGCCATTCAACGAtgtg                <  1:5733929/35‑1 (MQ=9)
                     aaataaaaaacaaaAAACCGCCATTCAACGAtgtgc               >  1:4593995‑M1/15‑36 (MQ=255)
                     aaaTAAACAACAAGAAACCTCCATTCAACGAtgtgc               <  1:167565/36‑1 (MQ=18)
                      aaaaaaaaacaaaAAACCGCCATTCAACGAtgtgcg              >  1:3267140‑M1/14‑36 (MQ=255)
                      aaTTAACAACCAGAAACCGCCATTCAACGAtgtgcg              >  1:343236/1‑36 (MQ=255)
                      aaTTAACAACCAGAAACCGCCATTCAACGAtgtgcg              >  1:2431886/1‑36 (MQ=255)
                      aaTTAACAACCAGAAACCGCCATTCAACGAtgtgcg              >  1:6551779/1‑36 (MQ=255)
                       aTTAACAACCAGAAACCGCCATTCAACGAtgtgcgt             >  1:4758528/1‑36 (MQ=255)
                       aTTAACAACCAGAAACCGCCATTCAACGAtgtgcgt             >  1:8161163/1‑36 (MQ=255)
                         tAACAACCAGAAACCGCCATTCAACGAtgtgcgtgt           >  1:3765166/1‑36 (MQ=255)
                         tAACAACCAGAAACCGCCATTCAACGAtgtgcgtgt           >  1:2560066/1‑36 (MQ=255)
                         tAACAACAAAAAACCGCCATACAACGAtgtgcgtgt           <  1:3045779/36‑1 (MQ=21)
                         aaacaacaaaAAACCGCCATTCAACGAtgtgcgtgt           >  1:3058301‑M1/11‑36 (MQ=255)
                         aaaaaacaaaAAACCGCCATTCAACGAtgtgcgtgt           >  1:6298793‑M1/11‑36 (MQ=255)
                          aaaaacaaaAAACCGCCATTCAACGAtgtgcgtgtg          >  1:8372387‑M1/10‑36 (MQ=255)
                            caacCAGAAACCGCCATTCAACGAtgtgcgtgtgcg        >  1:1308442/1‑36 (MQ=255)
                            caacAAAAACCCGCCATTCAACGGtgtgcgtgtgcg        <  1:8146561/36‑1 (MQ=19)
                             aacaagAAACCGCCATTCAACGATgtgcgtgtgcgt       >  1:7575555‑M1/7‑36 (MQ=255)
                             aacaaaAAACCGCCATTCAACGATgtgcgtgtgcgt       >  1:3885427‑M1/7‑36 (MQ=255)
                             aacCAGAAACCGCCATTCAACGATgtgcgtgtgcgt       >  1:6720912/1‑36 (MQ=255)
                              acCAGAAACCGCCATTCAACGATGTGCGTGTGCGTa      >  1:5100140/1‑36 (MQ=255)
                               ccAGAAACCGCCATTCAACGATGTGCGTGTGCGTAc     >  1:3923468/1‑36 (MQ=255)
                               ccAGAAACCGCCATTCAACGATGTGCGTGTGCGTAc     >  1:7627673/1‑36 (MQ=255)
                               ccAGAAACCGCCATTCAACGATGTGCGTGTGCGTAc     >  1:3847288/1‑36 (MQ=255)
                                 aaAAACCGCCATTCAACGATGTGCGTGTGCGTACCg   >  1:479846‑M1/3‑36 (MQ=255)
                                 aaAAACCGCCATTCAACGATGTGCGTGTGCGTACCg   >  1:4807118‑M1/3‑36 (MQ=255)
                                 aGAAACCGCCATTCAACGATGTGCGTGTGCGTACCg   >  1:6996025/1‑36 (MQ=255)
                                  gAAACCGCCATTCAACGATGTGCGTGTGCGTACCGc  >  1:287621/1‑36 (MQ=255)
                                  gAAACCGCCATTCAACGATGTGCGTGTGCGTACCGc  >  1:2819639/1‑36 (MQ=255)
                                  aAACCCGCCATTCAACGATGTGCGTGTGCGTACCGc  >  1:6288633‑M1/2‑36 (MQ=255)
                                  aAAACCGCCATTCAACGATGTGCGTGTGCGTACCGc  >  1:764896‑M1/2‑36 (MQ=255)
                                   |                                  
ACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGC  >  REL606/1299659‑1299728

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG