New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3805341 =154 (4.070)6 (0.160) 6/250 2.8 3.9% coding (610/798 nt) waaP kinase that phosphorylates core heptose of lipopolysaccharide
?NC_000913 3805362 = 144 (3.870)coding (589/798 nt) waaP kinase that phosphorylates core heptose of lipopolysaccharide
Rejected: Frequency below/above cutoff threshold.

ATACTGCAGGGATAAATCACCGCGACTGTTACATTTGCCACTTTTTGCTTCATTTGCCATTTACTGGCCGGGAAGATGAATTAAAAATTTCAGTTATCGATCTGCATCGGGCACAGATACGTGCAAAAGTACCGCGCCGCTGGCGCGATA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3805490‑3805341
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑taCTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCA  >  NC_000913/3805362‑3805493
                                                                                                                                                                                                                                                                                          
ATACTGCAGGGATAAATCACCGCGACTGTTACATTTGCCACTTTTTGCTTCATTTGCCATTTACTGGCCGGGAAGATGAATTAAAAATTTCAGTTATCGATCTGCATCGGGCACAGATACGTGCAAAAGTACCGCGCCGCTGGCGCGAT                                                                                                                                       <  6:32939/149‑1
                                                                                                                        GTGCAAAAGTACCGCGCCGCTGGCGCGATACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAA                                                                                     >  1:57665/1‑79
                                                                                                                        GTGCAAAAGTACCGCGCCGCTGGCGCGATACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAA                                                                                     <  2:57665/79‑1
                                                                                                                                GTACCGCGCCGCTGGCGCGATACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAA                                                                                            <  1:12192/64‑1
                                                                                                                                GTACCGCGCCGCTGGCGCGATACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAA                                                                                            >  2:12192/1‑64
                                                                                                                                      CGCCGCTGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCA  >  6:78414/1‑148
                                                                                                                                      CGCCGCTGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATACATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCA  <  5:78414/148‑1
                                                                                                                                                     ACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTA      >  3:1001/1‑129
                                                                                                                                                     ACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTA      <  4:1001/129‑1
                                                                                                                                                                                                                                                                                          
ATACTGCAGGGATAAATCACCGCGACTGTTACATTTGCCACTTTTTGCTTCATTTGCCATTTACTGGCCGGGAAGATGAATTAAAAATTTCAGTTATCGATCTGCATCGGGCACAGATACGTGCAAAAGTACCGCGCCGCTGGCGCGATA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3805490‑3805341
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑taCTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCA  >  NC_000913/3805362‑3805493

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.