New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4348362 =	NA (NA)	4 (0.250)	4/266	1.6	NA	coding (318/1422 nt)	PP_3820	group II intron‑encoding maturase
	?	NC_002947	4392757 =	NA (NA)					coding (302/1422 nt)	PP_3868	group II intron‑encoding maturase

TGATGGAGCGGGTGCTTGCACCGGCCAACCTCAGGCGTGCGTATCAACGCGTGGTTAGCAACAAGGGGGCTCCGGGTGCCGATGGCATGACGGTCGCTGACTTGGCGGGCTACGTGAAACAGTATTGGCCAACCCTCAAGGCCAGGTTGCTGGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/4348516‑4348362 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctggccTTGAGGGTTGGCCAATACTGTTTCACGTAGCCCGCCAAGTCAGCGA  >  NC_002947/4392757‑4392802                                                                                                                                                                                                           TGATGGAGCGGGTGCTTGCACCGGCCAACCTCAGGCGTGCGTATCAACGCGTGGTTAGCAACAAGGGGGCTCCGGGTGCCGATGGCATGACGGTCGCTGACTTGGCGGGCTACGTGAAACAGTATTGGCCAACCCTCAAGGCCAGGTTGC                                                     >  1:38313/1‑150      GAGCGGGTGCTTGCACCGGCCAACCTCAGGCGTGCGTATCAACGCGTGGTTAGCAACAAGGGGGCTCCGGGTGCCGATGGCATGACGGTCGCTGACTTGGCGGGCTACGTGAAACAGTATTGGCCAACCCTCAAGGCCAGGTTGCTGGCC                                                <  2:61601/150‑1                                                                                                                                       CTCAAGGCCAGGTTGCTGGCCTTGAGGGTTGGCCAATACTGTTTCACGTAGCCC               >  1:363693/1‑54                                                                                                                                       CTCAAGGCCAGGTTGCTGGCCTTGAGGGTTGGCCAATACTGTTTCACGTAGCCC               <  2:363693/54‑1                                                                                                                                        TCAAGGCCAGGTTGCTGGCCTTGAGGGTTGGCCAATACTGTTTCACGTAGCCCGCCAAGTCAGCGA  <  1:283826/66‑1                                                                                                                                        TCAAGGCCAGGTTGCTGGCCTTGAGGGTTGGCCAATACTGTTTCACGTAGCCCGCCAAGTCAGCGA  >  2:283826/1‑66                                                                                                                                                                                                           TGATGGAGCGGGTGCTTGCACCGGCCAACCTCAGGCGTGCGTATCAACGCGTGGTTAGCAACAAGGGGGCTCCGGGTGCCGATGGCATGACGGTCGCTGACTTGGCGGGCTACGTGAAACAGTATTGGCCAACCCTCAAGGCCAGGTTGCTGGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/4348516‑4348362 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctggccTTGAGGGTTGGCCAATACTGTTTCACGTAGCCCGCCAAGTCAGCGA  >  NC_002947/4392757‑4392802

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