New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4077584 =	12 (0.030)	49 (0.220) +CTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCA	27/182	12.4	39.6%	intergenic (‑204/+237)	PP_3586/tpx	transposase/2‑Cys peroxiredoxin
	?	NC_002947	5547596 =	216 (0.620)					intergenic (‑766/+111)	rpsF/rlmB	30S ribosomal protein S6/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Rejected: Coverage evenness skew score above cutoff.

CCTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATGTCACAGCCTTT  >  NC_002947/4077573‑4077733            |                                                                                                                                                      ccTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATg             >  1:835891/1‑150 (MQ=32) ccTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATg             >  2:3901523/1‑150 (MQ=32) ccTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATg             >  2:6527056/1‑150 (MQ=32)   tAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATGTc           >  1:7086001/1‑150 (MQ=34)   tAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATGTc           >  1:7714907/1‑150 (MQ=34)    aaCAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATGTca          >  1:3217882/1‑150 (MQ=35)      cAAAGCAAACGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATGTcaca        <  2:3457768/150‑1 (MQ=37)      cAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATGTcaca        <  2:1514125/150‑1 (MQ=37)      cAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATGTcaca        <  2:3457470/150‑1 (MQ=37)      cAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATGTcaca        <  2:5490795/150‑1 (MQ=37)      cAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCATTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATGTcaca        <  1:3998642/150‑1 (MQ=37)           cAAACGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATGTCACAGCCtt   <  2:1745053/150‑1 (MQ=255)            aaacGGCGGCTGGCCCTACCGGGCCAACCGCATTTTCTTTTACAAAAGGGCCGCGTCGCTGCGGAACTGAAGCCGAACACCTCATGTTGGCGGCGCCACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATGTCACAGCCttt  <  1:1591636/150‑1 (MQ=14)            aaaaGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATGTCACAGCCttt  >  1:2456055/1‑150 (MQ=255)            aaaaGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATGTCACAGCCttt  <  2:7798804/150‑1 (MQ=255)            aaaaCGCGGCTCGCCCTACCGCGCCAACCGCCTTTTCTTTTTCAAAATCACAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTCCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATGTCACAGCCttt  <  2:7628638/150‑1 (MQ=14)            |                                                                                                                                                      CCTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGGATGTCACAGCCTTT  >  NC_002947/4077573‑4077733

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