New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb = 45473723 (0.080)15 (0.500) 11/82 NT 82.3% intergenic (+57/‑425) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin
?W3110S.gb = 4547668 4 (0.130)intergenic (+353/‑129) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin

ATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAA  >  W3110S.gb/4547602‑4547734
                                                                  |                                                                  
aTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtg                                                                  >  1:1589259‑M2/1‑67 (MQ=255)
  aGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggc                                                                >  1:1652719‑M2/1‑65 (MQ=255)
        ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaat                                                          <  1:167749‑M2/69‑11 (MQ=255)
             ccAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtct                                                     >  1:1111791‑M2/1‑54 (MQ=255)
             ccAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtct                                                     >  1:2738734‑M2/1‑54 (MQ=255)
                   aaTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCtt                                               <  1:1418639/69‑1 (MQ=255)
                    aTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaa                                                           >  1:2359800‑M2/1‑47 (MQ=255)
                              aaaGTTTGTCCTCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttt                                     >  1:2748578‑M2/1‑37 (MQ=255)
                                 gTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgttt                                      <  1:1935853‑M2/64‑31 (MQ=255)
                                   ttGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaac                               <  1:55442‑M2/69‑38 (MQ=255)
                                    tGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaac                               <  1:636103‑M2/68‑38 (MQ=255)
                                                   ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                     <  1:67856/47‑1 (MQ=255)
                                                   ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                     <  1:743061/47‑1 (MQ=255)
                                                   ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                     <  1:3202886/47‑1 (MQ=255)
                                                    ctctATGAGTCAAAAtggccccaattgtcttgtatttatttg                                         <  1:1687999‑M2/42‑28 (MQ=255)
                                                    ctctATGAGTCAAAAtggccccaattgtcttgtatttatttg                                         <  1:2523806‑M2/42‑28 (MQ=255)
                                                     tctATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagtt              >  1:282027‑M2/1‑14 (MQ=255)
                                                       tATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaa            <  1:2547200‑M2/68‑57 (MQ=255)
                                                       tATGAGTCAAAAtggccccaattgtcttgtatttatttgt                                        <  1:3309985‑M2/40‑29 (MQ=255)
                                                           aGTCAAGAtggccccaattgtcttgtatttatttgtttttttaactttattatcaat                   <  1:476387‑M2/57‑50 (MQ=255)
                                                           aGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaa       >  1:1968567‑M2/1‑8 (MQ=255)
                                                           aGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatca        >  1:1688816‑M2/1‑8 (MQ=255)
                                                           aGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatca        >  1:2595029‑M2/1‑8 (MQ=255)
                                                               aaaaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                            >  1:828463/1‑44 (MQ=255)
                                                               aaaaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                            >  1:3520528/1‑44 (MQ=255)
                                                               aaaaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                            >  1:678158/1‑44 (MQ=255)
                                                               aaaaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAg                            >  1:605134/1‑44 (MQ=255)
                                                                 aaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaa  <  1:829817/68‑1 (MQ=255)
                                                                  |                                                                  
ATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAA  >  W3110S.gb/4547602‑4547734

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