New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3970088 =33 (0.960)4 (0.120) 4/268 3.8 10.5% intergenic (+11/‑45) wzzE/wecB Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein/UDP‑N‑acetyl glucosamine‑2‑epimerase
?NC_000913 3970114 = 36 (1.080)intergenic (+37/‑19) wzzE/wecB Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein/UDP‑N‑acetyl glucosamine‑2‑epimerase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTGACGCAAACTTTAGCCTCAAAAAAAGGATCTTTTGCCAACGCATGCACCAACGGCGCCATCTTGATGGCTTCCGGGCGCGTACCAAATACAGTCAGTACTTTCACATCGATTCTCTTCGAATAAGCGGCGAGCGCCTTTGCGCTCACCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3970239‑3970088
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccgcTTATTCGAAGAGAATCGATGTGAAAG  >  NC_000913/3970114‑3970139
                                                                                                                                                                                  
GTGACGCAAACTTTAGCCTCAAAAAAAGGATCTTTTGCCAACGCATGCACCAACGGCGCCATCTTGATGGCTTCCGGGCGCGTACCAAATACAGTCAGTACTTTCACATCGATTCTCTTCGAATAAGCGGCGAGCGCCTTTGCGCTCAC                               <  3:65339/149‑1
GTGACGCAAACTTTAGCCTCAAAAAAAGGATCTTTTGCCAACGCATGCACCAACGGCGCCATCTTGATGGCTTCCGGGCGCGTACCAAATACAGTCAGTACTTTCACATCGATTCTCTTCGAATAAGCGGCGAGCGCCTTTGCGCTCAC                               >  4:103625/1‑149
                                          GCATGCACCAACGGCGCCATCTTGATGGCTTCCGGGCGCGTACCAAATACAGTCAGTACTTTCACATCGATTCTCTTCGAATAAGCGGCGAGCGCCTTTGCGCTCACCGCTTAT                        >  7:111188/1‑114
                                          GCATGCACCAACGGCGCCATCTTGATGGCTTCCGGGCGCGTACCAAATACAGTCAGTACTTTCACATCGATTCTCTTCGAATAAGCGGCGAGCGCCTTTGCGCTCACCGCTTAT                        <  8:111188/114‑1
                                                                                                                     TTCGAATAAGCGGCGAGCGCCTTTGCGCTCACCGCTTATTCGAAGAGAATCGATGTGAAAG  >  7:58977/1‑61
                                                                                                                     TTCGAATAAGCGGCGAGCGCCTTTGCGCTCACCGCTTATTCGAAGAGAATCGATGTGAAAG  <  8:58977/61‑1
                                                                                                                                                                                  
GTGACGCAAACTTTAGCCTCAAAAAAAGGATCTTTTGCCAACGCATGCACCAACGGCGCCATCTTGATGGCTTCCGGGCGCGTACCAAATACAGTCAGTACTTTCACATCGATTCTCTTCGAATAAGCGGCGAGCGCCTTTGCGCTCACCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3970239‑3970088
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccgcTTATTCGAAGAGAATCGATGTGAAAG  >  NC_000913/3970114‑3970139

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