New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3512437 =17 (0.700)3 (0.130) 3/164 1.8 15.6% intergenic (‑14/+387) aroK/hofQ shikimate kinase AroK/DNA uptake porin HofQ
?NZ_CP009273 3645396 = NA (NA)pseudogene (746/747 nt) BW25113_RS18155 hypothetical protein
Rejected: Frequency below/above cutoff threshold.

GACTCTCGCAATATCTTATGAGGTTTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3512523‑3512437
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttagGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAAT  >  NZ_CP009273/3645396‑3645483
                                                                                                                                                                               
GACTCTCGCAATATCTTATGAGGTTTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAGGGA                                                                                       >  2:639336/1‑90
                   GAGGTTTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAATAAGCGGGGA                                                                    <  1:58103/90‑1
                        TTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAATAAGCGGGGAAATTC                                                               >  2:235613/1‑90
                                                                                   TTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAA    <  1:143637/90‑1
                                                                                    TAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAA   <  1:293914/90‑1
                                                                                     AGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAAT  <  1:165779/90‑1
                                                                                     AGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAAT  <  1:372448/90‑1
                                                                                     AGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAAT  <  1:437111/90‑1
                                                                                      GGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATT                                     <  1:311988/54‑1
                                                                                      GGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATT                                     >  2:311988/1‑54
                                                                                                                                                                               
GACTCTCGCAATATCTTATGAGGTTTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3512523‑3512437
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttagGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAAT  >  NZ_CP009273/3645396‑3645483

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