New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4532499 =3 (0.100)9 (0.320) 8/148 1.2 56.8% intergenic (+49/‑433) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 4532813 = 11 (0.390)intergenic (+363/‑119) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

ATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAT  >  NZ_CP009273/4532742‑4532901
                                                                       |                                                                                        
atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggggg                                                                        <  2:310746/90‑1 (MQ=255)
atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggggg                                                                        <  1:935561/90‑1 (MQ=255)
  gttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTggggg                                                                        >  2:845155‑M2/70‑88 (MQ=255)
  gttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTggggg                                                                        <  1:845155‑M2/19‑1 (MQ=255)
   ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                                                     >  2:691218/1‑90 (MQ=255)
         caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCAAAATGTTTTATCTTTTGGGGGAAAACtgtg                                                               >  2:897850/1‑90 (MQ=255)
          aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                                          <  2:8471/62‑1 (MQ=255)
          aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                                          >  1:8471/1‑62 (MQ=255)
          aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCtttt                                                                             >  1:351485/1‑75 (MQ=255)
          aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCtttt                                                                             <  2:351485/75‑1 (MQ=255)
          aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGc                                                              >  1:440781/1‑90 (MQ=255)
                    gcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCa                                                    <  1:571779‑M2/39‑1 (MQ=255)
                                          tatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa                              >  2:25734‑M2/30‑90 (MQ=255)
                                           atttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaaa                             <  2:509934‑M2/62‑1 (MQ=255)
                                            tttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                            <  1:352088‑M2/63‑1 (MQ=255)
                                            ttCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                            <  1:445603/90‑1 (MQ=255)
                                                        tCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGa                >  2:748699/1‑90 (MQ=255)
                                                          agtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACt              >  1:817492‑M2/14‑90 (MQ=255)
                                                           gtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTg             >  2:182775‑M2/13‑90 (MQ=255)
                                                           gtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTg             >  2:441499‑M2/13‑90 (MQ=255)
                                                               ggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCa         >  2:165569‑M2/9‑90 (MQ=255)
                                                                      aATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAt  <  1:12954‑M2/89‑1 (MQ=255)
                                                                       |                                                                                        
ATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAT  >  NZ_CP009273/4532742‑4532901

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