New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4532499 =4 (0.130)23 (0.860) 15/150 0.4 79.9% intergenic (+49/‑433) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 4532813 = 8 (0.300)intergenic (+363/‑119) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

ATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGA  >  NZ_CP009273/4532737‑4532900
                                                                            |                                                                                       
ttttgtcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCtttt                                                                            >  1:148152‑M2/77‑90 (MQ=255)
ttttgtcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCtttt                                                                            >  1:148153‑M2/77‑90 (MQ=255)
ttttgtcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCtttt                                                                            <  2:187269‑M2/14‑1 (MQ=255)
ttttgtcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCtttt                                                                            >  1:85835‑M2/77‑90 (MQ=255)
        ttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGaaa                                                                    >  1:777100‑M2/69‑90 (MQ=255)
        ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                                                    >  1:9208/1‑90 (MQ=255)
               aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGc                                                             >  2:498911/1‑90 (MQ=255)
               aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGc                                                             >  2:233215/1‑90 (MQ=255)
               aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGc                                                             >  1:92641/1‑90 (MQ=255)
                            tgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTc                                                >  1:324581‑M2/49‑90 (MQ=255)
                            tgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTc                                                >  1:630400‑M2/49‑90 (MQ=255)
                                     gTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGtt                                       <  1:498911/90‑1 (MQ=255)
                                        cgtaacttatttatgatatggacagtttggccccaaATGTTTCATCtttt                                                                            <  2:807975‑M2/14‑1 (MQ=255)
                                        cgtaacttatttatgatatggacagtttggccccaaATGTTTCATCtttt                                                                            >  1:807975‑M2/37‑50 (MQ=255)
                                            acttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaa                                <  1:864796‑M2/58‑1 (MQ=255)
                                              ttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaaca                              <  1:259462‑M2/60‑1 (MQ=255)
                                               tatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa                             <  2:324581‑M2/61‑1 (MQ=255)
                                                atttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaaa                            >  2:27487‑M2/29‑90 (MQ=255)
                                                      ctATGAGTCAAAATGGCCCCAAATGTTTCATCtttt                                                                            <  1:618844/36‑1 (MQ=255)
                                                      ctATGAGTCAAAATGGCCCCAAATGTTTCATCtttt                                                                            >  2:618844/1‑36 (MQ=255)
                                                       atatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACa                     >  1:731432‑M2/22‑90 (MQ=255)
                                                           aGTCGAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAAc                 >  2:749251/1‑90 (MQ=255)
                                                              cagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCGGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGAc              >  1:802899‑M2/15‑90 (MQ=255)
                                                              cagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGAc              <  1:449322‑M2/76‑1 (MQ=255)
                                                                gtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTg            >  1:192423‑M2/13‑90 (MQ=255)
                                                                gtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTg            >  2:533941‑M2/13‑90 (MQ=255)
                                                                gtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGGGCAGTGTTGGCAGTCAAACTCGTTTAAAAAAAAAAGTGTAAAGAACGACTg            >  1:110339‑M2/13‑90 (MQ=255)
                                                                gtttggccccaaAGGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTAAAAAACAAAGTGTACAGAACGACTg            >  1:7025‑M2/13‑90 (MQ=255)
                                                                  ttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGtt                                                       <  1:169775‑M2/35‑1 (MQ=255)
                                                                  ttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGtt                                                       >  2:169775‑M2/11‑45 (MQ=255)
                                                                  ttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGcc          >  2:685365‑M2/11‑90 (MQ=255)
                                                                         caaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCg   <  1:233215‑M2/87‑1 (MQ=255)
                                                                         caaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCg   <  1:27487‑M2/87‑1 (MQ=255)
                                                                          aaATGTTTCATCTTTTGGGGGAAAACTGTGCAGGGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGa  >  1:189532‑M2/3‑90 (MQ=255)
                                                                            |                                                                                       
ATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGA  >  NZ_CP009273/4532737‑4532900

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