New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb 652499 =74 (1.520)3 (0.070) 3/92 NT 4.7% coding (1042/1659 nt) citA sensory histidine kinase in two‑component regulatory system with citB
?W3110S.gb 652526 = 54 (1.210)coding (1069/1659 nt) citA sensory histidine kinase in two‑component regulatory system with citB
Rejected: Frequency below/above cutoff threshold.

CCATCGCCAGCACGCGATCATACTCTTTCATCTGCAACAGACCATTGAGCGTCGACATCCAATTGAGATGCTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  W3110S.gb/652571‑652499
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gctcAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCAGGGGGAGTCTCAGGCCCAG  >  W3110S.gb/652526‑652594
                                                                                                                                              
CCATCGCCAGCACGCGATCATACTCTTTCATCTGCAACAGACCATTGAGCGTCGACATCCAATTGAGATGC                                                                         >  1:448575/1‑71
  ATCGCCAGCACGCGATCATACTCTTTCATCTGCAACAGACCATTGAGCGTCGACATCCAATTGAGATGCTC                                                                       >  1:2897562/1‑71
              CGATCATACTCTTTCATCTGCAACAGACCATTGAGCGTCGACATCCAATTGAGATGCT                                                                        >  1:2064573/1‑58
                                             TGAGCGTCGACATCCAATTGAGATGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGA                            >  1:4299637/1‑71
                                                   TCGACATCCAATTGAGATGCTCAATGGTCTGTTGC                                                          <  1:4320866/35‑1
                                                     GACATCCAATTGAGATGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGC                                  >  1:4398618/1‑57
                                                                      CTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCAGGGGGAGTCTCAGGCCCA   <  1:1963224/71‑1
                                                                      CTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGAT                           <  1:482827/47‑1
                                                                        CAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCAGGGGGAGTCTCAGGCCCAG  >  1:3751920/1‑70
                                                                        CAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCAGGGGGAGTCTCAG        >  1:17600/1‑64
                                                                        CAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCAG                    >  1:3615802/1‑52
                                                                        CAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATG                          >  1:4429400/1‑46
                                                                        CAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGA                            <  1:3401939/44‑1
                                                                                                                                              
CCATCGCCAGCACGCGATCATACTCTTTCATCTGCAACAGACCATTGAGCGTCGACATCCAATTGAGATGCTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  W3110S.gb/652571‑652499
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gctcAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCAGGGGGAGTCTCAGGCCCAG  >  W3110S.gb/652526‑652594

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