New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4893598 =	18 (1.380)	4 (0.380)	4/220	0.9	25.4%	intergenic (+28/‑180)	pyk/PP_4302	pyruvate kinase/urea transporter
	?	NC_002947	4893677 =	9 (0.860)					intergenic (+107/‑101)	pyk/PP_4302	pyruvate kinase/urea transporter

TACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTC  >  NC_002947/4893534‑4893795                                                                                                                                                |                                                                                                                       ctctttgggcctcggttttccgcggtcgggtagagcggtattcagggagggctttgccctccaaTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAgg                                                                                                                  >  1:71798‑M1/65‑150 (MQ=255)           acGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGGCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGcccc                                                                                                        >  2:227199/1‑149 (MQ=255)              cTGCGGATCGAGACTATGTTCTGAGGTTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGCGCGGCCTTGTGTCGCGATTGGAGGGCAAAGccctccc                                                                                                     <  2:254174/150‑1 (MQ=255)                                                  cgcaacgcgcCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGATGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGc                                                                 <  2:264490/144‑1 (MQ=255)                                                     gTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGa                                                                    <  2:78405/144‑1 (MQ=255)                                                     gTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGa                                                                    >  1:78405/1‑144 (MQ=255)                                                          cctccaatcgcgacacaaggccgctcccacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAagag                                                         >  1:71798‑M2/87‑150 (MQ=255)                                                           ccccTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTACCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAgg                                                        <  2:139350/150‑1 (MQ=255)                                                                        ACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCaaa                                                            >  1:166867‑M1/1‑133 (MQ=255)                                                                        ACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCaaa                                                            <  2:166867‑M1/133‑1 (MQ=255)                                                                        ACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTcccc                                           >  2:99061‑M1/1‑150 (MQ=255)                                                                                  TCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTc                                 <  2:145033‑M1/150‑1 (MQ=255)                                                                                     cacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                              <  2:71798‑M2/91‑1 (MQ=255)                                                                                      acagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGccctcc                                                                                                      >  1:36192‑M2/59‑77 (MQ=255)                                                                                      acagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGccctcc                                                                                                      <  2:36192‑M2/19‑1 (MQ=255)                                                                                                    gCACAATGGAATGTAGTCCCTGTGGTAGCGGCCTTGTGTGCCTATTGTATTGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAACCCGATGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa               <  2:30731/150‑1 (MQ=255)                                                                                                    gCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa               <  1:145033/150‑1 (MQ=255)                                                                                                                 cAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTc  >  1:217813/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                       TACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTC  >  NC_002947/4893534‑4893795

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