New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4532499 =6 (0.320)10 (0.590) 7/152 0.8 55.0% intergenic (+49/‑433) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 4532813 = 11 (0.650)intergenic (+363/‑119) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

TAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTA  >  NZ_CP009273/4532411‑4532582
                                                                                        |                                                                                   
tcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                                    >  2:152181‑M1/89‑90 (MQ=255)
taatCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAtt                                                                                    <  1:195129/90‑1 (MQ=255)
  atCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTgg                                                                                  >  1:67503/1‑90 (MQ=255)
      cATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGcc                                                                              >  1:394844/1‑90 (MQ=255)
          aaCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCattt                                                                          <  2:258585/90‑4 (MQ=255)
                    tAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                                                <  2:394844/90‑1 (MQ=255)
                      acactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGc                                                                               >  2:445296‑M1/67‑73 (MQ=255)
                      acactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGc                                                                               <  1:445296‑M1/7‑1 (MQ=255)
                          tttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgggg                                                                                <  1:237564‑M1/6‑1 (MQ=255)
                          tttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgggg                                                                                >  2:237564‑M1/63‑68 (MQ=255)
                              ttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCa                                                                    <  2:81390‑M1/18‑1 (MQ=255)
                              ttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCa                                                                    <  2:262667‑M1/18‑1 (MQ=255)
                              ttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCa                                                                    >  1:81390‑M1/59‑76 (MQ=255)
                              ttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCa                                                                    >  1:262667‑M1/59‑76 (MQ=255)
                                                       cactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaaa                             >  2:405853‑M1/34‑90 (MQ=255)
                                                          tgcacagttttcccccaaaagatgaaacatTTGGGGc                                                                               <  1:387189‑M1/7‑1 (MQ=255)
                                                          tgcacagttttcccccaaaagatgaaacatTTGGGGc                                                                               >  2:387189‑M1/31‑37 (MQ=255)
                                                           gcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat                         <  2:261434‑M1/61‑1 (MQ=255)
                                                                   ttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCa                                                                    <  2:110339‑M1/18‑1 (MQ=255)
                                                                   ttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCa                                                                    >  1:110339‑M1/22‑39 (MQ=255)
                                                                          aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa          >  1:123135‑M1/15‑90 (MQ=255)
                                                                                cAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATc    >  1:140132/1‑90 (MQ=255)
                                                                                  aaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa  >  2:217460‑M1/7‑90 (MQ=255)
                                                                                        |                                                                                   
TAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTA  >  NZ_CP009273/4532411‑4532582

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