New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	1444956 =	NA (NA)	6 (0.570)	4/254	0.6	100%	intergenic (‑331/‑328)	yveA/PP_1260	aspartate‑proton symporter/transposase
	?	NC_002947	= 1849442	0 (0.000)					coding (343/2754 nt)	gacS	sensor protein GacS

AGCCCAGCAGGCGCTCGGCTTCGGCCGGCACGGCGTCGGTGGCCAGGTCGCGGTGGTGGCCGAAAACTGGCATCAGGTAGCGGGTGGCGTCGTTGCCACTGCGCTGCAGGAGTTGGGTGCCGGTACCGCCGCTGGGCGGCAGATTGAGCATGCTCGGGCCGGCGTGGGCCAGGCGCGTGCGGTCCGGGGCCAGAAACGCCACTGCGCGCACATCGGCCTGCTCCAGGGTTTGTGCGGCAATTCT  >  NC_002947/1849331‑1849574                                                                                                                |                                                                                                                                     ccccaagcccggccggccggggcttggctccgcaacgcccggaaactggggtatcaactggtggccgcttgagtgcttagaaaacaatgggttgcattttgtttgatgaGAGTTGGGTGCCGGTACCGCCGCTggg                                                                                                              >  2:195440‑M1/110‑136 (MQ=255)         ccggccggtcgggggttggcttcgcaaagcccggaaactggggtatcaactggtggccgcttgagtgcttagaaaacaatgggttgcattttgtttgatgaGAGTTGGGTGCCGGTACCGCCGCTGGGCGGCAGAt                                                                                                      >  1:227500‑M1/102‑136 (MQ=255)         agGCGCTCGGCTTCGGCCGGCACGGCGTCGGTGGCCAGGTCGCGGTGGTGGCCGAAAACTGGCATCAGGTAGCGGGTGGCGTCGTTGCCACTGCGCTGCAGGAGatcacccagcccttccgggctgcgctgtcgcg                                                                                                      >  1:287912‑M2/1‑104 (MQ=255)           gcgcTCGGCTTCGGCCGGCACGGCGTCGGTGGCCAGGTCGCGGTGGTGGCCGAAAACTGGCATCAGGTAGCGGGTGGCGTCGTTGCCACTGCGCTGCAGGAGctcacccagcccttccggggggggcggtcgcgcg                                                                                                    <  1:181456‑M2/136‑35 (MQ=255)           gcgcTCGGCTTCGGCCGGCACGGCGTCGGTGGCCAGGTCGCGGTGGTGGCCGAAAACTGGCATCAGGTAGCGGGTGGCGTCGTTGCCACTGCGCTGCAGGAGatcacccagcccttccgggctgcgctgtcgcgca                                                                                                    <  2:112851‑M2/136‑35 (MQ=255)                                          gaaactggggtatcaactggtggccgcttgagtgcttagaaaacaatgggttgcattttgtttgatgaGAGTTGGGTGCCGGTACCGCCGCTGGGCGGCAGATTGAGCATGCTCGGGCCGGCGTGGGCCAGGcgcg                                                                     >  2:68775‑M1/69‑136 (MQ=255)                                                    tatcaactggtggccgcttgagtgcttagaaaacaatgggttgcattttgtttgatgaGAGTTGGGTGCCGGTACCGCCGCTGGGCGGCAGATTGAGCATGCTCGGGCCGGCGTGGGCCAGGCGCGTGCGGTCCgg                                                           <  1:234142‑M1/78‑1 (MQ=255)                                                           tggtggccgcttgagtgcttagaaaacaatgggttgcattttgtttgatgaGAGTTGGGTGCCGGTACCGCCGCTGGGCGGCAGATTGAGCATGCTCGGGCCGGCGTGGGCCAGGCGCGTGCGGTCCGGGGCCAGa                                                    >  2:175886‑M1/52‑136 (MQ=255)                                                                             gTAGCGGGTGGCGTCGTTGCCACTGCGCTGCAGGAGatcacccagcccttccgggctgcgctgtcgcgcagcaaacatgaggtgtacggcttcagttccgcagcgacgctgacattttgaaaaagaaaaggcggtt                                  >  1:112851‑M2/1‑36 (MQ=255)                                                                                         gggttgcattttgtttgatgaGAGTTGGGTGCCGGTACCGCCGCTGGGCGGCAGATTGAGCATGCTCGGGCCGGCGTGGGCCAGGCGCGTGCGGTCCGGGGCCAGAAACGCCACTGCGCGCACATCGGCCTGCTcc                      <  2:235508‑M1/115‑1 (MQ=255)                                                                                                    ggcgctGCTGGTGaccccccagcccttccgggctgcggtgtcgcggagcaaacatgaggtgtacggcttcagttccgcagcgacggtgaaattttgaaaaagaaaaggcggtgggcccggtagggccagccgcctt           >  2:4643‑M2/2‑13 (MQ=255)                                                                                                        tgatgaGAGGTGGGTGCCGGTACCGCCGCTGGGCGGCAGATTGAGCATGCTCGGGCCGGCGGGGGCCAGGCGCGTGCTGTCCGGGGCCAGAAACGCCACTGCGCGCACATCGGCCTGCTCCAGGGTTTGTGCGGCa       >  1:205237‑M1/7‑136 (MQ=255)                                                                                                          gCAGGAGatcaccccgcccttccgggctgcgctgccgcgcagcaaacatgaggtgtacggcttcagttccgcagcgacgctgacattttgaaaaagaaacggcggtgggcccggtagggccagccgccttttgctt     >  1:219280‑M2/1‑7 (MQ=255)                                                                                                             aGAGTTGGGTGCCGGTACCGCCGCTGGGCGGCAGATTGAGCATGCTCGGGCCGGCGTGGGCCAGGCGCGTGCGGTCCGGGGCCAGAAACGCCACTGCGCGCACATCGGCCTGCTCCAGGGTTTGTGCGGCAATtct  <  1:141936‑M1/135‑1 (MQ=255)                                                                                                                |                                                                                                                                     AGCCCAGCAGGCGCTCGGCTTCGGCCGGCACGGCGTCGGTGGCCAGGTCGCGGTGGTGGCCGAAAACTGGCATCAGGTAGCGGGTGGCGTCGTTGCCACTGCGCTGCAGGAGTTGGGTGCCGGTACCGCCGCTGGGCGGCAGATTGAGCATGCTCGGGCCGGCGTGGGCCAGGCGCGTGCGGTCCGGGGCCAGAAACGCCACTGCGCGCACATCGGCCTGCTCCAGGGTTTGTGCGGCAATTCT  >  NC_002947/1849331‑1849574

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