New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4617850 =133 (1.230)8 (0.080) 6/170 NT 6.2% coding (528/780 nt) deoC 2‑deoxyribose‑5‑phosphate aldolase, NAD(P)‑linked
?NC_000913 4617871 = 115 (1.110)coding (549/780 nt) deoC 2‑deoxyribose‑5‑phosphate aldolase, NAD(P)‑linked
Rejected: Frequency below cutoff threshold.

TCCGCAGTACGCACGCCGCCCGCCGGTTTGAAACCAACGGTTTTTTCTACGCCCATATCACGGATCACTTCCATCATGATGCGCGCGCTTTCCGGCGTCGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4617951‑4617850
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgcgCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATA  >  NC_000913/4617871‑4617966
                                                                                                                                                                                                      
TCCTCAGTACTCACTCCGCCCGCCGGTTTGAAACCAACGGTTTTTTCTACGCCCATATCACGGATCACTTCCATCATGATGCGCGCGCTTTCCGGCGTCG                                                                                                    <  1:1592382/100‑1
TCCGCAGTACGCACGCCGCCCGCCGGTTTGAAACCAACGGTTTTTTCTACGCCCATATCACGGATCACTTCCATCATGATGCGCGCGCTTTCCGGCGTCG                                                                                                    >  2:2822667/1‑100
  CGCAGTACGCACGCCGCCCGCCGGTTTGAAACCAACGGTTTTTTCTACGCCCATATCACGGATCACTTCCATCATGATGCGCGCGCTTTCCGGCGTCGC                                                                                                   >  1:209089/1‑99
  CGCAGTACGCACACAACCCGCCGGTTTGAAACAAACGGCTTTTTCTACGCCCATATCACGGATCACTTCCATCATCATGCGCGCGCTTTCCGGCGTCGC                                                                                                   >  2:1130833/1‑99
   GCAGTACGCACGCCGCCCGCCGGTTTGAAACCAACGGTTTTTTCTACGCCCATATCACAGATCACTTCCATCATGATGCGCGCGCTTTCCGGCGTCGC                                                                                                   >  1:2768880/1‑98
                                                                    TTCCATCATGATGCGCGCGCTTTCCGGCGTCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGG                                >  1:686837/1‑100
                                                                     TCCATCATGATGCGCGCGCTTTCCGGCGTCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGC                               <  2:686833/100‑1
                                                                        ATCATGATGCGCGCGCTTTCCGGCGTCGCGCATCATGATGGAAGTGATCCGTGATATGGGAGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTG                            >  1:143205/1‑100
                                                                          CATGATGCGCGCGCTTTCCGGCGTCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCG                          >  2:2460849/1‑100
                                                                           ATGATGCGCGCGCTTTCCGGCGTCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGT                         >  1:1752339/1‑100
                                                                           ATGATGCGCGCGCTTTCCGGCGTCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGT                         >  1:2320251/1‑100
                                                                            TGATGCGCGCGCTTTCCGGCGTCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTA                        >  1:1350798/1‑100
                                                                            TGATGCGCGCGCTTTCCGGCGTCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTA                        >  1:2380439/1‑100
                                                                                                  ctCGCATCATGATAGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATA  >  1:2057247/3‑100
                                                                                                  CGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATA  <  2:2775174/100‑1
                                                                                                                                                                                                      
TCCGCAGTACGCACGCCGCCCGCCGGTTTGAAACCAACGGTTTTTTCTACGCCCATATCACGGATCACTTCCATCATGATGCGCGCGCTTTCCGGCGTCGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4617951‑4617850
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgcgCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATA  >  NC_000913/4617871‑4617966

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