New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4532499 =11 (0.110)11 (0.120) 9/148 2.8 17.9% intergenic (+49/‑433) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 4532813 = 91 (0.980)intergenic (+363/‑119) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
Rejected: Frequency below/above cutoff threshold.

TAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTT  >  NZ_CP009273/4532411‑4532588
                                                                                        |                                                                                         
tcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                                          >  2:429694‑M1/89‑90 (MQ=255)
tcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                                          >  2:305357‑M1/89‑90 (MQ=255)
tcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                                          >  2:1329905‑M1/89‑90 (MQ=255)
 aatCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTG                                                                                         >  1:1381902‑M1/1‑90 (MQ=255)
  atCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGG                                                                                        >  1:445093‑M1/1‑90 (MQ=255)
    catgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgggg                                                                                      <  2:2447607‑M1/6‑1 (MQ=255)
      cATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCC                                                                                    >  1:1176819‑M1/1‑90 (MQ=255)
      cATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCC                                                                                    >  1:2340425‑M1/1‑90 (MQ=255)
      cATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCC                                                                                    >  1:1857880‑M1/1‑90 (MQ=255)
          cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                >  2:2328417‑M1/79‑90 (MQ=255)
          cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                >  1:1154015‑M1/79‑90 (MQ=255)
          aaCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                                                <  2:604246/90‑1 (MQ=255)
          aaCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                                                <  2:110172/90‑1 (MQ=255)
            gtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTg                                                                              >  2:2581748‑M1/77‑90 (MQ=255)
              aaaaaTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGC                                                                                     >  1:2854860‑M1/1‑81 (MQ=255)
              aaaaaTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGC                                                                                     <  2:2854860‑M1/81‑1 (MQ=255)
                 tctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGc                                                                                     >  1:374823‑M1/72‑78 (MQ=255)
                 tctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGc                                                                                     <  2:374823‑M1/7‑1 (MQ=255)
                  aTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcat                                                                        <  2:1914962‑M1/90‑12 (MQ=255)
                   tgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCatat                                                                       >  2:2367567‑M1/70‑90 (MQ=255)
                         ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                >  1:1660397‑M1/64‑75 (MQ=255)
                         ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                <  2:1660397‑M1/12‑1 (MQ=255)
                          tttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaat                                                                <  1:1952071‑M1/28‑1 (MQ=255)
                            gaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaataa                                                              >  1:406124/1‑90 (MQ=255)
                            gaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaataa                                                              >  1:2724802/1‑90 (MQ=255)
                                 ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGC                                                                                     >  2:2640945‑M1/1‑62 (MQ=255)
                                 ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGC                                                                                     <  1:2640945‑M1/62‑1 (MQ=255)
                                   taaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGt                                                       <  2:1754646‑M1/37‑1 (MQ=255)
                                              ttATTGGTAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCa                                            <  2:1657353/90‑1 (MQ=255)
                                              ttATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCa                                            <  2:1537168/90‑1 (MQ=255)
                                                aTTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttca                                          >  2:690669‑M1/1‑49 (MQ=255)
                                                   ccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCAt                                       >  1:1389180‑M1/38‑90 (MQ=255)
                                                           aGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgtt                               <  2:263485‑M1/90‑53 (MQ=255)
                                                             ttAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttgg                             >  1:2564454‑M1/1‑36 (MQ=255)
                                                              tAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggc                            >  2:139113‑M1/1‑35 (MQ=255)
                                                                gttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaaaa                          <  2:519357‑M1/66‑1 (MQ=255)
                                                                aaaaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                          <  1:2940946‑M1/90‑58 (MQ=255)
                                                                 aaaaaCAAATAAATACAAGACAATTGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaaaaaa                        <  2:1381902/90‑1 (MQ=255)
                                                                  aaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaa                        <  1:97415‑M1/90‑60 (MQ=255)
                                                                       aaataaatACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattc                   <  1:690669‑M1/90‑65 (MQ=255)
                                                                       aaataaatACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattc                   <  1:483176‑M1/90‑65 (MQ=255)
                                                                          taaataCAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctat                <  2:2952440‑M1/90‑68 (MQ=255)
                                                                          taaataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa                <  2:2342419/90‑1 (MQ=255)
                                                                                cAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATc          >  2:132834/1‑90 (MQ=255)
                                                                                    acatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAAc      >  1:519357‑M1/5‑90 (MQ=255)
                                                                                       aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGt   >  1:603425/1‑90 (MQ=255)
                                                                                       aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGt   >  1:2275119/1‑90 (MQ=255)
                                                                                        ttGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttattt  >  2:2109471‑M1/1‑9 (MQ=255)
                                                                                        |                                                                                         
TAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTT  >  NZ_CP009273/4532411‑4532588

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