New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 290192	39 (1.670)	21 (0.900)	13/166	0.5	35.0%	coding (319/855 nt)	BW25113_RS01425	hypothetical protein
	?	NZ_CP009273	= 378046	NA (NA)					pseudogene (262/264 nt)	BW25113_RS23340	acetyltransferase

GCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAA  >  NZ_CP009273/377965‑378129                                                                                  |                                                                                    gctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTa                                                                             >  1:651118/1‑90 (MQ=38)   tcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATc                                                                           >  1:322343/1‑90 (MQ=40)   tcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATc                                                                           >  1:368550/1‑90 (MQ=40)                                       gggTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATa                                       >  1:43285/1‑90 (MQ=255)                                            aaGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTa                                                                             >  1:556453/1‑47 (MQ=255)                                            aaGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTa                                                                             <  2:556453/47‑1 (MQ=255)                                            aaGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGc                                  >  2:673565/1‑90 (MQ=255)                                               tGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCAcc                               <  1:673565/90‑1 (MQ=255)                                                 aTAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCAt                             >  1:424242/1‑90 (MQ=255)                                                  tAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTgcgc                                                              >  2:76795/1‑56 (MQ=255)                                                  tAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTgcgc                                                              <  1:76795/56‑1 (MQ=255)                                                       gATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACaaa                       >  2:78325/1‑90 (MQ=255)                                                        aTGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAAt                      >  2:122397/1‑90 (MQ=255)                                                          gTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATAc                    >  2:314531/1‑90 (MQ=255)                                                            cTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACgg                  >  2:582088/1‑90 (MQ=255)                                                                     taGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAccc                                             >  1:184101/1‑54 (MQ=255)                                                                     taGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAccc                                             <  2:184101/54‑1 (MQ=255)                                                                      aGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCgcgc        <  2:620851/90‑1 (MQ=255)                                                                      aGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCgcgc        <  2:121643/90‑1 (MQ=255)                                                                         ggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAcc                                              >  2:42928/1‑49 (MQ=255)                                                                         ggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAcc                                              <  1:42928/49‑1 (MQ=255)                                                                            aaaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaa  >  1:127955/1‑90 (MQ=255)                                                                                  |                                                                                    GCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAA  >  NZ_CP009273/377965‑378129

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