New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 5696567	15 (0.210)	5 (0.070)	5/584	5.4	32.7%	intergenic (+19/‑193)	pyrC'/hslV	dihydroorotase‑like protein/peptidase component of the ATP‑dependent HslVU protease
	?	AE015451	5696629 =	6 (0.090)					intergenic (+81/‑131)	pyrC'/hslV	dihydroorotase‑like protein/peptidase component of the ATP‑dependent HslVU protease
Rejected: Coverage evenness skew score above cutoff.

CTGGTGCAGGACGGCCTGCTCGACCTGCCGACCTTGCTGGCCCGCCTGAGCAGCGGCCCGGCTGCCGCCCTGCGCCTGCCGGCCGGTGAGCTGAAAGTGGGTGGCGCGGCCGACCTGGTGCTGTTCGACCCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAA  >  AE015451/5696291‑5696692                                                                                                                                                                                                                                                                                     |                                                                                                                              cttgTGCAGGACGGCCTGCGCGACCTGCCGTGCTTGCTGGCCCGGCTGCGCAGCTGCCCGGCTGCCGCCCGGCGCCTGCCGGCCGGGGAGCTGAAAGTGGGTGGCGCGGCCGACCTGGTTCTGTTCGACCCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCATGCCCGcgaatgccttcacatagt                                                                                                       <  3:75150/298‑18 (MQ=255)       cAGGACGGCCTGCTCGACCTGCCGACCTTGCTGGCCCGCCTGAGCAGCGGCCCGGCTGCCGCCCTGCGCCTGCCGGCCGGTGAGCTGAAAGTGGGTGGCGCGGCCGACCTGGTGCTGTTCGACCCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTcc                                                                                                 <  4:106451/301‑1 (MQ=255)        aGGACGGCCTGCTCGACCTGCCGACCTTGCTGGCCCGCCTGAGCAGCGGCCCGGCTGCCGCCCTGCGCCTGCCGGCCGGTGAGCTGAAAGTGGGTGGCGCGGCCGACCTGGTGCTGTTCGACCCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTTCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTccc                                                                                                >  1:289314/1‑300 (MQ=255)        aGGACGGCCTGCTCGACCTGCCGACCTTGCTGGCCCGCCTGAGCAGCGGCCCGGCTGCCGCCCTGCGCCTGCCGGCCGGTGAGCTGAAAGTGGGTGGCGCGGCCGACCTGGTGCTGTTCGACCCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCa                                                                                                >  4:422226/1‑301 (MQ=255)              gCCTGCTCGACCTGCCGACCTTGCTGGCCCGCCTGAGCAGCGGCCCGGCTGCCGCCCTGCGCCTGCCGGCCGGTGAGCTGAAAGTGGGTGGCGCGGCCGACCTGGTGCTGTTCGACCCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCAt                                                                                          <  3:289313/301‑1 (MQ=255)              cccTGCTCGCCCTGCCGACCTGTCTGGCCCTCCTTCGCAGCGGCCCGGCTGCCGCCCTGCTCCTGCCGTCCGGTTATCTGAAAGTGGGTGTCGCGGCCGACCTGGTGCTGTTCGACCCGCAAGCGTCCACCGTTGCTGGCGAACAATTGTTCTCCCGTGGCTAGAACTGCCCTTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCTATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCAt                                                                                          <  1:96846/300‑1 (MQ=255)              cccTGCTCGACCTGCCGACCTTGCTGGCCCGCCTGAGCAGCGGCCCGGCTGCCGCCCTGCGCCTGCCGGCCGGTGAGCTGAAAGTGGGTGGCGCGGCCGACCTGGTGCTGTTCGACCCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCAt                                                                                          <  3:10393/300‑1 (MQ=255)                  ccccgcccTGCCGCCCTTGCTGGCCCGGCTGAGCAGCGGCCCGGCGGCCGCCCTGCGCCTGCCGGCCGGTGAGCTGAAAGTGGGTGGCGCGGCCGACCTGGTGCTGTTCGACCCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTGCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTgg                                                                                      <  2:159030/295‑1 (MQ=255)                   cTCGACCTGCCGACCTTGCTGGCCCGCCTGAGCAGCGGCCCGGCTGCCGCCCTGCGCCTGCCGGCCGGTGAGCTGAAAGTGGGTGGCGCGGCCGACCTGGTGCTGTTCGACCCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGt                                                                                     >  1:632855/1‑301 (MQ=255)                       aCCTGCCGACCTTGCTGGCCCGCCTGAGCAGCGGCCCGGCTGCCGCCCTGCGCCTGCCGGCCGGTGAGCTGAAAGTGGGTGGCGCGGCCGACCTGGTGCTGTTCGACCCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCCCAATTCCAAGGCATGTGGTAAcc                                                                                 >  3:368042/1‑301 (MQ=255)                             cgcgcatgctggcgctcatgcgcagcGGCCAGGCTGCCGCTCTTCGCCTGCCGGCCGGTGGGCTGAAAGTGGGTGGCGCGGCCGACCTGGTGCTGTTCGACCCGCAAGCGGCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTgg                                                                           <  2:422117/280‑1 (MQ=255)                                     cTGGCCCGCCTGAGCAGCGGCCCGGCTGCCGCCCTGCGCCTGCCGGCCGGTGAGCTGAAAGTGGGTGGCGCGGCCGACCTGGTGCTGTTCGACCCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGTGCACGCCCGCTCCCCCAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGCGAGCGGGc                                                                   >  4:554006/1‑301 (MQ=255)                                                                           ccTGCCGGCCGGTGAGCTGAAAGTGGGTGGCGCGGCCGACCTGGTGCTGTTCGACCCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTACCCCCGCCAAgc                             >  1:35656/1‑301 (MQ=255)                                                                             tGCCGGCCGGTGAGCTGAAAGTGGGTGGCGCGGCCGACCTGGTGCTGTTCGACCCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGACCGCGAATGCCTTCACATAGTCCACCCGCCCAgccc                           >  3:213855/1‑299 (MQ=255)                                                                                         aGCTGAAAGTGGGTGGCGCGGCCGACCTGGTGCTGTTCGACCCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATcc               >  1:266708/1‑301 (MQ=255)                                                                                                      tggCGCGGCCGACCTGGTGCTGTTCGACCCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCACCCCGCCAAgccccaccctccaccctctccctctacaa  >  4:341818/1‑274 (MQ=255)                                                                                                                                                                                                                                                                                     |                                                                                                                              CTGGTGCAGGACGGCCTGCTCGACCTGCCGACCTTGCTGGCCCGCCTGAGCAGCGGCCCGGCTGCCGCCCTGCGCCTGCCGGCCGGTGAGCTGAAAGTGGGTGGCGCGGCCGACCTGGTGCTGTTCGACCCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAA  >  AE015451/5696291‑5696692

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