New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 12966098 (0.120)8 (0.190)
+CTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCA
8/174 4.8 27.0% intergenic (‑362/‑204) nhaA‑I/PP_1133 Na(+)/H(+) antiporter NhaA/transposase
?NC_002947 5545626 = 59 (0.890)intergenic (‑766/+111) rpsF/rlmB 30S ribosomal protein S6/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Rejected: Coverage evenness skew score above cutoff.

GGGACCATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGTTAGG  >  NC_002947/1296462‑1296620
                                                                                                                                                   |           
gggACCATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTg            >  7:100920/1‑149 (MQ=255)
     cATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCCGCAAACATGAGGTGTAAGGCTTAAGTTCCGCAGCGCCGCTGACATTTTGAAAAAGAAAAGGCGGTTTGCCCGGTAGGGCCAGCCGCCGTtttgtttg       >  5:268870/1‑149 (MQ=11)
     cATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCACACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGTCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTg       >  7:281357/1‑149 (MQ=255)
      aTTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGt      <  2:9325/149‑1 (MQ=37)
      aTTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGt      <  6:223085/149‑1 (MQ=37)
          ttCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGTTAgg  <  2:108088/149‑1 (MQ=33)
          ttCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGTTAgg  <  3:158527/149‑1 (MQ=33)
          ttCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGTTAgg  <  7:377831/149‑1 (MQ=33)
                                                                                                                                                   |           
GGGACCATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGTTAGG  >  NC_002947/1296462‑1296620

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