New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4832391 =	38 (0.710)	4 (0.080)	4/252	NT	10.0%	coding (800/1014 nt)	pvdY	hydroxyproline acetylase
	?	NC_002947	4832416 =	36 (0.700)					coding (775/1014 nt)	pvdY	hydroxyproline acetylase

TCCATATGCTGGTGGGAGAGGAGAGCCACCGAGGGCCGCACAAGGTGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/4832439‑4832391 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggtggcTCTCCTCTCCCACCAGCATATGGATACCACGGTCGTAGTCATCGGCCGGGTAGAAGGGGGCAATGCGGTCCTCCTTGGCCCAGTACGCCTCGAAGTAGGCAAACGGCTCGTCATCGAAACAACCGATCAACGTCAACGTATGCGGGTC  >  NC_002947/4832416‑4832563                                                                                                                                                                                                       TCCATATGCTGGTGGGAGAGGAGAGCCACCGAGGGCCGCACAAGGTGGC                                                                                                                                                      >  4:179556/1‑49                     GAGAGCCACCGAGGGCCGCACAAGGTGGCTCTCCTCTCCCACCAGCATATGGATACCACGGTCGTAGTCATCGGCCGGGTAGAAGGGGGCAATGCGGTCCTCCTTGG                                                                        >  5:304947/1‑107                     GAGAGCCACCGAGGGCCGCACAAGGTGGCTCTCCTCTCCCACCAGCATATGGATACCACGGTCGTAGTCATCGGCCGGGTAGAAGGGGGCAATGCGGTCCTCCTTGG                                                                        <  6:304947/107‑1                       GAGCCACCGAGGGCCGCACAAGGTGGCTCTCCTCTCCCACCAGCATATGGATACCACGGTCGTAGTCATCGGCCGGGTAGAAGGGGGCAAT                                                                                      >  1:170993/1‑91                       GAGCCACCGAGGGCCGCACAAGGTGGCTCTCCTCTCCCACCAGCATATGGATACCACGGTCGTAGTCATCGGCCGGGTAGAAGGGGGCAAT                                                                                      <  2:170993/91‑1                                                 CTCTCCTCTCCCACCAGCATATGGATACCACGGTCGTAGTCATCGGCCGGGTAGAAGGGGGCAATGCGGTCCTCCTTGGCCCAGTACGCCTCGAAGTAGGCAAACGGCTCGTCATCGAAACAACCGATCAACGTCAACGTATGCGGGTC  >  6:15358/1‑149                                                                                                                                                                                                       TCCATATGCTGGTGGGAGAGGAGAGCCACCGAGGGCCGCACAAGGTGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/4832439‑4832391 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggtggcTCTCCTCTCCCACCAGCATATGGATACCACGGTCGTAGTCATCGGCCGGGTAGAAGGGGGCAATGCGGTCCTCCTTGGCCCAGTACGCCTCGAAGTAGGCAAACGGCTCGTCATCGAAACAACCGATCAACGTCAACGTATGCGGGTC  >  NC_002947/4832416‑4832563

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