New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 454269024 (0.650)3 (0.090) 3/208 3.5 8.7% intergenic (+57/‑425) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
?NC_000913 = 4542986 41 (1.200)intergenic (+353/‑129) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAA  >  NC_000913/4542552‑4542828
                                                                                                                                          |                                                                                                                                          
cGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCa                                                                                                                                            >  1:580398/1‑139 (MQ=255)
  ccAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaaa                                                                                                                                          <  2:78447/139‑1 (MQ=255)
             gctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaaa                                                                                                                                          >  1:194225/1‑128 (MQ=255)
             gctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaaa                                                                                                                                          <  2:194225/128‑1 (MQ=255)
                cGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCata                                                                                                                            >  2:108855/1‑139 (MQ=255)
                                       aaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCa                                                                                                     <  2:580398/139‑1 (MQ=255)
                                                                            gTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTaa                                                                <  1:108855/139‑1 (MQ=255)
                                                                                           ttGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTaa                                                 >  1:268516/1‑139 (MQ=255)
                                                                                                aataaAGTTAAAAAAACAAATAAATACAATACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTATAAAACGACAAAAAGCATCTaa                                                                <  2:484872/119‑1 (MQ=255)
                                                                                                aataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTaa                                                                >  1:484872/1‑119 (MQ=255)
                                                                                                  taaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTAttt                                          >  2:771929/1‑139 (MQ=255)
                                                                                                                 aaataaatACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAaattaat                           <  2:268516/139‑1 (MQ=255)
                                                                                                                            aGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCAcc                <  1:653353/139‑1 (MQ=255)
                                                                                                                             gACAATTGGGGCCAAACTGTCCATATCATAAATATGTTACGTATTTTTtctc                                                                                                      <  2:461608/52‑1 (MQ=255)
                                                                                                                             gACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtctc                                                                                                      >  1:461608/1‑52 (MQ=255)
                                                                                                                             gACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGata                                                      <  1:448432/100‑1 (MQ=255)
                                                                                                                             gACAATTGGGGCCAAACTATCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGata                                                      >  2:448432/1‑100 (MQ=255)
                                                                                                                                aaTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCg            <  2:390575/139‑1 (MQ=255)
                                                                                                                                  ttGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGttt                                                          <  2:495581/91‑1 (MQ=255)
                                                                                                                                  ttGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCATCTCCGCt          >  2:221436/1‑139 (MQ=255)
                                                                                                                                  gtGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTATTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGttt                                                          >  1:495581/2‑91 (MQ=255)
                                                                                                                                       gCCAAACTGTCCATATCGTAAATAAGTTACGTATTTTTTCACAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATg     >  2:638722/1‑139 (MQ=255)
                                                                                                                                          aaaCTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACa                                                                            >  2:410768/1‑65 (MQ=255)
                                                                                                                                          aaaCTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACa                                                                            <  1:410768/65‑1 (MQ=255)
                                                                                                                                          aaaCTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTaa  >  2:750958/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                                          
CGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAA  >  NC_000913/4542552‑4542828

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