New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 65278252 (0.780)16 (0.280) 15/240 4.2 41.2% coding (1961/1962 nt) ctpL methyl‑accepting chemotaxis protein CtpL
?NC_002947 = 4808159 1 (0.020)intergenic (‑268/‑287) PP_4234/dsbD‑II bifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2
Rejected: Coverage evenness skew score above cutoff.

ACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCG  >  NC_002947/4808011‑4808198
                                                                                                                                                    |                                       
aCAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTgg                                         >  6:137442/1‑149 (MQ=255)
                tgtATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTACGCCGGGAGAGTTGGTGCGCCGGGGAGATCGAGcgccgca                         >  6:111025/1‑148 (MQ=38)
                      gggCATTGGGCGGATTTGGCGCGTTGGGGCGTAGGTGTTGGGCGGCTGGTTTGCTGCGCCTGGTAGGTGTTCGCCGGGTTTTTTTCATCGGCTGGGAGGTGTTCGGCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGccgcgcgggcg                   <  8:103105/149‑5 (MQ=255)
                       ggCAATGGGCGGATATCGCGCCTTGGCGGGTCTGTGTTTGGCGTCACGTTTTCAGTGCCTGGTAGGTTTTGGGCGGTTTTTTTTCAGCGGCTGGGAGGTGTTCGCCGGGAGTGTTTCAGCGCCTGGGAGATCGAgc                               <  6:133384/136‑1 (MQ=16)
                       ggCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCCGCGCCTCGGAGGTCTTCGCCGGGAGAGTTTCAGCGCCTGGGATATCGAgc                               >  5:133384/1‑136 (MQ=255)
                                    aTTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGCGATTTTCAGCGCCTGGGAGGTGTACGCCGGGACAGTTTCAGCGGCTCGGAGATCGAGCGCCGCACGCGCGGCGCttcgcgtgt      >  4:132149/1‑145 (MQ=25)
                                      tCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGA‑TTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCGTGGGAGTTCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaaa  >  5:51623/1‑144 (MQ=37)
                                       cgcgcgGTGGCGGGTAGGTGTTTGGCTGCAGGTTTTCTTTGGCTTGTAGGTTTTGGGCGGTTTTTTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGATTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaa   <  2:322287/148‑5 (MQ=12)
                                                                                                                                                    |                                       
ACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCG  >  NC_002947/4808011‑4808198

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