New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4583137145 (0.830)6 (0.040) 5/266 NT 3.8% intergenic (‑89/+112) hsdM/hsdR DNA methyltransferase M/endonuclease R Type I restriction enzyme
?NC_000913 = 4583224 163 (0.980)intergenic (‑176/+25) hsdM/hsdR DNA methyltransferase M/endonuclease R Type I restriction enzyme
Rejected: Frequency below/above cutoff threshold.

AACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4582983‑4583137
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttcaggCTGCCTCTGCGTTGGCTGCGCTCGTTCACCCCGGTCACGTACTTCTGTACGCTCCCGGGGATTCACTCACTTGCCGCCTTGATGCAACCTGAATGATTTTGTGTATATTACCCTCGGCAATTTCTTCTTCTGCGGCTCGATGAATTTGG  <  NC_000913/4583224‑4583076
                                                                                                                                                                                                                                                                                                                
AACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCAT                                                                                                                                                            >  2:1469078/1‑150
AACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCAT                                                                                                                                                            >  2:2358776/1‑150
AACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCAT                                                                                                                                                            <  1:1507279/150‑1
AACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCAT                                                                                                                                                            <  1:1507285/150‑1
AACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCAT                                                                                                                                                            >  2:2529048/1‑150
 ACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATT                                                                                                                                                           >  1:234389/1‑150
 ACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATT                                                                                                                                                           <  1:704578/150‑1
 ACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATT                                                                                                                                                           <  1:585983/150‑1
 ACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATT                                                                                                                                                           <  1:2821304/150‑1
 ACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATT                                                                                                                                                           <  1:2701399/150‑1
 ACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATT                                                                                                                                                           >  1:2401487/1‑150
 ACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATT                                                                                                                                                           <  2:282660/150‑1
 ACGCCGCCATCGCGCAGGGTGTCGCACAGCTTCCACAGCGTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAAGTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATT                                                                                                                                                           <  2:834826/150‑1
   GCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGACGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCA                                                                                                                                                         >  1:2715454/1‑150
   GCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCA                                                                                                                                                         <  1:1967844/150‑1
   GCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCA                                                                                                                                                         >  1:485394/1‑150
    CCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAG                                                                                                                                                        <  1:1785484/150‑1
                                   ACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTAATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACTTCATCCTTCAGGCTGCCTCTGCGTTGGCTGCGCTCGTTCACC                                                                                                                         >  1:1272100/1‑150
                                                                                         CGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGCTGCCTCTGCGTTGGCTGCGCTCGTTCACCCCGGTCACGTACTTCTGTACGCTCCCGGGGATTCACTCACTTGCCGCCTTGATG                                                                   <  1:1277763/150‑1
                                                                                             CCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACTTCATCCTTCAGGCTGCCTCTGCGTTGGCTGCGCTCGTTCACCCCGGTCACGTACTTCTGTACGCTCCCGGGGATTCACTCACTTGCCGCCTTGATGCAAC                                                               <  1:2207796/150‑1
                                                                                             CCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACTTCATCCTTCAGGCTGCCTCTGCGTTGGCTGCGCTCGTTCACCCCGGTCACGTACTTCTGTACGCTCCCGGGGATTCACTCACTTGCCGCCTTGATGCAAC                                                               <  1:1439080/150‑1
                                                                                               AAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGCTGCCTCTGCGTTGGCTGCGCTCGTTCACCCCGGTCACGTACTTCTGTACGCTCCCGGGGATTCACTCACTTGCCGCCTTGATGCAACCT                                                             <  2:993482/150‑1
                                                                                                                GAAGAAGAAATTGCCGAGGGTAATATACACTTCATCCTTCAGGCTGCCTCTGCGTTGGCTGCGCTCGTTCACCCCGGTCACGTACTTCTGTACGCTCCCGGGGATTCACTCACTTGCCGCCTTGATGCAACCTGAATGATTTTGTGTATA                                            <  2:1272100/150‑1
                                                                                                                                                      TCAGGCTGCCTCTGCGTTGGCTGCGCTCGTTCACCCCGGTCACGTACTTCTGTACGCTCCCGGGGATTCACTCACTTGCCGCCTTGATGCAACCTGAATGATTTTGTGTATATTACCCTCGGCAATTTCTTCTTCTGCGGCTCGATGAAT      <  2:1205544/150‑1
                                                                                                                                                      TCAGGCTGCCTCTGCGTTGGCTGCGCTCGTTCACCCCGGTCACGTACTTCTGTACGCTCCCGGGGATTCACTCACTTGCCGCCTTGATGCAACCTGAATGATTTTGTGTATATTACCCTCGGCAATTTCTTCTTCTGCGGCTCGATGAAT      <  1:1350123/150‑1
                                                                                                                                                          GCTGCCTCTGCGTTGGCTGCGCTCGTTCACCCCGGTCACGTACTTCTGTACGCTCCCGGGGATTCACTCACTTGCCGCCTTGATGCAACCTGAATGATTTTGTGTATATTACCCTCGGCAATTTCTTCTTCTGCGGCTCGATGAATTTGG  <  2:1336640/150‑1
                                                                                                                                                          GCTGCCTCTGCGTTGGCTGCGCTCGTTCACCCCGGTCACGTACTTCTGTACGCTCCCGGGGATTCACTCACTTGCCGCCTTGATGCAACCTGAATGATTTTGTGTATATTACCCTCGGCAATTTCTTCTTCTGCGGCTCGATGAATTTGG  <  1:935351/150‑1
                                                                                                                                                                                                                                                                                                                
AACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4582983‑4583137
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttcaggCTGCCTCTGCGTTGGCTGCGCTCGTTCACCCCGGTCACGTACTTCTGTACGCTCCCGGGGATTCACTCACTTGCCGCCTTGATGCAACCTGAATGATTTTGTGTATATTACCCTCGGCAATTTCTTCTTCTGCGGCTCGATGAATTTGG  <  NC_000913/4583224‑4583076

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