New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4557049 =129 (1.050)10 (0.090) 9/162 NT 8.4% coding (57/393 nt) iraD RpoS stabilzer after DNA damage, anti‑RssB factor
?NC_000913 4557077 = 100 (0.890)coding (85/393 nt) iraD RpoS stabilzer after DNA damage, anti‑RssB factor
Rejected: Frequency below cutoff threshold.

TATGCACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAGCGAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4557155‑4557049
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcagcgatCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATC  >  NC_000913/4557077‑4557173
                                                                                                                                                                                                            
TATGCACTTGCCACTCGGGCATGCCTTAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACG                                                                                                          <  2:2256076/100‑1
TATGCACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACG                                                                                                          <  2:557589/100‑1
TATGCACTTGCCACACGAGCATGCGTGAAACAAGTAATGAAGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGCTCGCTGCATACCGATTTACG                                                                                                          <  1:593784/100‑1
 ATGCACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGC                                                                                                         <  2:1181852/100‑1
  TGCACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCA                                                                                                        >  2:235801/1‑100
   GCACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAG                                                                                                       >  2:382053/1‑100
     ACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCA                                                                                                        >  1:2105645/1‑97
        TGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAAAGAAAAAGAGTGAGGAGATCGCTGCATACCGACTTACGCAGCGAT                                                                                                   >  1:2789472/1‑99
                                                           GAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTT                                               >  1:1340456/1‑100
                                                             ATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCA                                             >  2:1843761/1‑100
                                                              TTAAAAAGAGTGAGGCGATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAG                                            <  2:1340462/100‑1
                                                                       GTGAGGAGATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCG                                   >  2:1984273/1‑100
                                                                        TGAGGAGATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCC                                    >  1:2616050/1‑98
                                                                          AGGAGATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGCCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGT                                >  2:3225129/1‑100
                                                                            GAGATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGG                              >  2:2874177/1‑100
                                                                             AGATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGC                             <  1:1843750/100‑1
                                                                             AGATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGC                             >  2:1424433/1‑100
                                                                             AGATCGCTGCATACCGATTTACGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGC                             >  1:241639/1‑100
                                                                                                    CAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAA      <  2:528313/100‑1
                                                                                                    CAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAA      >  2:1214927/1‑100
                                                                                                    CAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAA       >  1:3027334/1‑99
                                                                                                     AGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAAC     <  2:401066/100‑1
                                                                                                      GCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACA    <  2:1398246/100‑1
                                                                                                      GCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACA    >  1:545679/1‑100
                                                                                                        GATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATn  <  1:187719/100‑2
                                                                                                        GATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATC  >  2:1537018/1‑100
                                                                                                        GATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGCGACAAACATC  <  1:1048712/100‑1
                                                                                                                                                                                                            
TATGCACTTGCCACTCGGGCATGCCTGAAACAAGTAATGACGGTAATGCCGAATGAGGAGAATTAAAAAGAGTGAGGAGATCGCTGCATACCGATTTACGCAGCGAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4557155‑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.