New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4636809141 (1.010)5 (0.040) 5/164 NT 3.5% coding (114/1425 nt) creC sensory histidine kinase in two‑component regulatory system with CreB or PhoB
?NC_000913 = 4636834 144 (1.120)coding (139/1425 nt) creC sensory histidine kinase in two‑component regulatory system with CreB or PhoB
Rejected: Frequency below cutoff threshold.

GGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4636705‑4636809
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agcaacgTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAAAATGGCCAGTACGAACCAGGCTGCCACCGCCACCA  <  NC_000913/4636834‑4636736
                                                                                                                                                                                                            
GGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAG                                                                                                          <  2:224944/100‑1
 GCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGC                                                                                                         <  2:3494694/100‑1
 GCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGC                                                                                                         >  2:2260197/1‑100
   ATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAA                                                                                                       >  2:2673300/1‑100
    gGCGGTTGTTGGTGGGCTATTTTTTACTGGTGGCGGGGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAAC                                                                                                      <  1:3531830/99‑1
    TGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTCTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAAC                                                                                                      <  1:2015378/100‑1
     GCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACG                                                                                                     >  1:2439544/1‑100
     GCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACG                                                                                                     <  1:2009119/100‑1
                                    GCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGC                                                                      <  2:1012516/100‑1
                                      GGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTC                                                                    <  1:3096629/100‑1
                                         GGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTC                                                                 <  1:2581031/100‑1
                                          GCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGTTGCGGTGCCGATCAACGTCCGCTCCGTTGCTCTTCG                                                                <  1:2472094/100‑1
                                                tcGTTCGTACTGGATAGTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGTTGCGGTGGCGATCAACGTCCCCTCCGTTGCT                                                                     <  1:26191/87‑1
                                                                                                  AGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAAAATGGCCAGTACGAACCAGGCTGCCACCG        <  2:1618662/100‑1
                                                                                                  AGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAAAATGGCCAGTACGAACCAGGCTGCCACCG        <  2:918155/100‑1
                                                                                                   GCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAAAATGGCCAGTACGAACCAGGCTGCCACCGC       >  1:651886/1‑100
                                                                                                   GCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAAAATGGCCAGTACGAACCAGGCTGCCACCGC       >  2:1749419/1‑100
                                                                                                     AACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAAAATGGCCAGTACGAACCAGGCTGCCACCGCCA     <  2:686745/100‑1
                                                                                                       CGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAAAATGGCCAGTACGAACCAGGCTGCCACCGCCACC   <  2:1951182/100‑1
                                                                                                       CGTTGCGGTGCCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAACATCGCCAGTACGAACCAGGCTGCCACCGCCACC   <  1:2500415/100‑1
                                                                                                        GTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAAAATGGCCAGTACGAACCAGGCTGCCACCGCCACCA  <  2:2622483/100‑1
                                                                                                        GTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAAAATGGCCAGTACGAACCAGGCTGCCACAGCCACCA  >  1:2830272/1‑100
                                                                                                                                                                                                            
GGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4636705‑4636809
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agcaacgTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAAAATGGCCAGTACGAACCAGGCTGCCACCGCCACCA  <  NC_000913/4636834‑4636736

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