New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3798325	21 (1.050)	8 (0.420)	4/126	1.6	26.7%	coding (988/1074 nt)	waaU	lipopolysaccharide core biosynthesis
	?	NC_000913	= 3798516	24 (1.260)					coding (797/1074 nt)	waaU	lipopolysaccharide core biosynthesis

TGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAAT  >  NC_000913/3798463‑3798575                                                      |                                                            tGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGAtt                                               >  2:421485/1‑68 (MQ=255) tGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGAtt                                               >  2:104241/1‑68 (MQ=255)           aCTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAgtattg                                     >  2:528057/1‑68 (MQ=255)               taCTTGTTCATCTCAACTGTATAAATAAACTCATCAAATGgtatccccaacatatacggtataagata                                 >  2:14336‑M2/1‑40 (MQ=255)               taCTTGGTCAACGCAACGTTATAATTAATCTCATCAAATTgtatccccaacatatacggtaaaagata                                 >  2:195890‑M2/1‑40 (MQ=255)               taCTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagata                                 >  2:358834‑M2/1‑40 (MQ=255)                aCTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGa                                                 <  2:487023/51‑1 (MQ=255)                aCTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGa                                                 >  1:487022/1‑51 (MQ=255)                aCTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAgtattggtatt                                >  2:253813/1‑68 (MQ=255)                 cTTGTTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagatatc                               <  1:195890‑M2/68‑31 (MQ=255)                 cTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAgtattggtattg                               >  2:441002/1‑68 (MQ=255)                  ttGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagatatcg                              <  1:14336‑M2/68‑32 (MQ=255)                  ttGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagata                                 >  1:183526‑M2/1‑37 (MQ=255)                  ttGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagata                                 <  2:183526‑M2/65‑29 (MQ=255)                  ttGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTTGTGTTTCGATTTCAAGTATTGGTATTGt                              <  1:216814/68‑1 (MQ=255)                   tGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagatatcga                             <  2:222668‑M2/68‑33 (MQ=255)                    ggTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCa                            <  1:331858/68‑1 (MQ=255)                     gTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCaa                           >  1:227455/1‑68 (MQ=255)                           gcaacgGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATc                     >  2:736183/1‑68 (MQ=255)                              acgGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTAAATAAATCTTg                  >  2:341445/1‑68 (MQ=255)                                        tataCTCATCAATTGTCAGTGTTTCGATTCCTTGTATTTGTTTTGTCTATAAATCTAGCTGTAACCCg        <  1:389385/65‑1 (MQ=255)                                        taagCTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCg        >  1:34638/1‑68 (MQ=255)                                        taaaCTCATCAAATGGTAGTGTTTCGATTTCAAgt                                         <  2:90917/35‑1 (MQ=255)                                        taaaCTCATCAAATGGTAGTGTTTCGATTTCAAgt                                         >  1:90917/1‑35 (MQ=255)                                        taaaCTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCg        >  2:116599/1‑68 (MQ=255)                                           acacatcaAAAGGAAGAGATTCGATATCAAGAAAAGGAATTGTCAATAAAACATGCGCAACCCCGGTa     >  2:256554/4‑68 (MQ=255)                                           aCTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTa     >  2:255527/1‑68 (MQ=255)                                            cTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTaa    >  1:633976/1‑68 (MQ=255)                                            cTCATCAAATGCTAGTGAATCGATTTCAAGTATTGGTATTGTCAAAAAATCTTGCGGTAACCCGGTaa    >  1:633979/1‑68 (MQ=255)                                             tcatcaAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTaaa   <  2:537492/68‑1 (MQ=255)                                                     tGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAat  >  1:466798/1‑61 (MQ=255)                                                     tGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAat  <  2:466799/61‑1 (MQ=255)                                                      |                                                            TGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAAT  >  NC_000913/3798463‑3798575

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