New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi = 332818119 (1.100)4 (0.250) 4/252 1.8 17.2% intergenic (+98/+58) yhbY/greA RNA binding protein associated with pre‑50S ribosomal subunits/transcript cleavage factor
?NC_000913_3_pae_pgi = 3328184 21 (1.310)intergenic (+101/+55) yhbY/greA RNA binding protein associated with pre‑50S ribosomal subunits/transcript cleavage factor
Rejected: Frequency below/above cutoff threshold.

TTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_pae_pgi/3328127‑3328181
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaaaggccgcTATGCGGCCTTTTATCAACGAACAGAGCGTGGCATTTTGCTCTCCTGCCTGCGGAAAACCCCTCGTTTTACACAGCAAATGTGTGTAACTTTAGGATAATCTTAGCGTGGCAGCGAGATTTTACGTTCTTTAGTTGGGCGATAA  <  NC_000913_3_pae_pgi/3328184‑3328041
                                                                                                                                                                                                       
TTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCTATGCGGCCTTTTATCAACGAACAGAGCGTGGCATT                                                                                                              >  1:162589/1‑91
TTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCTATGCGGCCTTTTATCAACGAACAGAGCGTGGCATT                                                                                                              <  2:162589/91‑1
                   CAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCTATGCGGCCTTTTATCAACGAACAGAGCGTGGCATTTTGCTCTCCTGCCTGCGGAAAACCCCTCG                                                                                 <  1:4531/101‑1
                   CAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCTATGCGGCCTTTTATCAACGAACAGAGCGTGGCATTTTGCTCTCCTGCCTGCGGAAAACCCCTCG                                                                                 >  2:4531/1‑101
                                                 GGCCGCTATGCGGCCTTTTATCAACGAACAGAGCGTGGCATTTTGCTCTCCTGCCTGCGGAAAACCCCTCGTTTTACACAGCAAATGTGTGTAACTTTAGGATAATCTTAGCGTGGCAGCGAGATTTTACGTTCTTTAGTTGGGCGATAA  <  1:140516/150‑1
                                                 GGCCGCTATGCGGCCTTTTATCAACGAACAGAGCGTGGCATTTTGCTCTCCTGCCTGCGGAAAACCCCTCGTTTTACACAGCAAATGTGTGTAACTTTAGGATAATCTTAGCGTGGCAGCGAGATTTTACGTTCTTTAGTTGGGCGATAA  <  4:82209/150‑1
                                                                                                                                                                                                       
TTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_pae_pgi/3328127‑3328181
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaaaggccgcTATGCGGCCTTTTATCAACGAACAGAGCGTGGCATTTTGCTCTCCTGCCTGCGGAAAACCCCTCGTTTTACACAGCAAATGTGTGTAACTTTAGGATAATCTTAGCGTGGCAGCGAGATTTTACGTTCTTTAGTTGGGCGATAA  <  NC_000913_3_pae_pgi/3328184‑3328041

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