New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4542682 =7 (0.230)6 (0.230) 5/114 1.9 34.3% intergenic (+49/‑433) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
?NC_000913 4542996 = 17 (0.640)intergenic (+363/‑119) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)

AAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAA  >  NC_000913/4542618‑4542745
                                                                |                                                               
actttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgg                                                              <  2:449931‑M1/4‑1 (MQ=255)
 ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatttgtg                                                             >  1:95346‑M1/64‑66 (MQ=255)
 ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTggg                                                             >  2:1143869‑M1/64‑68 (MQ=255)
   ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGc                                                           <  1:711868‑M1/7‑1 (MQ=255)
     gttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCa                                                         >  2:168067‑M1/60‑68 (MQ=255)
           cttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                   <  1:302715/67‑1 (MQ=255)
                 agtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCatat                                             >  2:259616‑M1/48‑68 (MQ=255)
                  tAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                            <  1:640259/68‑1 (MQ=255)
                  tAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                            <  2:33416/68‑1 (MQ=255)
                  gtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATc                                            >  2:1002575‑M1/47‑68 (MQ=255)
                   aaCTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcataga                                           <  2:980907‑M1/68‑15 (MQ=255)
                     tgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCata                                         >  1:697956‑M1/44‑68 (MQ=255)
                     cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagagg                                         >  2:1086374‑M1/1‑52 (MQ=255)
                            caacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGt                                  >  1:449931‑M1/37‑68 (MQ=255)
                            caacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGt                                  >  1:608361‑M1/37‑68 (MQ=255)
                                 ctgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGt                             <  1:868399‑M1/37‑1 (MQ=255)
                                                aataaataCAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagt              <  1:1086379‑M1/68‑44 (MQ=255)
                                                     ataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaaa         >  1:1145313/1‑68 (MQ=255)
                                                         aaGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat     >  2:312594/1‑68 (MQ=255)
                                                            aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaa  >  1:830701/1‑68 (MQ=255)
                                                            aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaa  >  1:1120150/1‑68 (MQ=255)
                                                                |                                                               
AAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAA  >  NC_000913/4542618‑4542745

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