New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2904081 =41 (0.630)5 (0.080) 5/246 NT 11.4% intergenic (+707/+666) ygcE/queE putative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
?NC_000913 2904099 = 39 (0.640)intergenic (+725/+648) ygcE/queE putative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis

CCGCGCCAGCGGGGATAAACCGAATAGCAATAGTCCATAGATTTGCGAAAACAGGTGTTCCCCGCGCCAGCGGGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/2904156‑2904081
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcggggaACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGC  >  NC_000913/2904099‑2904228
                                                                                                                                                                                                               
CCGCGCCAGCGGGGATAAACCGAATAGCAATAGTCCATAGATTTGCGAAAACAGGTGTTCCCCGCTGGCGCGGGGAACACCT                                                                                                                              >  1:189266/1‑82
CCGCGCCAGCGGGGATAAACCGAATAGCAATAGTCCATAGATTTGCGAAAACAGGTGTTCCCCGCTGGCGCGGGGAACACCT                                                                                                                              <  2:189266/82‑1
                                    ATAGATTTGCGAAAACAGGTGTTCCCCGCGCCAGCGGG                                                                                                                                      <  1:123533/38‑1
                                    ATAGATTTGCGAAAACAGGTGTTCCCCGCGCCAGCGGG                                                                                                                                      >  2:123533/1‑38
                                                         TTCCCCGCGCCAGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGC  >  1:433962/1‑149
                                                           CCCCGCGCCAGCGGGGAACACCTGTTTTCGCAAATCTATGGA                                                                                                           >  3:274111/1‑42
                                                           CCCCGCGCCAGCGGGGAACACCTGTTTTCGCAAATCTATGGA                                                                                                           <  4:274111/42‑1
                                                                           GAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAG                                                >  3:394545/1‑86
                                                                           GAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAG                                                <  4:394545/86‑1
                                                                                                                                                                                                               
CCGCGCCAGCGGGGATAAACCGAATAGCAATAGTCCATAGATTTGCGAAAACAGGTGTTCCCCGCGCCAGCGGGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/2904156‑2904081
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcggggaACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGC  >  NC_000913/2904099‑2904228

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