New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 3493871	23 (0.880)	3 (0.130)	3/138	1.9	13.5%	coding (603/1338 nt)	frlA	fructoselysine/psicoselysine transporter
	?	NZ_CP009273	= 3494083	18 (0.780)					coding (815/1338 nt)	frlA	fructoselysine/psicoselysine transporter
Rejected: Frequency below/above cutoff threshold.

TTCTGGTTTAAAGCGGAGAATTTTGCCGCCCCTACCACCACTGCGATTGGCGCAACGGGCAGCTTTATGGCGCTGCTGGCGGGGATCTCTGCCACCAGTTGGTCGTATACCGGCATGGCCTCTATCTGTTATATGACCGGCGAAATTAAAAACCCCGGAAAAACCATGCCACG  >  NZ_CP009273/3493785‑3493957                                                                                       |                                                                                       ttctGGTTTAAAGCGGAGAATTTTGCCGCCCCTACCACCACTGCGATTGGCGCAACGGGCAGCTTTATGGCGCTGCTGGCGGGGATCcct                                                                                     >  1:455484‑M1/1‑87 (MQ=255)    tGGTTTAAAGCGGAGAATTTTGCCGCCCCTACCACCACTGCGATTGGCGCAACGGGCAGCTTTATGGCGCTGCTGGCGGGGATCTCTGcc                                                                                  <  2:455484/90‑1 (MQ=255)           aaGCGGAGAATTTTGCCGCCCCTACCACCACTGCGATTGGCGCAACGGGCAGCTTTATGGCGCTGCTGGCGGGGATCTCTGCCACCAGtt                                                                           <  1:47595/90‑1 (MQ=255)           aaGCGGAGAATTTTGCCGCCCCTACCACCACTGCGATTGGCGCAACGGGCAGCTTTATGGCGCTGCTGGCGGGGATCTCTGCCACCAGtt                                                                           <  1:543939/90‑1 (MQ=255)                     ttttGCCGCCCCTACCACCACTGCGATTGGCGCAACGGGCAGCTTTATGGCGCTGCTGGCGGGGATCTCTGCCACCAGTTGGTCGTATAc                                                                 <  2:607719/90‑1 (MQ=255)                            gccCCTACCACCACTGCGATTGGCGCAACGGGCAGCTTTATGGCGCTGCTGGCGGGGATCTCTGCCACCAGTTGGTCGTATACCGGCATg                                                          <  1:408125/90‑1 (MQ=255)                                      ccacTGCGATTGGCGCAACGGGCAGCTTTATGGCGCTGCTGGCGGGGATCTCTGCCACCAGTTGGTCGTATACCGGCATGGCCTCTATCt                                                <  2:562001/90‑1 (MQ=255)                                         cTGCGATTGGCGCAACGGGCAGCTTTATGGCGCTGCTGGCGGGGGTCTCTGCCACCAGTTGGGCGTATACCGGCATGGCCTCTATCTGtt                                             >  1:477384/1‑90 (MQ=255)                                         cTGCGATTGGCGCAACGGGCAGCTTTATGGCGCTGCTGGCGGGGATCcaggtcagggcgtcg                                                                         >  1:303136‑M1/1‑47 (MQ=255)                                         cTGCGATTGGCGCAACGGGCAGCTTTATGGCGCTGCTGGCGGGGATCcaggtcagggcgtcg                                                                         <  2:303136‑M1/62‑16 (MQ=255)                                         cTGCGATTGGCGCAACGGGCAGCTTTATGGCGCTGCTGGCGGGGATCTCTGCCACCAGTTGGTCGTATACCGGCATGGCCTCTa                                                   <  1:415749/84‑1 (MQ=255)                                         cTGCGATTGGCGCAACGGGCAGCTTTATGGCGCTGCTGGCGGGGATCTCTGCCACCAGTTGGTCGTATACCGGCATGGCCTCTa                                                   >  2:415749/1‑84 (MQ=255)                                           gCGATTGGCGCAACGGGCAGCTTTATGGCGCTGCTGGCGGGGGTCTCTGCCCCCCGTTGGTCGTATACCGGCATGGGCTCTATCTGTtat                                           >  2:802725/1‑90 (MQ=255)                                                     cAACGGGCAGCTTTATGGCGCTGCTGGCGGGGATCTCTGCCACCAGtt                                                                           <  2:596541/48‑1 (MQ=255)                                                     cAACGGGCAGCTTTATGGCGCTGCTGGCGGGGATCTCTGCCACCAGtt                                                                           >  1:596541/1‑48 (MQ=255)                                                      aaCGGGCAGCTTTATGGCGCTGCTGGCGGGGATCTCTGCCACCAGTTGGTCGTATACCGGCATGGCCTCTATCTGTTATATGACCGGCGa                                <  1:870214/90‑1 (MQ=255)                                                      aaCGGGCAGCTTTATGGCGCTGCTGGCGGGGATCTCTGCCACCAGTTGGTCGTATACCGGCATGGCCTCTATCTGTTATATGACCGGCGa                                <  1:802725/90‑1 (MQ=255)                                                                 ttATGGCGCTGCTGGCGGGGATCTCTGCCACCAGTTGGTCGTATACCGGCATGGCCTCTATCTGTTATATGACCGGCGAAATTAAAAAcc                     >  1:115860/1‑90 (MQ=255)                                                                  tATGGCGCTGCTGGCGGGGATCTCTGCCACCAGTTGGTCGTATACCGGCATGGCCTCTATCTGTTATATGACCGGCGAAATTAAAAAccc                    <  1:419999/90‑1 (MQ=255)                                                                        gctgctGGCGGGGATCTCTGCCACCAGTTGGTCGTATACCGGCATgg                                                         >  1:65498/1‑47 (MQ=255)                                                                        gctgctGGCGGGGATCTCTGCCACCAGTTGGTCGTATACCGGCATgg                                                         <  2:65498/47‑1 (MQ=255)                                                                        gctgctGGCGGGGATCTCTGCCACCAGTTGGTCGTATACCGGCATGGCctc                                                     <  2:147007/51‑1 (MQ=255)                                                                        gctgctGGCGGGGATCTCTGCCACCAGTTGGTCGTATACCGGCATGGCctc                                                     >  1:147007/1‑51 (MQ=255)                                                                        gctgctGGCGGGGATCTCTGCCACCAGTTGGTCGTATACCGGCATGGCCTCTATCTGTTATATGACCGGCGAAATTAAAAACCCCGGaaa              <  2:443782/90‑1 (MQ=255)                                                                        gctgctGGCGGGGATCTCTGCCACCAGTTGGTCGTATACCGGCATGGCCTCTATCTGTTATATGACCGGCGAAATTAAAAACCCCGGaaa              >  1:443782/1‑90 (MQ=255)                                                                                 ggggATCcaggtcagggcgtcggaaatcggcgtttcagaattggcgagtttgtcgaagggcatcaggccggaaatcaccagcgccagcag     <  2:346415‑M1/90‑84 (MQ=255)                                                                                 ggggATCcaggtcagggcgtcggaaatcggcgtttcagaattggcgagtttgtcgaagggcatcaggccggaaatcaccagcgccagcag     <  2:516208‑M1/90‑84 (MQ=255)                                                                                 ggggATCcaggtcagggcgtcggaaatcggcgtttcagaattggcgagtttgtcgaagggcatcaggccggaaatcaccagcgccagcag     <  2:57550‑M1/90‑84 (MQ=255)                                                                                 ggggATCcaggtcagggcgtcggaaatcggcgtttcagaattggcgagtttgtcgaagggcatcaggccggaaatcaccagcgccagcag     <  1:282619‑M1/90‑84 (MQ=255)                                                                                 ggggATCcaggtcagggcgtcggaaatcggcgtttcagaattggcgagtttgtcgaagggcatcaggccggaaatcaccagcgccagcag     <  2:716206‑M1/90‑84 (MQ=255)                                                                                    gATCTCTGCCACCAGTTGGTCGTATACCGGCATGGCCTCTATCTGTTATATGACCGGCGAAATTAAAAACCCCGGAAAAACCATGCCACg  <  2:35179/90‑1 (MQ=255)                                                                                       |                                                                                       TTCTGGTTTAAAGCGGAGAATTTTGCCGCCCCTACCACCACTGCGATTGGCGCAACGGGCAGCTTTATGGCGCTGCTGGCGGGGATCTCTGCCACCAGTTGGTCGTATACCGGCATGGCCTCTATCTGTTATATGACCGGCGAAATTAAAAACCCCGGAAAAACCATGCCACG  >  NZ_CP009273/3493785‑3493957

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