New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 2040272	54 (0.810)	5 (0.080)	4/138	8.4	10.6%	coding (1112/1113 nt)	vioA	VioA, involved in dTDP‑N‑acetylviosamine synthesis
	?	REL606	= 2040459	34 (0.550)					coding (925/1113 nt)	vioA	VioA, involved in dTDP‑N‑acetylviosamine synthesis
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GCTAGCTGAAATCAATCCTCGATACATAGGCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACGACTCATATGAAAACCATCATCTATGAGTATAG  >  REL606/2040391‑2040530                                                                     |                                                                        gCTAGCTGAAATCAATCCTCGATACATAGGCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGAtatt                                                                   <  1:343490/75‑1 (MQ=255)  cTAGCTGAAATCAATCCTCGATACATAGGCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGAtatta                                                                  <  2:543045/75‑1 (MQ=255)  cTAGCTGAAATCAATCCTCGATACATAGGCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGAtatta                                                                  <  2:1307730/75‑1 (MQ=255)    aGCTGAAATCAATCCTCGATACATAGGCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGATattatt                                                                >  2:435725/1‑75 (MQ=255)     gCTGAAATCAATCCTCGATACATAGGCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGATATTAttt                                                               >  2:1949372/1‑75 (MQ=255)      cTGAAATCAATCCTCGATACATAGGCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGATATTAtttt                                                              <  1:1137423/75‑1 (MQ=255)           atcaatcCTCGATACATAGGCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGATATTATTTTTCTTa                                                         <  1:1304339/75‑1 (MQ=255)           atcaatcCTCGATACATAGGCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGATATTATTTTTCTTa                                                         <  1:947735/75‑1 (MQ=255)            tcaatcCTCGATACATAGGCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGATATTATTTTTCTTaa                                                        >  1:52395/1‑75 (MQ=255)              aatcCTCGATACATAGGCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGATATTATTTTTCTTAAGt                                                      <  2:482968/75‑1 (MQ=255)                     gATACATAGGCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCat                                               <  1:448673/75‑1 (MQ=255)                        aCATAGGCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATg                                            >  2:1646273/1‑75 (MQ=255)                          aTAGGCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGcc                                          >  2:767658/1‑75 (MQ=255)                          aTAGGCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGcc                                          >  1:621890/1‑75 (MQ=255)                            aGGCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGCCTg                                        <  1:881008/75‑1 (MQ=255)                              gCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGCCTGAt                                      <  2:90638/75‑1 (MQ=255)                                  atTACTGATTAATGGATAAAAATACTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACg                                  >  2:397145/1‑75 (MQ=255)                                   ttACTGATTAATGGATAAAAATACTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACGa                                 >  2:145124/1‑75 (MQ=255)                                        gATTAATGGATAAAAATACTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACGACTCat                            >  2:1827338/1‑75 (MQ=255)                                          ttAATGGATAAAAATACTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACGACTCatat                          >  1:1882332/1‑75 (MQ=255)                                             aTGGATAAAAATACTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACGACTCATATGaa                       >  1:292412/1‑75 (MQ=255)                                                        tACTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCATGACTCATATGAAAACCATCATCt            <  2:1232902/75‑1 (MQ=255)                                                        tACTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACGACTCATATGAAAACCATCATCt            >  1:1553347/1‑75 (MQ=255)                                                        tACTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACGACTCATATGAAAACCATCATCt            <  2:1232889/75‑1 (MQ=255)                                                        tACTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACGACTCATATGAAAACCATCATCt            >  1:387190/1‑75 (MQ=255)                                                         tcTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACGACTCATATGAAAACCATCATCTa           <  1:2151005/74‑1 (MQ=255)                                                          cTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACGACTCATATGAAAACCATCATCTAt          >  1:859742/1‑75 (MQ=255)                                                           tttCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACGACTCATATGAAAACCATCATCTATg         <  1:1681392/75‑1 (MQ=255)                                                              cGAGATAGGATATTATTTTTCTTAAGTAATTCATATGTCTGATCACGACTCATATGAAAACCATCATCTATGAGt      <  2:594960/75‑1 (MQ=255)                                                               gagaTAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACGACTCATATGAAAACCATCATCTATGAGta     >  2:312770/1‑75 (MQ=255)                                                                agaTAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACGACTCATATGAAAACCATCATCTATGAGtat    <  2:529317/75‑1 (MQ=255)                                                                 gaTAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACGACTCATATGAAAACCATCATCTATGAGtata   <  2:688015/75‑1 (MQ=255)                                                                 gaTAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACGACTCATATGAAAACCATCATCTATGAGtata   >  2:891377/1‑75 (MQ=255)                                                                  aTAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACGACTCATATGAAAACCATCATCTATGAGTATAg  <  1:155595/75‑1 (MQ=255)                                                                     |                                                                        GCTAGCTGAAATCAATCCTCGATACATAGGCATATTACTGATTAATGGATAAAAATACTTTCGAGATAGGATATTATTTTTCTTAAGTAATTCATATGCCTGATCACGACTCATATGAAAACCATCATCTATGAGTATAG  >  REL606/2040391‑2040530

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