New junction evidence
  seq id position reads (cov) reads (cov) score skew annotation gene product
* ? REL606 = 161953134 (0.780)10 (0.230) 6/56 4.1 intergenic (‑140/+77) rem/hokD predicted protein/small toxic polypeptide
?REL606 2353399 = 35 (0.640)coding (1545/1833 nt) yfbS predicted transporter
Rejected: Position hash score below cutoff.

GCGGTATTAGAGATAAATAGCCCAATAA.CCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCC  >  REL606/2353364‑2353434
                                    |                                   
gcgGTATTAGAGATAAATAGCCCACTAC.CCGCCgac                                     >  1:3378595/1‑36 (MQ=255)
gcgGTATTAGAGATAAATAGCCCAATAA.CCGCCgac                                     >  1:2941975/1‑36 (MQ=255)
gcgGTATTAGAGATAAATAGCCCAATAA.CCGCCgac                                     >  1:2886124/1‑36 (MQ=255)
 cgGTATTAGAGATAAATAGCCCAATAA.CCGCCgaca                                    >  1:1181281/1‑36 (MQ=255)
  ggTATTAGAGATAAATAGCCCAATAA.CCGCCgacaa                                   >  1:1088828/1‑36 (MQ=255)
  ggTATTAGAGATAAATAGCCCAATAA.CCACCgacag                                   >  1:6598743/1‑35 (MQ=255)
      ttAGAGATAAATAGCCCAATAA.CCGCCgaccagaca                               >  1:654592/1‑36 (MQ=255)
       tAGAGATAAATAGCCCAATAA.CCGCCgacaagacaa                              >  1:6785865/1‑36 (MQ=255)
       tAGAGATAAATAGCCCAATAA.CCGCCgacaagacaa                              >  1:6692489/1‑36 (MQ=255)
        agagATAAATAGCCCAATAA.CCGCCGACAAGacaaa                             >  1:212252/1‑36 (MQ=255)
          agaTAAATAGCCCAATAA.CCGCCGACAAGacaaaca                           >  1:6985198/1‑36 (MQ=255)
           gaTAAATAGCCCAATAA.CCGCCGACAAGacaaacaa                          >  1:3839327/1‑36 (MQ=255)
             tAAATAGCCCAATAA.CCGCCGACAAGACaaacaaac                        >  1:6225728/1‑36 (MQ=255)
              aaaTAGCCCAATAA.CCGCCGACAAGACAaacaaaca                       >  1:2701410/1‑36 (MQ=255)
                   acccAAAAACCCGCCGACAAGACAAACAAACAcccc                   <  1:6284851/35‑1 (MQ=21)
                    cccAATAA.CCGCCGACAAGACAAACAAACACCCCac                 >  1:4052582/1‑35 (MQ=255)
                     ccAATAA.CCGCCGACAAGACAAACAAACACCCCAGc                >  1:6961157/1‑36 (MQ=255)
                       aaTAA.CCGCCGACAAGCCAAACACCCACCCCACcac              >  1:2949257/1‑35 (MQ=18)
                       aaTAA.CCGCCGACAAGACAAACAAACACCCCAGcat              >  1:1994592/1‑36 (MQ=255)
                       aaTAA.CCGCCGACAAGACAAACAAACACCCCAGcat              >  1:1838703/1‑36 (MQ=255)
                        aTAA.CCGCCGACAAGACAAACAAACACCCCAGcatc             >  1:3153489/1‑36 (MQ=255)
                         tAA.CCGCCGACAAGACAAACAAACACCCCCGcatca            >  1:5897608/1‑36 (MQ=255)
                         tAA.CCGCCGACAAGACAAACAAACACCCCAGcatca            >  1:4929548/1‑36 (MQ=255)
                         tAA.CCGCCGACAAGACAAACAAACACCCCAGcatca            >  1:373835/1‑36 (MQ=255)
                         tAA.CCGCCGACAAGACAAACAAACACCCCAGcatca            >  1:2010629/1‑36 (MQ=255)
                          aa.CCGCCGACAAGACAAACAAACACCCCAGcatcat           <  1:3517215/36‑1 (MQ=255)
                          aa.CCGCCGACAAGACAAACAAACACCCCAGcatcat           <  1:1192713/36‑1 (MQ=255)
                           acccgccgACAAGACAAACAAACACCCCAGcatcat           <  1:1091491/34‑1 (MQ=255)
                           a.CCGCCGACAAGACAAACAAACACCCCAGcaccaca          >  1:1342598/1‑34 (MQ=25)
                           a.CCGCCGACAAGACAAACAAACACCCCAGCATCata          >  1:3967944/1‑36 (MQ=255)
                           a.CCGCCGACAAGACAAACAAACACCCCAGCATCata          >  1:7382238/1‑36 (MQ=255)
                            cccgccgACAAGACAAACAAACACCCCAGCATCata          <  1:5278856/35‑1 (MQ=255)
                             ccgccgACAAGACAAACAAACACCCCAGCATCatat         <  1:2681226/36‑1 (MQ=255)
                             ccgccgACAAGACAAACAAACACCCCAGCATCatat         >  1:7364381/1‑36 (MQ=255)
                             ccgccgACAAGACAAACAAACACCCCAGCAACatat         >  1:6913912/1‑36 (MQ=255)
                                ccgACAAGACAAACAAACACCCCAGCATCATAtgtg      >  1:2616797/1‑36 (MQ=255)
                                 cgACAAGACAAACAAACACCCCAGCATCATATGTgg     <  1:5436844/36‑1 (MQ=255)
                                  gacaagacaaACAAACACCCCAGCATCATATGCGGc    >  1:6769577/1‑36 (MQ=255)
                                    caagacaaACAAACACCCCAGCATCATATGTGGccc  >  1:1149587/1‑36 (MQ=255)
                                    |                                   
GCGGTATTAGAGATAAATAGCCCAATAA.CCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCC  >  REL606/2353364‑2353434

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 17 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG