New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3621534NA (NA)4 (0.050) 4/260 NT 100% coding (2343/4236 nt) rhsB rhs element protein RhsB
?NC_000913 = 3621592 0 (0.000)coding (2401/4236 nt) rhsB rhs element protein RhsB

CGGCTGACCGGTGAGCGTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTT  >  NC_000913/3621436‑3621592
                                                                                                  |                                                          
cGGCTGACCGGTGAGCGTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGC         >  1:106373‑M2/1‑150 (MQ=255)
  gCTGACCGGTGAGCGTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACAGCAG       <  1:965379‑M2/150‑1 (MQ=255)
   cTGACCGGTGAGCGTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGC      >  1:1594204‑M2/1‑150 (MQ=255)
   cTGACCGGTGAGCGTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGC      <  1:762135‑M2/150‑1 (MQ=255)
     gACCGGTGAGCGTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGG    >  1:547631‑M2/1‑150 (MQ=255)
     gACCGGTGAGCGTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGG    >  2:411730‑M2/1‑150 (MQ=255)
      aCCGGTGAGCGTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGT   >  1:1603352‑M2/1‑150 (MQ=255)
       ccGGTGAGCGTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTT  <  1:944018‑M2/150‑1 (MQ=255)
                                                                                                  |                                                          
CGGCTGACCGGTGAGCGTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTT  >  NC_000913/3621436‑3621592

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