New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4542682 =7 (0.230)6 (0.230) 5/114 1.9 34.3% intergenic (+49/‑433) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
?NC_000913 4542996 = 17 (0.640)intergenic (+363/‑119) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)

AGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGT  >  NC_000913/4542931‑4543059
                                                                 |                                                               
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)
             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)
                            acgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAg                                   <  1:868399‑M2/31‑1 (MQ=255)
                                 acttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTg                              >  1:608361‑M2/33‑68 (MQ=255)
                                 acttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTg                              >  1:449931‑M2/33‑68 (MQ=255)
                                 aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTg                              >  2:456195/1‑68 (MQ=255)
                                        tatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCa                       >  1:697956‑M2/26‑68 (MQ=255)
                                        ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                          <  2:1082474/49‑1 (MQ=255)
                                        ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                          >  1:1082479/1‑49 (MQ=255)
                                           gatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAAc                    >  2:1002575‑M2/23‑68 (MQ=255)
                                            atatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACt                   >  2:259616‑M2/22‑68 (MQ=255)
                                                   cAAAATGGCCCCTATTGTTTCATCTTTTGGGGGAAATCTGTGCAGTGTTGGCAGTCAAACTCGTTTac            <  2:948503/68‑1 (MQ=255)
                                                       aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaga        >  2:423226/1‑66 (MQ=255)
                                                        tggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaac       >  2:168067‑M2/10‑68 (MQ=255)
                                                          gccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa     <  1:711868‑M2/61‑1 (MQ=255)
                                                            cccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg   >  2:1143869‑M2/6‑68 (MQ=255)
                                                            cacaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg   >  1:95346‑M2/6‑68 (MQ=255)
                                                             ccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAgt  <  2:449931‑M2/64‑1 (MQ=255)
                                                                 |                                                               
AGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGT  >  NC_000913/4542931‑4543059

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