New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4532499 =1 (0.080)10 (0.860) 9/148 0.3 77.2% intergenic (+49/‑433) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 4532813 = 5 (0.430)intergenic (+363/‑119) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

AAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGA  >  NZ_CP009273/4532727‑4532900
                                                                                      |                                                                                       
gttagatgctttttgtcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGt                                                                                      >  1:69264‑M2/87‑90 (MQ=255)
                  ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                                                    >  1:240088/1‑90 (MQ=255)
                  ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                                                    >  1:369566/1‑90 (MQ=255)
                    aaGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                                                  >  2:377187/1‑90 (MQ=255)
                                                   gtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTaca                                   >  1:309292‑M2/36‑90 (MQ=255)
                                                         tatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa                             >  1:347946‑M2/30‑90 (MQ=255)
                                                           tttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                           <  2:254916‑M2/63‑1 (MQ=255)
                                                           gtCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCAAAAAGGGAATAACGCTGCACTTTTTGCACTCAAACTCGTTTGCACAACAAAg                           <  1:349121/89‑1 (MQ=255)
                                                                gatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                                           >  1:361211‑M2/23‑53 (MQ=255)
                                                                gatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                                           <  2:361211‑M2/31‑1 (MQ=255)
                                                                        cAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGAc              <  2:240088/90‑1 (MQ=255)
                                                                          gtttggccccaaATGTTTCATCTTTTGGGGGAAAACtgtg                                                              <  1:359747‑M2/28‑1 (MQ=255)
                                                                          gtttggccccaaATGTTTCATCTTTTGGGGGAAAACtgtg                                                              >  2:359747‑M2/13‑40 (MQ=255)
                                                                          gtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTgg                                                     >  2:250162‑M2/13‑49 (MQ=255)
                                                                          gtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTgg                                                     <  1:250162‑M2/37‑1 (MQ=255)
                                                                               gccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCAt       >  1:333279‑M2/8‑90 (MQ=255)
                                                                                 cccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGt     <  1:330143‑M2/85‑1 (MQ=255)
                                                                                    aaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGa  >  1:96134‑M2/3‑90 (MQ=255)
                                                                                      |                                                                                       
AAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGA  >  NZ_CP009273/4532727‑4532900

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