New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 2035140	NA (NA)	7 (0.350)	4/134	2.2	20.0%	coding (717/768 nt)	wbbD	putative glycosyltransferase WbbD
	?	REL606	= 2035556	28 (1.420)					coding (301/768 nt)	wbbD	putative glycosyltransferase WbbD

AAGAGCTAAGAATATGAACATCACCATTCTTAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCC  >  REL606/2035484‑2035629                                                                         |                                                                          aaGAGCTAAGAATATGAACATCACCATTCTTAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTtat                                                                        >  2:357080‑M2/1‑73 (MQ=255)  agagCTAAGAATATGAACATCACCATTCTTAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTtata                                                                       >  2:642574‑M2/1‑72 (MQ=255)        aaGAATATGAACATCACCATTCTTAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTTGTAtcatca                                                                 >  2:365922/1‑76 (MQ=255)         aGAATATGAACATCACCATTCTTAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTTATATCATCAg                                                                <  1:597545/76‑1 (MQ=255)          gAATATGAACATCACCATTCTTAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGa                                                               >  1:114807/1‑76 (MQ=255)            atatGAACATCACCATTCTTAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAAt                                                             <  1:579855/76‑1 (MQ=255)              atGAACATCACCATTCTTAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTtatatcatatcaagac                                                           >  2:664065‑M2/1‑60 (MQ=255)                 aaCATCACCATTCTTAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCa                                                        <  1:683243/76‑1 (MQ=255)                   cATCACCATTCTTAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACg                                                      >  1:546159/1‑76 (MQ=255)                     tCACCATTCTTAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTa                                                    >  2:400667/1‑76 (MQ=255)                     tCACCATTCTTAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTa                                                    >  2:328602/1‑76 (MQ=255)                     tCACCATTCTTAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTa                                                    >  2:328592/1‑76 (MQ=255)                         cATTCTTAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTtatatcatatcaagactttataataaa                                                >  2:95311‑M2/1‑49 (MQ=255)                         cATTCTTAGTCATAAAACCAAATATACATTCCAATCGATTGGTTCTATTtatgtcatatcaagactttataataaa                                                >  2:244526‑M2/1‑49 (MQ=255)                             cTTAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTATTACTaaa                                            >  2:513870/1‑76 (MQ=255)                             cTTAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTATTACTaaa                                            <  1:497952/76‑1 (MQ=255)                               tAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTtatatcatatcaagactttataataaatattat                                          >  2:422554‑M2/1‑43 (MQ=255)                               tAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTtatatcatatcaagactttataataaatattat                                          >  2:422553‑M2/1‑43 (MQ=255)                               tAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTATTACTAAATc                                          <  1:122875/76‑1 (MQ=255)                                   cATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGc                                      <  1:621631/76‑1 (MQ=255)                                      aaaaTCAAGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATa                                   >  2:597572/1‑76 (MQ=255)                                      aaaaTCAAGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCAAGCATa                                   >  2:711096/1‑76 (MQ=255)                                         aTCAAGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACAt                                <  1:132302/76‑1 (MQ=255)                                           taaatATACATTCGAATCGATTTGTTCTATTtatatcatatcaagacttaataataaatattatcacaaggatatt                              >  2:672473‑M2/2‑31 (MQ=255)                                           cAAGTATACATTCGAATCGATTTGTTCTATTtatatcatatcaagactttataataaatattatcacaaggatatt                              >  2:672460‑M2/1‑31 (MQ=255)                                            aaGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGa                             >  1:349075/1‑76 (MQ=255)                                            aaGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGa                             >  1:583795/1‑76 (MQ=255)                                            aaGTATACATTCGAATCCATTTGTTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGa                             >  1:170353/1‑76 (MQ=255)                                             aGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGAt                            <  2:561467/76‑1 (MQ=255)                                              gTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGAtt                           >  2:157678/1‑76 (MQ=255)                                                       cGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCaaa                  >  1:525755/1‑76 (MQ=255)                                                       cGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCaaa                  >  1:474955/1‑76 (MQ=255)                                                            cGATTTGTTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTa             >  1:183645/1‑76 (MQ=255)                                                               cttGTTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACg          >  1:372203/2‑76 (MQ=255)                                                                ttGTTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGc         >  2:376662/1‑76 (MQ=255)                                                                     cTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGccc    >  1:12333/1‑76 (MQ=255)                                                                       aTTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGccccc  <  2:414413/76‑1 (MQ=255)                                                                         |                                                                          AAGAGCTAAGAATATGAACATCACCATTCTTAGTCATAAAATCAAGTATACATTCGAATCGATTTGTTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCC  >  REL606/2035484‑2035629

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