New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3037860 =17 (0.920)3 (0.190) 3/110 1.5 17.1% coding (247/1734 nt) recJ ssDNA exonuclease, 5' ‑‑> 3'‑specific
?NC_000913 3037943 = 15 (0.970)coding (164/1734 nt) recJ ssDNA exonuclease, 5' ‑‑> 3'‑specific
Rejected: Frequency below/above cutoff threshold.

TCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCG  >  NC_000913/3037894‑3038010
                                                 |                                                                   
gcttttcgcgaaggaacgcggattattgtggtcggtgatttcgacgccgCTCAGTTGCTGCCAGGgca                                                   >  2:682293‑M2/50‑68 (MQ=255)
  cTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGgcagc                                                 <  2:30122/68‑1 (MQ=255)
  cTTCGCGAAAAGCGGTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGgcagc                                                 <  2:592501/68‑1 (MQ=255)
      gcgAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATAc                                             <  2:52189/68‑1 (MQ=255)
        gAACAGCTTAGTATACGATCACTACGGCCTTTTCGACGCCGCGATGCTGCTGCCAGGGCAGCATACCt                                           <  2:128990/68‑1 (MQ=255)
        gAAAAGCGTTGTAACGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCt                                           <  2:442355/68‑1 (MQ=255)
        gAAAAGCGTTGTAAAGGATCTCACCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCt                                           >  1:114465/1‑68 (MQ=255)
        gAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCt                                           <  2:93051/68‑1 (MQ=255)
        gAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCt                                           <  2:568181/68‑1 (MQ=255)
        gAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCt                                           <  2:545899/68‑1 (MQ=255)
        gAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCt                                           <  2:31525/68‑1 (MQ=255)
        gAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCt                                           <  2:287244/68‑1 (MQ=255)
          aaggaacgcggattattgtggtcggtgatttcgacgccgCTCAGTTGCTGCCAGGGCAGCATACCttt                                         <  1:682297‑M2/29‑1 (MQ=255)
                    aaaGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGt                               >  2:48700/1‑68 (MQ=255)
                         aTCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCt                                           <  1:574750/51‑1 (MQ=255)
                         aTCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCt                                           >  2:574747/1‑51 (MQ=255)
                               tcggtgatttcgacgccgCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTg                    >  2:333139‑M2/19‑68 (MQ=255)
                                     ttttCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACt              <  1:302918/68‑1 (MQ=255)
                                      tttcgacgccgCTCAGTAGCTGCCAGGTCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTg             <  1:52112‑M2/57‑1 (MQ=255)
                                                 cTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCg  <  1:640238/68‑1 (MQ=255)
                                                 |                                                                   
TCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCG  >  NC_000913/3037894‑3038010

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