New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	2103478 =	18 (1.080)	7 (0.430) +TT	4/282	1.6	31.4%	intergenic (+49/+27)	msrC/PP_1878	free methionine‑(R)‑sulfoxide reductase/hypothetical protein
	?	NC_002947	4898402 =	13 (0.780)					intergenic (+68/‑140)	PP_4304/sbp‑I	VIC family cation transporter/sulfate ABC transporter substrate‑binding protein

GACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACG  >  NC_002947/4898256‑4898542                                                                                                                                                   |                                                                                                                                             gACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCGAAAGGCGCGGCCCGGGCGCGGTCCCCCGTAg                                                                                                                                           >  2:326274/1‑150 (MQ=255)         acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTAGCGGGTAAATCCCCACCAACCGGCGCGGACCCGGGGCGGGCCCCCGGAGGAGCGGGt                                                                                                                                   >  2:43896/1‑150 (MQ=255)                       gcatggccccggcagaagacgagcacaacgaattgaccgatgggctcaagaccgatgaggaagaaacctggaagaagaaggtgatctacgactgattggttaacctgtgccggcctcttcgcttGTAGGAGCGGGTTTACCCGCGAagag                                                                                                                     <  1:379965‑M2/26‑1 (MQ=255)                           aGGGGGCGGCGTTCTGTTCTCGCTGTGGGAAGGGGTTGTTTCCCAAAGCGTGATTGCGCCTGCGCCGGCCTCTTCGCGGTTAAATCCGCTCCTACCGGCGCGGACCGTGCGCGGTCCCCCGTCGGAGCGGGTTTACCCGCGAAGAGGCCg                                                                                                                 <  1:349465/150‑1 (MQ=39)                            gagtgcggcgATCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCgg                                                                                                                <  2:119464/146‑1 (MQ=255)                                        ctgttctCGCTGTGGGATTGGTTTGTTTCCCAAAGCGTGATAGCGCCTGCGCGGGCCTCTTCGCGGTTAAATCCGCCCCTACCGGCGCGGACCGTGCGCGGTCCCCCGTCGGAGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAg                                                                                                    <  2:81243/150‑1 (MQ=38)                                                           gggTTGTTCCCAAATGGGTTATTGGGCCTGCGCGGGCCTCTTCGCGGTTAAACCCGCTCCAACTGGGGGGGACGGTGGGGGGGCCCCCTTAGGTGGGGTTTTACCCGCTAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAtttttt                                                                                 <  1:337138/150‑1 (MQ=255)                                                                            gatgaggaagaaacctggaagaagaaggtgatctacgactgattggttaacctgtgccggcctcttcgcttGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                                   >  2:12987‑M2/72‑99 (MQ=255)                                                                            gatgaggaagaaacctggaagaagaaggtgatctacgactgattggttaacctgtgccggcctcttcgcttGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                                   >  2:12993‑M2/72‑99 (MQ=255)                                                                            gatgaggaagaaacctggaagaagaaggtgatctacgactgattggttaacctgtgccggcctcttcgcttGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                                   <  1:12987‑M2/28‑1 (MQ=255)                                                                            gatgaggaagaaacctggaagaagaaggtgatctacgactgattggttaacctgtgccggcctcttcgcttGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                                   <  1:12993‑M2/28‑1 (MQ=255)                                                                                                     cgcgGGTAAATCCGCTCCTACAGGCGTGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCa                                       >  2:250257/1‑150 (MQ=255)                                                                                                        gggTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAcc                                                                      >  1:62178/1‑116 (MQ=255)                                                                                                        gggTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAcc                                                                      <  2:62178/116‑1 (MQ=255)                                                                                                                    tCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACAcc                        >  1:285579/1‑150 (MQ=255)                                                                                                                           ggCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAg                 <  2:149469/150‑1 (MQ=255)                                                                                                                                    cGTGCGCGGTCCCCCGTTGTAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGAc        <  1:285430/150‑1 (MQ=255)                                                                                                                                         gtccgcgcctGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAac                           >  2:212209‑M2/11‑126 (MQ=255)                                                                                                                                         gtccgcgcctGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAac                           <  1:212209‑M2/116‑1 (MQ=255)                                                                                                                                          cgGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACg  >  2:231339/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                             GACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACG  >  NC_002947/4898256‑4898542

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