New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 598256	7 (0.520)	7 (0.540)	6/256	NT	42.7%	coding (778/1224 nt)	cusB	copper/silver export system membrane fusion protein
	?	NC_000913	636166 =	12 (0.920)					coding (404/1221 nt)	ybdN	putative PAPS reductase/DUF3440 domain‑containing protein YbdN

TCACCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTT  >  NC_000913/598131‑598528                                                                                                                              |                                                                                                                                                                                                                                                                                 tCACCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATg                                                                                                                                                                                                                                                                               >  2:56896/1‑129 (MQ=255)    ccGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCt                                                                                                                                                                                                                                                                            <  1:56896/129‑1 (MQ=255)                                    gCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGcagaaaaagtcagggtcggttatcgcatcttggggaggttgacgtacccattcgacatcaggttcccagcactgccattcgggttggtattgtgaaagggaattttgcgtcgtaagcgggagcgcaacccagtaaaactcttcgatgacatcggtgtacaactcgcgcagggactgaacatag                                                                                           >  1:74947‑M1/1‑91 (MQ=255)                                    gCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGcagaaaaagtcagggtcggttatcgcatcttggggaggttgacgtacccattcgacatcaggttcccagcactgccattcgggttggtattgtgaaagggaattttgcgtcgtaagcgggagcgcaacccagtaaaactcttcgatg                                                                                                                               >  1:139490‑M1/1‑91 (MQ=255)                                                                             ggTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGcagaaaaagtcagggtcggttatcgcatcttggggaggttgacg                                                                                                                                                                                                                                      >  2:242668‑M1/1‑50 (MQ=255)                                                                                gtgGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGcagaaaaagtcagggtcggttatcgcatcttggggaggttgacgtac                                                                                                                                                                                                                                   <  1:242669‑M1/94‑48 (MQ=255)                                                                                    gTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCaaa                                                                                                                                                                                                                    >  2:188881/1‑105 (MQ=255)                                                                                       aCTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATgg                                                                                                                                                                                                                 <  1:188882/105‑1 (MQ=255)                                                                                         tgctgcGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCaaa                                                                                                                                                                                                                    >  2:82158/1‑100 (MQ=255)                                                                                            tgcGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATgg                                                                                                                                                                                                                 <  1:82159/100‑1 (MQ=255)                                                                                               gATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGa                                                                                                                      >  1:96490/1‑188 (MQ=255)                                                                                                 tCCCGGAGTCTATCGCCTGGCTGGTGAAAGcagaaaaagtcagggtcggttatcgcatcttggggaggttgacgtacccattcgacatcaggttcccagcactgccattc                                                                                                                                                                                                  >  1:25974‑M1/1‑30 (MQ=255)                                                                                                  cccGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACg                                                                                                                   <  2:96489/188‑1 (MQ=255)                                                                                                    cGGAGTCTATCGCCTGGCTGGTGAAAGcagaaaaagtcagggtcggttatcgcatcttggggaggttgacgtacccattcgacatcaggttcccagcactgccattcgggttggtattgtgaaagggaattttgcgtcgtaagcgggagcgcaacccagtaaaactcttcgatgacatcggtgtacaactcgcgcagggactgaaaatagttaatggtgcaagagaattgcgcctcccagtcgataaacaaaacacagatttttttgcccattt                           >  1:52309‑M1/1‑27 (MQ=255)                                                                                                    cGGAGTCTATCGCCTGGCTGGTGAAAGcagaaaaagtcagggtcggttatcgcatcttggggaggttgacgtacccattcgacatcaggttcccagcactgccattcggg                                                                                                                                                                                               <  2:25974‑M1/110‑84 (MQ=255)                                                                                                                            aaGcagaaaaagtcagggtcggttatcgcatcttggggaggttgacgtacccattcgacatcaggttcccagcactgc                                                                                                                                                                                                       >  1:192662‑M1/1‑3 (MQ=255)                                                                                                                             aGcggaaaaattcagcgtcggttttcgcttcttgtgtgtgttgacggacccattcgacatcaggttccccgcactgccattcgggttggtagtgtgaaagggaattttgcgtcgtaagcgggagcgcaacccagtaaaactcttcgatgacatcggtgtacaactcgcgcagggactgaacatagttaatggtgcaagagaattgcgcctcccagtcgataaacaaaacgcagatttttttgcccatttgtcgggcaagttcggctgtcagatg  <  2:52309‑M1/274‑273 (MQ=255)                                                                                                                              |                                                                                                                                                                                                                                                                                 TCACCGCGTTTGATCTGCGCGCGGGAATGAATATCGCCAAAGATAACGTGGTAGCGAAAATTCAGGGTATGGACCCGGTGTGGGTCACTGCTGCGATCCCGGAGTCTATCGCCTGGCTGGTGAAAGATGCCTCGCAGTTTACCCTCACCGTTCCGGCGCGACCGGATAAAACACTCACCATCCGCAAATGGACGCTGCTACCTGGCGTGGATGCCGCGACCCGCACGCTGCAGCTGCGTCTGGAAGTCGACAACGCCGACGAGGCGCTAAAACCGGGCATGAACGCCTGGTTGCAACTCAACACCGCCAGCGAACCGATGCTGCTCATTCCGTCACAAGCGCTGATTGATACCGGCAGCGAACAGCGGGTGATTACCGTTGATGCCGACGGGCGCTTT  >  NC_000913/598131‑598528

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