New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4601572 =95 (1.100)4 (0.050) 4/162 NT 4.6% intergenic (‑55/+52) dnaT/yjjB DNA biosynthesis protein (primosomal protein I)/DUF3815 family inner membrane protein
?NC_000913 4601608 = 79 (1.000)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
                                                                                                                                                                                                             
ATCGGTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGTCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTC                                                                                                           <  1:500998/100‑1
 TCGGTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCC                                                                                                          <  1:185550/100‑1
  CGGTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCG                                                                                                         <  2:504464/100‑1
   GGTCTTTCCATTCCTGGATTATGGATGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGATACAGCCTCCTCTCCCCCG                                                                                                         >  2:775086/1‑99
    GTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTG                                                                                                       <  2:20865/100‑1
    GTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTG                                                                                                       <  2:2215410/100‑1
      CTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTG                                                                                                     <  1:1862464/100‑1
                                  GCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGA                                                                           >  1:532862/1‑98
                                                                CTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGTTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAAT                                           >  2:1694365/1‑100
                                                                    AGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAA                                       >  1:1310504/1‑100
                                                                      CACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAG                                     >  1:786236/1‑100
                                                                                                   CCGTGTGCTTAAGCGAGGCAATTTTATACGCGCGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATG        >  2:812892/1‑100
                                                                                                   CCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATG        >  2:2273030/1‑100
                                                                                                   CCGTGTGCTTAAGCGAGGCAATTATATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATG        >  2:494886/1‑100
                                                                                                    CGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGA       >  2:2036408/1‑100
                                                                                                      TGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAG     <  2:424396/100‑1
                                                                                                       GTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGC    >  2:1882432/1‑100
                                                                                                       GTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAG     <  1:479674/99‑1
                                                                                                         GCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGCTG  <  1:2171964/100‑1
                                                                                                         GCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCAAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGCTG  >  2:861130/1‑100
                                                                                                         GCTTAAGCGAGGCAATTTTATACACAAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGAGAACGCACCAACAATCGAAGAAGCTG  >  2:669724/1‑100
                                                                                                         GCTGAAGCGAGGCAATATTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATCGAAAGACCGCTGGATAACGCACCAACAATCGATGAAGC    >  2:228455/1‑98
                                                                                                                                                                                                             
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 < 12 ≤ ATCG/ATCG < 22 ≤ 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.