New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2574614 =75 (1.380)4 (0.090) 4/126 5.3 6.3% intergenic (+456/+362) sseA/sseB 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
?REL606 2574646 = 57 (1.260)intergenic (+488/+330) sseA/sseB 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTTTGCGAGACGGCGGAAAGGGCCGCAGTGGATAATACACTAACAGAACAATAATAACTCTCTTTTGCTTGACAAAAAGAGAGTTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2574699‑2574614
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaaagagagttaTTATTGTTCTGTTAGTGTATTATCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAG  >  REL606/2574646‑2574719
                                                                                                                                                                
GTTTGCGAGACGGCGGAAAGGGCCGCAGTGGATAATACACTAACAGAACAATAATAACTCTCTTTTGCTTGACAAA                                                                                      <  1:1593450/76‑1
  TTGCGAGACGGCGGAAAGGGCCGCAGTGGATAATACACTAACAGAACAATAATAACTCTCTTTTGCTTGACAAAAA                                                                                    >  1:1102483/1‑76
   gGCGAGACGGCGGAAAGGNCCGCAGTGGATAATACACTAACAGAACAATAATAACTCTCTTTTGCTTGACAAAAAG                                                                                   <  2:1573615/75‑1
   TGCGAGACGGCGGAAAGGGCCGCAGTGGATAATACACTAACAGAACAATAATAACTCTCTTTTGCTTGACAAAAAG                                                                                   <  2:875895/76‑1
       AGACGGCGGAAAGGGCCGCAGTGGATAATACACTAACAGAACAATAATAACTCTCTTTTGCTTGACAAAAAGAGAG                                                                               <  1:1334927/76‑1
        GACGGCGGAAAGGGCCGCAGTGGATAATACACTAACAGAACAATAATAACTCTCTTTTGCTTGACAAAAAGAGAGT                                                                              >  2:512804/1‑76
        GACGGCGGAAAGGGCCGCAGTGGATAATACACTAACAGAACAATAATAACTCTCTTTTGCTTGACAAAAAGAGAGT                                                                              >  1:318246/1‑76
          CGGCGGAAAGGGCCGCAGTGGATAATACACTAACAGAACAATAATAACTCTCTTTTGCTTGACAAAAAGAGAGTTA                                                                            <  2:896838/76‑1
                            TGGATAATACATTAACAGAACAATAATAACTCTCTTTTGCTTGACAAAAAGAGAGTTATTATTGTTCTGTTAGTGT                                                          >  1:572543/1‑76
                                 AATACACTAACAGAACAATAATAACTCTCTTTTGCTTGACAAAAAGAGAGTTATTATTGTTCTGTTAGTGTATT                                                       <  2:1006877/74‑1
                                 AATACACTAACAGAACAATAATAACTCTCTTTTGCTTGACAAAAAGAGAGTTATTATTGTTCTGTTAGTGTATT                                                       >  1:1006877/1‑74
                                    ACACTAACAGAACAATAATAACTCTCTTTTGCTTGACAAAAAGAGAGTTATTATTGTTCTGTTAGTGTATTATCCA                                                  <  2:572543/76‑1
                                                                          AAAAGAGAGTTATTATTGTTCTGTTAGTGTATTATCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGC            >  2:736754/1‑76
                                                                          AAAAGAGAGTTATTATTGTTCTGTTAGTGTATTATCCACTGCGGCCCTTTCCGCCGACTCGCAAACGGGCGCAGGC            >  1:1270582/1‑76
                                                                             AGAGAGTTATTATTGTTCTGTTAGTTTATTATCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTT         <  2:341876/76‑1
                                                                             AGAGAGTTATTATTGTTCTGTTAGTGTATTATCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTT         <  2:1152663/76‑1
                                                                             AGAGAGTTATTATTGTTCTGTTAGTGTATTATCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTT         <  2:808392/76‑1
                                                                              GAGAGTTATTATTGTTCTGTTAGTGTATTATCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTA        <  2:1387933/76‑1
                                                                                  GTTTTTATTGTTCTGTTAGTGTATTATCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAA    <  1:1425799/76‑1
                                                                                    TATTATTGTTCTGTTAGTGTATTATCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAG  >  1:1413032/1‑76
                                                                                                                                                                
GTTTGCGAGACGGCGGAAAGGGCCGCAGTGGATAATACACTAACAGAACAATAATAACTCTCTTTTGCTTGACAAAAAGAGAGTTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2574699‑2574614
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaaagagagttaTTATTGTTCTGTTAGTGTATTATCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAG  >  REL606/2574646‑2574719

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.