New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 333707453 (1.520)3 (0.090) 3/164 2.6 5.4% coding (336/558 nt) lptA lipopolysaccharide ABC transporter substrate‑binding protein LptA
?NZ_CP009273 3337087 = 53 (1.550)coding (349/558 nt) lptA lipopolysaccharide ABC transporter substrate‑binding protein LptA
Rejected: Frequency below/above cutoff threshold.

TAAAGAAGTGATTGACGGCTACGGTAAACCGGCAACGTTCTACCAGATGCAGGACAACGGTAAACCCGTTGAAGGTCACGCTTCCCAGA‑‑‑‑‑TG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/3336984‑3337074
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgGCAAAAGATTTTGTCGTTCTGACGGGTAATGCTTATCTGCAGCAGGTCGATAGCAACATTAAGGGCGATAAGATCACTTACCTGGTGAA  >  NZ_CP009273/3337087‑3337175
                                                                                                                                                                                         
TAAAGAAGTGATTGACGGCTACGGTAAACCGGCAACGTTCTACCAGATGCAGGACAACGGTAAACCCGTTGAAGGTCACGCTTCCCAGA‑‑‑‑‑T                                                                                            <  1:58622/90‑1
 AAAGAAGTGATTGACGGCTACGGTAAACCGGCAACGTTCTACCAGATGCAGGACAACGGTAAACCCGTTGAAGGTCACGCTTCCCAGA‑‑‑‑‑TG                                                                                           <  2:356150/90‑1
                                              ATGCAGGACAACGGTAAACCCGTTGAAGGTCACGCTTCCCAGA‑‑‑‑‑TGGCAAAAGATTTTGTCGTTCTGACGGGTAATGCTTATCTGCAGCAG                                              >  1:568434/1‑90
                                                   GGACAACGGTAAACCCGTTGAAGGTCACGCTTCCCAGATGCACTGGCAAAAGATTTTGTCGTTCTGACGGGTAATGCTTATCTGCAGCAG                                              >  2:411899/1‑90
                                                        ACGGTAAACCCGTTGAAGGTCACGCTTCCCAGA‑‑‑‑‑TGGCAAAAGATTTTGTCGTTCTGACGGGTAATGCTTATCTGCAGCAGGTCGATAGCA                                    <  2:568434/90‑1
                                                                                               GGCAAAAGATTTTGTCGTTCTGACGGGTAATGCTTATCTGCAGCAGGTCGATAGCAACATTAAGGGCGATAAGATCACTTACCTGGTGAA  <  2:573525/90‑1
                                                                                                                                                                                         
TAAAGAAGTGATTGACGGCTACGGTAAACCGGCAACGTTCTACCAGATGCAGGACAACGGTAAACCCGTTGAAGGTCACGCTTCCCAGA‑‑‑‑‑TG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/3336984‑3337074
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgGCAAAAGATTTTGTCGTTCTGACGGGTAATGCTTATCTGCAGCAGGTCGATAGCAACATTAAGGGCGATAAGATCACTTACCTGGTGAA  >  NZ_CP009273/3337087‑3337175

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