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

GCTGTTAACGGCATGGGCGGCCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACG  >  NC_000913/355435‑355548                                                    |                                                               gctgTTAACGGCATGGGCGGCCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCaa                                                >  1:541663/1‑68 (MQ=255)    gTTAACGGCATGGGCGGCCTGGACGCAGTAAAAGCGGGCGTTCCCGCACAACCGTTAGATTTCAATGt                                             <  2:157104/68‑1 (MQ=255)        aCGGCATGGGCGGCCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTcgcg                                         >  1:309948/1‑68 (MQ=255)        aCGGCATGGGCGGCCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTcgcg                                         >  1:316497/1‑68 (MQ=255)         cGGCATGGGCGGCCTGGACGCATTAAAAGCGGTCGTTCCCGCACtgataaatgaccccacaaccagcg                                        >  2:255549‑M1/1‑44 (MQ=255)                ggcggcCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGAtt                                                    <  1:610462/49‑1 (MQ=255)                ggcggcCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGAtt                                                    >  2:610458/1‑49 (MQ=255)                 gcggcCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCTATGTCgcgctggcgctg                                >  1:545650/1‑68 (MQ=255)                   ggcCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGt                              <  1:558934/68‑1 (MQ=255)                   ggcCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGt                              <  2:472437/68‑1 (MQ=255)                            gCATTAAAAGCGGTCGTTCCCGCACtgataaatgaccccacaaccagcgccagcgcgacattgaaatc                     <  1:255550‑M1/68‑44 (MQ=255)                            gCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTc                     <  1:598444/68‑1 (MQ=255)                            gCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTc                     <  1:313536/68‑1 (MQ=255)                              aTTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCAt                   <  1:219740/68‑1 (MQ=255)                              aTTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCAt                   <  2:240220/68‑1 (MQ=255)                              aTTAAAAGCGGGCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCAt                   <  2:51347/68‑1 (MQ=255)                                 aaaaGCGGTCGTTCCCGCACtgataaatgaccccacaaccagcgccagcgcgacattgaaatctaacg                <  1:277210‑M1/68‑49 (MQ=255)                                   aaGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATc              <  2:294024/68‑1 (MQ=255)                                    aGCGGTCGTTCCCGCGCAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATcg             >  2:295981/1‑67 (MQ=255)                                    aGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCa             >  1:519903/1‑68 (MQ=255)                                          cGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAgtgcgg       <  2:537638‑M2/68‑7 (MQ=255)                                               ccGCACtgataaatgaccccacaaccagcgccagcgcgacattgaaatctaacggttgtgcgggaacg  <  2:537638‑M1/68‑63 (MQ=255)                                                cGCACtgataaatgaccccacaaccagcgccagcgcgacattgaaatc                     >  2:459950‑M1/1‑5 (MQ=255)                                                cGCACtgataaatgaccccacaaccagcgccagcgcgacattgaaatc                     <  1:459952‑M1/48‑44 (MQ=255)                                                    |                                                               GCTGTTAACGGCATGGGCGGCCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACG  >  NC_000913/355435‑355548

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