New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2515407 =	9 (0.170)	3 (0.060)	3/484	4.1	6.8%	noncoding (1135/1345 nt)	IS186	repeat region
	?	NC_000913	2963185 =	73 (1.420)					coding (532/564 nt)	ppdB	prepilin‑type N‑terminal cleavage/methylation domain‑containing protein PpdB
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AATTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTT‑CAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCT  >  NC_000913/2515374‑2515652                                  |                                                                                                                                                                                                                                                       aaTTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTT‑CAACCCTTAAGTTAGCGCTTATGGGATAATTcc                              <  1:298486/251‑1 (MQ=32)          cATAATCCAGCCATCGCTGGATTTCCCCCCCATAAGTGCCTGATCCGAAAAGAAGAACTAACTCGTTGTGGATAATAACAAAAATGTTCATCTGGAGCTTACAGTTGGCCATTCGTGGGACAGTATCCCTGACAGCCTGCAAAACGCAATTGAATAGCGCTAGGCATCGTCGTAACGATGTACCGATGCGTCGCATTCTGCAG‑TGGTTTCTGCCCTTAAGTTATCGCTGATGGGATAATTCCCCGGttttt                     <  1:65694/251‑1 (MQ=11)          cATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGACGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTT‑CAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGttttt                     <  1:209169/251‑1 (MQ=25)             aaTCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTT‑CAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTa                  <  2:111675/251‑1 (MQ=37)              aTCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTGGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATGGAAGAACGCGAGGCAGCGTCTTAACGAGGCACCGAGGCGTCGCTTTCTTAAGATGGTT‑CAACCCTTAAGTTAGCGCTTAGGGGTTAATTCCCCGTTTTTTTTAc                 >  2:437216/1‑251 (MQ=11)                     cATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTT‑CAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGtt          >  2:70503/1‑251 (MQ=40)                       tCGCTGGATTTCCCCCCCAGAAGAGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTACCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTT‑CAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTa         >  2:272388/1‑250 (MQ=255)                        cGCTGGAGTTCCCACCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTT‑CAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAAt       <  2:62987/251‑1 (MQ=255)                             gATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATGCTTCAGATGGTT‑CAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAAGCAGCt  <  1:383862/251‑1 (MQ=255)                                  |                                                                                                                                                                                                                                                       AATTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTT‑CAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCT  >  NC_000913/2515374‑2515652

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