New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4541922 =101 (1.110)8 (0.090) 6/170 NT 7.9% intergenic (+363/‑115) fimB/fimE tyrosine recombinase/inversion of on/off regulator of fimA/tyrosine recombinase/inversion of on/off regulator of fimA
?NC_000913 4541952 = 91 (1.040)intergenic (+393/‑85) fimB/fimE tyrosine recombinase/inversion of on/off regulator of fimA/tyrosine recombinase/inversion of on/off regulator of fimA
Rejected: Frequency below cutoff threshold.

CCTAACTGAGATATATTTGAATGAATACCTATAGGAAACCTCAATCGGTCAAAATTAGCCCAATAGATAAAAAGATAACCACACACGGAGTGATGTGGCTATC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4542024‑4541922
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tatcTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCG  >  NC_000913/4541952‑4542050
                                                                                                                                                                                                          
CCTAACTGAGATATATTTGAATGAATACCTATAGGAAACCTCAATCGGTCAAAATTAGCCCAATAGATAAAAAGATAACCACACACGGAGTGATGTGGCT                                                                                                        >  1:855158/1‑100
CCTAACTGAGATATATTTGAATGAATACCTATAGGAAACCTCAATCGGTCAAAATTAGCCCAATAGATAAAAAGATAACCACACACGGAGTGATGTGGCT                                                                                                        >  2:1345348/1‑100
CCTAACTGAGATATATTTGAATGAATACCTATAGGAAACCTCAATCGGTCAAAATTAGCCCAATAGATAAAAAGATAACCACACACGGAGTGATGTGGCT                                                                                                        >  2:64574/1‑100
  TAACTGAGATATATTTGAATGAATACCTATAGGAAACCTCAATCGGTCAAAATTAGCCCAATAGATAAAAAGATAACCACACACGGAGTGATGTGGCTAT                                                                                                      <  1:1648802/100‑1
   AACTGAGATATATTTGAATGAATACCTATAGGAAACCTCAATCGGTCAAAATTAGCCCAATAGATAAAAAGATAACCACACACGGAGTGATGTGGCTATC                                                                                                     <  2:856162/100‑1
         GATATATTTGAATGAATACCTATAGGAAACCTCAATCGGTCAAAATTATCACAATAGATAAAAAGATAACCACACACGGAGTGATGTGGTTATCTTTTT                                                                                                >  2:1045938/1‑99
                                AGGAAACCTCAATCGGTCAAAATTAGTCCAATAGATAAAAAGATAACCACACACGGAGTGATGTGGATATCTTTTTATCTATTGGGCTAATTTTGACGGA                                                                        >  2:672824/1‑100
                                                          CCCAATAGATAAAAAGATAACCACACACGGAGTGATGTGGCTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTC                                              >  1:102274/1‑100
                                                               TAGATAAAAAGATAACCACACACGGAGTGATGTGGCTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATA                                         >  1:568657/1‑100
                                                                AGATAAAAAGATAACCACACACGGAGTGATGTGGCTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATAT                                        >  1:107654/1‑100
                                                                 GATAAAAAGATAACCACACACGGAGTGATGTGGCTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATA                                       >  1:752083/1‑100
                                                                 GATAAAAAGATAACCACACACGGAGTGATGTGGCTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATA                                       >  2:191232/1‑100
                                                                         GATAACCACACACGGAGTGATGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTT                               >  2:56822/1‑100
                                                                         GATAACCACACACGGAGTGATGTGGCTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTAATTCAAATATATCTCAGTT                               >  2:1918511/1‑100
                                                                                                     TCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTC   >  2:937879/1‑100
                                                                                                      CTTTTTATCTATTGGGTTAATTTTGACCGATTGAGGTTTCCTATAGGTATTTATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCG  <  2:2288602/100‑1
                                                                                                                                                                                                          
CCTAACTGAGATATATTTGAATGAATACCTATAGGAAACCTCAATCGGTCAAAATTAGCCCAATAGATAAAAAGATAACCACACACGGAGTGATGTGGCTATC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4542024‑4541922
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tatcTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCG  >  NC_000913/4541952‑4542050

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