New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3953036 =	10 (0.600)	3 (0.220)	3/230	1.5	31.5%	intergenic (‑34/+161)	PP_3486/PP_3487	cytochrome c/hypothetical protein
	?	NC_002947	3953152 =	5 (0.370)					intergenic (‑150/+45)	PP_3486/PP_3487	cytochrome c/hypothetical protein

GACCCGGGCTTGCAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCC  >  NC_002947/3952889‑3953179                                                                                                                                                    |                                                                                                                                                gatgttgcaggtaaacgctggagcaagattcgcccggatgcaccgccagcggcgattgcccagcgcatgcagctgttgacccagccgctggccggtcggtgagtatcatgggccgcctggggaattggcctctgtggcatccgcttttgt                                                                                                                                               >  1:331691‑M1/148‑150 (MQ=255) gatgttgcaggtaaacgctggagcaagattcgcccggatgcaccgccagcggcgattgcccagcgcatgcagctgttgacccagccgctggccggtcggtgagtatcatgggccgcctggggaattggcctctgtggcatccgcttttgt                                                                                                                                               >  1:331719‑M1/148‑150 (MQ=255) gACCCGGGCTTGCAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCtgt                                                                                                                                               >  2:269965/1‑150 (MQ=255)       gcaggtaaacgctggagcaagattcgcccggatgcaccgccagcggcgattgcccagcgcatgcagctgttgacccagccgctggccggtcggtgagtatcatgggccgcctggggacttggcctctgtggcatccgctttTGTGGGAGc                                                                                                                                         <  1:216972‑M1/9‑1 (MQ=255)               aTTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGccgggg                                                                                                                                 >  1:90275/1‑146 (MQ=255)                         gCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAg                                                                                                                       <  1:5260/150‑1 (MQ=255)                             gcctggaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCt                                                                                                                   <  1:175986/145‑1 (MQ=255)                                            cgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCgg                                                                                                    <  2:276409/150‑1 (MQ=255)                                                                      gTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAagag                                                                          <  2:203598/150‑1 (MQ=255)                                                                                  gTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTTCGGCGCGAAAGGGCAGCGCAGCAGCCCCCCCAGTGGCAGCTCAAGCGCAGGAAGAGGGGGGCCCTggg                                                              >  2:155300/1‑148 (MQ=255)                                                                                            cgggggggtgtttatcttgggccgcctgggtatttggcctctggggcatccgctttTGGGGGCGCGGCCTTGCGTCGCCAACGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTc                                                    >  1:322922‑M1/57‑150 (MQ=255)                                                                                             gggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCtctct                                                   <  1:269965/150‑2 (MQ=255)                                                                                             cggtcggtgagtatcatgggccgcctggggaattggcctctgtggcatccgctttTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTcg                                                   >  2:216972‑M1/56‑150 (MQ=255)                                                                                                tcggtgagtatcatgggccgcctggggaattggcctctgtggcatccgctttTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGa                                                >  1:283723‑M1/53‑150 (MQ=255)                                                                                                                                         tgcGCACGCCGGGTGGGTGGGGCCTTGCGGCGCCATGGGGGGGCGCCGCCCCCCCGCCTGTGGCAGCTCACGCGAATAATGTGGGTCGGGGGGCGGCCGCCTTCCGCGACGAATGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCaa       <  2:277829/148‑1 (MQ=255)                                                                                                                                              ccgCCCTGTGGGCGCGGCCTTGCCTCGCGAACGGGCTGCGCCGCAGCCCCGGCCGTGGCAGCACACGCGCAGAAAGCGGGGAGCGCTGCGCGCCCCTTTCGCGACGCATGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTccc  <  2:377027/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                GACCCGGGCTTGCAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCC  >  NC_002947/3952889‑3953179

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