New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 3981024	32 (0.800)	4 (0.110)	4/320	NT	11.0%	coding (678/1068 nt)	cobW	Protein CobW
	?	AE015451	= 3981044	34 (0.890)					coding (658/1068 nt)	cobW	Protein CobW

GCCGCATCGGCTTGCCCGGGATGGCGGCGAAGCCTTTGGCGCGCAGGATACCGAACTCCACCACCAACTGGGTCAATGCATCGAGCAGCAGGCTTTCGTCGGCTTCAGGCAGGTCGATGGAGATGGAGTCGAAAGCGTCGTGGTCATGATCGTCATGATCGTCGCCATCGTGGTGCGAATCGTGGTGGGTTCGGCGGCCATCGATATGCGCCTCGGACTCGGCACCCACGCCCAGCAGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/3980786‑3981024 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccagcagcGGCAAGCTGCCACTGGAGGTGCTGCTGGGCGTGGGTGCCGAGTCCGAGGCGCATATCGATGGCCGCCGAACCCACCACGATTCGCACCACGATGGCGACGATCATGACGATCATGACC  <  AE015451/3981044‑3980927                                                                                                                                                                                                                                                                                                                                                                       GCCGCATCGGCTTGCCCGGGATGGCGGCGAAGCCTTTGGCGCGCAGGATACCGAACTCCACCACCAACTGGGTCAATGCATCGAGCAGCAGGCTTTCGTCGGCTTCAGGCAGGTCGATGGAGATGGAGTCGAAAGCGTCGTGGTCATGATCGTCATGATCGTCGCCATCGTGGTGCGAATCGTGGTGGGTTCGGCGGCCATCGATATGCGCCTCGGACTCGGCACCCACGCCCAG                                                                                                                            <  1:722486/235‑1                                 CCTTTGGCGCGCAGGATACCGAACTCCACCACCAACTGGGTCAATGCATCGAGCAGCAGGCTTTCGTCGGCTTCAGGCAGGTCGATGGAGATGGAGTCGAAAGCGTCGTGGTCATGATCGTCATGATCGTCGCCATCGTGGTGCGAATCGTGGTGGGTTCGGCGGCCATCGATATGCGCCTCGGACTCGGCACCCACGCCC                                                                                                                              <  1:230731/201‑1                                                                      GGGTCAATGCATCGAGCAGCAGGCTTTCGTCGGCTTCAGGCAGGTCGATGGAGATGGAGTCGAAAGCGTCGTGGTCATGATCGTCATGATCGTCGCCATCGTGGTGCGAATCGTGGTGGGTTCGGCGGCCATCGATATGCGCCTCGGACTCGGCACCCACGCCCAG                                                                                                                            >  1:606690/1‑166                                                                      GGGTCAATGCATCGAGCAGCAGGCTTTCGTCGGCTTCAGGCAGGTCGATGGAGATGGAGTCGAAAGCGTCGTGGTCATGATCGTCATGATCGTCGCCATCGTGGTGCGAATCGTGGTGGGTTCGGCGGCCATCGATATGCGCCTCGGACTCGGCACCCACGCCCAG                                                                                                                            <  2:606689/166‑1                                                                                         CAGGCTTTCGTCGGCTTCAGGCAGGTCGATGGAGATGGAGTCGAAAGCGTCGTGGTCATGATCGTCATGATCGTCGCCATCGTGGTGCGAATCGTGGTGGGTTCGGCGGCCATCGATATGCGCCTCGGACTCGGCACCCACGCCCAGCAGCGGCAAG                                                                                                                  <  2:527063/157‑1                                                                                                                                                                                     CGTGGTGGGTTCGGCGGCCATCGATATGCGCCTCGGACTCGGCACCCACGCCCAGCAGCGGCAAGCTGCCACTGGAGGTGCTGCTGGGCGTGGGTGCCGAGTCAGAGGCGCATATCGATGGCCGCCGAACCCACCAC                                          >  1:723182/1‑137                                                                                                                                                                                                                GCGCCTCGGACTCGGCACCCACGCCCAGCAGCGGCAAGCTGCCACTGGAGGTGCTGCTGGGCGTGGGTGCCGAGTCCGAGGCGCATATCGATGGCCGCCGAACCCACCACGATTCGCACCACGATGGCGACGATCATGACGATCATGACC  <  2:723181/150‑1                                                                                                                                                                                                                      CGGACTCGGCACCCACGCCCAGCAGCGGCAAGCTGCCACTGGAGGTGCTGCTGGG                                                                                           >  1:555552/1‑55                                                                                                                                                                                                                                                                                                                                                                       GCCGCATCGGCTTGCCCGGGATGGCGGCGAAGCCTTTGGCGCGCAGGATACCGAACTCCACCACCAACTGGGTCAATGCATCGAGCAGCAGGCTTTCGTCGGCTTCAGGCAGGTCGATGGAGATGGAGTCGAAAGCGTCGTGGTCATGATCGTCATGATCGTCGCCATCGTGGTGCGAATCGTGGTGGGTTCGGCGGCCATCGATATGCGCCTCGGACTCGGCACCCACGCCCAGCAGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/3980786‑3981024 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccagcagcGGCAAGCTGCCACTGGAGGTGCTGCTGGGCGTGGGTGCCGAGTCCGAGGCGCATATCGATGGCCGCCGAACCCACCACGATTCGCACCACGATGGCGACGATCATGACGATCATGACC  <  AE015451/3981044‑3980927

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