New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	5650169 =	90 (1.260)	3 (0.040)	3/330	NT	3.4%	coding (1024/2127 nt)	PP_4959	diguanylate cyclase
	?	AE015451	5650216 =	84 (1.250)					coding (977/2127 nt)	PP_4959	diguanylate cyclase
Rejected: Frequency below/above cutoff threshold.

TGGACATGTCGGCAGTGCGTAACTTGCACCAGCAGCTGGAATACCAGGCCATTACTGACCCGCTAACCGGCCTGCTCAACCGCCGTGGTTTCTACCAGGCGGCTGAAGGTGCACTGTTGCGCAACGAGCGCTCTGACAAGGCGCAGGCGCTGATGTACATGGACCTGGATGGCTTCAAGCGCATCAACGACTCGCTCGGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCCGAGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/5650408‑5650169 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtggccgagcGAGTCGTTGATGCGCTTGAAGCCATCCAGGTCCATGTACATCAGCGCCTGCGCCTTGTCAGAGCGCTCGTTGCGCAACAGTGCACCTTCAGCCGCCTGGTAGAAACCACGGCGGTTGAGCAGGCCGGTTAGCGGGTCAGTAATGGCCTGGTATTCCAGCTGCTGGTGCAAGTTACGCACTGCCGACATGTCCAGCACGGTCACCACCATGGCCTGTT  >  AE015451/5650216‑5650432                                                                                                                                                                                                                                                                                                                                                                                                                                                                           TGGACATGTCGGCAGTGCGTAACTTGCACCAGCAGCTGGAATACCAGGCCATTACTGACCCGCTAACCGGCCTGCTCAACCGCCGTGGTTTCTACCAGGCGGCTGAAGGTGCACTGTTGCGCAACGAGCGCTCTGACAAGGCGCAGGCGCTGATGTACATGGACCTGGATGGCTTCAAGCGCATCAACGACTCGCTCGGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCC                                                                                                                                                                                                                               >  1:680969/1‑236   GACATGTCGGCAGTGCGTAACTTGCACCAGCAGCTGGAATACCAGGCCATTACTGACCCGCTAACCGGCCTGCTCAACCGCCGTGGTTTCTACCAGGCGGCTGAAGGTGCACTGTTGCGCAACGAGCGCTCTGACAAGGCGCAGGCGCTGATGTACATGGACCTGGATGGCTTCAAGCGCATCAACGACTCGCTCGGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCCG                                                                                                                                                                                                                              <  1:968541/235‑1       TGTCGGCAGTGCGTAACTTGCACCAGCAGCTGGAATACCAGGCCATTACTGACCCGCTAACCGGCCTGCTCAACCGCCGTGGTTTCTACCAGGCGGCTGAAGGTGCACTGTTGCGCAACGAGCGCTCTGACAAGGCGCAGGCGCTGATGTACATGGACCTGGATGGCTTCAAGCGCATCAACGACTCGCTCGGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCCG                                                                                                                                                                                                                              <  2:815161/231‑1            GCAGTGCGTAACTTGCACCAGCAGCTGGAATACCAGGCCATTACTGACCCGCTAACCGGCCTGCTCAACCGCCGTGGTTTCTACCAGGCGGCTGAAGGTGCACTGTTGCGCAACGAGCGCTCTGACAAGGCGCAGGCGCTGATGTACATGGACCTGGATGGCTTCAAGCGCATCAACGACTCGCTCGGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCCG                                                                                                                                                                                                                              >  2:137476/1‑226                                            CCAGGCCATTACTGACCCGCTAACCGGCCTGCTCAACCGCCGTGGTTTCTACCAGGCGGCTGAAGGTGCACTGTTGCGCAACGAGCGCTCTGACAAGGCGCAGGCGCTGATGTACATGGACCTGGATGGCTTCAAGCGCATCAACGACTCGCTCGGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCCGAG                                                                                                                                                                                                                            >  1:996500/1‑196                                            CCAGGCCATTACTGACCCGCTAACCGGCCTGCTCAACCGCCGTGGTTTCTACCAGGCGGCTGAAGGTGCACTGTTGCGCAACGAGCGCTCTGACAAGGCGCAGGCGCTGATGTACATGGACCTGGATGGCTTCAAGCGCATCAACGACTCGCTCGGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCCGAG                                                                                                                                                                                                                            <  2:996495/196‑1                                                   ATTACTGACCCGCTAACCGGCCTGCTCAACCGCCGTGGTTTCTACCAGGCGGCTGAAGGTGCACTGTTGCGCAACGAGCGCTCTGACAAGGCGCAGGCGCTGATGTACATGGACCTGGATGGCTTCAAGCGCATCAACGACTCGCTCGGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCC                                                                                                                                                                                                                               <  1:41113/186‑1                                                                  ACCGGCCTGCTCAACCGCCGTGGTTTCTACCAGGCGGCTGAAGGTGCACTGTTGCGCAACGAGCGCTCTGACAAGGCGCAGGCGCTGATGTACATGGACCTGGATGGCTTCAAGCGCATCAACGACTCGCTCGGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCCGAGC                                                                                                                                                                                                                           >  2:531902/1‑175                                                                                CCGCCGTGGTTTCTACCAGGCGGCTGAAGGTGCACTGTTGCGCAACGAGCGCTCTGACAAGGCGCAGGCGCTGATGTACATGGACCTGGATGGCTTCAAGCGCATCAACGACTCGCTCGGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCC                                                                                                                                                                                                                               <  2:352425/157‑1                                                                                 CGCCGTGGTTTCTACCAGGCGGCTGAAGGTGCACTGTTGCGCAACGAGCGCTCTGACAAGGCGCAGGCGCTGATGTACATGGACCTGGATGGCTTCAAGCGCATCAACGACTCGCTCGGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCCG                                                                                                                                                                                                                              <  1:365600/157‑1                                                                                                                                  GCTCTGACAAGGCGCAGGCGCTGATGTACATGGACCTGGATGGCTTCAAGCGCATCAACGACTCGCTCGGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCCG                                                                                                                                                                                                                              <  1:373693/108‑1                                                                                                                                  GCTCTGACAAGGCGCAGGCGCTGATGTACATGGACCTGGATGGCTTCAAGCGCATCAACGACTCGCTCGGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCCG                                                                                                                                                                                                                              >  2:373688/1‑108                                                                                                                                                    CGCTGCTGTACATGGACCTGGATGGCTTCAAGCGCATCAACGACTCGCTCGGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCCGAGCGAGTCGTTGATGCGCTTGAAGCCATCCAGGTCCATGTACATCAGCGCCT                                                                                                                                                                          <  1:212050/142‑1                                                                                                                                                                           GGCTTCAAGCGCATCAACGACTCGCTCGGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCCG                                                                                                                                                                                                                              >  1:780537/1‑67                                                                                                                                                                           GGCTTCAAGCGCATCAACGACTCGCTCGGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCCG                                                                                                                                                                                                                              <  2:780532/67‑1                                                                                                                                                                                                      GGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCCGAGCGAGTCGTTGATGCGCTTGAAGCCATCCAGGTCCATGTACATCAGCGCCTGCGCCTTGTCAGAGCGCTCGTTGCGCAACAGTGCACCTTCAGCCGCCTGGTAGAAACCACGG                                                                                                            >  2:1156792/1‑154                                                                                                                                                                                                      GGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCCGAGCGAGTCGTTGATGCGCTTGAAGCCATCCAGGTCCATGTACATCAGCGCCTGCGCCTTGTCAGAGCGCTCGTTGCGCAACAGTGCACCTTCAGCCGCCTGGTAGAAACCACGG                                                                                                            <  1:1156797/154‑1                                                                                                                                                                                                                                       GTGGCCGAGCGAGTCGTTGATGCGCTTGAAGCCATCCAGGTCCATGTACATCAGCGCCTGCGCCTTGTCAGAGCGCTCGTTGCGCAACAGTGCACCTTCAGCCGCCTGGTAGAAACCACGGCGGTTGAGCAGGCCGGT                                                                                           <  2:510535/138‑1                                                                                                                                                                                                                                         GGCCGAGCGAGTCGTTGATGCGCTTGAAGCCATCCAGGTCCATGTACATCAGCGCCTGCGCCTTGTCAGAGCGCTCGTTGCGCAACAGTGCACCTTCAGCCGCCTGGTAGAAACCACGGCGGTTGAGCAGGCCGGTTAGCGGGTCAGTAATGGCCTGGTATTCCAGCTGCTGGTGCAAGTTACGCACTGCCGACATGTCCAGCACGGTCACCACCATGGCCTGTT  <  1:1057656/225‑1                                                                                                                                                                                                                                                                                                                                                                                                                                                                           TGGACATGTCGGCAGTGCGTAACTTGCACCAGCAGCTGGAATACCAGGCCATTACTGACCCGCTAACCGGCCTGCTCAACCGCCGTGGTTTCTACCAGGCGGCTGAAGGTGCACTGTTGCGCAACGAGCGCTCTGACAAGGCGCAGGCGCTGATGTACATGGACCTGGATGGCTTCAAGCGCATCAACGACTCGCTCGGCCACGACGCGGGCGACCGGGTGTTGCGCTGGGTGGCCGAGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/5650408‑5650169 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtggccgagcGAGTCGTTGATGCGCTTGAAGCCATCCAGGTCCATGTACATCAGCGCCTGCGCCTTGTCAGAGCGCTCGTTGCGCAACAGTGCACCTTCAGCCGCCTGGTAGAAACCACGGCGGTTGAGCAGGCCGGTTAGCGGGTCAGTAATGGCCTGGTATTCCAGCTGCTGGTGCAAGTTACGCACTGCCGACATGTCCAGCACGGTCACCACCATGGCCTGTT  >  AE015451/5650216‑5650432

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