New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4638146108 (0.880)9 (0.080) 8/156 NT 7.7% intergenic (+26/‑32) creC/creD sensory histidine kinase in two‑component regulatory system with CreB or PhoB/inner membrane protein
?NC_000913 = 4638170 122 (1.130)intergenic (+50/‑8) creC/creD sensory histidine kinase in two‑component regulatory system with CreB or PhoB/inner membrane protein
Rejected: Frequency below cutoff threshold.

TGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4638038‑4638146
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acatagtcttcTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCACCTTCC  <  NC_000913/4638170‑4638072
                                                                                                                                                                                                                
TGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCAC                                                                                                              >  2:1875816/1‑100
 GTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACA                                                                                                             >  1:1080953/1‑100
 GTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACA                                                                                                             >  1:1081070/1‑100
 GTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACA                                                                                                             >  2:3257130/1‑100
  TTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACAT                                                                                                            <  1:1531947/100‑1
  TTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACAT                                                                                                            >  1:1750695/1‑100
  TTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACAT                                                                                                            <  1:2326135/100‑1
   TTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATA                                                                                                           >  1:477638/1‑100
   TTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATA                                                                                                           >  2:1466256/1‑100
   TTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATA                                                                                                           >  1:2348853/1‑100
   TTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATA                                                                                                           >  1:142919/1‑100
   TTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATA                                                                                                           >  1:139847/1‑100
     AACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGT                                                                                                         <  1:286506/100‑1
     AACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGT                                                                                                         <  2:24853/100‑1
      ACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTC                                                                                                        >  1:99819/1‑100
       CGGTGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGGCACTTCACATAGCTTCAAATTCTTCCCACATAGT                                                                                                         >  1:2899669/1‑98
       CGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCT                                                                                                       >  1:2806963/1‑100
        GGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTT                                                                                                      <  1:1490198/100‑1
        GGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTT                                                                                                      >  2:2189696/1‑100
        GGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTT                                                                                                      >  1:849661/1‑100
        GGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTT                                                                                                      >  1:3281176/1‑100
         GCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTC                                                                                                     <  1:1060030/100‑1
                                           CGTGCTGGCCTCGCTTCGCCTTCAGCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATG                                                                   <  2:3002544/100‑1
                                            GTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGT                                                                  <  1:762107/100‑1
                                             TGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTG                                                                 <  2:1910235/100‑1
                                             TGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATGCTTCCACATAGTCTTCTCCTGTGCAATGGCAGCAGGATACGAAGACTATGTGG                                                                <  1:2991003/100‑1
                                              GCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGG                                                                <  2:616385/100‑1
                                               CTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGG                                                               <  2:3414700/100‑1
                                                TGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGA                                                              <  1:2760539/100‑1
                                                 GGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAA                                                             <  2:2205308/100‑1
                                                      CGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATT                                                        <  1:1466250/100‑1
                                                                                                  ACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCAC            <  1:2238920/100‑1
                                                                                                      AGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCA        <  1:2795642/100‑1
                                                                                                       GTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCAC       <  2:2691641/100‑1
                                                                                                         CTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCACCT     <  2:1343903/100‑1
                                                                                                          TTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCACCTT    >  2:2939995/1‑100
                                                                                                          TTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCACCTT    >  2:3177402/1‑100
                                                                                                            CTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCACCTTCC  <  2:3252542/100‑1
                                                                                                                                                                                                                
TGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4638038‑4638146
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acatagtcttcTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCACCTTCC  <  NC_000913/4638170‑4638072

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