New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4601572 =175 (1.160)7 (0.050) 7/162 NT 4.5% intergenic (‑55/+52) dnaT/yjjB DNA biosynthesis protein (primosomal protein I)/DUF3815 family inner membrane protein
?NC_000913 4601608 = 135 (0.990)intergenic (‑91/+16) dnaT/yjjB DNA biosynthesis protein (primosomal protein I)/DUF3815 family inner membrane protein
Rejected: Frequency below cutoff threshold.

CATCGGTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4601678‑4601572
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tccgtgtgCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGCTG  >  NC_000913/4601608‑4601706
                                                                                                                                                                                                              
CATCGGTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCT                                                                                                            >  1:2512801/1‑100
CATCGGTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCACTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCT                                                                                                            >  1:224829/1‑100
  TCGGTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCC                                                                                                          <  1:646430/100‑1
    GGTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGT                                                                                                        >  2:1700501/1‑100
     GTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTG                                                                                                       <  1:244138/100‑1
     GTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTG                                                                                                       >  1:611268/1‑100
     GTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTG                                                                                                       <  1:4032230/100‑1
       CTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTG                                                                                                     >  1:764123/1‑100
       CTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTG                                                                                                     <  1:1898226/100‑1
       CTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTG                                                                                                     <  1:2919591/100‑1
                                                                GCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAA                                            >  2:3611341/1‑100
                                                                 CTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAAT                                           >  1:2162132/1‑100
                                                                   TAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGG                                         >  1:446590/1‑100
                                                                     AGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAA                                       >  1:391832/1‑100
                                                                      GCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAA                                      >  2:1543532/1‑100
                                                                       CACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAG                                     >  2:3142825/1‑100
                                                                         CACGGATGAGAGACAGCCTCCTCTCCTCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGAC                                   >  1:1249886/1‑100
                                                                                                   TCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGAT         >  2:3315637/1‑100
                                                                                                   TCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAAC                <  2:1969156/93‑1
                                                                                                    CCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGCATGGAAAGACCGATGGATAACGCACCAACAATCGATG        >  2:2814371/1‑100
                                                                                                    CCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATG        <  1:3893834/100‑1
                                                                                                    CCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATG        >  2:236605/1‑100
                                                                                                     CGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGA       <  2:1896053/100‑1
                                                                                                     CGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGA       <  1:1613610/100‑1
                                                                                                      GTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAA      <  2:2235467/100‑1
                                                                                                      GTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAA      >  1:1272051/1‑100
                                                                                                       TGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAG     <  2:2833711/100‑1
                                                                                                        GTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGC    >  2:3139121/1‑100
                                                                                                        GTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGC    <  2:973770/100‑1
                                                                                                        GTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGC    <  1:399730/100‑1
                                                                                                        GTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGC    <  1:3801368/100‑1
                                                                                                        GTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGC    >  1:3640662/1‑100
                                                                                                        GTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGC    <  1:1547159/100‑1
                                                                                                        GTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAG     >  2:3822177/1‑99
                                                                                                         TGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGCT   >  1:1894149/1‑100
                                                                                                         TGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGAGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAG     >  2:2227696/1‑98
                                                                                                          GCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGCTG  >  2:2531643/1‑100
                                                                                                          GCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGCTG  <  2:1172881/100‑1
                                                                                                          GCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGCTG  >  1:680295/1‑100
                                                                                                          GCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGATGGATAACGCACCAACAATCGATGAAGCTG  <  2:3989325/100‑1
                                                                                                                                                                                                              
CATCGGTCTTTCCATTCCTGGATTATGGTTGTACCGCAAGCGCCCTCGCGTATAAAATTGCCTCGCTTAAGCACACGGATGAGAGACAGCCTCCTCTCCTCCGTGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4601678‑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 < 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.