New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4601057 =173 (1.270)9 (0.070) 8/168 NT 5.3% coding (461/540 nt) dnaT DNA biosynthesis protein (primosomal protein I)
?NC_000913 4601089 = 157 (1.220)coding (429/540 nt) dnaT DNA biosynthesis protein (primosomal protein I)
Rejected: Frequency below cutoff threshold.

GAACTGGCCTCATTCATTGCCTACTGGCAGGCGGAAGGTAAAGTCTTTCACCATGTGCAGTGGCAACAAAAACTGGCGCGCAGCCTGCAAATCGGTCGTGCCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4601160‑4601057
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccagTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGCCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCG  >  NC_000913/4601089‑4601186
                                                                                                                                                                                                          
GAACTGGCCTCATTCATTGCCTACTGGCAGGCGGAAGGTAAAGTCTTTCACCATGTGCAGTGGCAACAAAAACTGGCGCGCAGCCTGCAAATCGGTCGTG                                                                                                        <  1:664239/100‑1
GAACTGGCCTCATTCATTGCCTACTGGCAGGCGGAAGGTAAAGTCTTTCACCATGTGCAGTGGCAACAAAAACTGGCGCGCAGCCTGCAAATCGGTCGTG                                                                                                        >  1:3006648/1‑100
 AACTGGCCTCATTCATTGCCTACTGGCAGGCGGAAGGTAAAGTCTTTCACCATGTGCAGTGGCAACAAAAACTGGCGCGCAGCCTGCAAATCGGTCGTGC                                                                                                       >  2:1410999/1‑100
  ACTGGCCTCATTCATTGCCTACTGGCAGGCGGAAGGTAAAGTCTTTCACCATGTGCAGTGGCAACAAAAACTTGCGCGCAGCCTGCAAATCGGTCGTGCC                                                                                                      <  1:81733/100‑1
  ACTGGCCTCATTCATTGCCTACTGGCAGGCGGAAGGTAAAGTCTTTCACCATGTGCAGTGGCAACAAAAACTGGCGCGCAGCCTGCAAATCGGTCGTGCC                                                                                                      <  1:471203/100‑1
   CTGGCCTCATTCATTGCCTACTGGCAGGCGGAAGGTAAAGTCTTTCACCATGTGCAGTGGCAACAAAAACTGGCGCGCAGCCTGCAAATCGGTCGTGCCA                                                                                                     >  1:1084273/1‑100
   CTGGCCTCATTCATTGCCTACTGGCAGGCGGAAGGTAAAGTCTTTCACCATGTGCAGTGACAACAAAAAATGGCGCGCAGACTGCAAATCGGTAGTGCCA                                                                                                     >  2:2268699/1‑100
    TGGCCTCATTCATTGCCTACTGGCAGGCGGAAGGTAAAGTCTTTCACCATGTGCAGTGGCAACAAAAACTGGCGCGCAGCCTGCAAATCGGTCGTGCCAG                                                                                                    <  2:657244/100‑1
    TGGCCTCATTCATTGCCTACTGGCAGGCGGAAGGTAAAGTCTTTCACCATGTGCAGTGGCAACAAAAACTGGCGCGCAGCCTGCAAATCGGTCGTGCCAG                                                                                                    <  2:3489774/100‑1
    TGGCCTCATTCATTGCCTACTGGCAGGCGGAAGGTAAAGTCTTTCACCATGTGCAGTGGCAACAAAAACTGGCGCGCAGCCTGCAAATCGGTCGTGCCAG                                                                                                    <  2:3347418/100‑1
                                            CTTTCACCATGTGCAGTGGCAACAAAAACTGGCGCGCAGCCTGCAAATCGGTCGTGCCAGTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCC                                                            >  1:639659/1‑100
                                                             GGCAACAAAAACTGGCGCGCAGCCTGCAAATCGGTCGTGCCAGTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAAT                                           >  1:3753828/1‑100
                                                               CAACAAAAACTGGCGCGCAGCCTGCAAATCGGTCGTGCCAGTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGA                                         >  1:753522/1‑100
                                                                AACAAAAACTGGCGCGCAGCCTGCAAATCGGTCGTGCCAGTTTTTTTTTTCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGA                                         >  2:279616/1‑99
                                                                 ACAAAAACTGGCGCGCAGCCTGCAAATCGGTCGTGCCAGTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAAT                                       >  2:3455236/1‑100
                                                                 ACAAAAACTGGCGCGCAGCCTGCAAATCGGTCGTGCCAGTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAAT                                       >  2:3062318/1‑100
                                                                  CAAAAACTGGCGCGCAGCCTGCAAATCGGTCGTGCCAGTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAAT                                       >  1:887767/1‑99
                                                                    AAAACTGGCGCGCAGCCTGCAAATCGGTCGTGCCAGTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAG                                    >  1:504017/1‑100
                                                                          GGCGCGCAGCCTGCAAATCGGTCGTGCCAGTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGCCAGT                              >  1:627493/1‑100
                                                                                                cgcgcagGTTATTGTGGCCACTGCACATGGTGAAAGACTTTACCTTCCGACTGCAAGTAGGCAATGAATGAGGCCAGTTCTTCGGTGGTCACCGGCTCTC        <  2:1022033/93‑1
                                                                                                   GCCAGTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCGGTAGGCAATGAATGAGCCCAGTTCTTCGGTGGTCACGGGCTCGCTTA     >  2:2834868/1‑100
                                                                                                   GCCAGTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGCCAGTTCTTCGGTGGTCACCGGCTCTCTTA     >  1:2737867/1‑100
                                                                                                   GCCAGTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGCCAGTTCTTCGGTGGTCACCGGCTCTCTTA     >  2:2184043/1‑100
                                                                                                    CCAGTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGCCAGTTCTTCGGTGGTCACCGGCTCTCTTAG    <  1:2135517/100‑1
                                                                                                     CAGTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGCCAGTTCTTCGGTGGTCACCGGCTCTCTTAGC   <  1:1070564/100‑1
                                                                                                      AGTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGCCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCG  >  1:2074100/1‑100
                                                                                                                                                                                                          
GAACTGGCCTCATTCATTGCCTACTGGCAGGCGGAAGGTAAAGTCTTTCACCATGTGCAGTGGCAACAAAAACTGGCGCGCAGCCTGCAAATCGGTCGTGCCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4601160‑4601057
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccagTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGCCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCG  >  NC_000913/4601089‑4601186

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