New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3587686 =	87 (0.940)	5 (0.060)	4/272	NT	5.4%	coding (428/744 nt)	ugpQ	glycerophosphodiester phosphodiesterase
	?	NC_000913	3587731 =	89 (0.980)					coding (383/744 nt)	ugpQ	glycerophosphodiester phosphodiesterase

GGCAAAATGGTGGCGCTGGCGGCACGCGAGCTGTGGGCCGGTATGACGCCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3587738‑3587686 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCA  >  NC_000913/3587731‑3587874                                                                                                                                                                                                       GGCAAAATGGTGGCGCTGGCGGCACGCGAGCTGTGGGCCGGTATGACGCCGCCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCA                                                 >  1:753663/1‑150  GCAAAATGGTGGCGCTGGCGGCACGCGAGCTGTGGGCCGGTATGACGCCGCCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAA                                                >  2:1032420/1‑150   CAAAATGGTGGCGCTGGCGGCACGCGAGCTGTGGGCCGGTATGACGCCGCCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAAC                                               >  2:322530/1‑150     AAATGGTGGCGCTGGCGGCACGCGAGCTGTGGGCCGGTATGACGCCGCCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGG                                             >  2:13913/1‑150     AAATGGTGGCGCTGGCGGCACGCGAGCTGTGGGCCGGTATGACGCCGCCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGG                                             >  2:8063/1‑150                                                GCCGCCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCA  <  1:322530/150‑1                                                                                                                                                                                                       GGCAAAATGGTGGCGCTGGCGGCACGCGAGCTGTGGGCCGGTATGACGCCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3587738‑3587686 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGGCTCACCTTTAAACATTTTGCTGTACCAACTGCCCGCA  >  NC_000913/3587731‑3587874

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