New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 355486	17 (0.920)	3 (0.180)	3/118	1.6	21.3%	coding (565/1260 nt)	codB	cytosine transporter
	?	NC_000913	= 355537	7 (0.420)					coding (616/1260 nt)	codB	cytosine transporter

AGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGACTTTGTCCGGTTTGGTCGCAATGCCAAACTGGCGGTGCTGGTGGC  >  NC_000913/355470‑355604                                                                    |                                                                    aGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCa                                                                     >  1:519903/1‑68 (MQ=255)       cGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAgtgcgg                                                               <  2:537638‑M2/68‑7 (MQ=255)                     cGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAgtgcgggaacgaccgctttt                                                 <  1:277210‑M2/68‑21 (MQ=255)                        tAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGAc                                              >  1:64240/1‑68 (MQ=255)                          gATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAgtgcgggaacgaccgcttttaatgc                                            <  1:255550‑M2/68‑26 (MQ=255)                          gATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAgtgcg                                                                <  1:459952‑M2/48‑6 (MQ=255)                          gATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAgtgcg                                                                >  2:459950‑M2/1‑43 (MQ=255)                          gATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGACtt                                            >  1:364725/1‑68 (MQ=255)                                   gTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGACTTTGTCCGGtt                                   >  1:311319/1‑68 (MQ=255)                                             cgctggTTGTGGGGTCATTTATCAgtgcgggaacgaccgcttttaatgcgtccaggccgcccatgccg                         >  2:255549‑M2/1‑24 (MQ=255)                                             cgctggGTGTGTGGTCATTTATCAGTGCGGGTACGCTCACCGCTGACTTTGTCCGGCTTGGTCGCAAt                         <  2:659973/68‑1 (MQ=255)                                                                 aTCAGTGCGGGTACGCTCACCGCTGACTTTGTCCGGTTTGGTCGCAATGCCAAACTGGCGGTGCtggt     >  2:243595/1‑68 (MQ=255)                                                                 aTCAGTGCGGGTACGCTCACAGCTGACTTTGTCCGGTTTGGTCGCAATGCCAAACTGGCGGTGCtggt     >  1:17363/1‑68 (MQ=255)                                                                    aGTGCGGGTACGCTCACCGCTGACTTTGTCCGGTTTGGTCGCAATGCCAAACTGGCGGTGCTGGTGGc  >  2:686289/1‑68 (MQ=255)                                                                    |                                                                    AGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGACTTTGTCCGGTTTGGTCGCAATGCCAAACTGGCGGTGCTGGTGGC  >  NC_000913/355470‑355604

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