New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3620569 =NA (NA)22 (0.180) 8/166 NT 100% coding (1378/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 3620582 = 0 (0.000)coding (1391/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

TCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCA  >  NC_000913/3620483‑3620681
                                                                                                   |                                                                                                   
caattcatcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcc                                                                                                     <  2:3029870‑M2/1‑1 (MQ=255)
 aattcatcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCg                                                                                                    <  2:2314783‑M2/2‑1 (MQ=255)
 aattcatcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCg                                                                                                    >  2:2253096‑M2/99‑100 (MQ=255)
  attcatcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCgg                                                                                                   <  1:1243870‑M2/3‑1 (MQ=255)
  attcatcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCgg                                                                                                   <  2:2343717‑M2/3‑1 (MQ=255)
  attcatcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCgg                                                                                                   <  2:646004‑M2/3‑1 (MQ=255)
   ttcatcatattcccggcgcaattccagcccgtcaggcccggtgtctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCGGc                                                                                                  <  1:2274260‑M2/4‑1 (MQ=255)
   ttcatcatattcccggcgcaattccagcccgtcaggcccggtggctgccgttacctggttgtgggggtttcagtaaaacgccgctcccctgccatcCGGc                                                                                                  >  2:1208546‑M2/97‑100 (MQ=255)
   ttcatcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCGGc                                                                                                  >  2:2075618‑M2/97‑100 (MQ=255)
    tcatcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCGGCg                                                                                                 >  1:3122530‑M2/96‑100 (MQ=255)
    tcatcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCGGCg                                                                                                 >  1:937939‑M2/96‑100 (MQ=255)
    ncatcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCGGCg                                                                                                 >  1:2571290‑M2/96‑100 (MQ=255)
     catcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCGGCg                                                                                                 >  1:2688777‑M2/95‑99 (MQ=255)
     catcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCGGCGt                                                                                                >  2:3535366‑M2/95‑100 (MQ=255)
     catcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCGGCGt                                                                                                >  2:1370508‑M2/95‑100 (MQ=255)
     catcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCGGCGt                                                                                                >  2:123041‑M2/95‑100 (MQ=255)
      atcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCGGCGt                                                                                                <  2:2926169‑M2/6‑1 (MQ=255)
                                                                  gtggttatagtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGgctggaaatacg                                   >  2:1690539‑M2/34‑95 (MQ=255)
                                                                    ggttatagtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGcc                                 >  2:1403409‑M2/32‑100 (MQ=255)
                                                                      ttatagtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGACATACGCCgg                               >  1:794678‑M2/30‑100 (MQ=255)
                                                                         tagtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGa                             >  1:1378533‑M2/27‑99 (MQ=255)
                                                                          agtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAat                           >  2:3209292‑M2/26‑100 (MQ=255)
                                                                          agtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAat                           >  2:2674746‑M2/26‑100 (MQ=255)
                                                                          agtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAat                           >  1:1903268‑M2/26‑100 (MQ=255)
                                                                          agtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGGCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGaa                            >  1:1545968‑M2/26‑99 (MQ=255)
                                                                           gtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAata                          >  1:2151214‑M2/25‑100 (MQ=255)
                                                                           gtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAata                          >  1:745081‑M2/25‑100 (MQ=255)
                                                                           gtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAata                          >  2:3395073‑M2/25‑100 (MQ=255)
                                                                           gtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAata                          >  2:3180079‑M2/25‑100 (MQ=255)
                                                                           gtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAata                          >  2:3050336‑M2/25‑100 (MQ=255)
                                                                           gtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAata                          >  1:990428‑M2/25‑100 (MQ=255)
                                                                            taaaacgccgatgccctgccatcCGGCGTTTTACTATAACCAACACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAatat                         >  1:2979807‑M2/24‑100 (MQ=255)
                                                                             aaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatg                        >  2:378206‑M2/23‑100 (MQ=255)
                                                                             aaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatg                        >  2:1884105‑M2/23‑100 (MQ=255)
                                                                             aaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatg                        >  2:669545‑M2/23‑100 (MQ=255)
                                                                              aaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatga                       >  1:3320131‑M2/22‑100 (MQ=255)
                                                                               aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCCGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatgat                      >  2:2829330‑M2/21‑100 (MQ=255)
                                                                                acgccgatgccctgccatcCGGCGTTTTACTATAACCGCCACAACGAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatgatg                     <  1:378204‑M2/81‑1 (MQ=255)
                                                                                 cgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATAtgatga                    >  1:1813030‑M2/19‑100 (MQ=255)
                                                                                                   CGGCGTTTTACTATAACCACCACATCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGAtt    <  1:2556104‑M2/98‑1 (MQ=255)
                                                                                                   CGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCa  <  2:776367‑M2/100‑1 (MQ=255)
                                                                                                   CGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCa  >  2:1939263‑M2/1‑100 (MQ=255)
                                                                                                   CGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCa  <  1:3110095‑M2/100‑1 (MQ=255)
                                                                                                   CGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCa  <  1:1244206‑M2/100‑1 (MQ=255)
                                                                                                   CGGCGTTTTACTATAACCACCACAACCAGTTAACGCCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGACTTGGGCCGTCTGATTCa  >  1:635205‑M2/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCA  >  NC_000913/3620483‑3620681

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