New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4306577 =	30 (1.450)	3 (0.150)	3/250	NT	9.9%	coding (21/1986 nt)	yjcS	alkyl sulfatase
	?	NC_000913	4306596 =	26 (1.330)					coding (2/1986 nt)	yjcS	alkyl sulfatase

CGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGC  >  NC_000913/4306435‑4306712                                                                                                                                               |                                                                                                                                        aatgcctaataaaatggcgcgttccgctcacacttccacttattaaagtcgcaatttccccattatgataacgccagctttaaatgacgtattttgcatggattacatatttacgacataaatgataaatggagatttatatg                                                                                                                                         <  2:104552‑M1/1‑1 (MQ=255)     cctaataaaatggcgcgttccgctcacacttccacttattaaagtcgcaatttcccaattatgataacgccagctttaaatgccgtattttgcatggattacatatttacgacataaattataactggagatttatatGAATa                                                                                                                                     <  1:4696‑M1/5‑1 (MQ=255)     cctaataaaatggcgcgttccgctcacacttccacttattaaagtcgcaatttcccaattatgataacgccagctttaaatgacgtatcttgcatggattacatatttacgacataaatgataaatggagacttatatGAATa                                                                                                                                     <  2:14777‑M1/5‑1 (MQ=255)      aaCGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGTCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAGCGGAATaa                                                                                                                                    >  1:14777/1‑143 (MQ=255)       taataaaatggcgcgttccgctcacacttccacttattaaagtcgcaatttcccacttctgataacgccagctttaaatgccgtattttccatggcttacatatttacgacataaatgataaatggtgatttctatgaatctc                                                                                                                                   <  2:381587‑M1/7‑4 (MQ=255)            aaTGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCgaga                                                                                                                              >  1:297920/1‑143 (MQ=255)                aTTGGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGACGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGACTAACCGAGAGTTa                                                                                                                          >  2:347240/1‑143 (MQ=255)                  ttgttgTAATTAGGATTAACCTGGTGGGTATACTGAGTGGCGGCGTTCGCTTCCCCTGTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTAtt                                                                                                                        <  1:10689/138‑1 (MQ=255)                         cgctcacacttccacttattaaagtcgcaatttcccaattatgataacgccagctttaaatgacgtatcttgcatgtattacatatttacgacataaatgataaatggagatttatatGAATAACCGAGAGTTATTCATATaa                                                                                                                 >  2:215263‑M1/119‑143 (MQ=255)                          aTTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATaaa                                                                                                                >  1:96279/1‑143 (MQ=255)                             cttaTTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAAtct                                                                                                             <  2:199655/140‑1 (MQ=255)                                    aCCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATCACCGAGAGTTAATCATATAAGTCTCcaatta                                                                                                      >  2:133470/1‑143 (MQ=255)                                       cttattaaagtagcagtttcccaatgatgataacgctagcgttaaatgacgtattttgcatggattacatatttacgacataaatgataaatggagatttatatGAATAACCGAGAGTTATTCATATAAATCTCcatttatca                                                                                                   >  2:360811‑M1/105‑143 (MQ=255)                                          tGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatttatcattt                                                                                                <  1:351079/143‑1 (MQ=255)                                               atatTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGATAAAACGGAATAACGGAGAGTTATTCATATAAATCTCTATTTATAATTTATGTc                                                                                           >  1:65537/1‑143 (MQ=255)                                                 atTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATCACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATCTAAATCTCCATTTATCATTTATGTCgg                                                                                         >  1:104552/1‑142 (MQ=255)                                                      gTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGACCGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTAACATTTATGTGATAAata                                                                                    >  2:352287/1‑143 (MQ=255)                                                                tatgataacgccagctttaaatgacgtattttgcatggattacatatttacgacataaatgataaatggagatttatatGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCAt                                                                          >  1:359158‑M1/80‑143 (MQ=255)                                                                tCGCTGCCTCGCTCGCGTTAGCGGTATCTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCAt                                                                          <  2:21728/143‑1 (MQ=255)                                                                 cGCTTCCTCTTGCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATg                                                                         <  2:20714/143‑1 (MQ=255)                                                                   cTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCa                                                                       >  1:300245/1‑143 (MQ=255)                                                                    ttCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCaa                                                                      >  1:214726/1‑143 (MQ=255)                                                                    ttCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAAAAATCCGGCTCGAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCAGTGATGTCGTAAATATGTAATCCAGGCaa                                                                      >  2:154451/1‑143 (MQ=255)                                                                               cggtaccGGTATTTACCATCATGCTGCTGGCGGCAGTTAACTCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCAGTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCAtt                                                           <  2:288580/137‑1 (MQ=255)                                                                               cgTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCAtt                                                           >  1:353322/1‑143 (MQ=255)                                                                                gTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCTCTTTATCGTTTATGTCGTAAATATGTAATCCATGCAAAATACGTCAttt                                                          >  1:376256/1‑143 (MQ=255)                                                                                    gCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGGATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAg                                                      >  2:66366/1‑143 (MQ=255)                                                                                    gCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTGTTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAg                                                      >  1:208580/1‑143 (MQ=255)                                                                                         aTTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGc                                                 <  1:297068/143‑1 (MQ=255)                                                                                         aTTTACCAGCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGc                                                 <  1:209414/143‑1 (MQ=255)                                                                                                    tGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATaac                                      >  2:310660/1‑142 (MQ=255)                                                                                                    tGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAAt                                      >  1:58260/1‑143 (MQ=255)                                                                                                               gcAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAACTCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAATTGc                           >  1:322018/1‑143 (MQ=255)                                                                                                                cAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATGTATCATTTATGTCGTAAATATGTAATCCATCCAAAATACGTCATTTAAACCTGGCGTGATCATCATTTGGAAATTGCg                          >  2:87975/1‑143 (MQ=255)                                                                                                                                aaTCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAATTGCGACTTGAAAAAGTGGaa          >  2:128130/1‑143 (MQ=255)                                                                                                                                aaTCCTGCGCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGaa          >  2:179375/1‑143 (MQ=255)                                                                                                                                        tCAAACGGAAGAACCGAGAGTTATTCATATAAATCTCCATGTATCATTTATGTCGTAAATATGTAAGCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGc  <  1:133470/143‑1 (MQ=255)                                                                                                                                               |                                                                                                                                        CGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGC  >  NC_000913/4306435‑4306712

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