New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 45426907 (0.230)3 (0.110) 3/114 2.5 20.7% intergenic (+57/‑425) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
?NC_000913 = 4542986 17 (0.640)intergenic (+353/‑129) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)

AATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAA  >  NC_000913/4542919‑4543053
                                                                   |                                                                   
aaTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCaaaa                                                                     <  1:522211/68‑1 (MQ=255)
   aGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtgg                                                                  <  1:758263‑M2/68‑4 (MQ=255)
   aGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtgg                                                                  <  1:640332‑M2/68‑4 (MQ=255)
    gAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggc                                                                 >  2:640331‑M2/1‑64 (MQ=255)
      atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccc                                                               <  1:627193‑M2/68‑7 (MQ=255)
      atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccc                                                               <  2:237297‑M2/68‑7 (MQ=255)
        atTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggcccca                                                             <  1:167197‑M2/68‑9 (MQ=255)
        atTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggcccca                                                             <  2:799809‑M2/68‑9 (MQ=255)
            aGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATg                                                         <  2:924373/68‑1 (MQ=255)
             gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCTAAATGGCCCCAAATGt                                                        >  1:592410/1‑68 (MQ=255)
             gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGt                                                        >  1:140717/1‑68 (MQ=255)
                         acagAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGATTCATCTATTgg                                            >  2:1124994/4‑68 (MQ=255)
                         aCTGAAAAAGTTTGTCCTAGATGCTTTACTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTATgg                                            >  1:457652/1‑68 (MQ=255)
                         aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatt                                            <  1:1086379‑M2/68‑26 (MQ=255)
                         aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggcccca                                                             >  1:540548‑M2/1‑43 (MQ=255)
                         aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggcccca                                                             <  2:540548‑M2/51‑9 (MQ=255)
                         aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTgg                                            >  1:998812/1‑68 (MQ=255)
                         aCTGAAAAAGTTTGTCAGCGATGCTTTCCTCTAAGAGTCAAAAAGGCAACAAAAGTTTAATCTTTTgg                                            >  1:161040/1‑68 (MQ=255)
                          cTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                           <  2:627192/68‑1 (MQ=255)
                               aaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                      >  2:20752/1‑68 (MQ=255)
                                 aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                    >  2:201117/1‑68 (MQ=255)
                                       tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGcg                              >  2:1147081/1‑67 (MQ=255)
                                        ccGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                             <  1:799812/68‑1 (MQ=255)
                                             aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTg                        >  2:456195/1‑68 (MQ=255)
                                                    ccTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataag                 >  2:1086374‑M2/1‑16 (MQ=255)
                                                    ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                    >  1:1082479/1‑49 (MQ=255)
                                                    ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                    <  2:1082474/49‑1 (MQ=255)
                                                      tctATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagtt               <  2:980907‑M2/68‑55 (MQ=255)
                                                               cAAAATGGCCCCTATTGTTTCATCTTTTGGGGGAAATCTGTGCAGTGTTGGCAGTCAAACTCGTTTac      <  2:948503/68‑1 (MQ=255)
                                                                   aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaga  >  2:423226/1‑66 (MQ=255)
                                                                   |                                                                   
AATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAA  >  NC_000913/4542919‑4543053

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