New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4557056114 (0.810)14 (0.110) 8/162 NT 11.3% coding (64/393 nt) iraD RpoS stabilzer after DNA damage, anti‑RssB factor
?NC_000913 = 4557068 116 (0.910)coding (76/393 nt) iraD RpoS stabilzer after DNA damage, anti‑RssB factor
Rejected: Frequency below cutoff threshold.

CAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4556950‑4557056
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atcgctgcATACCGATTTACGCAGCGATGACGTTTTATTGCCTGAAATTTCAGGTAAAACAGCCTGAAGTGATTGTCGCATCATTTTGCGCACTCCTGACGTTTAGC  <  NC_000913/4557068‑4556970
                                                                                                                                                                                                              
CAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCA                                                                                                            >  2:2931925/1‑100
CAGGGCAAGCCGCATACGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCA                                                                                                            <  2:1835315/100‑1
 AGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCAT                                                                                                           >  1:680023/1‑100
 AGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCAT                                                                                                           >  1:1349776/1‑100
  GGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATC                                                                                                          >  2:926613/1‑100
  GGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATC                                                                                                          >  2:3631566/1‑100
  GGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATC                                                                                                          >  1:1279262/1‑100
   GGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCG                                                                                                         <  2:1245522/100‑1
      AAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTG                                                                                                      <  1:2640197/100‑1
       AGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTCCAGGCAATAAAACGTCATCGCTGC                                                                                                     >  2:2714308/1‑100
                        AACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCATACCGATTTACGCAtc                                                                                    <  2:2632423/100‑3
                        AACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCATACCGATTTACGCAGC                                                                                    <  1:2100952/100‑1
                           GTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCATACCGATTAACGCAGCGAT                                                                                 <  1:2640850/100‑1
                           GTCAGGAGTACACAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCATACCGATTTACGCAGCGAT                                                                                 <  2:2584775/100‑1
                            TCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCATACCGATTTACGCAGCGATG                                                                                <  2:406782/100‑1
                             CAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGGTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCAAACCGATTTACGCAGCGAT                                                                                 <  1:1102422/98‑1
                             CAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCATACCGATTTACGCAGCGATGA                                                                               <  2:2889752/100‑1
                             CAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCATACCGATTTACGCAGCGATGA                                                                               <  1:3360645/100‑1
                              AGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCATACCGATTTACGCAGCGATGAC                                                                              <  2:132791/100‑1
                              AGGAGTGCGAAAAAGGATGCGACAATCACTTCATGCTATTTTACCTGAAAGTTCAGGCAATAATACGTCATCGCTGCATACCGATTTACGCAGCGATG                                                                                <  1:328108/98‑1
                               GGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCATACCGATTTACGCAGCGATGACG                                                                             <  2:237855/100‑1
                                 AGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCATACCGATTTACGCAGCGATGACGTT                                                                           <  1:1059551/100‑1
                                 AGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGACAATAAAACGTCATCGCTGCATACCGATTTACGCAGCGATGACG                                                                             <  1:1137247/98‑1
                                  GTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCATACCGATTTACGCAGCGATGACGTTT                                                                          <  2:2129235/100‑1
                                                                                                  aATCGCTGCATACCGATTTACGCAGCGATGACGTTTTATTGCCTGAAATTTCAGGTAAAACAGCCTGAAGTGATTGTCGCATCATTTTGAGAACTCCTGA          >  1:1342560/2‑100
                                                                                                    TCGCTGCATACCGATTTACGCAGCGATGACGTTTTATTGCCTGAAATTTCAGGTAAAACAGCCTGAAGTGATTGTCGCATCATTTTGCGCACTCCTGACG        <  1:30583/100‑1
                                                                                                       CTGCATACCGATTTACGCAGCGATGACGTTTTATTGCCTGAAATTTCAGGTAAAACAGCCTGAAGTGATTGTCGCATCATTTTGCGCACTCCTGACGTTT     >  1:3784314/1‑100
                                                                                                       CTGCATACCGATTTACGCAGCGATGACGTTTTATTGCCTGAAATTTCAGGTAAAACAGCCTGAAGTGATTGTCGCATCATTTTGCGCACTCCTGACGTTT     <  1:3166048/100‑1
                                                                                                       CTGCATACCGATTTACGCAGCGATGACGTTTTATTGCCTGAAATTTCAGGTAAAACAGCCTGAAGTGATTGTCGCATCATTTTGCGCACTCCTGACGTTT     >  2:3123821/1‑100
                                                                                                        TGCATACCGATTTACGCAGCGATGACGTTTTATTGCCTGAAATTTCAGGTAAAACAGCCTGAAGTGATTGTCGCATCATTTTGCGCACTCCTGACGTTTA    <  2:2386515/100‑1
                                                                                                        TGCATACCGATTTACGCAGCGATGACGTTTTATTGCCTGAAATTTCAGGTAAAACAGCCTGAAGTGATTGTCGCATCATTTTGCGCACTCCTGACGTTTA    >  1:3578998/1‑100
                                                                                                         GCATACCGATTTACGCAGCGATGACGTTTTATTGCCTGAAATTTCAGGTAAAACAGCCTGAAGTGATTGTCGCATCATTTTGCGCACTCCTGACGTTTAG   <  1:482882/100‑1
                                                                                                         GCATACCGATTTACGCAGCGATGACGTTTTATTGCCTGAAATTTCAGGTAAAACAGCCTGAAGTGATTGTCGCATCATTTTGCGCACTCCTGACGTTTAG   <  2:2898823/100‑1
                                                                                                         GCATACCGATTTACGCAGCGATGACGTTTTATTGCCTGAAATTTCAGGTAAAACAGCCTGAAGTGATTGTCGCATCATTTTGCGCACTCCTGACGTTTAG   >  1:289816/1‑100
                                                                                                          CATACCGATTTACGCAGCGATGACGTTTTATTGCCTGAAATTTCAGGTAAAACAGCCTGAAGTGATTGTCGCATCATTTTGCGCACTCCTGACGTTTAGC  >  1:1658050/1‑100
                                                                                                          CATACCGATTTACGCAGCGATGACGTTTTATTGCCTGAAATTTCAGGTAAAACAGCCTGAAGTGATTGTCGCATCATTTTGCGCACTCCTGACGTTTAGC  <  1:79074/100‑1
                                                                                                                                                                                                              
CAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4556950‑4557056
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atcgctgcATACCGATTTACGCAGCGATGACGTTTTATTGCCTGAAATTTCAGGTAAAACAGCCTGAAGTGATTGTCGCATCATTTTGCGCACTCCTGACGTTTAGC  <  NC_000913/4557068‑4556970

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.