New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 461993882 (0.950)6 (0.070) 5/166 NT 7.6% coding (336/1224 nt) deoB phosphopentomutase
?NC_000913 = 4619964 70 (0.870)coding (362/1224 nt) deoB phosphopentomutase
Rejected: Frequency below cutoff threshold.

GCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4619835‑4619938
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgtttTCGTGATCGGAGAAATATCCCCACTCAAACAGAACCGGGACACCGGCAATTTCCCAGTGACCAGACGGGGTATCTTTACCGGATGACATTTCGTGCGCC  <  NC_000913/4619964‑4619866
                                                                                                                                                                                                           
GCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCT                                                                                                         >  1:1919322/1‑100
   accCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCAATGGGAAATTGCCGGTGTCCCGGTTCTGTT                                                                                                      <  2:1165405/97‑1
    CGCTGAAGTTATCGGCGCGTACGCATTGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTT                                                                                                     <  2:176959/100‑1
    CGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTT                                                                                                     >  1:1308473/1‑100
    CGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTT                                                                                                     <  2:341559/100‑1
    CGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTT                                                                                                     <  1:1808615/100‑1
    CGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTT                                                                                                     >  2:1732372/1‑100
            TTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATGCCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTTCGTGATC                                                                                             <  1:1260519/100‑1
                               GGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTTCGTGATCGGAGAAATATCCCCACTCA                                                                          <  1:895942/100‑1
                                    ACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTTCGTGATCGGAGAAATATCCCCACTCAAACAG                                                                     <  1:58866/100‑1
                                       AAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTTCGTGATCGGAGAAATATCCCCACTCAAACAGAAn                                                                  <  1:416628/100‑2
                                       AAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTTCGTGATCGGAGAAATATCCCCACTCAAACAGAAC                                                                  <  2:607793/100‑1
                                        AATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTTCGTGATCGGAGAAATATCCCCACTCAAACAGAACC                                                                 <  1:1641859/100‑1
                                                                                                  CTGTTTTCGTGATCGGAGAAATATCCCCACTCAAACAGAACCGGGACACCGGCAATTTCCCAGTGACCAGACGGGGTATCTTTACCGGATGACATTTCGT       <  2:357867/100‑1
                                                                                                   TGTTTTCGTGATCGGAGAAATATCCCCACTCAAACAGAACCGGGACACCGGCAATTTCCCAGTGACCAGACGGGGTATCTTTACCGGATGACATTTCGTG      >  1:894922/1‑100
                                                                                                    GTTTTCGTGATCGGAGAAATATCCCCACTCAAACAGAACCGGGACACCGGCAATTTCCCAGTGACCAGACGGGGTATCTTTACCGGATGACATTTCGTGC     <  1:1194414/100‑1
                                                                                                       TTCGTGATCGGAGAAATATCCCCACTCAAACAGAACCGGGACACCGGCAATTTCCCAGTGACCAGACGGGGTATCTTTACCGGATGACATTTCGTGCGCC  <  1:13974/100‑1
                                                                                                                                                                                                           
GCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4619835‑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.