New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4037450	NA (NA)	84 (0.460)	24/182	6.6	24.2%	intergenic (+115/‑69)	alaT/rrlA	tRNA‑Ala/23S ribosomal RNA of rrnA operon
	?	NC_000913	= 4168568	284 (1.450)					intergenic (+121/‑73)	gltT/rrlB	tRNA‑Glu/23S ribosomal RNA of rrnB operon
Rejected: Coverage evenness skew score above cutoff.

AATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATC  >  NC_000913/4037358‑4037494                                                                                             |                                             aagaCTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAgttgtt                                        <  2:81954/96‑1 (MQ=1)                            tCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAacac             <  1:2202274/99‑1 (MQ=32)                            tCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAacac             >  1:3350096/1‑99 (MQ=32)                              aCTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACga           >  1:1564106/1‑99 (MQ=34)                              aCTGCGAAGTTTTGCTCTTTAAAAAGCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCTTGAGTCTCTCGAATTTTCGCAACACga           <  2:3537797/99‑1 (MQ=255)                                tGCGAAGTTTTGCTCTTTGAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg         <  2:4289289/99‑1 (MQ=18)                                tGCGAAGTTTTGCTCTTTAAAAATCTGGGTCAAGCTTAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGCGTCTCTCAAATTTTCGCAACACgatg         >  2:1029445/1‑99 (MQ=255)                                tGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCACCACgatg         >  2:3409927/1‑99 (MQ=18)                                tGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg         >  2:3410905/1‑99 (MQ=34)                                tGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg         >  2:1821279/1‑99 (MQ=34)                                tGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg         >  1:1094417/1‑99 (MQ=33)                                 gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGCAACACTGAACAACGAAAGTTGTTCGTTAGTCTCTCAAATTTTCGCATCGCgatga        <  2:4261339/99‑1 (MQ=11)                                 gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatga        >  2:307176/1‑99 (MQ=34)                                  cGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACATCGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatgat       <  2:2845769/99‑1 (MQ=18)                                  cGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAACTTTTCGCAACACgatgat       >  2:4869169/1‑99 (MQ=18)                                  cGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatgat       >  2:1826554/1‑99 (MQ=34)                                  cGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGAAACACgatgat       <  1:4824888/99‑1 (MQ=21)                                   gAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGatgatg      >  1:3804557/1‑99 (MQ=34)                                   gAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGatgatg      <  2:2030121/99‑1 (MQ=34)                                   gAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGatgatg      >  1:2227061/1‑99 (MQ=33)                                    aaGTTTTGCTCTTTGAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGAtgatga     <  2:1547668/99‑1 (MQ=21)                                    aaGTTTTGCTCTTTAAAAATCTGGATCCAGCTGAAAATTGAAACAATGAACAACGAAAGTTGTGGGTGAGTCTATCAAATTTTCGCAACACGAtgatga     <  1:5018106/99‑1 (MQ=11)                                     aGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGaa    <  1:4448598/99‑1 (MQ=33)                                       ttttGCTCTTTAAAAATCTGGATCGAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATc  <  1:5161369/99‑1 (MQ=14)                                       ttttGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATc  <  1:4992943/99‑1 (MQ=31)                                       ttttGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATc  >  1:2429289/1‑99 (MQ=31)                                       ttttGCTCTTTAAAAATCTGCATCAAGCTGAAAATTGAAACACTCAACAACGAAAGTTGTTCCTGAGTCTCTCAAATTTTCGCAACACGATGATGAATc  <  2:847567/99‑1 (MQ=255)                                                                                             |                                             AATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATC  >  NC_000913/4037358‑4037494

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