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

CGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4636704‑4636809
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agcaacgTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAAAATGGCCAGTACGAACCAGGCTGCCACCGCCACCA  <  NC_000913/4636834‑4636736
                                                                                                                                                                                                             
CGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGA                                                                                                           <  2:2423547/100‑1
CGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGA                                                                                                           >  1:2535582/1‑100
     TGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAAC                                                                                                      >  2:1588148/1‑100
      GCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGGAGTTAAACCGGGCGTGCGAAGAGCAACG                                                                                                     >  1:145272/1‑100
      GCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACG                                                                                                     >  2:1657248/1‑100
      GCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGGGCTGGGCATTTTTGTCGAAGAAGTTAAACCGGGCGTGCGAAGAGCAACG                                                                                                     >  1:253389/1‑100
                                      CGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCT                                                                     <  1:1214213/100‑1
                                        GTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCT                                                                   <  2:371718/100‑1
                                          GGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTC                                                                 <  2:2834837/100‑1
                                           GCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCG                                                                <  2:241878/100‑1
                                            CAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGC                                                               <  1:1179275/100‑1
                                                                                                   AGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAAAATGGTTTGTACGAACCAGGCGGCCACCG        <  1:838525/100‑1
                                                                                                   AGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAAAATGGCCAGTACGAACCAGGCTGCCACCG        <  2:2734063/100‑1
                                                                                                   AGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAAAATGGCCAGTACGAACCAGGCTGCCACCG        <  2:2734640/100‑1
                                                                                                   AGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGCTTAACTTCTTCGACAAAAATGGCCAGTACGAACCAGGCTGCCACC         >  1:68579/1‑99
                                                                                                        CGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAAAATGGCCAGTACGAACCAGGCTGCCACCGCCACC   <  2:692724/100‑1
                                                                                                         GTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAAAATGGCCAGTACGAACCAGGCTGCCACCGCCACCA  <  2:2244661/100‑1
                                                                                                                                                                                                             
CGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4636704‑4636809
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agcaacgTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTGACAAAAATGGCCAGTACGAACCAGGCTGCCACCGCCACCA  <  NC_000913/4636834‑4636736

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.