New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	3423817 =	0 (0.000)	11 (0.280)	7/274	4.2	100%	intergenic (‑30/+63)	rrfD/rrlD	5S ribosomal RNA of rrnD operon/23S ribosomal RNA of rrnD operon
	?	NC_000913_3_pae_tpiA	= 4040499	NA (NA)					intergenic (+76/‑18)	rrlA/rrfA	23S ribosomal RNA of rrnA operon/5S ribosomal RNA of rrnA operon
Rejected: Coverage evenness skew score above cutoff.

AGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAAT  >  NC_000913_3_pae_tpiA/4040465‑4040641                                   |                                                                                                                                               aGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGc                              <  2:281610/149‑1 (MQ=32)     tGATACAGATTAAATCATAACGCAGAAGCGGGCTGATAAAACAGAATTTTCCTGGCGGCAGTAGCGCGGTGGCCCCACCAGACCCCATCAAGCCCTCAGACGTCAACCGCCGTATCCATCATCGTAGTGGCGCGTCTGCTCATGCgaga                          <  2:430514/149‑1 (MQ=2)          cAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAgg                     <  2:594839/149‑1 (MQ=32)           aGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAAATCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGGGGGGTCTACTCATGCGAGCGTAggg                    >  1:292732/1‑149 (MQ=255)             aTTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGaa                  >  1:555643/1‑149 (MQ=32)                 aagcAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGc              <  2:555643/146‑1 (MQ=17)                             gAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAAt  <  2:405150/149‑1 (MQ=32)                                   gTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACt                <  1:99902/129‑1 (MQ=33)                                   gTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACt                >  2:99902/1‑129 (MQ=32)                                   |                                                                                                                                               AGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAAT  >  NC_000913_3_pae_tpiA/4040465‑4040641

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