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

TCCGTACAGACGCTCCCGTCACCCGGCCCCCGCGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/4160708‑4160673
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgttTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCAATGCTCAAGGTAACCGCCTGGGGCGTAGGCACAGTCTTGTCGGCGTACTGCTGGAACTGCTCACGCAGCCCACTCCCCAGTTCCATGACCCGTTCACTGCTCGCGCCACTGAGCAAGATGACGAACTCATCCCCGCC  >  AE015451/4160817‑4161033
                                                                                                                                                                                                                                                             
TCCGTACAGACGCTCCCGTCACCCGGCCCCCGCGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCC                                                                                                                                                                               >  1:205416/1‑80
TCCGTACAGACGCTCCCGTCACCCGGCCCCCGCGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCC                                                                                                                                                                               <  2:205412/80‑1
                 GTCACCCGGCCCCCGCGTTTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCAATGCTCAAGGTAACCGCCTGGGGCGTAGGCACAGTCTTGTCGGCGTACTGCTGGAACTGCTCACGCAGCCCACTCCCCAGTTCCATGACCCGTTCACTGCTCGCGCCACTGAGCAAGATGACGAACTCATCCCCGCC  >  1:48005/1‑236
                             CCGCGTTTGGACTGGGACAACGCCACATCCCCCGGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCAATGCTCAAGGTAACCGCCTGGGGCGTAGGCACAGTCTTGTCGG                                                                                                <  2:33926/130‑1
                                                                                                                                                                                                                                                             
TCCGTACAGACGCTCCCGTCACCCGGCCCCCGCGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/4160708‑4160673
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgttTGGACTGGTACAACGCCACATCCCCCTGCTCGATCAACGCCGCCAGGCTCGCCGGCGGCTGATCGAACACGTTGGCACCAATGCTCAAGGTAACCGCCTGGGGCGTAGGCACAGTCTTGTCGGCGTACTGCTGGAACTGCTCACGCAGCCCACTCCCCAGTTCCATGACCCGTTCACTGCTCGCGCCACTGAGCAAGATGACGAACTCATCCCCGCC  >  AE015451/4160817‑4161033

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