New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 526508NA (NA)3 (0.040) 3/260 NT 38.7% coding (3248/4281 nt) rhsD protein RhsD
?NC_000913 = 526599 5 (0.060)coding (3339/4281 nt) rhsD protein RhsD

GGTGACGTGGTATGGCTGGGACGGAGACAGGCTGACGACGGTGCAGACTGACACCACACGTATCCAGACGGTATACGAGCCGGGAAGCTTCACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGcagcctg‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/526362‑526508
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CAGCCTGTCCAGCAGCCGCACCAGTTCAGCCGGGAACACCACGCCGTGGCCGTTCTCACTCCCTTCCTGCTGGAGCGTCTCTGCCAGGCTGCGCCGCTGCGCTTTTTCCCGCTCGCCGTTCTCTGTCTCGACCCGGATGAGCGGCGTGAAGC  <  NC_000913/526599‑526448
                                                                                                                                                                                                                                                                                                           
GGTGACGTGGTATGGCTGGGACGGAGACAGGCTGACGACGGTGCAGACTGACACCACACGTATCCAGACGGTATACGAGCCGGGAAGCTTCACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAG                                                                                                                                                       >  2:1358372/1‑150
GGTGACGTGGTATGGCTGGGACGGAGACAGGCTGACGACGGTGCAGACTGACACCACACGTATCCAGACGGTATACGAGCCGGGAAGCTTCACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAG                                                                                                                                                       <  2:657321/150‑1
GGTGACGTGGTATGGCTGGGACGGAGACAGGCTGACGACGGTGCAGACTGACACCACACGTATCCAGACGGTATACGAGCCGGGAAGCTTCACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAG                                                                                                                                                       >  2:883212/1‑150
  TGACGTGGTATGGCTGGGACGGAGACAGGCTGACGACGGTGCAGACTGACACCACACGTATCCAGACGGTATACGAGCCGGGAAGCTTCACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCC                                                                                                                                                     <  2:1303366/150‑1
   GACGTGGTATGGCTGGGACGGAGACAGGCTGACGACGGTGCAGACTGACACCACACGTATCCAGACGGTATACGAGCCGGGAAGCTTCACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCT                                                                                                                                                    >  2:1398703/1‑150
                                               CTGACACCACACGTATCCAGACGGTATACGAGCCGGGAAGCTTCACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTG                                                                                                                                                   <  1:1345117/107‑1
                                               CTGACACCACACGTATCCAGACGGTATACGAGCCGGGAAGCTTCACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTG                                                                                                                                                   >  2:1345117/1‑107
                                                        CACGTATCCAGACGGTATACGAGCCGGGAAGCTTCACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGTCCAGCAGCCGCACCAGTTCAGCCGGGAACACCACGCCGTGGCCGTTCTCAC                                                                                               >  2:1167966/1‑150
                                                                                          CACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGTCCAGCAGCCGCACCAGTTCAGCCGGGAACACCACGCCGTGGCCGTTCTCACTCCCTTCCTGCTGGAGCGTCTCTGCCAGGCTGCG                                                             <  1:1167966/150‑1
                                                                                                        GGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGTCCAGCAGCCGCACCAGTTCAGCCGGGAACACCACGCCGTGGCCGTTCTCACTCCCTTCCTGCTGGAGCGTCTCTGCCAGGCTGCGCCGCTGCGCTTTTT                                               <  2:118349/150‑1
                                                                                                                                                     GCCTGTCCAGCAGCCGCACCAGTTCAGCCGGGAACACCACGCCGTGGCCGTTCTCACTCCCTTCCTGCTGGAGCGTCTCTGCCAGGCTGCGCCGCTGCGCTTTTTCCCGCTCGCCGTTCTCTGTCTCGACCCGGATGAGCGGCGTGAAGC  <  2:216449/150‑1
                                                                                                                                                                                                                                                                                                           
GGTGACGTGGTATGGCTGGGACGGAGACAGGCTGACGACGGTGCAGACTGACACCACACGTATCCAGACGGTATACGAGCCGGGAAGCTTCACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGcagcctg‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/526362‑526508
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CAGCCTGTCCAGCAGCCGCACCAGTTCAGCCGGGAACACCACGCCGTGGCCGTTCTCACTCCCTTCCTGCTGGAGCGTCTCTGCCAGGCTGCGCCGCTGCGCTTTTTCCCGCTCGCCGTTCTCTGTCTCGACCCGGATGAGCGGCGTGAAGC  <  NC_000913/526599‑526448

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