New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3423376 =	16 (0.170)	11 (0.120)	5/196	NT	48.2%	intergenic (+182/+47)	yhdZ/rrfF	putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
	?	NC_000913	= 4213162	8 (0.090)					intergenic (+3/‑72)	rrfE/yjaA	5S ribosomal RNA of rrnE operon/stress‑induced protein

AGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGC  >  NC_000913/4213051‑4213201                                                                                                                |                                        aGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGACATCAAataaaacaaaaggcccagtcttccgactgagccttt     >  3:47975‑M2/1‑112 (MQ=255)        cgGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaa                                    <  4:1034728‑M2/110‑6 (MQ=255)        cgGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaa                                    >  3:1034728‑M2/1‑105 (MQ=255)           tggtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                        >  1:343197‑M2/1‑102 (MQ=255)           tggtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                        <  2:343197‑M2/103‑2 (MQ=255)             gtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaa                                    <  1:265316‑M2/105‑6 (MQ=255)             gtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaa                                    >  2:265316‑M2/1‑100 (MQ=255)               cccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGACTTCAAataaaacaaaaggcccagtcttccgactgagcc        >  4:1038098‑M2/1‑98 (MQ=255)               cccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGACATCAActaaaacaaaaggcccagtcttccgactgagcc        <  3:1038098‑M2/131‑34 (MQ=255)                ccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCgg            >  1:128625/1‑126 (MQ=255)                ccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCgg            <  2:128625/126‑1 (MQ=255)                      gACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaa                                    >  1:119114‑M2/1‑91 (MQ=255)                      gACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaa                                    <  2:119114‑M2/96‑6 (MQ=255)                          ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                        <  1:70762‑M2/88‑2 (MQ=255)                          ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                        >  2:70762‑M2/1‑87 (MQ=255)                           cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggcccagtcttccgactga           <  2:491071‑M2/116‑31 (MQ=255)                           cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggcccagtcttccgactga           >  1:491071‑M2/1‑86 (MQ=255)                           cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                        <  1:1103330‑M2/87‑2 (MQ=255)                           cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                        >  2:1103330‑M2/1‑86 (MQ=255)                           cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGt                        <  2:680080/103‑1 (MQ=255)                           cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGt                        >  1:680080/1‑103 (MQ=255)                                  aaCTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATa                   >  2:918964/1‑101 (MQ=255)                                  aaCTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATa                   <  1:918964/101‑1 (MQ=255)                                        gAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGACATCAAataaaacaaaaggcccagtcttccgactg            <  4:881736‑M2/102‑30 (MQ=255)                                        gAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGACATCAAataaaacaaaaggcccagtcttccgactg            >  3:881736‑M2/1‑73 (MQ=255)                                               aaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                        <  4:965349‑M2/67‑2 (MQ=255)                                               aaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAa                                        >  3:965349‑M2/1‑66 (MQ=255)                                                 cGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggcccagtctttcgactgagcct       >  1:358880‑M2/1‑64 (MQ=255)                                                 cGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggcccagtctttcgactgagcct       <  2:358880‑M2/98‑35 (MQ=255)                                                 cGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaa                                    >  1:905020‑M2/1‑64 (MQ=255)                                                 cGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAataaa                                    <  2:905020‑M2/69‑6 (MQ=255)                                                                             ctcCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGc  >  3:80187/1‑75 (MQ=255)                                                                             ctcCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGc  <  4:80187/75‑1 (MQ=255)                                                                               cctcATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggcccagtcttccgactgagcc        >  3:926911‑M2/1‑34 (MQ=255)                                                                               cctcATGCGAGAGTAGGGAACTGCCAGGCATCAAataaaacgaaaggcccagtcttccgactgagcc        <  4:926911‑M2/67‑34 (MQ=255)                                                                                                                |                                        AGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGC  >  NC_000913/4213051‑4213201

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