New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 302638	11 (0.940)	3 (0.260)	3/304	1.0	19.4%	intergenic (‑65/‑353)	paoA/yagU	aldehyde dehydrogenase, Fe‑S subunit/DUF1440 domain‑containing inner membrane protein YagU
	?	NC_000913	= 1924823	14 (1.200)					coding (147/375 nt)	yobA	copper‑binding protein
Rejected: Frequency below/above cutoff threshold.

GGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTC‑ATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGC‑TTTTTTAAACGTTAAGCATAGTCGGCGACAGCAGGAGGATTATTCTTAAAACGCGAATGCAGCTATGAGCCCGGCTAATAAATTCATGTTCGCCGGGATGTTGATGATGATGGGAGCTGGTTTATTGCTTTAGTTGTACGATGCAAAAACCAATAAGGAAACCTGTGATTTTCAGCTCTACATCACCCTGCAAATCTCTGTCACTTCTAATATAAAAATAGGGAGAA  >  NC_000913/302401‑302858                                                                                                                                                                                                                                                |                                                                                                                                                                                                                             ggCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGGGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGTTTGCTC‑ATACCAGGCTTCCGGTATTCGAAAGGGGAAAATAAGTGTGCGTTGCGAGAGTGCTGATGC‑TTTTTTAAACGTTAAGCa                                                                                                                                                                                                                   <  1:157676/249‑1 (MQ=255)        gtgtAGCTGCCGGAACGCTTGCCGCCAGCGTAGAAAGAGGATAAACCACGGCGGTCGACGCCGGTGCGGGGGGCCCCTTTATCAAATAGCCACGGGGCAAAACTAAATAGTGGGGGTCGGGGTTCCCAACCCCAATATATTCGGGGGATTCGACTTGGTTG‑TCGCTAACAGGGCTCCGGGATTGTAAAGGGGAAACAAAGAGAGGGTGGCGAGAGGGCAGGTGG‑TTTTTTaaaagaaaagcaaagtagggg                                                                                                                                                                                                          >  1:92347/1‑233 (MQ=255)          gtAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTC‑ATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGC‑TTTTTTAAACGTTAAGCATAGTCGGCGNc                                                                                                                                                                                                        <  1:215817/251‑1 (MQ=255)                                                                                                gcgACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTC‑ATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGC‑TTTTTTAACCGTTACGCATAGTCGGCGACAGCAGGAGGATTATTCTTAAAACGCGAATGCAGCTATGAGCCCGGCTAATAAATTCATGt                                                                                                                                            >  1:37369/1‑225 (MQ=255)                                                                                                       gTCAAACTTAAATCGTGCGGCTCGTGATTCACAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTC‑ATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGC‑TTTTTTAAACGTTACGCATAGTCGGCGACAGCAGGAGGATTATTCTTAAAACGCGAATGCAGCTATGAGCCCGGCTAATAAATTCATGTTCGCCGGGATGTTGATGATGATGGGAGCTGGtt                                                                                                           >  1:234302/1‑251 (MQ=255)                                                                                                             cTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTC‑ATACGAGGCCTCCTGTATTCTAAAGGGGAAAATAAGATTGCGTTGCGAGAGTGCTGATGC‑TTTTTTAAACGTTAAGCATAGTCGGCGAc                                                                                                                                                                                                        <  1:151245/152‑1 (MQ=255)                                                                                                                     gTGCGGCTCGTGCTTCCCAACCCGATTATCTGCGGGGTATTCGCCGTGGTTGCTC‑ATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCGTTGTTCTAACGGTAAGCATAGTCGGCTACAGCAGGAGGATTATTCTTAAAACGCGAATGCAGCTATGAGCCCGGCGAATAAATTCATGTGCGCCGGGATGTTGATGATGATGGGAGCTTGTTTATTGCTTTAGtt                                                                                              <  1:145627/251‑1 (MQ=255)                                                                                                                                     ccAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTC‑ATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGC‑TTTTTTAtttaaacttctcggaaggtgttgaaacaggattcagtggtgcaaaaatcacggggccaaaaaacgaaaatattaaaacattacc                                                                                                                                          <  1:64182‑M1/190‑85 (MQ=255)                                                                                                                                     ccAACCCGATTATCTTCGGGGGATTCNCCTTGGNNGCNC‑ATACCAGGCCTCCGGTATNCTNNNGNNNANCATAANAGNNCGTTGCGAGNGTGCTGATGC‑TTTTTTAtttaaacttctcggaaggtgttgaaacaggattcagnggtgcaaaanncacggggccaaaaaacgaaaanattnnannantacn                                                                                                                                          >  2:64182‑M1/1‑106 (MQ=255)                                                                                                                                                           aTTCGCCTTGGTGGCTC‑ATATCAGGGCTCCGGTATTCTAAAGGGGAAAATAAGAGTTCGTTTCGAGAGTGCTTATGC‑TTTTGTAAACGTTAAGCATAGTCGACGTGAGCAGGAGGATTATTCTTAATACGTGATGGCAGCTAGGAGTCCGTCTAATAAATTCATGTTCGCTGTGATGTTGATGATTATTGGATCTGGTTTATTGCTTTAGTTGTACGATGCAAAAACCAATAAGGAAACCTGTGATTTTc                                                        <  1:281507/250‑1 (MQ=255)                                                                                                                                                                   tGGTTGCTC‑ATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGC‑TTTTTTAAACGTTAAGCATAGTCGGCGACAGCAGGAGGATTATTCTTAAAACGCGAATGCAGCTATGAGCCCGGCTAATAAATTCATg                                                                                                                                             <  1:126501/157‑1 (MQ=255)                                                                                                                                                                         cTC‑ATACCAGGCCTCCGGTATTCTCAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGC‑TTTTTTAtac                                                                                                                                                                                                                           >  1:133046‑M1/1‑70 (MQ=255)                                                                                                                                                                                                                  gTGCGTTGCGAGAGTGCTGACGC‑TTTTTTAACCGTTAAGCATAGTCGGCGACAGCAGGAGGATTATTCTTAAAACGCGAATGCAGCTATGCGCCCGGATAATAAATTCATGTTCGCCGGGATGTTGATGATGATGGGAGCTGGTTTATTGATTTAGATGTACGATGCAAAAACCAATAAGGAAACCTGTGATTTTCAGCTCTACATCACCCTGCAAATCTCTGTCACTTCTAATATAAAAATAGGGAGaa  >  1:100647/1‑250 (MQ=255)                                                                                                                                                                                                                                  cTGATNC‑TNNTTNNACANTTTAGCATAGTCGGCGACAGCAGGAGGATTATTCtn                                                                                                                                                                                      >  2:154850/1‑53 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                             GGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTC‑ATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGC‑TTTTTTAAACGTTAAGCATAGTCGGCGACAGCAGGAGGATTATTCTTAAAACGCGAATGCAGCTATGAGCCCGGCTAATAAATTCATGTTCGCCGGGATGTTGATGATGATGGGAGCTGGTTTATTGCTTTAGTTGTACGATGCAAAAACCAATAAGGAAACCTGTGATTTTCAGCTCTACATCACCCTGCAAATCTCTGTCACTTCTAATATAAAAATAGGGAGAA  >  NC_000913/302401‑302858

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