New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 45325070 (0.000)14 (0.470) 10/146 1.0 93.3% intergenic (+57/‑425) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 = 4532803 2 (0.070)intergenic (+353/‑129) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

CATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGA  >  NZ_CP009273/4532716‑4532887
                                                                                       |                                                                                    
cATCCATAAATGAAGATAACAATAGAATATTAAGCCAAAAAATAAACTGAAAAAGTTTGTCCGCGATGCATTCCTCTATGAGTCAAAAtg                                                                                    >  2:503155‑M2/1‑88 (MQ=255)
   cTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggcc                                                                                 <  1:745617‑M2/90‑6 (MQ=255)
     ataaataaAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggcccc                                                                               <  2:945322‑M2/90‑8 (MQ=255)
     ataaataaAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggcccc                                                                               <  2:61693‑M2/90‑8 (MQ=255)
        aataaaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTTCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaat                                                                            >  2:833408‑M2/1‑80 (MQ=255)
                        gAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttattt                                                            <  2:4838‑M2/90‑27 (MQ=255)
                          atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatt                                                                 <  1:892080‑M2/83‑22 (MQ=255)
                          atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatt                                                                 >  2:892080‑M2/1‑62 (MQ=255)
                             ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttt                                                       <  1:107956‑M2/90‑32 (MQ=255)
                             ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                                       >  1:203913/1‑90 (MQ=255)
                               atgccAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaat                                                                            <  1:335106‑M2/65‑11 (MQ=255)
                               atgccAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaat                                                                            >  2:335106‑M2/3‑57 (MQ=255)
                               aaGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttt                                                     <  2:436244‑M2/90‑34 (MQ=255)
                                  ccAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaa                                                  <  1:670634‑M2/90‑37 (MQ=255)
                                    aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaact                                                <  2:984302‑M2/90‑39 (MQ=255)
                                    aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaat                                                                            <  2:119321‑M2/62‑11 (MQ=255)
                                    aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaat                                                                            >  1:119321‑M2/1‑52 (MQ=255)
                                                gAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaata                                    <  2:513823‑M2/90‑51 (MQ=255)
                                                                     tttCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaaa               <  1:945322/90‑1 (MQ=255)
                                                                                  tCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttc  <  2:704294‑M2/90‑85 (MQ=255)
                                                                                       |                                                                                    
CATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGA  >  NZ_CP009273/4532716‑4532887

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