New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 362270937 (0.730)6 (0.120) 6/246 5.4 7.2% coding (3518/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 = 3767498 119 (2.450)pseudogene (278/945 nt) rhsJ Rhs‑family protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATA  >  NC_000913/3622580‑3622841
                                                                                                                                 |                                                                                                                                    
aaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCtgaatgaa                                                                                                                   <  4:85196/149‑1 (MQ=32)
aaaaCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCtgaatgaa                                                                                                                   <  6:32081/149‑1 (MQ=32)
        aTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAAccc                                                                                                           <  3:14924/149‑1 (MQ=34)
                         gCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTgcagccg                                                                                           >  3:69807/1‑146 (MQ=37)
                                 gAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCtgaatgaa                                                                                                                   >  5:87578/1‑116 (MQ=33)
                                 gAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCtgaatgaa                                                                                                                   <  6:87578/116‑1 (MQ=33)
                                 gAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAag                                                                                                                  <  1:266097/117‑1 (MQ=34)
                                 gAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAag                                                                                                                  >  2:266097/1‑117 (MQ=34)
                                        ccACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCgg                                                                           >  2:224308/1‑149 (MQ=255)
                                            cTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCgggg                                                                         >  8:99338/1‑147 (MQ=255)
                                            cTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCgggg                                                                         <  7:99338/147‑1 (MQ=255)
                                            cTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCCTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCgggg                                                                         >  4:56466/1‑147 (MQ=255)
                                                  cACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTgcagca                                                                                            <  2:55594/122‑1 (MQ=38)
                                                  cACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTgcagca                                                                                            >  1:55594/1‑122 (MQ=38)
                                                    cTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTat                                                                >  7:186225/1‑148 (MQ=255)
                                                         aCCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATgag                                                          <  3:56466/149‑1 (MQ=255)
                                                                 ggCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCgg                                                  >  3:78871/1‑149 (MQ=255)
                                                                         gCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTAtt                                          <  4:78871/149‑1 (MQ=255)
                                                                                 tcatcaGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGc                                  >  2:113780/1‑149 (MQ=255)
                                                                                            gAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTAtt                                          >  8:199371/1‑130 (MQ=40)
                                                                                            gAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTAtt                                          <  7:199371/130‑1 (MQ=40)
                                                                                             aaGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCtattat                      <  8:186225/149‑1 (MQ=39)
                                                                                               ggggCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGa                    >  8:193191/1‑149 (MQ=39)
                                                                                                gggCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGAc                   <  1:113780/149‑1 (MQ=39)
                                                                                                gggCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGAc                   >  3:143503/1‑149 (MQ=39)
                                                                                                    aaCAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATc                                                                                                      >  5:120900/1‑62 (MQ=40)
                                                                                                    aaCAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATc                                                                                                      <  6:120900/62‑1 (MQ=40)
                                                                                                      cAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACAGCCCCCGCTATTATGACCCGCTg             >  8:171136/1‑149 (MQ=255)
                                                                                                         cGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAg          <  7:193191/149‑1 (MQ=35)
                                                                                                          gtggtgCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCagcagc                                                                                          >  2:95713/1‑68 (MQ=37)
                                                                                                          gtggtgCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCagcagc                                                                                          <  1:95713/68‑1 (MQ=37)
                                                                                                            ggtgCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAgggg       <  7:171136/149‑1 (MQ=34)
                                                                                                                cgcAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGat   <  1:224308/149‑1 (MQ=35)
                                                                                                                 gcAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata  >  3:201248/1‑149 (MQ=37)
                                                                                                                 gcAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata  <  4:143503/149‑1 (MQ=37)
                                                                                                                            tgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTc                                                     >  7:70555/1‑87 (MQ=33)
                                                                                                                            tgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTc                                                     <  8:70555/87‑1 (MQ=33)
                                                                                                                                 |                                                                                                                                    
AAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGGCAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATA  >  NC_000913/3622580‑3622841

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