New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AM260479 = 302546611 (0.430)4 (0.170) 4/194 1.8 32.7% coding (556/579 nt) h16_A2798 conserved hypothetical protein
?AM260479 = 3025492 6 (0.250)intergenic (+3/+13) h16_A2798/h16_A2799 conserved hypothetical protein/Dinucleotide‑utilizing enzymes involved in molybdopterin and thiamine biosynthesis family 1

CGCGGATTATCCGGAGCTGATCTTCTTCGACAACCTTGCACCCGGCGACGCCGAGATCGACGGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AM260479/3025402‑3025466
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgggcCGCTCAGGCGCTCAGCGCGCGCAGCAGCCCGTCGATCTCGGCGTCGCCGGGTGCAAGGTTGTCGAAGAAGATCAGCTCCGGATAATCCGCGGGCGCCGGCATGAAGCGCCGGGTGCGGTAGGCTTCCCAGT  <  AM260479/3025492‑3025362
                                                                                                                                                                                                    
CGCGGATTATCCGGAGCTGATCTTCTTCGACAACCTTGCACCCGGCGACGCCGAGATCGACGGGCCG                                                                                                                                   <  1:611096/67‑1
CGCGGATTATCCGGAGCTGATCTTCTTCGACAACCTTGCACCCGGCGACGCCGAGATCGACGGGCCG                                                                                                                                   >  2:611096/1‑67
                                       ACCCGGCGACGCCGAGATCGACGGGCCGCTCAGGCGCT                                                                                                                         <  1:626104/38‑1
                                       ACCCGGCGACGCCGAGATCGACGGGCCGCTCAGGCGCT                                                                                                                         >  2:626104/1‑38
                                                              GGCCGCTCAGGCGCTCAGCGCGCGCAGCAGCCCGTCGATCTCGGCGTCGCCGGGTGCAAGGTTGTCGAAGAAGATCAGCTCCGGATAATCCGCGGGCGCCGGCATGAAGCGCCGGGTGCGGTAGGCTTCCCAGT  <  2:454711/134‑1
                                                               GCCGCTCAGGCGCTCAGCGCGCGCAGCAGCCCGTCGATCTCGGCGTCGCCGGGTGCAAGGTTGT                                                                       >  1:251783/1‑64
                                                               GCCGCTCAGGCGCTCAGCGCGCGCAGCAGCCCGTCGATCTCGGCGTCGCCGGGTGCAAGGTTGT                                                                       <  2:251783/64‑1
                                                               GCCGCTCAGGCGCTCAGCGCGCGCAGCAGCCCGTCGATCTCGGCGT                                                                                         >  1:117868/1‑46
                                                               GCCGCTCAGGCGCTCAGCGCGCGCAGCAGCCCGTCGATCTCGGCGT                                                                                         <  2:117868/46‑1
                                                                                                                                                                                                    
CGCGGATTATCCGGAGCTGATCTTCTTCGACAACCTTGCACCCGGCGACGCCGAGATCGACGGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AM260479/3025402‑3025466
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgggcCGCTCAGGCGCTCAGCGCGCGCAGCAGCCCGTCGATCTCGGCGTCGCCGGGTGCAAGGTTGTCGAAGAAGATCAGCTCCGGATAATCCGCGGGCGCCGGCATGAAGCGCCGGGTGCGGTAGGCTTCCCAGT  <  AM260479/3025492‑3025362

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