New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4619938160 (1.040)4 (0.030) 4/166 NT 2.6% coding (336/1224 nt) deoB phosphopentomutase
?NC_000913 = 4619964 154 (1.080)coding (362/1224 nt) deoB phosphopentomutase
Rejected: Frequency below cutoff threshold.

GGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4619834‑4619938
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgtttTCGTGATCGGAGAAATATCCCCACTCAAACAGAACCGGGACACCGGCAATTTCCCAGTGACCAGACGGGGTATCTTTACCGGATGACATTTCGTGCGCC  <  NC_000913/4619964‑4619866
                                                                                                                                                                                                            
GGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTTGGAAATTGCCGGTGTCCCGGTTC                                                                                                          <  2:1811695/100‑1
GGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTC                                                                                                          <  2:2360859/100‑1
GGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTC                                                                                                          <  1:2834348/100‑1
 GCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCT                                                                                                         <  1:1787337/100‑1
 GCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCT                                                                                                         >  1:2057729/1‑100
 GCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCT                                                                                                         >  1:2969255/1‑100
 GCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCT                                                                                                         >  2:1344126/1‑100
  CAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTG                                                                                                        <  1:242924/100‑1
    ACGCTGAAGTTACGGGCGCGTTCGCATGGGCGCAGGAAATGTCAGTCGGTAAAGATACCCCTTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGT                                                                                                       <  1:1741221/99‑1
     CGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTT                                                                                                     >  1:1903567/1‑100
                                      CGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTTCGTGATCGGAGAAATATCCCCACTCAAACAGA                                                                    <  1:1726550/100‑1
                                         AATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTTCGTGATCGGAGAAATATCCCCACTCAAACAGAACC                                                                 <  2:397054/100‑1
                                          ATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGACATTGCCGGTGTCCCGGTTCTGTTTTCGTGATCGGAGAAATATCCCCACTCAATCAGAACCG                                                                <  1:1216882/100‑1
                                            GTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTTCGTGATCGGAGAAATATCCCCACTCAAACAGAACCGGG                                                              <  1:1799278/100‑1
                                                                                                   CTGTTTTCGTGATCGGAGAAATATCCCCACTCAAACAGAACCGGGACACCGGCAATTTCCCAGTGACCAGACGGGGTATCTTTACCGGATGACATTTCG        >  1:1212913/1‑99
                                                                                                    TGTTTTCGTGATCGGAGAAATATCCCCACTCAAACAGAACCGGGACACCGGCAATTTCCCAGTGACCAGACGGGGTATCTTTACCGGATGACATTTCGTG      >  2:375890/1‑100
                                                                                                        TTCGTGATCGGAGAAATATCCCCACTCAAACAGAACCGGGACACCGGCAATTTCCCAGTGACCAGACGGGGTATCTTTACCGGATGACATTTCGTGCGCC  >  2:2057573/1‑100
                                                                                                                                                                                                            
GGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4619834‑4619938
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgtttTCGTGATCGGAGAAATATCCCCACTCAAACAGAACCGGGACACCGGCAATTTCCCAGTGACCAGACGGGGTATCTTTACCGGATGACATTTCGTGCGCC  <  NC_000913/4619964‑4619866

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.