New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4165491 =	62 (1.470)	3 (0.080)	3/250	NT	5.1%	coding (64/858 nt)	murI	glutamate racemase
	?	NC_000913	4165527 =	52 (1.300)					coding (100/858 nt)	murI	glutamate racemase

CTCAACAATAAACGCTTCGCTTTTTTCGCCATACGGGAAAGCGACGTTATCGAAAGCATAAATGTAATGGAGATCCGGTAAGAGATGCCGGATCTCGTCATAGACCGACAACCCACCGACGCCGGAGTCAAACACCAGCACGGTGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4165637‑4165491 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggtgggTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGAAAAAAGCGAAGCGTTTATTGTTGAGCGAGTGGTGGCAATTGTCACCGCGGTGCA  >  NC_000913/4165527‑4165666                                                                                                                                                                                                                                                                                                 CTCAACAATAAACGCTTCGCTTTTTTCGCCATACGGGAAAGCGACGTTATCGAAAGCATACATGTAATGGAGATCCGGTAAGAGATGCCGGATCTCGTCATAGACCGACAACCCACCGCCGCAGGAGGCAAACACCAGCACGG                                                                                                                                                  >  1:670986/1‑143    AACAATAAACGCTTCGCTTTTTTCGCCATACGGGAAAGCGACGTTATCGAAAGCATAAATGTAATGGAGATCCGGTAAGAGATGCCGGATCTCGTCATAGACCGACAACCCACCGACGCCGGAGTCAAACACCAGCACGGTGG                                                                                                                                               >  2:565838/1‑143    AACAATAAACGCTTCGCTTTTTTCGCCATACGGGAAAGCGACGTTATCGAAAGCATAAATGTAATCGAGATCCGGTAAGAGATGCCGGATCTCGTCATAGACCGACAACCCAGCGACGCCCGAGGCAAAGACCAGCACagagg                                                                                                                                               >  2:822059/1‑138                                                                                                              AACCCACCGACGCCGGAGTCAAACACCAGCACGGTGGGTTGTCGTTCTATGACGAGATCCGGCATCTCTTACCCTATCTCCATTACATTTATTCTTTCGATAACGTCCCTTTCCCGTATGGCGAAAAAAGCGAAGCGTTTATT                                     >  2:208975/1‑143                                                                                                                CCCACCGACGCCGGAGTCAAACACCAGCACGGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGAGCTCCATTACCTTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGAAAAAAGCGGAGCGTGTATTGT                                   >  1:362810/1‑143                                                                                                                 CCACCGACGCCGGAGTCAAACACCAGCAGGGATGGTTGACGGGCTAGGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTGATGCTTTCGATAACGTCGCTTTCCCGTATGGCGACAAAATCGAAGCGTGTATTGct                                  >  1:310786/1‑141                                                                                                                                             gatcGGGTTG‑CGGCCTTTGACGAGAGCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCGTTCCCGTAGGGCGAAAAAAGCGAAGCGTTTATTGTTGAGCGAGTGGTGGCAAGTGTCACCGCGGT     <  2:303141/139‑1                                                                                                                                              GGTGGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGAAAAAAGCGAAGCGTTTATTGTTGAGCGAGTGGTGGCAATTGTCACCGCGGT     <  1:325621/143‑1                                                                                                                                                 GGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGAAAAAAGCGAAGCGTTTATTGTTGAGCGAGTGGTGGCAATTGTCACCGCGGTGCA  <  1:164981/143‑1                                                                                                                                                 GGGTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGAAAAAAGCGAAGCGTTTATTGTTGAGCGAGTGGTGGCAATTGTCACCGCGGTGCA  <  2:203018/143‑1                                                                                                                                                                                                                                                                                                 CTCAACAATAAACGCTTCGCTTTTTTCGCCATACGGGAAAGCGACGTTATCGAAAGCATAAATGTAATGGAGATCCGGTAAGAGATGCCGGATCTCGTCATAGACCGACAACCCACCGACGCCGGAGTCAAACACCAGCACGGTGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4165637‑4165491 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggtgggTTGTCGGTCTATGACGAGATCCGGCATCTCTTACCGGATCTCCATTACATTTATGCTTTCGATAACGTCGCTTTCCCGTATGGCGAAAAAAGCGAAGCGTTTATTGTTGAGCGAGTGGTGGCAATTGTCACCGCGGTGCA  >  NC_000913/4165527‑4165666

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