New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1299498	0 (0.000)	20 (1.100)	17/124	0.1	100%	intergenic (+253/‑1684)	ychE/oppA	UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein
	?	NC_000913	1300698 =	0 (0.000)					intergenic (+1453/‑484)	ychE/oppA	UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein

TGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCAT  >  NC_000913/1299431‑1299565                                                                    |                                                                    tGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAA                                                                     <  2:75218‑M1/68‑1 (MQ=255)    gCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAata                                                                  >  2:365768‑M1/1‑65 (MQ=255)                     tGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaat                                                 >  2:475714‑M1/1‑48 (MQ=255)                                tattaATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatt                                      <  1:197671‑M1/68‑32 (MQ=255)                                   taATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttca                                   >  2:648417‑M1/1‑34 (MQ=255)                                   taATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttca                                   >  1:294492‑M1/1‑34 (MQ=255)                                     aTTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcact                                 >  2:552315‑M1/1‑32 (MQ=255)                                       tCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactac                               >  2:447587‑M1/1‑30 (MQ=255)                                       tCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaatagttaaattattgtgcatttcactac                               >  2:446488‑M1/1‑30 (MQ=255)                                         cACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactg                             >  1:126592‑M1/1‑28 (MQ=255)                                          aCTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactgg                            >  1:345868‑M1/1‑27 (MQ=255)                                           cTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactgga                           <  2:182596‑M1/68‑43 (MQ=255)                                             gCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaac                         <  1:365768‑M1/68‑45 (MQ=255)                                                tttGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgt                      <  1:102996‑M1/68‑48 (MQ=255)                                                    gTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatc                  >  2:494909‑M1/1‑17 (MQ=255)                                                    gTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatc                  >  2:57012‑M1/1‑17 (MQ=255)                                                     tAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatca                 >  1:172998‑M1/1‑16 (MQ=255)                                                        agGAAAGTGCTAAataataatcaattgttaaattaatgtgcatttcactactggatctgtaatcagaa              <  2:700148‑M1/68‑56 (MQ=255)                                                         ggAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaa             >  1:334128‑M1/1‑12 (MQ=255)                                                           aaaGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaag           >  1:194293‑M1/1‑10 (MQ=255)                                                                gCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagataga      >  1:710916‑M1/1‑5 (MQ=255)                                                                    aataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatg  <  2:294492‑M1/68‑68 (MQ=255)                                                                    |                                                                    TGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCAT  >  NC_000913/1299431‑1299565

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