New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	= 3946053	8 (0.040)	45 (0.230)	27/250	15.8	85.3%	noncoding (2350/2904 nt)	rrlC	23S ribosomal RNA of rrnC operon
	?	NC_000913_3_pae_pgi	= 3946069	NA (NA)					noncoding (2366/2904 nt)	rrlC	23S ribosomal RNA of rrnC operon
Rejected: Coverage evenness skew score above cutoff.

TGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTGGGGCGGTCTCCTCCTAAAGAGTAACGGAGGAGCACGAAGGTTGGCTAATCCTGGTCGGACAT‑CAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGAC  >  NC_000913_3_pae_pgi/3945905‑3946062                                                                                                                                                      |          tGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTGGGGCGGTCTCCTCCTAAAGAGTAACGGAGGAGCACGAAGGTTGGCTAATCCTGGTCGGACAT‑CAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCg          <  3:1407938/150‑1 (MQ=33) tGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTGGGGCGGTCTCCTCCTAAAGAGTAACGGAGGAGCACGAAGGTTGGCTAATCCTGGTCGGACAT‑CAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCg          >  3:243207/1‑150 (MQ=33)  ggACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTGGGGCGGTCTCCTCCTAAAGAGTAACGGAGGAGCACGAAGGTTGGCTAATCCTGGTCGGACAT‑CAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGa         <  3:1103192/150‑1 (MQ=33)  attcccGTGATCCGGGTTGCGGACAGTGTCTGGTTGGTAGTTTGAAGGGGGCGGTCTCCTCCTGTAGTGTAACGGAGGAGCACGAAGGTTGGCTCACCCTGGGC‑ATCATGCAGGAGTTTAGTGCTATTGCATAAGCCAGCTTGACTGcaa         <  4:519658/147‑3 (MQ=2)      ccGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTGGGGCGGTCTCCTCCTAAAGAGTAACGGAGGAGCACGAAGGTTGGCTAATCCTGGTCGGACAT‑CAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGctcgcgc     >  4:889642/1‑144 (MQ=17)      ccGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTGGGGCGGTCTCCTCCTAAAGAGTAACGGAGGAGCACGAAGGTTGGCTAATCCTGGTCGGACAT‑CAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGt     <  1:1836481/150‑1 (MQ=32)      ccGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTGGGGCGGTCTCCTCCTAAAGAGTAACGGAGGAGCACGAAGGTTGGCTAATCCTGGTCGGACAT‑CAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGt     <  1:265649/150‑1 (MQ=32)      ccGTGATACGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTGGGGCGGTCTCCTCCTAAAGAGTAACGGAGGAGCAAGAAGGTTGGCTAATCCTGGTCGGACAT‑CAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGt     >  1:831947/1‑150 (MQ=17)        gTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTGGGGCGGTCGCCTCCTAAAGAGTAACGGAGGAGCACGAAGGTTGGCTAATCCTGGTCGGACAT‑CAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGa   <  1:1279015/150‑1 (MQ=17)         tGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACAGGGGCGGTCTCCTCCTAAAGAGTAACGGAGGAGCACGAAGGTTGGCTAATCCTGGTCGGACAT‑CAGGAGGTTAGTTCAATGGCATAAGCCAGCTTGACTGCGAGCGTGAc  <  3:1257408/150‑1 (MQ=14)                                                                                                                                                      |          TGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTGGGGCGGTCTCCTCCTAAAGAGTAACGGAGGAGCACGAAGGTTGGCTAATCCTGGTCGGACAT‑CAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGAC  >  NC_000913_3_pae_pgi/3945905‑3946062

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