New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4542683 =	34 (0.790)	3 (0.070)	3/264	6.3	7.2%	intergenic (+49/‑433)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	NC_000913	4542997 =	46 (1.130)					intergenic (+363/‑119)	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.

TGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTA  >  NC_000913/4542543‑4542822                                                                                                                                             |                                                                                                                                            tgCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGACAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATAGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTg                                                                                                                                           >  1:40896/1‑143 (MQ=255)   cGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTggg                                                                                                                                         <  1:478925/143‑1 (MQ=255)       ataCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCa                                                                                                                                     >  1:238642/1‑143 (MQ=255)            gCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGCGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCattttg                                                                                                                                >  2:166661/1‑138 (MQ=255)                        gctCGTTTTGCCGGATGATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGAGTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataa                                                                                                                    <  1:60805/143‑1 (MQ=255)                            gtcgcgCCGTATTAGCGCAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAATAGGTTTTATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaataa                                                                                                                <  2:668770/138‑1 (MQ=255)                            gTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaacaa                                                                                                                <  2:213877/143‑1 (MQ=255)                             ttttGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAACAACAAAAAAATAAAAGACAATGGGGGCCAAACTGTCCATATCATAAATAAg                                                                                                               >  1:496432/1‑143 (MQ=255)                                tgacgaaTTATGCGAAAGAAATAAACGCAAAAACGAAAAATTAAAAAGAGAAGAAGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTa                                                                                                            <  1:573285/137‑1 (MQ=255)                                      aTTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCAAAACTGTCCATATCATAAATAAGTTACGTAtt                                                                                                      >  1:277160/1‑143 (MQ=255)                                         aTGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTAAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTAttttt                                                                                                   <  1:17550/143‑1 (MQ=255)                                           gggAAAGAAATAATATCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtc                                                                                                 <  1:683598/143‑1 (MQ=255)                                           gggAAAGAAAGAAGCTCAGAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTAGTGATAATAAAGATAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTAATTTGtc                                                                                                 <  1:309532/143‑1 (MQ=255)                                             gaaagaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtctc                                                                                               >  2:16292/1‑143 (MQ=255)                                              aaagaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCa                                                                                              >  2:188966/1‑143 (MQ=255)                                               aagaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCaa                                                                                             <  2:699154/143‑1 (MQ=255)                                                                   aaaaaTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAc                                                                         >  2:509828/1‑143 (MQ=255)                                                                    aaaaTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACg                                                                        <  1:146963/143‑1 (MQ=255)                                                                          aaaaaGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa                                                                  >  2:292729/1‑143 (MQ=255)                                                                                        tttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTTAAGCATAAAAATATTCAAAAACAACAAAAAGCATCTAACTAtt                                                    >  1:118600/1‑143 (MQ=255)                                                                                        tttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAACAAACGACAAAAAGCATCTAACTGtt                                                    >  2:332008/1‑143 (MQ=255)                                                                                         ttGATTTAACTGATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGttt                                                   <  2:316724/143‑1 (MQ=255)                                                                                                       tGATAGTACAGTAAGAAAAACAAATAAATAGAAGACAATTGGGGCCAAACTGTACATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAattat                                     <  1:105521/143‑1 (MQ=255)                                                                                                         ataataAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTAttt                                   <  2:40896/143‑1 (MQ=255)                                                                                                         ataataAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCAATATCAAAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAattatta                                   >  1:461480/1‑142 (MQ=255)                                                                                                             taaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTAt                               >  2:213504/1‑143 (MQ=255)                                                                                                             taaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTAt                               <  2:238642/143‑1 (MQ=255)                                                                                                             taaAGTTAAAAAAACAAATAAATACAAAACAATAGGGGCCAAACTGTCCATATCATAAATAAGTAACATATTTTTTATCAAGCATAAAAATATTAAAAAACGAAAAAAAGCATCTAACTCTTTGAGATATAAattattaatat                               >  1:110756/1‑138 (MQ=255)                                                                                                                    aaaaaaaCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAaat                        <  2:556863/143‑1 (MQ=255)                                                                                                                       aaaaCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTGTGCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATGGTaactaaa                     >  1:245919/1‑138 (MQ=255)                                                                                                                        aaaCAAATAAATACAAGACAATTGGGGCCAAACTGACCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATAAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAaattaat                    <  1:111341/143‑1 (MQ=255)                                                                                                                          aCAAAAAAATATAAGACAATTGCGGCCAAACTGTCCATATCATAAAAAAGTGACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAAttt                  <  1:16292/143‑1 (MQ=255)                                                                                                                              ataaataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCACGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTTTTAGATATGTAAATTATTTCTATTGTAAATTAATTTcaca              >  1:656703/1‑143 (MQ=255)                                                                                                                                     tgaaaCATTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTc       <  2:446447/139‑1 (MQ=255)                                                                                                                                          cAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCta  <  2:227250/143‑1 (MQ=255)                                                                                                                                             |                                                                                                                                            TGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTA  >  NC_000913/4542543‑4542822

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