New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 413620	12 (0.720)	3 (0.190)	3/272	1.8	20.4%	coding (2571/2934 nt)	glnE	glutamate‑ammonia‑ligase adenylyltransferase
	?	NC_002947	= 413635	12 (0.760)					coding (2586/2934 nt)	glnE	glutamate‑ammonia‑ligase adenylyltransferase

CGTGTGCTGGTGGGTTGCAAGCAGGTGGGGGCGGCATTCGAAGGTGTGCGCGCCAAGGTACTGGGCCAGGCCCGCGACCTGGAAAAGCTGCGCGGCGAGGTGAGTGAAATGCGCGCCAAGATGCGTGACAACCTCGGTACCAAGGCGACGGCTGCCGGTACCGCGGCGAACGCCTTCGACGCCGGTGTACCCTTCGATATCAAGCAGGATGCCGGCGGTATCGTCGATATCGAATTTATGGTGCAATACGCCGCTTTGGCCTGGTCCCATGACCATCCGGCCATACTTCG  >  NC_002947/413480‑413769                                                                                                                                             |                                                                                                                                                      cGTGTGCTGGTGGGTTGCAAGCAGGTGGGGGCGGCATTCGAAGGTGTGCGCGCCAAGGTACTGGGCCAGGCCCGCGACCTGGAAAAGCTGCGCGGCGAGGTGAGTGAAATGCGCGCCAAGATGCGTGACAACCTCGGTACCAAGGCGACg                                                                                                                                              <  1:31469/150‑1 (MQ=255)        tggtggGTTGCAAGCAGGTGGGGGCGGCATTCGAAGGTGTGCGCGCCAAGGTACTGGGCCAGGCCCGCGACCTGGAAAAGCTGCGCGGCGAGGTGAGTGAAATGCGCGCCAAGATGCGTGACAACCTCGGTACCggcagccgtcgccttg                                                                                                                                       <  1:27781‑M1/150‑17 (MQ=255)                               gcggcATTCGAAGGTGTGCGCGCCAAGGTCCTGGGCCAGGCCCGCGACCTGGAAAAGCTGCGCGGCGAGGTGAGTGAAATGCGCGCCAAGATGCGTGACAACCTCGGTACCAAGGCGACGGCTGCCGGTACCGCGGCGAACGCCTTCGAc                                                                                                                <  2:361658/150‑1 (MQ=255)                                                        aGGTACTGGGCCAGGCCCGCGACCTGGAAAAGCTGCGCGGCGAGGTGAGTGAAATGCGCGCCAAGATGCGTGACAACCTCGGTACCAAGGCGACGGCTGCCGGTACCGCGGCGAACGCCTTCGACGCCGGTGTACCCTTCGATATCAAGc                                                                                       >  1:99898/1‑150 (MQ=255)                                                         ggTACTGGGCCGGGCCCGCGACCTGGAAAAGCTGCGCGGCGAGGTGAGTGAAATGCGCGCCAAGATGCGTGACAACCTCGGTACCAAGGCGACGGCTGCCGGTACCGCGGCGAACGCCTTCGACGCCGGTGTACCCTTCGATATCAAGCa                                                                                      <  1:20392/150‑1 (MQ=255)                                                              tGGGCCAGGCCCGCGACCTGGAAAAGCTGCGCGGCGAGGTGAGTGAAATGCGCGCCAAGATGCGTGACAACCTCGGTACCggcagccgtcgccttggtac                                                                                                                                   >  1:277622‑M1/1‑80 (MQ=255)                                                              tGGGCCAGGCCCGCGACCTGGAAAAGCTGCGCGGCGAGGTGAGTGAAATGCGCGCCAAGATGCGTGACAACCTCGGTACCggcagccgtcgccttggtac                                                                                                                                   <  2:277622‑M1/100‑21 (MQ=255)                                                              tGGGCCAGGCCCGCGACCTGGAAAAGCTGCGCGGCGAGGTGAGTGAAATGCGCGCCAAGATGCGTGACAACCTCGGTACCAAGGCGACGGCTGCCGGTACCGCGGCGAACGCCTTCGACGCCGGTGTACCCTTCGATATCAAGCAGGATg                                                                                 >  2:247478/1‑150 (MQ=255)                                                                         cgcgACCTGGAAAAGCTGCGCGGCGAGGTGAGTGAAATGCGCGCCAAGATGCGTGACAACCTCGGTACCAAGGCGACGGCTGCCGGTACCGCGGCGAACGCCTTCGACGCCGGTGTACCCTTCGATATCAAGCAGGATGCCGGCGGTAtc                                                                      <  1:247478/150‑1 (MQ=255)                                                                            gACCTGGAAAAGCTGCGCGGCGAGGTGAGTGAAATGCGCGCCAAGATGCGTGACAACCTCGGTACCAAGGCGACGGCTGCCGGTACCGCGGCGAACGCCTTCGACGCCGGTGTACCCTTCGATATCAAGCAGGATGCCGGCGGTAtcgtc                                                                   <  2:354990/150‑1 (MQ=255)                                                                            gACCTGGAAAAGCTGCGCGGCGAGGTGAGTGAAATGCGCGCCAAGATGCGTGACAACCTCGGTACCAAGGCGACGGCTGCCGGTACCGCGGCGAACGCCTTCGACGCCGGTGTACCCTTCGATATCAAGCAGGATGCCGGCGGTAtcgtc                                                                   <  2:240281/150‑1 (MQ=255)                                                                                           cgcgGCGAGGTGAGTGAACTGCGCGCCAAGAGGCGGGACAACCTCGGTCCCAAGGCGACGGCTGCCGGTACCGCGGCGAACGCCTTCGACGCCGGTGTACCCTTCGATATCACGCAGGATGCCGGCGGTATCGTCGATATCGAATTTCTg                                                    <  1:66821/150‑1 (MQ=255)                                                                                                           aaaTGCGCGCCAAGATGCGTGACAACCTCGGTACCAAGGCGACGGCTGCCGGTACCGCGGCGAACGCCTTCGACGCCGGTGTACCCTTCGATATCAAGCAGGATGCCGGCGGTATCGTCGATATCGAATTTATGGTGCAATACGCCGCtt                                    >  1:71583/1‑150 (MQ=255)                                                                                                                cgcgcCAAGATGCGTGACAACCTCGGTACCAAGGCGACGGCTGCCGGTACCGCGGCGAACGCCTACGACGCCGGTGTACCCTTCGATAACAAGCAGGATGCCGGCGGTATCGTCGATATCGAATTTATCGTGCAATACGCCGCTTTGccc                               <  2:226074/150‑1 (MQ=255)                                                                                                                         aTGCGTGACAACCTCGGTACCAAGGCGACGGCTGCCGGTACCGCGGCGAACGCCTTCGACGCCGGTGTACCCTTCGATATCAAGCAGGATGCCGGCGGTATCGTCGATATCGAATTTATGGTGCAATACGCCGCTTTGGCCTGGTCCCAt                      <  1:230126/150‑1 (MQ=255)                                                                                                                                         gTACCAAGGCGACGGCTGCCggtaccgaggttgtcacgcatcttggcgcgcatttcactcacctcgccgcgcagcttttccaggtcgcgggcctggccca                                                        <  2:277622‑M2/100‑81 (MQ=255)                                                                                                                                         gTACCAAGGCGACGGCTGCCggtaccgaggttgtcacgcatcttggcgcgcatttcactcacctcgccgcgcagcttttccaggtcgcgggcctggccca                                                        >  1:277622‑M2/1‑20 (MQ=255)                                                                                                                                             cAAGGCGACGGCTGCCggtaccgaggttgtcacgcatcttggcgcgcatttcactcacctcgccgcgcagcttttccaggtcgcgggcctggcccagtaccttggcgcgcacaccttcgaatgccgcccccacctgcttgcaacccacca  <  1:27781‑M2/150‑135 (MQ=255)                                                                                                                                             |                                                                                                                                                      CGTGTGCTGGTGGGTTGCAAGCAGGTGGGGGCGGCATTCGAAGGTGTGCGCGCCAAGGTACTGGGCCAGGCCCGCGACCTGGAAAAGCTGCGCGGCGAGGTGAGTGAAATGCGCGCCAAGATGCGTGACAACCTCGGTACCAAGGCGACGGCTGCCGGTACCGCGGCGAACGCCTTCGACGCCGGTGTACCCTTCGATATCAAGCAGGATGCCGGCGGTATCGTCGATATCGAATTTATGGTGCAATACGCCGCTTTGGCCTGGTCCCATGACCATCCGGCCATACTTCG  >  NC_002947/413480‑413769

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