New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 21373017 (0.620)43 (2.170)
+AGGACTTCTTGAA
4/70 1.7 75.3% coding (203/255 nt) rof modulator of Rho‑dependent transcription termination
?NC_000913 3958648 = 22 (0.800)coding (679/1476 nt) ilvC ketol‑acid reductoisomerase (NADP(+))

GTTGCGGAAGTGAAATCTGACCTGATGGGCGAGCAAACCATCCTGTGCGGTATGTTGCAGGCTGGCTCTCTGCTGTGCTTCGACAAGCTGGTGGAAG  >  NC_000913/3958600‑3958696
                                                |                                                
gTTGCGGAAGTGAAATCTGACCTGATGGGCGAGCAAACCATCCTGTGCg                                                  >  1:692954/1‑49 (MQ=255)
gTTGCGGAAGTGAAATCTGACCTGATGGGCGAGCAAACCATCCTGTGCg                                                  <  1:2078981/49‑1 (MQ=255)
gTTGCGGAAGTGAAATCTGACCTGATGGGCGAGCAAACCATCCTGTGCg                                                  >  1:794031/1‑49 (MQ=255)
 ttGCGGAAGTGAAATCTGACCTGATGGGCGAGCAAACCATCCTGTGCgg                                                 >  1:457964/1‑49 (MQ=255)
   gCGGAAGTGAAATCTGACCTGATGGGCGAGCAAACCATCCTGTGCGGTa                                               >  1:1981972/1‑49 (MQ=255)
          tGAAATCTGACCTGATGGGCGAGCAAACCATCCTGTGCGGTATGTTg                                          >  1:335407/1‑47 (MQ=255)
          tGAAATCTGACCTGATGGGCGAGCAAACCATCCTGTGCGGTATGTTg                                          >  1:867376/1‑47 (MQ=255)
            aaaTCTGACCTGATGGGCGAGCAAACCATCCTGTGCGGTATGTTGCAgg                                      >  1:1182836/1‑49 (MQ=255)
            aaaTCTGACCTGATGGGCGAGCAAACCATCCTGTGCGGTATGTTGCAgg                                      <  1:1935764/49‑1 (MQ=255)
                 tatggctaaagctggtaaaggacttcttgaaGGTACGTTGCAggctgg                                  >  1:1791002‑M2/32‑48 (MQ=255)
                 tatggctaaagctggtaaaggacttcttgaaGGTACGTTGCAggctgg                                  >  1:551606‑M2/32‑48 (MQ=255)
                 tatggctaaagctggtaaaggacttcttgaaGGTACGTTGCAggctgg                                  >  1:711153‑M2/32‑48 (MQ=255)
                 tatggctaaagctggtaaaggacttcttgaaGGTACGTTGCAggctgg                                  >  1:1691891‑M2/32‑48 (MQ=255)
                 tatggctaaagctggtaaaggacttcttgaaGGTACGTTGCAggctgg                                  >  1:1235302‑M2/32‑48 (MQ=255)
                 tatggctaaagctagtaaaggacttcttgaaGGTACGTTGCAggctgg                                  >  1:793328‑M2/32‑48 (MQ=255)
                  gACCTGATGGGCGAGCAAACCATCCTGTGCGGTATGTTGCAggctggct                                >  1:319246/1‑49 (MQ=255)
                  atggctaaagctggtaaaggacttcttgaaGGTACGTTGCaggctggca                                >  1:1462623‑M2/31‑48 (MQ=255)
                  atggctaaagctggtaaaggacttcttgaaGGTACGTTGCaggctggca                                >  1:1472983‑M2/31‑48 (MQ=255)
                   tggctaaagctggtaaaggacttcttgaaGGTACGTTGcaggctggcac                               <  1:1050710‑M2/20‑3 (MQ=255)
                   tggctaaagctggtaaaggacttcttgaaGGTACGTTGcaggctggcac                               >  1:2123531‑M2/30‑47 (MQ=255)
                   tggctaaagctggtaaaggacttcttgaaGGTACGTTGcaggctggcac                               <  1:2014973‑M2/20‑3 (MQ=255)
                   tggctaaagctggtaaaggacttcttgaaGGTACGTTGcaggctggcac                               >  1:1551477‑M2/30‑47 (MQ=255)
                   tggctaaagctggtaaaggacttcttgaaGGTACGTTGcaggctggcac                               >  1:1177417‑M2/30‑47 (MQ=255)
                   tggctaaagctggtaaaggacttcttgaaGGTACGTTGCaggctggca                                <  1:873063‑M2/19‑2 (MQ=255)
                   tggctaaagctggtaaaggacttcttgaaGGTACGTTGCaggctggca                                <  1:1778699‑M2/19‑2 (MQ=255)
                   tggctaaagctggtaaaggacttcttgaaGGTACGTTGCaggctggca                                <  1:186992‑M2/19‑2 (MQ=255)
                   tggctaaagctggtaaaggacttcttgaaGGTACGTTGCaggctggca                                <  1:75617‑M2/19‑2 (MQ=255)
                   tggctaaagctggtaaaggacttcttgaaGGTACGTTGCaggctggca                                <  1:1329121‑M2/19‑2 (MQ=255)
                   tggctaaagctggtaaaggacttcttgaaGGTACGTTGCaggctggca                                <  1:225096‑M2/19‑2 (MQ=255)
                   tggctaaagctggtaaaggacttcttgaaGGTACGTTGCAggctggct                                <  1:1551517‑M2/19‑1 (MQ=255)
                   tggctaaagctggtaaaggacttcttgaaGGTACGTTGCAggctggc                                 <  1:1285579‑M2/18‑1 (MQ=255)
                   tggctaaagctggtaaaggacttcttgaaGGTACGTTGCAggctggc                                 <  1:1318636‑M2/18‑1 (MQ=255)
                    ggctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTGGcact                              >  1:311653‑M2/29‑49 (MQ=255)
                    ggctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTGGcact                              >  1:1399351‑M2/29‑49 (MQ=255)
                    ggctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTGGcact                              >  1:140713‑M2/29‑49 (MQ=255)
                    ggctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTGGcact                              >  1:1586595‑M2/29‑49 (MQ=255)
                     gctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTGgcactt                             <  1:2116350‑M2/22‑2 (MQ=255)
                     gctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTGgcactt                             >  1:598452‑M2/28‑48 (MQ=255)
                     gctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTGgcactt                             >  1:1447647‑M2/28‑48 (MQ=255)
                     gctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTGgcactt                             <  1:448770‑M2/22‑2 (MQ=255)
                     gctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTGgcactt                             <  1:1318930‑M2/22‑2 (MQ=255)
                     gctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTGGcact                              <  1:62218‑M2/21‑1 (MQ=255)
                     gctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTGGcact                              <  1:721474‑M2/21‑1 (MQ=255)
                     gctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTGGcact                              <  1:513884‑M2/21‑1 (MQ=255)
                      ctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTggcacttc                            >  1:1347075‑M2/27‑47 (MQ=255)
                      ctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTggcacttc                            <  1:281719‑M2/23‑3 (MQ=255)
                      ctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTggcacttc                            >  1:833763‑M2/27‑47 (MQ=255)
                      ctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTggcacttc                            >  1:430259‑M2/27‑47 (MQ=255)
                      ctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTggcacttc                            >  1:703213‑M2/27‑47 (MQ=255)
                      ctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTggcacttc                            <  1:479034‑M2/23‑3 (MQ=255)
                      ctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTggcacttc                            >  1:962386‑M2/27‑47 (MQ=255)
                      ctaaagctggtaaaggacttcttgaaGGTACGTTGCAGGCTggcactta                            <  1:1132001‑M2/23‑3 (MQ=255)
                      ctaaagctggtaaaggacttcttgaaGGGACGTTGCAGGCTggcacttc                            >  1:70696‑M2/27‑47 (MQ=255)
                        aTGGGCGAGCAAACCATCCTGTGCGGTATGTTGCAGGCTGGCTCTctg                           <  1:1670721/48‑1 (MQ=255)
                         tGGGCGAGCAAACCATCCTGTGCGGTATGTTGCAGGCTGGCTATctgct                         <  1:192052/49‑1 (MQ=255)
                           ggCGAGCAAACCATCCTGTGCGGTATGTTGCAGGCTGGCTCTCTGCtgt                       >  1:2064101/1‑49 (MQ=255)
                           ggCGAGCAAACCATCCTGTGCGGTATGTTGCAGGCTGGCTCTCTGCtgt                       <  1:1465962/49‑1 (MQ=255)
                                 cAAACCATCCTGTGCGGTATGTTGCAGGCTGGCTCTCTGCTGTGCTTCg                 <  1:473266/49‑1 (MQ=255)
                                  aaaCCATCCTGTGCGGTATGTTGCAGGCTGGCTCTCTGCTGTGCTTCGa                <  1:1446833/49‑1 (MQ=255)
                                      cATCCTGTGCGGTATGTTGCAGGCTGGCTCTCTGCTGTGCTTCGACAAg            >  1:1669545/1‑49 (MQ=255)
                                       aTCCTGTGCGGTATGTTGCAGGCTGGCTCTCTGCTGTGCTTCGACAAGc           >  1:528297/1‑49 (MQ=255)
                                        tCCTGTGCGGTATGTTGCAGGCTGGCTCTCTGCTGTGCTTCGACAAGCt          <  1:385431/49‑1 (MQ=255)
                                         ccTGTGCGGTATGTTGCAGGCTGGCTCTCTGCTGTGCTTCGACAAGCtg         >  1:2074273/1‑49 (MQ=255)
                                          cTGTGCGGTATGTTGCAGGCTGGCTCTCTGCTGTGCTTCGACAAGCtgg        >  1:253999/1‑49 (MQ=255)
                                          cTGTGCGGTATGTTGCAGGCTGGCTCTCTGCTGTGCTTCGACAAGCtg         <  1:410567/48‑1 (MQ=255)
                                                ggTATGTTGCAGGCTGGCTCTCTGCTGTGCTTCGACAAGCtggtgg     >  1:1396103/1‑46 (MQ=255)
                                                ggTATGTTGCAGGCTGGCTCTCTGCTGTGCTTCGACAAGCTGGTGgaag  <  1:1318640/49‑1 (MQ=255)
                                                |                                                
GTTGCGGAAGTGAAATCTGACCTGATGGGCGAGCAAACCATCCTGTGCGGTATGTTGCAGGCTGGCTCTCTGCTGTGCTTCGACAAGCTGGTGGAAG  >  NC_000913/3958600‑3958696

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