New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 454269015 (0.580)7 (0.310) 6/114 1.3 32.6% intergenic (+57/‑425) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
?NC_000913 = 4542986 16 (0.720)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:904043/1‑68 (MQ=255)
   aGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtgg                                                                  <  2:300256‑M2/68‑4 (MQ=255)
     tatatTAAGCCAACAATTAAACTGATATAGTTTGTCCGCGATGCTTTCTTCTATGAGTCAAAATGGcc                                                                <  2:210420/67‑1 (MQ=255)
     aaTATTAAGCCAACAAATAAACTGAAAAAGTTTGGCCGCGATGCTTTCCTCTATGAGTCAAAAtggcc                                                                >  2:210418‑M2/1‑63 (MQ=255)
      atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccc                                                               <  1:198489‑M2/68‑7 (MQ=255)
      atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccc                                                               <  1:723342‑M2/68‑7 (MQ=255)
      atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccc                                                               <  1:718549‑M2/68‑7 (MQ=255)
      atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccc                                                               <  1:66740‑M2/68‑7 (MQ=255)
      atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccc                                                               <  1:664330‑M2/68‑7 (MQ=255)
      atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccc                                                               <  1:410263‑M2/68‑7 (MQ=255)
      atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGATATCCTCTATGAGTCAAAAtggccc                                                               <  1:723343‑M2/68‑7 (MQ=255)
        atTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggcccca                                                             <  2:440474‑M2/68‑9 (MQ=255)
        atTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggcccca                                                             <  1:701654‑M2/68‑9 (MQ=255)
         ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaa                                                            <  1:116841‑M2/68‑10 (MQ=255)
             gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgt                                                        <  2:348433‑M2/68‑14 (MQ=255)
             gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGt                                                        >  2:60788/1‑68 (MQ=255)
             gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGt                                                        >  2:60791/1‑68 (MQ=255)
             gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGt                                                        >  1:777365/1‑68 (MQ=255)
             gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGt                                                        >  1:4417/1‑68 (MQ=255)
                         aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtat                                                 <  2:722361‑M2/63‑21 (MQ=255)
                         aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtat                                                 >  1:722363‑M2/1‑43 (MQ=255)
                         aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTgg                                            >  1:668475/1‑68 (MQ=255)
                          cTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                           <  1:670086/68‑1 (MQ=255)
                                    ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                                 <  2:410952/68‑1 (MQ=255)
                                         cgcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgc                            <  2:701652/68‑2 (MQ=255)
                                          gcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                           <  2:4328/68‑1 (MQ=255)
                                             aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTg                        >  2:157709/1‑68 (MQ=255)
                                             aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTg                        >  1:50475/1‑68 (MQ=255)
                                               gCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatca                      >  1:640533‑M2/1‑21 (MQ=255)
                                                   tCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTc                  <  2:66740/68‑1 (MQ=255)
                                                    ccTCTATGAGTCAAAAtggccccaattttcttgtatttttttgtttttttaactttattatcaatatg                 >  1:610149‑M2/1‑16 (MQ=255)
                                                    ccTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataag                 <  2:44718‑M2/68‑53 (MQ=255)
                                                    ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCa                 >  1:932978/1‑68 (MQ=255)
                                                      tctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCaaa               >  1:583462/1‑68 (MQ=255)
                                                      tctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCaaa               <  1:358907/68‑1 (MQ=255)
                                                         aTGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaa            <  1:839612‑M2/68‑58 (MQ=255)
                                                                   aTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaa  <  2:4417/68‑1 (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.