New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 3006451 =87 (0.850)9 (0.100) 4/44 NT 11.8% coding (461/819 nt) yafZ conserved protein
?REL606 4115747 = 55 (0.580)coding (1806/2181 nt) katG catalase/hydroperoxidase HPI(I)

GCTTTTCCTCAATCCGGTCAAAAACCCCCACCACCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/3006486‑3006451
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctGACGCTGACCGCGCCGGAAATGACTGCGCTGGTG  >  REL606/4115747‑4115780
                                                                      
GCTTTTCCTCAATCCGGTCAAAAACCCCCACCACC                                     >  1:4133989/1‑35
GCTTTTCCTCAATCCGGTCAAAAACCCCCACCACC                                     >  1:18387520/1‑35
GCTTTTCCTCAATCCGGTCAAAAACCCCCACCACC                                     >  1:12916049/1‑35
 aTTTTCCTCAATCCGGTCAAAACCCCCCACCACCT                                    <  1:12343930/34‑1
 CTTTTCCTCAATCCGGTCAACAACCCCCACCACCT                                    >  1:17330895/1‑35
 CTTTTCCTCAATCCGGTCAAAAACCCCCCCCACCT                                    <  1:19477543/35‑1
 CTTTTCCTCAATCCGGTCAAAAACCCCCACCACCT                                    >  1:4480962/1‑35
 CTTTTCCTCAATCCGGTCAAAAACCCCCACCACCT                                    >  1:15198347/1‑35
 CTTTTCCTCAATCCGGTCAAAAACCCCCACCACCT                                    >  1:13360686/1‑35
 CTTTTCCTCAATCCGGTCAAAAACCCCCACCACCT                                    <  1:12425915/35‑1
 CTTTTCCTCAATCCGGTCAAAAACCCCCACCACCT                                    <  1:10350784/35‑1
 CTTTTCCTCAATCCGGTCAAAAACCCCCACCACCT                                    >  1:10874635/1‑35
                     AAACCCCCACCACCTGACGCTGACCGCGCCGGAAA                <  1:11358228/35‑1
                     AAACCCCCACCACCTGACGCTGACCGCGCCGGAAA                <  1:2578947/35‑1
                      AACCCCCACCACCTGACGCTGACCGCGCAGGAAAT               <  1:482906/35‑1
                       ACCCCCACCCCCTGACGCTGACCGCGCCGGAAATG              <  1:17304241/35‑1
                           caACCACCTGACGCTGACCGCGCCGGAAATGACTG          <  1:17535526/33‑1
                           aCACCACCTGACGCTGACCGCGCCGGAAATGACTG          <  1:9018580/34‑1
                           CCACCACCTGACGCTGACCGCGCCGGAAATGACTG          <  1:7306515/35‑1
                           CCACCACCTGACGCTGACCGCGCCGGAAATGACTG          <  1:18268953/35‑1
                           CCACCACCTGACGCTGACCGCGCCGGAAATGACTG          <  1:4182288/35‑1
                           CCACCACCTGACGCTGACCGCGCCGGAAATGACTG          <  1:2521851/35‑1
                               ccCCTGACGCTGACCGCGCCGGAAATGACTGCGCT      <  1:10425245/33‑1
                               CACCTGACGCTGACCGCGCCGGAAATGACTGCGCT      <  1:9598456/35‑1
                               CACCTGACGCTGACCGCGCCGGAAATGACTGCGCT      <  1:8940918/35‑1
                                cCCTGACGCTGCCCGCGCCGGAAATGACTGCGCTG     <  1:13207224/34‑1
                                cCCTGACGCTGCCCGCGCCGGAAATGACTGCGCTG     <  1:14739271/34‑1
                                ACCTGACGCTGACCGCGCCGGAAATGACTGCGCTG     <  1:7888284/35‑1
                                ACCTGACGCTGACCGCGCCGGAAATGACTGCGCTG     <  1:15020124/35‑1
                                 CCTGACGCTGACCGCGCCGGAAATGACTGCGCTGG    <  1:6144945/35‑1
                                 CCTGACGCTGACCGCGCCGGAAATGACTGCGCTGG    <  1:20055334/35‑1
                                  CTGACGCTGACCGCGCCGGAAATGACTGCGCTGGT   >  1:18971399/1‑35
                                  CTGACGCTGACCGCGCCGGAAATGACTGCGCTGGT   <  1:18933139/35‑1
                                  CTGACGCTGACCGCGCCGGAAATGACTGCGCTGGT   <  1:12507966/35‑1
                                  CTGACGCTGACCGCGCCGGAAATGACTGCGCTGGT   >  1:2915549/1‑35
                                  CTGACGCTGACCGCGCCGGAAATGACTGCGCTGGT   >  1:78794/1‑35
                                  CTGACGCTGACCGCGCCGGAAATGACTGCGCTGGT   <  1:20101299/35‑1
                                   TGACGCTGACCGCGCCGGAAATGACTGCGCTGGTG  <  1:7156370/35‑1
                                   TGACGCTGACCGCGCCGGAAATGACTGCGCTGGTG  <  1:8134000/35‑1
                                   TGACGCTGACCGCGCCGGAAATGACTGCGCTGGTG  <  1:14838829/35‑1
                                   TGACGCTGACCGCGCCGGAAATGACTGCGCTGGTG  >  1:20473329/1‑35
                                   TGACGCTGACCGCGCCGGAAATGACTGCGCTGGTG  >  1:2265940/1‑35
                                   TGACGCTGACCGCGCCGGAAATGACTGCACTGGTG  >  1:3714628/1‑35
                                                                      
GCTTTTCCTCAATCCGGTCAAAAACCCCCACCACCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/3006486‑3006451
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctGACGCTGACCGCGCCGGAAATGACTGCGCTGGTG  >  REL606/4115747‑4115780

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ 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.