New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AM260480	= 2637763	15 (0.680)	4 (0.200)	4/188	1.6	23.6%	coding (1347/1575 nt)	h16_B2327	GGDEF domain protein
	?	AM260480	= 2637780	12 (0.590)					coding (1364/1575 nt)	h16_B2327	GGDEF domain protein

CTGCATGGCCACGTGGCCGGCGACGACTGCCTGCGGCAGATTGCCGGCGCCCTGGCCCGCACGCTCTCCGGCGAGCAGTTCTGCGCCCGGCTGGGCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCACGCGTGAAGAGACCCGTGACGTGGCGCGCCGCATGGTGGCCGTTGTGCGCGATCTCGGCATCGCGCACAACGGCACGCCGGTGCAGCCGGTGGTCACCATCAGCGTGGGCGCGGTC  >  AM260480/2637632‑2637877                                                                                                                                    |                                                                                                                   cTGCATGGCCACGTGGCCGGCGACGACTGCCTGCGGCAGATTGCCGGCGCCCTGGCCCGCACGCTCTCCGGCGAGCAGTTCTGCGCCCGGCTGGGCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCAcgcg                                                                                                                  <  1:707912/134‑1 (MQ=255)                  cGGCGACGACTGCCTGCGGCAGATTGCCGGCGCCCTGGCCCGCACGCTCTCCGGCGAGCAGTTCTGCGCCCGGCTGGGCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCACGCGTGAag                                                                                                             >  1:366477/1‑122 (MQ=255)                  cGGCGACGACTGCCTGCGGCAGATTGCCGGCGCCCTGGCCCGCACGCTCTCCGGCGAGCAGTTCTGCGCCCGGCTGGGCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCACGCGTGAag                                                                                                             <  2:366477/122‑1 (MQ=255)                                             cGGCGCCCTGGCCCGCACGCTCTCCGGCGAGCAGTTCTGCGCCCGGCTGGGCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCACGCGTGAAGAGACCCGTGACGTGGCGCGCCGCATGGTGGCCGTTGTgc                                                                      <  2:827890/134‑1 (MQ=255)                                                                     cGGCGAGCAGTTCTGCGCCCGGCTGGGCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCACGCGTGAAGAGACCCGTGAcgtggc                                                                                             >  2:594957‑M2/1‑81 (MQ=255)                                                                     cGGCGAGCAGTTCTGCGCCCGGCTGGGCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCACGCGTGAAGAGACCCGTGAcgtggc                                                                                             <  1:594957‑M2/87‑7 (MQ=255)                                                                     cGGCGAGCAGTTCTGCGCCCGGCTGGGCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCACGCGTGAAGAGACCCGTGACGTGGCGCGCCGCATGGTGGCCGTTGTGCGCGATCTCGGCATCGCGCACAACg                                              >  1:722871/1‑134 (MQ=255)                                                                              gTTCTGCGCCCGGCTGGGCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCACGCGTGAag                                                                                                             >  2:592941/1‑62 (MQ=255)                                                                              gTTCTGCGCCCGGCTGGGCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCACGCGTGAag                                                                                                             <  1:592941/62‑1 (MQ=255)                                                                               ttCTGCGCCCGGCTGGGCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCACGtcacgggtctcttcacgcgtggc                                                                                             >  1:83127‑M1/1‑54 (MQ=255)                                                                               ttCTGCGCCCGGCTGGGCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCACGtcacgggtctcttcacgcgtggc                                                                                             <  2:83127‑M1/77‑24 (MQ=255)                                                                                        cGGCTGGGCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCACGCGTGAAGAGACCCGTGACGTGGCGCGCCGCAtggtg                                                                                >  1:721732/1‑81 (MQ=255)                                                                                        cGGCTGGGCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCACGCGTGAAGAGACCCGTGACGTGGCGCGCCGCAtggtg                                                                                <  2:721732/81‑1 (MQ=255)                                                                                             gggCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCACGtcacgggtctcttcacgcgt                                                                                                >  1:1026993‑M1/1‑40 (MQ=255)                                                                                             gggCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCACGtcacgggtctcttcacgcgt                                                                                                <  2:1026993‑M1/60‑21 (MQ=255)                                                                                             gggCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCACGCGTGAAGAGACCCGTGACGTGGCGCGCCGCATGGTGGCCGTTGTGCGCGATCTCGGCATCGCGCACAACGGCACGCCGGTGCAGCCGGTGGTCa                      >  1:1062237/1‑134 (MQ=255)                                                                                                         aTTCGCGGCGGTGCTGCCCGGCGCCACGCGTGAAGAGACCCGTGACGTGGCGCGCCGCATGGTGGCCGTTGTGCGCGATCTCGGCAt                                                         >  1:578714/1‑87 (MQ=255)                                                                                                         aTTCGCGGCGGTGCTGCCCGGCGCCACGCGTGAAGAGACCCGTGACGTGGCGCGCCGCATGGTGGCCGTTGTGCGCGATCTCGGCAt                                                         <  2:578714/87‑1 (MQ=255)                                                                                                         aTTCGCGGCGGTGCTGCCCGGCGCCACGCGTGAAGAGACCCGTGACGTGGCGCGCCGCATGGTGGCCGTTGTGCGCGATCTCGGCATCGCGCACAACGGCACGCCGGTGCAGCCGGTGGTCACCATCAGCGTg           >  2:142522/1‑133 (MQ=255)                                                                                                         aTTCGCGGCGGTGCTGCCCGGCGCCACGCGGGAAGAGACCCGTGACGTGGCGCGCCGCATGGTGGCCGTTGTGCGCGATCTCGGCATCGCGCACAACGGCACGCCGGTGCAGCCGGTGGTCACCATCAGCGTg           <  1:142522/133‑1 (MQ=255)                                                                                                                 gggTGCGCCCCGGCGCCACGCGTGAAGAGCCCCGTGCGGGGGCGCGCCGCATGGTGGCCTTTGTGCGCGATCTCGGCATCGCGCACAACGGCACGCCGGTGCAGCCGGTGGTCACCATCAGCGTGGGCGCGGTc  <  2:722871/133‑1 (MQ=255)                                                                                                                               gcCACGtcacgggtctcttcacgcgtggcgccgggcagcaccgccgcgaattcctcgccgcccagccgggcgcagaactgctcgccg                                   <  1:594957‑M1/87‑82 (MQ=255)                                                                                                                               gcCACGtcacgggtctcttcacgcgtggcgccgggcagcaccgccgcgaattcctcgccgcccagccgggcgcagaactgctcgccg                                   >  2:594957‑M1/1‑6 (MQ=255)                                                                                                                               gcCACGCGTGAAGAGACCCGTGAcgtggcgccgggcagcaccgccgcgaattcctcgccgcccagccgggcgcagaa                                             >  1:83127‑M2/1‑23 (MQ=255)                                                                                                                               gcCACGCGTGAAGAGACCCGTGAcgtggcgccgggcagcaccgccgcgaattcctcgccgcccagccgggcgcagaa                                             <  2:83127‑M2/77‑55 (MQ=255)                                                                                                                                  aCGCGTGAAGAGACCCGTGAcgtggcgccgggcagcaccgccgcgaattcctcgccgccc                                                           <  2:1026993‑M2/60‑41 (MQ=255)                                                                                                                                  aCGCGTGAAGAGACCCGTGAcgtggcgccgggcagcaccgccgcgaattcctcgccgccc                                                           >  1:1026993‑M2/1‑20 (MQ=255)                                                                                                                                    |                                                                                                                   CTGCATGGCCACGTGGCCGGCGACGACTGCCTGCGGCAGATTGCCGGCGCCCTGGCCCGCACGCTCTCCGGCGAGCAGTTCTGCGCCCGGCTGGGCGGCGAGGAATTCGCGGCGGTGCTGCCCGGCGCCACGCGTGAAGAGACCCGTGACGTGGCGCGCCGCATGGTGGCCGTTGTGCGCGATCTCGGCATCGCGCACAACGGCACGCCGGTGCAGCCGGTGGTCACCATCAGCGTGGGCGCGGTC  >  AM260480/2637632‑2637877

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