New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4348540 =	NA (NA)	5 (0.060)	4/238	NT	NA	coding (140/1422 nt)	PP_3820	group II intron‑encoding maturase
	?	NC_002947	4392901 =	NA (NA)					coding (158/1422 nt)	PP_3868	group II intron‑encoding maturase

CCGAACAGGTTATGCCGCCAGTAGGCGTCGCCGTCTCGTTGGTAACCGTAATGCAGAAATTTCCCACAGCGGAAACTGTCATCCCGAATCCCGGCCAGAAGCCGAGGGTGATGCCTGACAGTGCAAAGGTGCCGGCGGCGTCAGCGACGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/4348690‑4348540 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgtgTCCGGCTCCGCGTTCGTCCACGTC  >  NC_002947/4392901‑4392924                                                                                                                                                                                 CCGAACAGGTTATGCCGCCAGTAGGCGTCGCCGTCTCGTTGGTAACCGTAATGCAGAAATTTCCCACAGCGGAAACTGTCATCCCGAATCCCGGCCAGAAGCCGAGGGTGATGCCTGACAGTGCAAAGGTGCCGGCGGCGTCAGCGACGT                           <  2:1218898/150‑1                          CGTCGCCGTCTCGTTGGTAACCGTAATGCAGAAATTTCCCACAGCGGAAACTGTCATCCCGAATCCCGGCCAGAAGCCGAGGGTGATGCCTGACAGTGCAAAGGTGCCGGCGGCGTCAGCGACGTGTCCGGCTCCGCGTTCGTCCACGTC  <  1:361027/150‑1                                 GTCTCGTTGGTAACCGTAATGCAGAAATTTCCCACAGCGGAAACTGTCATCCCGAATCCCGGCCAGAAGCCGAGGGTGATGCCTGACAGTGCAAAGGTGCCGGCGGCGTCAGCGACGTG                          >  1:834499/1‑119                                 GTCTCGTTGGTAACCGTAATGCAGAAATTTCCCACAGCGGAAACTGTCATCCCGAATCCCGGCCAGAAGCCGAGGGTGATGCCTGACAGTGCAAAGGTGCCGGCGGCGTCAGCGACGTG                          <  2:834499/119‑1                                               CGTAATGCAGAAATTTCCCACAGCGGAAACTGTCATCCCGAATCCCGGCCAGAAGCCGAGGGTGATGCCTGACAGTGCAAAGGTGCCGGCGGCGTCAGCGAC                             <  1:1933651/102‑1                                               CGTAATGCAGAAATTTCCCACAGCGGAAACTGTCATCCCGAATCCCGGCCAGAAGCCGAGGGTGATGCCTGACAGTGCAAAGGTGCCGGCGGCGTCAGCGAC                             >  2:1933651/1‑102                                                                       GGAAACTGTCATCCCGAATCCCGGCCAGAAGCCGAGGGTGATGCCTGACAGTGCAAAGGTGCCGGCGGCGTCAGCGACGTGTCCGGCTCCGCGTTCGTCCACGTC  >  1:1183197/1‑105                                                                       GGAAACTGTCATCCCGAATCCCGGCCAGAAGCCGAGGGTGATGCCTGACAGTGCAAAGGTGCCGGCGGCGTCAGCGACGTGTCCGGCTCCGCGTTCGTCCACGTC  <  2:1183197/105‑1                                                                                                                 GCCTGACAGTGCAAAGGTGCCGGCGGCGTCAGCGACGTGTCCGGCTCCGCGT             >  1:1912924/1‑52                                                                                                                 GCCTGACAGTGCAAAGGTGCCGGCGGCGTCAGCGACGTGTCCGGCTCCGCGT             <  2:1912924/52‑1                                                                                                                                                                                 CCGAACAGGTTATGCCGCCAGTAGGCGTCGCCGTCTCGTTGGTAACCGTAATGCAGAAATTTCCCACAGCGGAAACTGTCATCCCGAATCCCGGCCAGAAGCCGAGGGTGATGCCTGACAGTGCAAAGGTGCCGGCGGCGTCAGCGACGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/4348690‑4348540 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgtgTCCGGCTCCGCGTTCGTCCACGTC  >  NC_002947/4392901‑4392924

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.