New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? CP009273 4532499 =15 (0.780)3 (0.170) 3/258 NT 18.7% intergenic (+49/‑433) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
?CP009273 4532813 = 12 (0.670)intergenic (+363/‑119) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)

TTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAAT  >  CP009273/4532674‑4532954
                                                                                                                                           |                                                                                                                                             
ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCATCTATGAGTAAACATGGCCCCaaagagatcatc                                                                                                                                      >  2:259473/1‑140 (MQ=255)
ttCCTCATTCTTCTCTCCAAAAACCACCTCATACAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGTCCCCAAATGTTACATc                                                                                                                                      >  1:10071/1‑149 (MQ=255)
 tCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCt                                                                                                                                     >  1:261979/1‑149 (MQ=255)
     cATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGGGTCAAAATGGCCCCAAATGTTTCATCTTTTg                                                                                                                                 >  2:260055/1‑149 (MQ=255)
      ggtgatgtgaaattaatttacaatagaaataatttacatatcaaacagttagatgctttttgtcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTgg                                                                                                                                <  2:328778‑M2/16‑1 (MQ=255)
                ccAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAACGTTTGTCCTCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtg                                                                                                                      >  1:332519/1‑149 (MQ=255)
                                                    taaaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGAGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGCAAACTGTGCAGTGTTGGCAGTCAAACCCGTTTacaaaaccaa                                                                                  >  2:176435/1‑149 (MQ=255)
                                                      ttagatgctttttgtcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAgt                                                                                <  2:182861‑M2/64‑1 (MQ=255)
                                                        gATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTATATCCGCGATGCTATCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAgtgt                                                                              >  1:322482/1‑149 (MQ=255)
                                                                       ttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAGCAAAGTGTACAGAACGACTGccc                                                               >  1:25140‑M2/69‑149 (MQ=255)
                                                                              aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTATACAAAACAAAGTGGACAGAACGACTggcactttcga                                                        >  2:3407/1‑139 (MQ=255)
                                                                              aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATAGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGa                                                        >  1:110146/1‑149 (MQ=255)
                                                                                   aTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCAAAATGTTTCATCTTTTGGGGGAAAACTGTGAAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGCCTGCCCATGTCGATTTAg                                                   >  1:55186/1‑149 (MQ=255)
                                                                                                           aTGCTTTCCTCGATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAagca                           <  1:34020/149‑1 (MQ=255)
                                                                                                                  ccTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACTACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGaa                    <  1:310016/149‑1 (MQ=255)
                                                                                                                   ctctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTGGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGACTAAGAAATAGTTTTGTGAAAGGAAAGAAGCATGaaa                   >  2:287676/1‑149 (MQ=255)
                                                                                                                                    gccccaaATGTTTCAACTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTTGATTTATAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAAt  <  1:262853‑M2/142‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                             
TTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAAT  >  CP009273/4532674‑4532954

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