New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 5222887 =NA (NA)7 (0.160)
+TGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG
6/184 5.9 72.7% intergenic (‑298/+81) PP_4602/PP_4603 AraC family transcriptional regulator/transposase
?NC_002947 = 5453138 4 (0.060)intergenic (+151/‑204) lnt/PP_4791 apolipoprotein N‑acyltransferase/transposase
Rejected: Coverage evenness skew score above cutoff.

GAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5222930‑5222887
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTT  <  NC_002947/5453138‑5453097
                                            ||||||||||||||||||||||||||||||||||||||||||||||||                                          
GAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTT  >  1:68345/1‑134
GAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTT  <  2:68345/134‑1
GAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGCTGGCCCTACCGGGCCAA                 <  7:296512/119‑1
GAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGCTGGCCCTACCGGGCCAA                 >  8:296512/1‑119
  GATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTT     <  1:195237/129‑1
              CTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTT    <  3:177080/118‑1
              CTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTT    >  4:177080/1‑118
                                                                          GCTGACATTTTGAAAAAGAAAAGGCGGCTGGCCCTACCGGGCCAACC               >  3:332277/1‑47
                                                                          GCTGACATTTTGAAAAAGAAAAGGCGGCTGGCCCTACCGGGCCAACC               <  4:332277/47‑1
                                            ||||||||||||||||||||||||||||||||||||||||||||||||                                          
GAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5222930‑5222887
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTT  <  NC_002947/5453138‑5453097

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