New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4601572 =119 (1.270)9 (0.110) 9/162 NT 8.1% intergenic (‑55/+52) dnaT/yjjB DNA biosynthesis protein (primosomal protein I)/DUF3815 family inner membrane protein
?NC_000913 4601608 = 95 (1.120)intergenic (‑91/+16) dnaT/yjjB DNA biosynthesis protein (primosomal protein I)/DUF3815 family inner membrane protein
Rejected: Frequency below cutoff threshold.

ATCGGTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4601677‑4601572
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tccgtgtgCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGCTG  >  NC_000913/4601608‑4601706
                                                                                                                                                                                                             
ATCGGTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTC                                                                                                           >  1:99984/1‑100
 TCGGTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACATGGATGAGAGACAGCCTCCTCTCCTCC                                                                                                          <  1:1584453/100‑1
   GGTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGT                                                                                                        >  1:1363639/1‑100
   GGTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGT                                                                                                        >  2:2250650/1‑100
   GGTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGT                                                                                                        >  2:1464577/1‑100
   GGTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGT                                                                                                        >  1:450499/1‑100
    GTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCTCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTG                                                                                                       <  2:1241232/100‑1
    GTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTG                                                                                                       <  1:605242/100‑1
    GTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTG                                                                                                       <  1:1610870/100‑1
    GTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTG                                                                                                       <  1:1639898/100‑1
    GTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTG                                                                                                       <  2:1909513/100‑1
    GTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTG                                                                                                       <  1:257558/100‑1
     TCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGACTCCTCTCCTCCGTG                                                                                                       >  1:105146/1‑99
      CTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTG                                                                                                     <  1:857732/100‑1
      CTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTG                                                                                                     >  1:2251808/1‑100
                        TGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTACGCGAGGCAATTT                                                                                    >  2:633855/1‑99
                                                           CCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCA                                                >  1:1598091/1‑100
                                                                CTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAAT                                           >  1:1194847/1‑100
                                                                   AAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGA                                        >  1:617251/1‑100
                                                                    AGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAA                                       >  2:829985/1‑100
                                                                     GCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAA                                       >  1:478385/1‑99
                                                                      CACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAG                                     >  2:2088699/1‑100
                                                                        CACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGAC                                   >  1:1256778/1‑100
                                                                          CGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCG                                 >  1:1932751/1‑100
                                                                                                   CCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGc          <  2:1031079/98‑2
                                                                                                    CGTGTGCTTAAGCGAGGCAATTTTCTACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGA       <  2:478381/100‑1
                                                                                                     GTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAA      >  2:175404/1‑100
                                                                                                     GTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAA      >  2:657754/1‑100
                                                                                                       GTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGC    <  1:949529/100‑1
                                                                                                       GTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGC    <  1:2115046/100‑1
                                                                                                       GTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGC    <  2:2047855/100‑1
                                                                                                       GTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGC    <  2:283641/100‑1
                                                                                                         GCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGCTG  >  1:491271/1‑100
                                                                                                         GCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGCTG  >  1:1611746/1‑100
                                                                                                                                                                                                             
ATCGGTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4601677‑4601572
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tccgtgtgCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGCTG  >  NC_000913/4601608‑4601706

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