New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 624576 =147 (1.060)4 (0.040)
+TGTCTGAGCACTTTTGTACTTTGTCA		 4/174 NT 5.6% intergenic (‑66/‑309) fepB/entC ferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1
?NC_000913 3525419 = 28 (0.200)coding (2549/2553 nt) mrcA penicillin‑binding protein 1a, murein transglycosylase and transpeptidase

CCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/624612‑624576
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TCTGATTAAAAAGGCGCTTCGGCGCCTTT  >  NC_000913/3525419‑3525447
                                     ||||||||||||||||||||||||||                             
CCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCTTCGGCGCCTT   >  1:619813/1‑91
CCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCTTCGGCGCCTT   <  2:619813/91‑1
                TCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCTTCGGCGCCTTT  <  3:143959/76‑1
                TCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCTTCGGCGCCTTT  >  4:143959/1‑76
                                              ACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCTTCGGCGCCTT   >  3:1415708/1‑45
                                              ACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCTTCGGCGCCTT   <  4:1415708/45‑1
                                     ||||||||||||||||||||||||||                             
CCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/624612‑624576
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TCTGATTAAAAAGGCGCTTCGGCGCCTTT  >  NC_000913/3525419‑3525447

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