New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1458931 =14 (0.540)3 (0.120) 3/124 2.1 17.1% coding (668/789 nt) paaG 1,2‑epoxyphenylacetyl‑CoA isomerase, oxepin‑CoA‑forming
?NC_000913 1458971 = 16 (0.660)coding (708/789 nt) paaG 1,2‑epoxyphenylacetyl‑CoA isomerase, oxepin‑CoA‑forming
Rejected: Frequency below/above cutoff threshold.

AAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTA  >  NC_000913/1458904‑1459023
                                                                   |                                                    
aagcGATAAATAGCGCTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGa                                                      >  2:737599/1‑68 (MQ=255)
 gataatctgcgctccgtccggcaagtcgctgatagtcacgttccagatccagttgcgtatcgagcgaa                                                     >  1:584482‑M2/67‑67 (MQ=255)
   cGATAAATAGCGCTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACtt                                                   <  1:439792/68‑1 (MQ=255)
     aTAAATAGCGCTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGc                                                 <  1:159440/68‑1 (MQ=255)
         aTAGCGCTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCcgga                                             >  2:768402/1‑68 (MQ=255)
         aTAGCGCTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCcgga                                             >  2:386525/1‑68 (MQ=255)
              gcTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAg                                        >  2:916718/1‑68 (MQ=255)
                        aTACGCTCGTTACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTAt                              <  1:794530/68‑1 (MQ=255)
                         ttCGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATc                             >  1:584482‑M1/3‑68 (MQ=255)
                                 aTACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGgt                     <  1:404188/68‑1 (MQ=255)
                                      acgttccagatccagttgcgtatcgagcgACTTGCCGGACGGAGCGCAGATTATCGTGAAg                       >  1:270256‑M2/30‑61 (MQ=255)
                                      acgttccagatccagttgcgtatcgagcgACTTGCCGGACGGAGCGCAGATTATCGTGAAg                       <  2:270258‑M2/32‑1 (MQ=255)
                                       aaCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGt               <  1:205322/68‑1 (MQ=255)
                                       aaCTGGATCTGGAACGCGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGt               <  2:134846/68‑1 (MQ=255)
                                          tGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCg            >  2:37198/1‑68 (MQ=255)
                                           ggATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGt           <  1:883335/68‑1 (MQ=255)
                                            gATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGtt          <  1:958411/68‑1 (MQ=255)
                                             aTCTGGGACGAGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGATc         >  1:538666/1‑68 (MQ=255)
                                                   aaCGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCt   >  2:194675/1‑68 (MQ=255)
                                                    aCGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTa  <  1:916714/68‑1 (MQ=255)
                                                                   |                                                    
AAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTA  >  NC_000913/1458904‑1459023

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