New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 404794	39 (0.780)	4 (0.090)	3/264	7.4	8.7%	coding (1092/1116 nt)	yaiC	diguanylate cyclase, cellulose regualtor
	?	NC_000913	= 404817	47 (1.000)					coding (1115/1116 nt)	yaiC	diguanylate cyclase, cellulose regualtor
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/404645‑404794 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccggacgtCAGGCCGCCACTTCGGTGCGGTTACGTCCGG  <  NC_000913/404817‑404787                                                                                                                                                                                       GGGCTAAATACATTACGTTTGCCGACTACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCAGCTTTCCCCACAAATGACTCAATATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGC                                        <  2:504065/143‑1  GGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCC                                       >  2:38362/1‑143       aatgcgTTACGCTTGCCGAATACGCCACAGGTAACTTTACGGATTAGGGTGGGGGTTGCGCCGCTGAACCCACAAAGGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACG                                  <  2:172020/137‑1          ACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAACTGAGTCACTATCGTGAGTGGTTGAAATCGGCCGCGTTGGCGCTTTACAAAGCCAAGAAAGCCGGACTTCA                               >  2:51789/1‑143                                       AACGTTACGCGTGAGTGGGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGATTGGTTGAAAGCGGCAGATGTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTCAGGCCGCCACTTCGGTGCGGTTACGTCCGG  <  1:749587/143‑1                                                                              ACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTCAGGCCGCCACTTCGGTGCGGTTACGTC     <  1:765826/101‑1                                                                              ACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTCAGGCCGCCACTTCGGTGCGGTTACGTC     >  2:765826/1‑101                                                                                                                                                                                       GGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/404645‑404794 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccggacgtCAGGCCGCCACTTCGGTGCGGTTACGTCCGG  <  NC_000913/404817‑404787

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