New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 45426907 (0.230)3 (0.110) 3/114 2.5 20.7% intergenic (+57/‑425) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
?NC_000913 = 4542986 17 (0.640)intergenic (+353/‑129) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)

AAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAA  >  NC_000913/4542623‑4542755
                                                                   |                                                                 
gttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCa                                                                   >  2:168067‑M1/60‑68 (MQ=255)
      cttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                             <  1:302715/67‑1 (MQ=255)
            agtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCatat                                                       >  2:259616‑M1/48‑68 (MQ=255)
             tAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                                      <  2:33416/68‑1 (MQ=255)
             tAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                                      <  1:640259/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)
                                                            tGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaat       <  1:167197‑M1/68‑61 (MQ=255)
                                                            tGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaat       <  2:799809‑M1/68‑61 (MQ=255)
                                                            tGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagt                        >  1:540548‑M1/1‑8 (MQ=255)
                                                            tGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagt                        <  2:540548‑M1/51‑44 (MQ=255)
                                                              gggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatat     <  2:237297‑M1/68‑63 (MQ=255)
                                                              gggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatat     <  1:627193‑M1/68‑63 (MQ=255)
                                                                gCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattc   >  2:640331‑M1/1‑4 (MQ=255)
                                                                 ccAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattct  <  1:758263‑M1/68‑66 (MQ=255)
                                                                 ccAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattct  <  1:640332‑M1/68‑66 (MQ=255)
                                                                   |                                                                 
AAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAA  >  NC_000913/4542623‑4542755

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.