New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4639713169 (1.080)6 (0.040) 4/172 NT 3.4% coding (594/717 nt) arcA response regulator in two‑component regulatory system with ArcB or CpxA
?NC_000913 = 4639734 176 (1.160)coding (573/717 nt) arcA response regulator in two‑component regulatory system with ArcB or CpxA
Rejected: Frequency below cutoff threshold.

AGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4639613‑4639713
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggCCGTGAGCTGAAACCGCACGACCGTACTGTAGACGTGACGATCCGCCGTATTCGTAAACATTTCGAATCTACGCCGGATACGCCGGAAATCATCGCCAC  <  NC_000913/4639734‑4639636
                                                                                                                                                                                                        
AGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTAC                                                                                                       >  1:872751/1‑99
                        TGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGCCGTGAGCTGAAACCGCACGACC                                                                              <  1:2206350/100‑1
                          GCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGCCGTGAGCTGAAACCGCACGACCGT                                                                            <  1:2119739/100‑1
                          GCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGCCGTGAGCTGAAACCGCACGACCGT                                                                            <  1:2119798/100‑1
                            GATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGCCGTGAGCTGAAACCGCACGACCGTAC                                                                          <  2:3245718/100‑1
                            GATGATTTCCGGCGTATCCGGCGTAGATTCGAAAGGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGCCGTGAGCTGAAACCGCACGACCGTAC                                                                          <  1:2069215/100‑1
                             ATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTATGAATGCGGCGGATCGTCACGTCTTCAGTACGGCCGTGAGCTGAAACCGCACGACCGTACT                                                                         <  1:1412276/100‑1
                                                                                                  CGGCCGTGAGCTGAAACCGCACGACCGTACTGTAGACGTGACGATCCGCCGTATTCGTAAACATTTCGAATCTACGCCGGATACGCCGGAAATCATCGCC    <  1:434351/100‑1
                                                                                                  CGGCCGTGAGCTGAAACCGCACGACCGTACTGTAGACGTGACGATCCGCCGTATTCGTAAACATTTCGAATCTACGCCGGATACGCCGGAAATCATCGCC    <  2:162962/100‑1
                                                                                                  CGGCCGTGAGCTGAAACCGCACGACCGTACTGTAGACGTGACGATCCGCCGTATTCGTAAACATTTCGAATCTACGCCGGATACGCCGGAAATCATCGCC    <  2:2120083/100‑1
                                                                                                  CGGCCGTGAGCTGAAACCGCACGACCGTACTGTAGACGTGACGATCCGCCGTATTCGTAAACATTTCGAATCTACGCCGGATACGCCGGAAATCATCGCC    <  2:2814803/100‑1
                                                                                                  CGGCCGTGAGCTGAAACCGCACGACCGTACTGTAGACGTGACGATCCGCCGTATTCGTAAACATTTCGAATCTACGCCGGATACGCCGGAAATCATCGCC    >  2:2937770/1‑100
                                                                                                  CGGCCGTGAGCTGAAACCGCACGACCGTACTGTAGACGTGACGATCCGCCGTATTCGTAAACATTTCGAATCTACGCCGGATACGCCGGAAATCATCG      <  1:3310835/98‑1
                                                                                                   GGCCGTGAGCTGAAACCGCACGACCGTACTGTAGACGTGACGATCCGCCGTATTCGTACACATTTCGAATCTACGCCGGATACGCCGGAAATCATCGCC    >  1:4310187/1‑99
                                                                                                   GGCCGTGAGCTGAAACCGCACGACCGTACTGTAGACGTGACGATCCGCCGTATTCGTAAACATTTCGAATCTACGCCGGATACGCCGGAAATCATCGCCA   >  2:2037755/1‑100
                                                                                                    GCCGTGAGCTGAAACCGCACGACCGTACTGTAGACGTGACGATCCGCCGTATTCGTAAACATTTCGAATCTACGCCGGATACGCCGGAAATCATCGCCAC  >  2:3442270/1‑100
                                                                                                    GCCGTGAGCTGAAACCGCACGACCGTACTGTAGACGTGACGATCCGCCGTATTCGTAAACATTTCGAATCTACGCCGGATACGCCGGAAATCATCGCCAC  >  2:4294990/1‑100
                                                                                                                                                                                                        
AGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4639613‑4639713
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggCCGTGAGCTGAAACCGCACGACCGTACTGTAGACGTGACGATCCGCCGTATTCGTAAACATTTCGAATCTACGCCGGATACGCCGGAAATCATCGCCAC  <  NC_000913/4639734‑4639636

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