New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 416438	11 (0.820)	4 (0.310)	4/250	NT	31.5%	coding (515/1203 nt)	sbcD	ATP dependent, structure specific DNA nuclease ‑ SbcD subunit
	?	NC_000913	= 416786	7 (0.550)					coding (167/1203 nt)	sbcD	ATP dependent, structure specific DNA nuclease ‑ SbcD subunit

TTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAAT  >  NC_000913/416254‑416568                                                                                                                                                                                         |                                                                                                                                   ttCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGtt                                                                                                                              >  1:145169/1‑191 (MQ=255)    cATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGc                                                                                                                           <  2:145168/191‑1 (MQ=255)                                                                                                cAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATc                                                                                                                    >  1:244714/1‑106 (MQ=255)                                                                                                   tataAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGc                                                                                                                 <  2:244713/106‑1 (MQ=255)                                                                                                          gTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATc                                                                                                                    >  2:266467/1‑96 (MQ=255)                                                                                                            cacgcacgGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCccagttacgcccgcacgttatacaaccgttttgttgtcaatttacagcaaactggct                                                                           >  1:611‑M1/1‑78 (MQ=255)                                                                                                             acgcacgGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGc                                                                                                                 <  1:266468/96‑1 (MQ=255)                                                                                                               gcacgGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCccagttacgcccgcacgttatacaaccgttttgttgtcaatttacagcaaactggctgtc                                                                        <  2:611‑M1/135‑61 (MQ=255)                                                                                                                                    cccccACGGTCGTTAAATGTCCCGTGGCGAGGATGGGGAGGGGCTGTTCGCCGCGCAGTTTGCAGGCCTCGGCATAGTGTTGTTGGTAATAATCGGTa                                                                                        >  1:259055/1‑98 (MQ=255)                                                                                                                                       ccACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATc                                                                                     <  2:259054/98‑1 (MQ=255)                                                                                                                                                                          agagGCTGATCGCCGCgcagttacgcccgcacgttatacaaccgttttgttgtcaatttacagcaaactggctgtcatctggtggtactggcaggaaaccatgactcggtcgcc                                  >  2:200252‑M1/1‑16 (MQ=255)                                                                                                                                                                             ggCTGATCGCCGCgcagttacgcccgcacgttatacaaccgttttgttgtcaatttacagcaaactggctgtcatctggtggtactggcaggaaaccatgactcggtcgccacg                               <  1:200253‑M1/114‑102 (MQ=255)                                                                                                                                                                                  aTCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGtttt                                                             >  1:258816/1‑79 (MQ=255)                                                                                                                                                                                  aTCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCct                                       >  2:13208/1‑101 (MQ=255)                                                                                                                                                                                  aTCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACgg     >  2:122572/1‑135 (MQ=255)                                                                                                                                                                                     gccgcGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCa                                                          <  2:258815/79‑1 (MQ=255)                                                                                                                                                                                     gccgcGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCctggc                                    <  1:13208/101‑1 (MQ=255)                                                                                                                                                                                     gccgcGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAAt  <  1:122573/135‑1 (MQ=255)                                                                                                                                                                                         |                                                                                                                                   TTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAAT  >  NC_000913/416254‑416568

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