New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3620131 =NA (NA)3 (0.080) 3/258 NT 100% coding (940/4236 nt) rhsB rhs element protein RhsB
?NC_000913 3620147 = 0 (0.000)coding (956/4236 nt) rhsB rhs element protein RhsB

GTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGG  >  NC_000913/3620125‑3620289
                      |                                                                                                                                              
cgccgcgtggcgtccagccataGCGGCGAACTGGCGGTGGTGTATGACCGTAGTgg                                                                                                               <  1:766357‑M2/34‑1 (MQ=255)
cgccgcgtggcgtccagccataGCGGCGAACTGGCGGTGGTGTATGACCGTAGTgg                                                                                                               >  2:766357‑M2/23‑56 (MQ=255)
   ccgcgtggcgtccaggcacGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAGACCGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCAGCGTAACACGGGCCGGCAggaaacccgggacag                     >  2:389047‑M2/20‑137 (MQ=255)
                      GCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGAcgggcgg  <  1:155231‑M2/143‑1 (MQ=255)
                      |                                                                                                                                              
GTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGG  >  NC_000913/3620125‑3620289

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