New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3530659 =	65 (1.680)	3 (0.080)	3/256	NT	4.8%	intergenic (+9/+54)	hslO/yhgE	heat shock protein Hsp33/DUF4153 family putative inner membrane protein
	?	NC_000913	3530683 =	59 (1.640)					intergenic (+33/+30)	hslO/yhgE	heat shock protein Hsp33/DUF4153 family putative inner membrane protein
Rejected: Frequency below/above cutoff threshold.

TTACTCACAGCTGCGAAGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGAATGTGAATCTGAATGAGTGACAAAGCGCATCAGGTGTATTAGTCGTAAAATCCGGCAGAGCCCTCTCTGCCGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3530814‑3530659 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctctgccggaTTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTT  >  NC_000913/3530683‑3530827                                                                                                                                                                                                                                                                                                               TTACTCACAGCTGCGAAGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGAATGTGAATCTGAATGAGTGACAAAGCGCATCAGGTGTATTAGTCGTAAAATCCGGCAGAGCCCTCTC                                                                                                                                                          <  2:624098/149‑1     TCACAGCTGCGAAGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGAATGTGAATCTGAATGAGTGACAAAGCGCATCAGGTGTATTAGTCGTAAAATCCGGCAGAGCCCTCTCTGCC                                                                                                                                                      <  2:74703/149‑1                                                                                                       CAAAGCGCATCAGGTGTATTAGTCGTAAAATCCGGCAGAGAGGGCTCTGCCGGATTTTACG                                                                                                                                            <  1:225146/61‑1                                                                                                       CAAAGCGCATCAGGTGTATTAGTCGTAAAATCCGGCAGAGAGGGCTCTGCCGGATTTTACG                                                                                                                                            >  2:225146/1‑61                                                                                                                                    ATCCGGCAGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCT                       >  2:535270/1‑149                                                                                                                                                      TGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTC     <  1:39494/149‑1                                                                                                                                                      TGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGCCTTCGTCCCCAGCTTCCCCTCATTCCGAGCTGTGAGTAATTTTTCTTTC     >  1:265054/1‑149                                                                                                                                                         CGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCCCGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTT  <  2:171339/149‑1                                                                                                                                                         CGGATTTTACGACTAATACACCTGATGCGCTTAGTCGTATAATCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCCGCTGTGAGTAATTTTTCTTTCGTT  >  1:166576/1‑149                                                                                                                                                         CGGATTTGACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTT  <  1:374348/149‑1                                                                                                                                                                                                                                                                                                               TTACTCACAGCTGCGAAGGATGGAAAACTGGGGACGAAGCCTAAAGCGTGGCGCGATCTTGTAGTGGATGGTGAAAGGCTGAATGTGAATCTGAATGAGTGACAAAGCGCATCAGGTGTATTAGTCGTAAAATCCGGCAGAGCCCTCTCTGCCGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3530814‑3530659 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctctgccggaTTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTT  >  NC_000913/3530683‑3530827

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