New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AE015451 = 456525835 (0.880)4 (0.110) 4/320 NT 11.1% coding (501/1560 nt) glgA glycogen synthase
?AE015451 = 4565305 31 (0.810)coding (548/1560 nt) glgA glycogen synthase

TGCCCAGAACTGTACCAGCGCGACGGCACACCCTATGGTGCCAACAATGGCCGCGACTGGCCCGACAACCACATCCGCTTCGCCCGCCTGGGCCTGGCTGCCGCCGAAATCGCGGCCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/4565139‑4565258
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggccggcCAGTCATGCGCGTGCACCACTTCCGGCGTCCAGTGGATCTTGCCTTCGCCGGCCGCGATTTCGGCGGCAGCCAGGCCCAGGCGGGCGAAGCGGATGTGGTTGTCGGGCCAGTCGCGGCCATTGTTGGCACCATAGGGTGTGCCGTCGCGCTGGTACAGTTCTGGGCAGATCAACACATAGATGACCAGGCCGTCGGCCATGTCC  <  AE015451/4565305‑4565102
                                                                                                                                                                                                                                                                                                                                    
TGCCCAGAACTGTACCAGCGCGACGGCAGACCCTATGGTGCCAACAATGGCCGCGACTGGCCCGACAACCACCTCCGCTTCGCCCGCCTGGGCCTGGCTGCCGCCGAAATCGCGGCC                                                                                                                                                                                                                 >  2:540905/1‑117
                                                                                 GCCCGCCTGGGCCTGGCTGCCGCCGAAATCGCGGCCGGCCAGTCATGCGCGTGCACCACTTCCGGCGTCCAGTGGATCTTGCCTTCGCCGGCCGCGATTTCGGCGGCAGCCAGGCCCAGGCGGGCGAAGCGGATGTGGTTGTCGGGCCAGTCGCGGCCATTGTTGGCACCATAGGGTGTGCCGTCGCGCTGGTACAGTTCTGGGCAGATCAACACATAGATGACCctgtctctta          >  1:547025/1‑225
                                                                                  CCCGCCTGGGCCTGGCTGCCGCCGAAATCGCGGCCGGCCAGTCATGCGCGTGCACCACTTCCGGCGTCCAGTGGATCTTGCCTTCGCC                                                                                                                                                            <  1:394282/88‑1
                                                                                  CCCGCCTGGGCCTGGCTGCCGCCGAAATCGCGGCCGGCCAGTCATGCGCGTGCACCACTTCCGGCGTCCAGTGGATCTTGCCTTCGCC                                                                                                                                                            >  2:394281/1‑88
                                                                                         GGGCCTGGCTGCCGCCGAAATCGCGGCCGGCCAGTCATGCGCGTGCACCACTTCCGGCGTCCAGTGGATCTTGCCTTCGCCGGCCGCGATTTCGGCGGCAGCCAGGCCCAGGCGGGC                                                                                                                        >  2:581337/1‑117
                                                                                                              CGCGGCCGGCCAGTCATGCGCGTGCACCACTTCCGGGGTCCGGTGGATCTTGCCTTCGCCGGCCGCGATTTCGGCGGCAGCCAGGCCCAGGCGGGCGAAGCGGATGTGGTTGTCGGGCCAGTCGCGGCCATTGTTGGCACCATAGGGTGTGCCGTCGCGCTGGTACAGTTCTGGGCAGATCA                                  <  2:547024/182‑1
                                                                                                                     GGCCAGTCATGCGCGTGCACCACTTCCGGCGTCCAGTGGATCTTGCCTTCGCCGGCCGCGATTTCGGCGGCAGCCAGGCCCAGGCGGGCGAAGCGGATGTGGTTGTCGGGCCAGTCGCGGCCATTGTTGGCACCATAGGGTGTGCCGTCGCGCTGGTACAGTTCTGGGCAGATCAACACATAGATGACCAGGCCGTCGGCCATGTCC  >  1:425465/1‑207
                                                                                                                     GGCCAGTCATGCGCGTGCACCACTTCCGGCGTCCAGTGGATCTTGCCTTCGCCGGCCGCGATTTCGGCGGCAGCCAGGCCCACGCGGGCGAAGCGGATGTGGTT                                                                                                         >  2:20807/1‑104
                                                                                                                     GGCCAGTCATGCGCGTGCACCACTTCCGGCGTCCAGTGGATCTTGCCTTCGCCGGCCGCGATTTCG                                                                                                                                               >  1:161759/1‑66
                                                                                                                     GGCCAGTCATGCGCGTGCACCACTTCCGGCGTCCAGTGGATCTTGCCTTCGCCGGCCGCGATTTCG                                                                                                                                               <  2:161759/66‑1
                                                                                                                                                                                                                                                                                                                                    
TGCCCAGAACTGTACCAGCGCGACGGCACACCCTATGGTGCCAACAATGGCCGCGACTGGCCCGACAACCACATCCGCTTCGCCCGCCTGGGCCTGGCTGCCGCCGAAATCGCGGCCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/4565139‑4565258
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggccggcCAGTCATGCGCGTGCACCACTTCCGGCGTCCAGTGGATCTTGCCTTCGCCGGCCGCGATTTCGGCGGCAGCCAGGCCCAGGCGGGCGAAGCGGATGTGGTTGTCGGGCCAGTCGCGGCCATTGTTGGCACCATAGGGTGTGCCGTCGCGCTGGTACAGTTCTGGGCAGATCAACACATAGATGACCAGGCCGTCGGCCATGTCC  <  AE015451/4565305‑4565102

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