New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4532499 =14 (0.180)13 (0.190) 9/148 2.5 30.8% intergenic (+49/‑433) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 4532813 = 46 (0.680)intergenic (+363/‑119) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

TAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGT  >  NZ_CP009273/4532411‑4532587
                                                                                        |                                                                                        
tcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                                         >  2:1824938‑M1/89‑90 (MQ=255)
tcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                                         >  2:134957‑M1/89‑90 (MQ=255)
tcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                                         >  1:1818901‑M1/89‑90 (MQ=255)
taatCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATT                                                                                         <  2:1354933‑M1/90‑1 (MQ=255)
taatCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATT                                                                                         <  1:880302‑M1/90‑1 (MQ=255)
     tcATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGC                                                                                    <  1:1650667‑M1/90‑1 (MQ=255)
     tcATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGC                                                                                    <  1:1651348‑M1/90‑1 (MQ=255)
     tcATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGC                                                                                    <  1:810299‑M1/90‑1 (MQ=255)
      cATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCC                                                                                   >  2:417913‑M1/1‑90 (MQ=255)
          cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                               >  1:1835973‑M1/79‑90 (MQ=255)
          aaCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                                               <  1:1403508/90‑1 (MQ=255)
          aaCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                                               <  2:724913/90‑1 (MQ=255)
          aaCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                                               <  1:742332/90‑1 (MQ=255)
          aaCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                                               <  1:312748/90‑1 (MQ=255)
           aCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCA                                                                                  >  1:1482362‑M1/1‑86 (MQ=255)
           aCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCA                                                                                  <  2:1482362‑M1/86‑1 (MQ=255)
            gtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTg                                                                             <  2:231707‑M1/14‑1 (MQ=255)
            gtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTg                                                                             <  2:207563‑M1/14‑1 (MQ=255)
                  aTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcat                                                                       <  1:2113394‑M1/90‑12 (MQ=255)
                           agaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaata                                                              <  1:1285412/90‑1 (MQ=255)
                                 ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgc                                                        <  1:417913‑M1/90‑27 (MQ=255)
                                              ttATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACCATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCa                                           <  1:2036898/90‑1 (MQ=255)
                                               tATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttc                                          <  2:466249‑M1/90‑41 (MQ=255)
                                                 ttGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGc                                        >  1:964824/1‑90 (MQ=255)
                                                        actgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaaaa                                 <  1:1865191‑M1/58‑1 (MQ=255)
                                                          tgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                               <  2:835529‑M1/12‑1 (MQ=255)
                                                          tgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                               >  1:835529‑M1/31‑42 (MQ=255)
                                                          tgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAata                               <  2:918863‑M1/60‑1 (MQ=255)
                                                           aGTTAAAAAAACAAATAAATACAAGACTATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgtt                              <  1:1977851‑M1/90‑53 (MQ=255)
                                                            cacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATAtt                             <  1:2027160‑M1/62‑1 (MQ=255)
                                                                gttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCat                                                                        >  1:1914644‑M1/25‑43 (MQ=255)
                                                                gttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCat                                                                        <  2:1914644‑M1/19‑1 (MQ=255)
                                                                gttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaata                                                              <  1:2143956‑M1/29‑1 (MQ=255)
                                                                gttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaata                                                              >  2:2143956‑M1/25‑53 (MQ=255)
                                                                  aaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaa                       <  2:1850003‑M1/90‑60 (MQ=255)
                                                                          taaataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa               <  1:1949269/90‑1 (MQ=255)
                                                                          aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa               >  2:2027160‑M1/15‑90 (MQ=255)
                                                                                  aGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaa                       >  1:1926822‑M1/1‑15 (MQ=255)
                                                                                  aGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaa                       <  2:1926822‑M1/74‑60 (MQ=255)
                                                                                  aGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                         <  1:713517‑M1/72‑58 (MQ=255)
                                                                                  aGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                         >  2:713517‑M1/1‑15 (MQ=255)
                                                                                  aGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa       >  1:1684625/1‑90 (MQ=255)
                                                                                     cAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCa                                           >  1:1507091/1‑51 (MQ=255)
                                                                                     cAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCa                                           <  2:1507091/51‑1 (MQ=255)
                                                                                       aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCa                                           >  1:1551378/1‑49 (MQ=255)
                                                                                       aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGt  >  1:2284214/1‑90 (MQ=255)
                                                                                       aTTGGGGCCAAACTGCCCATATCATAAATAAGTTACGTATTTTTTCTCa                                           <  2:1551378/49‑1 (MQ=255)
                                                                                        |                                                                                        
TAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGT  >  NZ_CP009273/4532411‑4532587

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