New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4638146117 (0.910)6 (0.050) 6/156 NT 5.2% intergenic (+26/‑32) creC/creD sensory histidine kinase in two‑component regulatory system with CreB or PhoB/inner membrane protein
?NC_000913 = 4638170 116 (1.040)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                                                                                                              <  1:3065505/100‑1
  TTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCTTCACTTCACATAGCTTCAAATTCTTCCCACAT                                                                                                            <  2:3282351/100‑1
  TTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACAT                                                                                                            <  2:1512225/100‑1
  TTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACAT                                                                                                            >  2:3207520/1‑100
   TTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATA                                                                                                           >  1:322938/1‑100
    TAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAG                                                                                                          <  2:1237680/100‑1
    TAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAG                                                                                                          >  1:1379572/1‑100
    TAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAG                                                                                                          >  2:1675819/1‑100
    TAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAG                                                                                                          >  2:165122/1‑100
     AACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGT                                                                                                         >  1:2629490/1‑100
       CGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCT                                                                                                       >  1:2252266/1‑100
       CGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCT                                                                                                       >  2:1854550/1‑100
        GGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTT                                                                                                      >  1:2277548/1‑100
                                         GGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTA                                                                     <  1:592571/100‑1
                                          GCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTAT                                                                    <  2:1229729/100‑1
                                              GCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGG                                                                <  1:1093731/100‑1
                                                       GCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTT                                                       <  1:1230531/100‑1
                                                        CTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTG                                                      <  1:739958/100‑1
                                                                                  GCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTC                            <  2:803132/100‑1
                                                                                                  ACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCAC            <  1:721049/100‑1
                                                                                                  ACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCAC            >  2:216250/1‑100
                                                                                                    ATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGC          <  1:3310708/100‑1
                                                                                                     TAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCC         >  1:2153585/1‑100
                                                                                                     TAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCC         <  2:682462/100‑1
                                                                                                      AGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCA        <  1:2730119/100‑1
                                                                                                       GTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATTTGAAGTGACGGTGAAGTCAAAGCGAGGCCAGCACACCAC       >  2:2760121/1‑100
                                                                                                       GTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCAC       >  1:2555779/1‑100
                                                                                                       GTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCAC       >  2:2064802/1‑100
                                                                                                       GTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCAC       >  2:527116/1‑100
                                                                                                       GTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCAC       >  1:1809627/1‑100
                                                                                                         CTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCACCT     >  1:353384/1‑100
                                                                                                          TTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCACCTT    >  2:3290446/1‑100
                                                                                                           TCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCACCTTC   >  2:2637040/1‑100
                                                                                                            CTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCACCTTCC  <  1:2839322/100‑1
                                                                                                            CTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCACCTTCC  <  1:2007856/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.