New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	2282398 =	NA (NA)	17 (1.050)	8/168	0.6	100%	pseudogene (3/2532 nt)	yejO	autotransporter YejO
	?	NZ_CP009273	3512437 =	0 (0.000)					intergenic (‑14/+387)	aroK/hofQ	shikimate kinase AroK/DNA uptake porin HofQ

AGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCC  >  NZ_CP009273/2282315‑2282479                                                                                    |                                                                                  aGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGt                                                                             >  2:22293/1‑90 (MQ=255) aGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGt                                                                             >  1:395239/1‑90 (MQ=255)                   aGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAt                                                           >  2:312047/1‑90 (MQ=255)                         acacCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATtcttct                                                     >  2:144314/1‑90 (MQ=255)                         acacCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATtcttct                                                     >  1:440147/1‑90 (MQ=255)                                    tGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCa                                          >  1:161447/1‑90 (MQ=255)                                          aCCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatt                                    <  1:144314/90‑1 (MQ=255)                                                       gTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATg                       >  2:148888/1‑90 (MQ=255)                                                                 aTGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGa             >  2:18448/1‑90 (MQ=255)                                                                 aTGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGa             >  1:206166/1‑90 (MQ=255)                                                                       tGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttcttc                          <  2:105271/71‑1 (MQ=255)                                                                       tGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttcttc                          >  1:105271/1‑71 (MQ=255)                                                                       tGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAtttttt       <  2:93142/90‑1 (MQ=255)                                                                            gAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTccc  >  2:260285/1‑90 (MQ=37)                                                                            gAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTccc  >  2:455756/1‑90 (MQ=37)                                                                            gAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTcc   <  1:425290/89‑1 (MQ=37)                                                                            gAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTcc   >  2:425290/1‑89 (MQ=37)                                                                                    |                                                                                  AGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCC  >  NZ_CP009273/2282315‑2282479

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