New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3653230NA (NA)6 (0.200) 6/260 2.4 10.2% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 = 3906171 53 (1.800)coding (397/837 nt) bglG transcriptional antiterminator of the bgl operon
Rejected: Frequency below/above cutoff threshold.

AATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3653128‑3653230
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTAGGGGAAGAAGCGTTAACCATTATTGATAAACGGTTGGGCGTGCAGTTACCGAAAGATGAAGTGGGCTTTATTGCCATGCATCTGGTCAGTGCCCAAATGAGCGGAAATATGGAGGATGTTGCAGGTGTC  <  NC_000913/3906171‑3906040
                                                                                                                                                                                                                                           
AATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGGGGAAGAAGCGTTAACCATT                                                                                                              <  3:106115/127‑1
AATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGGGGAAGAAGCGTTAACCATT                                                                                                              >  4:106115/1‑127
         TCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGGGGAAGAAGCGTTAACCATTATTGATAAACGGTTGGGCGTGCAGTTACCG                                                                                >  6:70447/1‑148
                     AGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGGGGAAGAAGCGTTAACCATTATTGATAAACGGTTGGGCGTGCAGTTACCGAAAGATGAAGTG                                                                    <  5:70447/148‑1
                                               TCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGGGGAAGAAGCGTTAACCATTATTGATAAACGGTTGGGCGTGCAGTTACCGAAAGATGAAGTGGGCTTTATTGCCATGCATCTGGTCAGT                                         >  1:5606/1‑149
                                                                                      ACTTGTTCGCACCTTCCCTAGGGGAAGAAGCGTTAACCATTATTGATAAACGGTTGGGCGTGCAGTTACCGAAAGATGAAGTGGGCTTTATTGCCATGCATCTGGTCAGTGCCCAAATGAGCGGAAATATGGAGGATGTTGCAGGTGTC  <  2:5606/149‑1
                                                                                                                                                                                                                                           
AATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3653128‑3653230
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTAGGGGAAGAAGCGTTAACCATTATTGATAAACGGTTGGGCGTGCAGTTACCGAAAGATGAAGTGGGCTTTATTGCCATGCATCTGGTCAGTGCCCAAATGAGCGGAAATATGGAGGATGTTGCAGGTGTC  <  NC_000913/3906171‑3906040

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