New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3567519 =47 (1.260)3 (0.080) 3/272 5.1 6.2% coding (515/1434 nt) glgA glycogen synthase
?NC_000913 3567558 = 46 (1.290)coding (476/1434 nt) glgA glycogen synthase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GGCTTGACCCATTCTGGCGTCCTGATGTGGTGCATGCGCACGACTGGCATGCAGGCCTTGCGCCTGCGTATCTGGCGGCGCGCGGGCGTCCGGCGAAGTCGGTGTTTACTGTGCACAACCTGGCCTATCAAGGCATGTTTTATGCACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3567666‑3567519
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgcacaGTAAACACCGACTTCGCCGGACGCCCGCGCGCCGCCAGATACGCAGGCGCAAGGCCTGCATGCCAGTCGTGCGCATGCACCACATCAGGACGCCAGAATGGGTC  >  NC_000913/3567558‑3567661
                                                                                                                                                                                                                                                            
GGCTTGCCCCATTCTGGCGTCCTTATGTGGTGCATGCGCACGACTGGCATGCAGGCCTTGCGCCTGCGTATCTGGCGGCGCGCGGGCGTCCGGCGAAGTCGGTGTTTACTCTGCACAACCTGGCCTATCAAGGCATGTTTTAT                                                                                                               <  2:351690/143‑1
                                                                                                            CTGTGCACAACCTGGCCTATCAAGGCATGTTTTATGCACAGTAAACACCGACTTCGCCGGACGCCCGCGCGCCGCCAGATACGCAGGC                                                          >  1:661051/1‑88
                                                                                                            CTGTGCACAACCTGGCCTATCAAGGCATGTTTTATGCACAGTAAACACCGACTTCGCCGGACGCCCGCGCGCCGCCAGATACGCAGGC                                                          <  2:661051/88‑1
                                                                                                             TGTGCACAACCTGGCCTATCAAGGCATGTTTTATGCACAGTAAACACCGACTTCGCCGGACGCCCGCGCGCCGCCAGATACGCAGGCGCAAGGCCTGCATGCCAGTCGTGCGCATGCACCACATCAGGACGCCAGGATGGGTC  >  1:324473/1‑143
                                                                                                                                                                                                                                                            
GGCTTGACCCATTCTGGCGTCCTGATGTGGTGCATGCGCACGACTGGCATGCAGGCCTTGCGCCTGCGTATCTGGCGGCGCGCGGGCGTCCGGCGAAGTCGGTGTTTACTGTGCACAACCTGGCCTATCAAGGCATGTTTTATGCACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3567666‑3567519
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgcacaGTAAACACCGACTTCGCCGGACGCCCGCGCGCCGCCAGATACGCAGGCGCAAGGCCTGCATGCCAGTCGTGCGCATGCACCACATCAGGACGCCAGAATGGGTC  >  NC_000913/3567558‑3567661

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