New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 4893579	22 (0.690)	3 (0.120)	3/236	4.0	16.3%	intergenic (+54/‑154)	pyk/PP_4302	pyruvate kinase/urea transporter
	?	NC_002947	= 4893604	13 (0.500)					intergenic (+79/‑129)	pyk/PP_4302	pyruvate kinase/urea transporter
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ACCGCTGGCGTGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTC  >  NC_002947/4893460‑4893750                                                                                                                                                 |                                                                                                                                                   aCCGCTGGCGTGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCGGGCACAAGGGAATGCAGTc                                                                                                                                                   >  4:124383/1‑146 (MQ=255)      tGGCGTGCCTTTCGGCAAGCCGGGTACGACTAACCCGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTGGCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTgg                                                                                                                                            <  7:115989/148‑1 (MQ=255)      tGGCGTGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTgg                                                                                                                                            <  7:24763/148‑1 (MQ=255)         cGTGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAg                                                                                                                                         <  2:8911/148‑1 (MQ=255)           tGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCg                                                                                                                                       <  5:144833/148‑1 (MQ=255)                   ggCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTtgt                                                                                                                                <  7:163716/147‑1 (MQ=255)                                cGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCCCGCTGGCACAAGGGAATGCAGTCCCTGTGGGCGCGGCCTTGTGTAGCGATTGGacg                                                                                                                  >  8:148850/1‑146 (MQ=255)                                    tAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACCCAAGGCCGCTCCCCCAGGGCCCCGCAGGCACAAGGGGATGCCGGCCCTGGGGGAGCCGCCGTGTGTCGCGCTTGGCGGGcca                                                                                                              >  3:171220/1‑146 (MQ=255)                                      acacGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAACGCCGCACCCACAAGGCCACGCTGGCACAAGGGAATGCACTCACAGTGGGAGCGGCCATGTGTCGCGAATTGAGGGCAAAg                                                                                                            >  2:3438/1‑148 (MQ=255)                                                    gagaCTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCCCAGGGCCACGCTGGCCCAAGGGCATGCAGTCCCAGTGGGAGCGGCCGTGTGTCGCGCTTGGAGGGCCAAGCCCTCCCTGAATAc                                                                                              >  1:23442/1‑148 (MQ=38)                                                         tttGATCTGAGATTGCCGGTGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCACCCACAGCGCCACCCTGGCACAAGGGAATGCCGTCCCTGTGGGCGCGGGCGTGGGTCGCGCTTGGAGGGGAAACCCCACCCGGAATAACGct                                                                                          >  6:136378/1‑147 (MQ=25)                                                                                gcgTGGCGCCCTTTGCGCCACACACGGCCGCTCCCACCGGGCCCCGCTGGCAAATGGTAATGCTGCCCGTGTGGGCGGGGCCTTGTGGCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAgg                                                                  <  6:8084/148‑1 (MQ=255)                                                                                       gcCCCTTTCGCGACCAAAACCCCCTCCCACAGGGCCCCCCGGGCACCAGGGAATGCAGTCCCTGTGGGAGCCGCCGTGTGTAGCGCATTGAGGGCAAAACCAACCAGAAAAAACACTCGACCCGCCCCCCGAAAAACGCGGcccca                                                             >  2:3308/1‑144 (MQ=255)                                                                                                       aaGGCCGCTCCCACAGGGCCACGCTGGCACAAGGTAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTcccc                                           <  2:23442/148‑1 (MQ=255)                                                                                                       aaGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTccc                                            >  2:179349/1‑147 (MQ=255)                                                                                                                        gCCACGCTGGCACATGGTAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCactgac                          <  1:3438/148‑1 (MQ=255)                                                                                                                        gCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGACCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCAAAAGAGGCCCACCGCTCCCCCATCAAACTACactgac                          >  2:97723/1‑148 (MQ=255)                                                                                                                                   acaAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGCCCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa               >  5:25179/1‑148 (MQ=255)                                                                                                                                                gTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTc  <  1:179349/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   ACCGCTGGCGTGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTC  >  NC_002947/4893460‑4893750

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