New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4557049 =170 (1.140)14 (0.100) 7/162 NT 8.5% coding (57/393 nt) iraD RpoS stabilzer after DNA damage, anti‑RssB factor
?NC_000913 4557077 = 147 (1.090)coding (85/393 nt) iraD RpoS stabilzer after DNA damage, anti‑RssB factor
Rejected: Frequency below cutoff threshold.

ATGCACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAGCGAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4557154‑4557049
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcagcgatCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATC  >  NC_000913/4557077‑4557173
                                                                                                                                                                                                           
ATGCACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGC                                                                                                         <  1:1955551/100‑1
ATGCACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGC                                                                                                         <  2:51712/100‑1
 TGCACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCA                                                                                                        >  2:939674/1‑100
 TGCACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCA                                                                                                        >  2:939468/1‑100
 TGCACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCA                                                                                                        >  1:3242404/1‑100
  GCACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAG                                                                                                       <  1:2240868/100‑1
   CACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTACTGCCGCATGAGGCGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAGC                                                                                                      >  2:243069/1‑100
   CACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAGC                                                                                                      >  2:1360520/1‑100
   CACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAGC                                                                                                      >  1:477540/1‑100
   CACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAGC                                                                                                      <  1:3340886/100‑1
    ACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAGCG                                                                                                     >  1:3092475/1‑100
    ACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCA                                                                                                        >  2:2420944/1‑97
     CTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAGCGA                                                                                                    >  1:2383479/1‑100
     CTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAGC                                                                                                      <  1:852192/98‑1
      TTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAGCGAT                                                                                                   <  2:1690463/100‑1
          CACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAGCGA                                                                                                    >  2:592116/1‑95
           ACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAGCGATCTCCT                                                                                              >  1:976004/1‑100
                             CAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCC                                                                            >  2:3932982/1‑100
                                                                          GGAGATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTG                               >  2:3313650/1‑100
                                                                          GGAGATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTG                               >  1:1633583/1‑100
                                                                          GGAGATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTCTTTCAGGCATGCCCcactg                               >  1:4179411/1‑95
                                                                           GAGATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCCGGCATGCCC                                    >  2:255369/1‑94
                                                                            AGATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGC                             >  1:3711455/1‑100
                                                                            AGATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGC                             >  1:1525105/1‑100
                                                                             GATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCA                            >  1:1261221/1‑100
                                                                             GATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCA                            >  1:793581/1‑100
                                                                             GATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGC                             >  2:4165837/1‑99
                                                                              ATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAA                           >  1:648533/1‑100
                                                                               TCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAG                          >  2:2199124/1‑100
                                                                               TCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAG                          >  2:1108447/1‑100
                                                                               TCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTAATTGTTTCAGGCATTCCCGATTCGCAAG                          >  2:4008478/1‑100
                                                                                                 tgtAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACA        <  1:923316/97‑1
                                                                                                  GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCACGTGCATAACCCGAGCGACAA       >  2:1485507/1‑100
                                                                                                  GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAA       >  2:474374/1‑100
                                                                                                   CAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAA      >  2:2466719/1‑100
                                                                                                    AGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAAC     <  2:1808633/100‑1
                                                                                                    AGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAAC     <  1:4056586/100‑1
                                                                                                    AGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAAC     >  1:3840201/1‑100
                                                                                                     GCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACA    <  2:309336/100‑1
                                                                                                       GATCTCCTCACTCTTTTTAATTCTCCTCATTCGCCATTACCGTCATTACTTGTTTTAGGCATGCCCGAGTAGCAAGTGCATAACCCGAGCGACAAACATC  >  2:2524765/1‑100
                                                                                                                                                                                                           
ATGCACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAGCGAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4557154‑4557049
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcagcgatCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATC  >  NC_000913/4557077‑4557173

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.