New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 225281	NA (NA)	10 (0.150)	8/272	7.6	NA	noncoding (1511/1542 nt)	rrsH	16S ribosomal RNA of rrnH operon
	?	NC_000913	= 225293	NA (NA)					noncoding (1523/1542 nt)	rrsH	16S ribosomal RNA of rrnH operon
Rejected: Coverage evenness skew score above cutoff.

GTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAA  >  NC_000913/225205‑225419                                                                                         |                                                                                                                               gtcggtagCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACAGa                                                                          <  1:254804/140‑1 (MQ=1)         cgtatccTTAGGGAGGTCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACAGATTGTCtga                                                                  <  2:630836/138‑1 (MQ=0)                ttggggAGGGCGCTTACCACTGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACAGATTGTCTGATGAAAAt                                                           <  2:1063947/140‑1 (MQ=1)                            ctgaCCACTTTGTGATTCATGCCTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTGCTTTGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAAc                                               <  1:1133268/140‑1 (MQ=0)                                       atgcgtcgtGACTGCGGTCAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCCGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAACCCTCTACAGGCt                                    <  2:1060593/135‑1 (MQ=1)                                                      ggggAAGGCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAggtggt                     <  1:32112/140‑1 (MQ=1)                                                                  gcagggTAACCGGCGGGGCACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACCGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCAccc         <  2:255371/139‑1 (MQ=1)                                                                         cccctatgggcagccTGCGGTTGGGTCCCCTCCTTACTGTAAAGAAGCGTACTTTGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATaa  <  2:565345/130‑1 (MQ=1)                                                                                         |                                                                                                                               GTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAA  >  NC_000913/225205‑225419

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