New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	2282400 =	NA (NA)	21 (1.330)	14/166	0.2	100%	pseudogene (1/2532 nt)	yejO	autotransporter YejO
	?	NZ_CP009273	= 3626937	0 (0.000)					pseudogene (882/1607 nt)	BW25113_RS18065	YhiJ/YhiL family protein

CATAGAATTTAGCAAAGTAAGGATATATTTATCACCAAAAATAGTACTAAATCTATCGCCTGTATGATCCCCAAGGAAGATACAAGGGGTTATCGCATCGTCCTTTAAAAAAACATGTTGATGCAACTCATTCAAGATCATGTCTCTCGTATTTAATAAATCGTAGAGCACTTTA  >  NZ_CP009273/3626850‑3627024                                                                                        |                                                                                        cATAGAATTTAGCAAAGTAAGGATATATTTATCACCAAAAATAGTACTAAATCTATCGCCTGTATGATCCCCAAGGAAGATACAAGGGaa                                                                                       <  1:240471‑M2/90‑3 (MQ=255)     gAATTTAGCAAAGTAAGGATATATTTATCACCAAAAATAGTACTAAATCTATCGCCTGTATGATCCCCAAGGAAGATACAAGGGaaggtg                                                                                   <  1:286628‑M2/90‑7 (MQ=255)              aaaGTAAGGATATATTTATCACCAAAAATAGTACTAAATCTATCGCCTGTATGATCCCCAAGGAAGATACAAGGGaaggtgcgaataagc                                                                          >  1:115992‑M2/1‑75 (MQ=255)               aaGTAAGGATATATTTATCACCAAAAATAGTACTAAATCTATCGCCTGTATGATCCCCAAGGAAGATACAAGGGaaggtgcgaataagcg                                                                         <  1:296812‑M2/90‑17 (MQ=255)                           atTTATCACCAAAAATAGTACTAAATCTATCGCCTGTATGATCCCCAAGGAAGATACAAGGGaaggtgcgaataagcggggaaattcttc                                                             <  2:352120‑M2/90‑29 (MQ=255)                              tATCACCAAAAATAGTACTAAATCTATCGCCTGTATGATCCCCAAGGAAGATACAAGGGaaggtgcgaataagcggggaaattcttctcg                                                          <  2:376112‑M2/90‑32 (MQ=255)                               aTCACCAAAAATAGTACTAAATCTATCGCCTGTATGATCCCCAAGGAAGATACAAGGGaaggtgcgaataagcggggaaattcttctcgg                                                         <  1:137145‑M2/90‑33 (MQ=255)                               aTCACCAAAAATAGTACTAAATCTATCGCCTGTATGATCCCCAAGGAAGATACAAGGGaaggtgcgaataagcggggaaattcttctcgg                                                         <  2:159280‑M2/90‑33 (MQ=255)                                tCACCAAAAATAGTACTAAATCTATCGCCTGTATGATCCCCAAGGAAGATACAAGGGaaggtgcgaataagcggggaaattcttctcggc                                                        <  1:211451‑M2/90‑34 (MQ=255)                                     aaaaaTAGTACTAAATCTATCGCCTGTATGATCCCCAAGGAAGATACAAGGGaaggtgcgaataagcggggaaattcttctcggctgact                                                   <  1:369147‑M2/90‑39 (MQ=255)                                                 aaaTCTATCGCCTGTATGATCCCCAAGGAAGATACAAGGGaaggtgcgaataagcggggaaattcttctcggctgactcagtcatt                                           >  2:455996‑M2/1‑40 (MQ=255)                                                 aaaTCTATCGCCTGTATGATCCCCAAGGAAGATACAAGGGaaggtgcgaataagcggggaaattcttctcggctgactcagtcatt                                           <  1:455996‑M2/86‑47 (MQ=255)                                                              gTATGATCCCCAAGGAAGATACAAGGGaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttg                          <  1:436303‑M2/90‑64 (MQ=255)                                                              gTATGATCCCCAAGGAAGATACAAGGGaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttg                          <  1:362498‑M2/90‑64 (MQ=255)                                                              gTATGATCCCCAAGGAAGATACAAGGGaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttg                          <  1:197666‑M2/90‑64 (MQ=255)                                                                     ccccAAGGAAGATACAAGGGaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgat                   <  2:127726‑M2/90‑71 (MQ=255)                                                                     ccccAAGGAAGATACAAGGGaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgat                   <  2:271331‑M2/90‑71 (MQ=255)                                                                     ccccAAGGAAGATACAAGGGaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgat                   <  1:226150‑M2/90‑71 (MQ=255)                                                                                    aaGGGaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaa    >  1:159280‑M2/1‑5 (MQ=255)                                                                                      gggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatg  <  1:280063‑M2/90‑88 (MQ=255)                                                                                      gggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatg  <  2:168448‑M2/90‑88 (MQ=255)                                                                                        |                                                                                        CATAGAATTTAGCAAAGTAAGGATATATTTATCACCAAAAATAGTACTAAATCTATCGCCTGTATGATCCCCAAGGAAGATACAAGGGGTTATCGCATCGTCCTTTAAAAAAACATGTTGATGCAACTCATTCAAGATCATGTCTCTCGTATTTAATAAATCGTAGAGCACTTTA  >  NZ_CP009273/3626850‑3627024

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