New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4578923 =437 (0.900)5 (0.010) 3/266 NT 1.2% coding (415/1380 nt) mcrB 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB
?NC_000913 4579019 = 377 (0.810)coding (319/1380 nt) mcrB 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB
Rejected: Frequency below/above cutoff threshold.

CCACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGCCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4579076‑4578923
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttgcctGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCC  >  NC_000913/4579019‑4579167
                                                                                                                                                                                                                                                                                                               
CCACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTT                                                                                                                                                           >  2:2483170/1‑150
CCACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTT                                                                                                                                                           >  1:2689400/1‑150
CCACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTT                                                                                                                                                           <  1:4152861/150‑1
CCACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTT                                                                                                                                                           <  1:4153385/150‑1
CCACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTT                                                                                                                                                           <  1:7909059/150‑1
CCACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTT                                                                                                                                                           >  2:6103022/1‑150
CCACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTT                                                                                                                                                           >  2:4318717/1‑150
CCACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGAGTACACCCGATTT                                                                                                                                                           >  2:3086071/1‑150
 CACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTG                                                                                                                                                          >  2:7127118/1‑150
  ACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGC                                                                                                                                                         <  1:6301671/150‑1
  ACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGC                                                                                                                                                         <  2:3623878/150‑1
  ACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGC                                                                                                                                                         <  1:942766/150‑1
  ACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGC                                                                                                                                                         <  2:410131/150‑1
  ACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGC                                                                                                                                                         <  1:7647546/150‑1
  ACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGC                                                                                                                                                         >  1:7066449/1‑150
  ACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGC                                                                                                                                                         <  1:6154659/150‑1
  ACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGC                                                                                                                                                         <  2:4558161/150‑1
  ACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGC                                                                                                                                                         >  1:469996/1‑150
  ACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGC                                                                                                                                                         <  2:6547210/150‑1
  ACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATCTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGC                                                                                                                                                         <  2:4558154/150‑1
  ACATCCCCAATCGCAGTACTCTTCAGACATACCTAAAACAATCGCAGAGTAATTTCAGGCAACATCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGC                                                                                                                                                         <  2:3828920/150‑1
    ATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGCCT                                                                                                                                                       >  1:7169241/1‑150
                                                     TTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTG                                                                                                      >  1:1516746/1‑150
                                                     TTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTG                                                                                                      >  1:1516351/1‑150
                                                            GCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGT                                                                                               >  1:2481565/1‑150
                                                                                                               TCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCT                                            <  2:1516746/150‑1
                                                                                                               TCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCT                                            <  2:1516351/150‑1
                                                                                                                                                    TTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAG       <  2:2244073/150‑1
                                                                                                                                                     TGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGTATATATCCCGTTAttctc      <  2:3215532/150‑6
                                                                                                                                                       CCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTT    <  2:1859087/150‑1
                                                                                                                                                       CCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTT    <  2:7185154/150‑1
                                                                                                                                                       CCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTT    <  1:492720/150‑1
                                                                                                                                                       CCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTT    <  2:347588/150‑1
                                                                                                                                                       CCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTT    <  2:347597/150‑1
                                                                                                                                                       CCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTT    <  1:7656420/150‑1
                                                                                                                                                        CTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTCGAAGCTTC   >  2:7033915/1‑150
                                                                                                                                                        CTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTC   <  2:2921781/150‑1
                                                                                                                                                        CTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTC   <  2:7349408/150‑1
                                                                                                                                                        CTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTC   >  1:5750726/1‑150
                                                                                                                                                        CTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTC   >  1:6836634/1‑150
                                                                                                                                                        CTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTC   <  1:7274588/150‑1
                                                                                                                                                        CTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTC   <  2:3288756/150‑1
                                                                                                                                                        CTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTC   >  1:1343267/1‑150
                                                                                                                                                        CTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTC   >  2:301949/1‑150
                                                                                                                                                        CTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCAGGGGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTC   >  2:1956751/1‑150
                                                                                                                                                         TGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAATCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCC  >  2:3143902/1‑150
                                                                                                                                                         TGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCC  >  1:7946429/1‑150
                                                                                                                                                         TGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCC  >  1:4477523/1‑150
                                                                                                                                                         TGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACAGTTAGAAGCTTCC  >  1:4656723/1‑150
                                                                                                                                                                                                                                                                                                               
CCACATGCCCAATGGCAGTTCTCTTCAGACATACCTAAAACAATCGCAGAGTATTTTCAGGCAACTTCGGGTGTATATCCTAAAAAATACGGACAGTCCTATTACGCCTGTTCCCAAAAAGTCTCACAGGGTATTGATTACACCCGATTTGCCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4579076‑4578923
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttgcctGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCC  >  NC_000913/4579019‑4579167

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