New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	4160673 =	10 (0.390)	3 (0.120)	3/310	NT	16.3%	intergenic (+39/+82)	PP_3662/PP_3663	conserved protein of unknown function/GGDEF domain protein
	?	AE015451	4160817 =	21 (0.840)					coding (1240/1302 nt)	PP_3663	GGDEF domain protein

CCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCAATGCTCAAGGTAACCGCCTGGGGCGTAGGCACAGTCTTGTCGGCGTACTGCTGGAACTGCTCACGCAGCCCACTCCCCAGTTCCATGACCCGTTCACTGCTCGCGCCACTGAGCAAGATGACGAACTCATCCCCGCC  >  AE015451/4160611‑4161033                                                                                                                                                                                                               |                                                                                                                                                                                                                         ccGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTa                                                                                                                                                                                                                 <  1:396045/216‑1 (MQ=255)                                ctTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTg                                                                                                                                                                                                                         >  1:439858/1‑177 (MQ=255)                                ctTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTg                                                                                                                                                                                                                         <  2:439854/177‑1 (MQ=255)                                           aCGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAaca                                                                                                                                                      >  2:168529/1‑233 (MQ=255)                                                      tGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACa                                                                                                                                                                                                               <  2:165126/165‑1 (MQ=255)                                                            tCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAg                                                                                                                                                                             >  2:507976/1‑193 (MQ=255)                                                              gAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCTg                                                                                                                                                                                               >  2:443945/1‑173 (MQ=255)                                                                         ccGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGATGAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCAATGCTCAAGGTAACCGCCTggg                                                                                                                     >  1:253111/1‑236 (MQ=255)                                                                                         gcTGTACCGATCAACCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGAt                                                                                                                                                                                         <  1:383844/151‑1 (MQ=255)                                                                                                               gTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGACCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCAATg                                                                                                                                        >  1:355072/1‑179 (MQ=255)                                                                                                                    aTGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATc                                                                                                                                                           <  1:349427/155‑1 (MQ=255)                                                                                                                    aTGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATc                                                                                                                                                           >  2:349423/1‑155 (MQ=255)                                                                                                                    aTGAAAGCTCCCCCACCTAGGTGACGCGCTCACTGCCGGAGTTCAACCACCGCCCGCCCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCAATGCTCAAGGTAACCGCCTGGGGCGTAGGcaca                                                                                                          >  2:469629/1‑204 (MQ=255)                                                                                                                        aaGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACCCGTTGGCACCa                                                                                                                                           >  1:241441/1‑167 (MQ=255)                                                                                                                        aaGCTCCCCCACCTAGGTGACTCGCTCACTGCCAGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCAATGCTCAAGGTAACCACCTgggg                                                                                                                    >  2:485383/1‑190 (MQ=255)                                                                                                                         aGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGCGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCAATGCTCAAGGTAACCGcc                                                                                                                         <  2:420557/184‑1 (MQ=255)                                                                                                                           cTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCGGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCAATGCTCAAGGTAACCGCCTGGGGCGTAGGCACAGTCTTGTCGGCGTACTGCTGGAACTGCTc                                                                             <  2:253107/226‑1 (MQ=255)                                                                                                                                 cACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCa                                                                                                                                           >  2:384550/1‑158 (MQ=255)                                                                                                                                 cACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCa                                                                                                                                           <  1:384554/158‑1 (MQ=255)                                                                                                                                                cTCACTGCCGGAGTTCAACAACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGCCTGGTCCAACGCCCCATCCCCCTGCCCGCTCAACCCCGCCAGGCTCGCCggcggc                                                                                                                                                                >  2:92258/1‑122 (MQ=255)                                                                                                                                                                           tccgtacagacgctcccgtcacccggcccccgcgttTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACgccgcc                                                                                                                                                                               <  2:205412‑M2/44‑1 (MQ=255)                                                                                                                                                                           tccgtacagacgctcccgtcacccggcccccgcgttTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACgccgcc                                                                                                                                                                               >  1:205416‑M2/37‑80 (MQ=255)                                                                                                                                                                               gAAAACCCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCCCATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCAATGCTCAAGGTAACCGCCTGGGGCGTAGGCACAGTCTTGTCGGCGTACTGCTGGAACTGCTCACGCa                                                                        <  1:65902/179‑1 (MQ=255)                                                                                                                                                                                  aaaCCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCAATGCTCAAGGTAACCGCCTGGGGCGTAGGCACAGTCTTGTCGGCGTACTGCTGGAACTGCTCACGCAGc                                                                      <  1:405478/178‑1 (MQ=255)                                                                                                                                                                                            gtcacccggcccccgcgttTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCAATGCTCAAGGTAACCGCCTGGGGCGTAGGCACAGTCTTGTCGGCGTACTGCTGGAACTGCTCACGCAGCCCACTCCCCAGTTCCATGACCCGTTCACTGCTCGCGCCACTGAGCAAGATGACGAACTCATCCCCGcc  >  1:48005‑M2/20‑236 (MQ=255)                                                                                                                                                                                               tCGCGCCCCCCACGTTGGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCAATGCTCAAGGTAACCGCCTGGGGCGTAGGCACAGTCTTGTCGGCGTACTGCTGGAACTGCTCACGCAGCCCACTCCCCAGTTCCATGACCCGt                                              <  1:95168/189‑1 (MQ=255)                                                                                                                                                                                                        ccgcgttTGGACTGGGACAACGCCACATCCCCCGGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCAATGCTCAAGGTAACCGCCTGGGGCGTAGGCACAGTCTTGTCgg                                                                                                <  2:33926‑M2/123‑1 (MQ=255)                                                                                                                                                                                                               |                                                                                                                                                                                                                         CCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCCGCTCGAAAAACCGAATGCTGTCGCGCCCCCCACGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCAATGCTCAAGGTAACCGCCTGGGGCGTAGGCACAGTCTTGTCGGCGTACTGCTGGAACTGCTCACGCAGCCCACTCCCCAGTTCCATGACCCGTTCACTGCTCGCGCCACTGAGCAAGATGACGAACTCATCCCCGCC  >  AE015451/4160611‑4161033

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