New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 4542683 =13 (0.900)8 (0.600) 8/242 0.6 40.9% intergenic (+49/‑433) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
?NC_000913_3_bme_pgi 4542997 = 11 (0.820)intergenic (+363/‑119) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)

ATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATA  >  NC_000913_3_bme_pgi/4542537‑4542824
                                                                                                                                                  |                                                                                                                                             
aTACTGTGCGTTATACCGCCATGAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGATGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTgg                                                                                                                                            <  2:100457/150‑1 (MQ=255)
 tACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTggg                                                                                                                                           <  1:1994/150‑1 (MQ=255)
    attgccagagttttaattttcatgctgctttcctttcaaaaaactatttctaaatcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGcc                                                                                                                                        <  2:113417‑M1/8‑1 (MQ=255)
          ttATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTCGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCattttg                                                                                                                                  >  2:119861/1‑145 (MQ=255)
           tataCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCattttga                                                                                                                                 <  1:119861/150‑7 (MQ=255)
           aataCCGCAAGTAATGCTGCTCGTTATGACGGATTATGGAAAAGAAATAATCAAATAAACGAAAAATGAAAAAGAGAAGAGGTTTGATTTAACTTATAGATAAAAAAATTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                                                                                                 <  2:17968/149‑1 (MQ=255)
                    aGTAATGCTGCACGATTTGACGGATTACGGGAACGAAATAATCTCATAAACGAAAAATTAAAAAGAAAAGAGGATAGAGTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGATAATTATGGCCAAACTGTACAtattat                                                                                                                        >  2:125447/1‑150 (MQ=255)
                          tgctgctttcctttcaaaaaactatttctaaatcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaataa                                                                                                                  >  1:113417‑M1/121‑150 (MQ=255)
                              gctttcctttcaaaaaactatttctaaatcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTa                                                                                                              >  2:152576‑M1/117‑150 (MQ=255)
                                               tGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGc                                                                                             >  3:53097/1‑150 (MQ=255)
                                                                    aaCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGAc                                                                        >  3:51733/1‑150 (MQ=255)
                                                                                                     aaCTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAattat                                       >  4:6473/1‑150 (MQ=255)
                                                                                                              ataataAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAATTTACGTATTTTTTCTCAAGCATAAAAAAAAAAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGt                              <  2:95828/150‑1 (MQ=255)
                                                                                                               aacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTa                             <  3:61702‑M1/115‑1 (MQ=255)
                                                                                                                          gttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGt                                                                                                                <  3:24175‑M1/32‑1 (MQ=255)
                                                                                                                          gttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGt                                                                                                                >  4:24175‑M1/25‑56 (MQ=255)
                                                                                                                           aaaaaCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTcac                 <  4:53097/150‑1 (MQ=255)
                                                                                                                                ccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTg                                                    >  1:64162‑M1/19‑110 (MQ=255)
                                                                                                                                ccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTg                                                    <  2:64162‑M1/92‑1 (MQ=255)
                                                                                                                                    taaataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTcc        >  1:114712/1‑150 (MQ=255)
                                                                                                                                        gatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCta    >  4:61702‑M1/11‑150 (MQ=255)
                                                                                                                                          cAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCtata  <  4:92763/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
ATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATA  >  NC_000913_3_bme_pgi/4542537‑4542824

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