New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 4542691	32 (0.700)	6 (0.140)	6/266	5.2	16.9%	intergenic (+57/‑425)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	NC_000913_3_bme_pgi	= 4542987	29 (0.670)					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.

GCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGT  >  NC_000913_3_bme_pgi/4542554‑4542838                                                                                                                                          |                                                                                                                                                    gCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTTTCCata                                                                                                                                          >  1:346126/1‑149 (MQ=255)  ccAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAACTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACTAGAAAATTAGGGCCACACTGCCCatat                                                                                                                                         >  2:807779/1‑149 (MQ=255)        aTGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaa                                                                                                                                   >  1:6990/1‑149 (MQ=255)             gctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGt                                                                                                                              <  1:233287/149‑1 (MQ=255)               tCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTa                                                                                                                            <  1:259700/149‑1 (MQ=255)                 gTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGCGAACGAGACCAAAATAACTTCTCTAGTATAAAGTTAAAAAAACAAATAATCACATGACAATTAGGGCCAAACTGTCCATATCATAAATAAGTTACg                                                                                                                          >  1:796436/1‑149 (MQ=255)                        cGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATGAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACGAATAAATACAAGGCGATTGCGGCCATACTGTCTATATCATTAATAAGTTACGTAttttt                                                                                                                   >  1:216285/1‑149 (MQ=255)                              aTGTGAAAGAAATAAACTCACAAACTAAAAATTAAAAAGAGAACAAGTTTGATTTAACTTATTTATAATAAAGCTAAAAAAACAAATAAATACAAGTCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCaa                                                                                                             <  2:719254/149‑1 (MQ=255)                               tGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAAAAAAGGTAAAAAAAAAAATAAATCCAAGACACTACAACCAAAACTGTCCATATCAAAACAAACATACGTATTTTTtctccac                                                                                                            >  2:310151/1‑146 (MQ=255)                                      gaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaaa                                                                                                     >  1:117689/1‑149 (MQ=255)                                        aataatCTCAGAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAat                                                                                                   <  1:48566/149‑1 (MQ=255)                                          taatCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat                                                                                                 <  1:61851/149‑1 (MQ=255)                                                 gttagagAAAAATGGAAAAGGAAACAGGTTAGATTTAATTGATAAATAAGAAAGTTAAAAAAACAAATAAATATAATACAATATGTTCCAAACGATCCAGGGCAAAAACAAGTTCCGTATCTTTTCTCAAGCATAACAATATTataaaa                                                                                          <  2:777905/143‑1 (MQ=255)                                                  tAAACGAACAATTAAAAAGTGAAGAGGTTTGATTTAACTTATAGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAc                                                                                         <  1:3414/149‑1 (MQ=255)                                                               aaaaaGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATc                                                                            <  2:117689/149‑1 (MQ=255)                                                                 aaaGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTTTCAAACTTAAAAATATTAAAAAACGACAAAAAGCAtata                                                                          >  1:755331/1‑149 (MQ=255)                                                                       aagagGTTTGAATTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGtt                                                                    <  1:86099/149‑1 (MQ=255)                                                                                 aTTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGAAAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTCAAATCCTACAGAAAGCATCAAAATGTTTGATATGTaa                                                          >  1:413433/1‑149 (MQ=255)                                                                                      aCTTATTGATAATAAAGTTAAAAAAACAAATACCCACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAattat                                                     >  2:378994/1‑149 (MQ=255)                                                                                                aataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTa                                           >  2:197269/1‑149 (MQ=255)                                                                                                      gTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCa                                                                                                                                            >  2:13022/1‑46 (MQ=255)                                                                                                      gTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCa                                                                                                                                            <  1:13022/46‑1 (MQ=255)                                                                                                       ttAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATAAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAacttaat                                    >  2:152641/1‑149 (MQ=255)                                                                                                           aaaaaCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTca                                <  1:351985/149‑1 (MQ=255)                                                                                                               aCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATc                            >  2:527794/1‑149 (MQ=255)                                                                                                                       ataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTTTTTTGTTTAAAAATATTAAAAAAAGACAAAATGTATATAACTGTTTGATATATAAATCATTCGTCTTTTAATTTCATTTCACATAACTTTCGc                    >  2:371296/1‑149 (MQ=255)                                                                                                                       ataCAAGACAATTGGGGCCAAACTGTCCATACCATAAATAAGTCACGTATTTTTTCTCAAGCATAAAAATAATAAAAAACAACAAAACGCATCTAACTCTTTCATATATAACTAATTTCTATCAAAAAATCATTACACCACACCTccac                    >  2:371335/1‑147 (MQ=255)                                                                                                                            aGACAATTGGGGCCAAACTGTCCATATCAAAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTatat               <  2:6990/149‑1 (MQ=255)                                                                                                                             gACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATg              <  1:378994/149‑1 (MQ=255)                                                                                                                                     gggCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCtat                 <  2:247358/138‑1 (MQ=255)                                                                                                                                     gggCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCtat                 >  1:247358/1‑138 (MQ=255)                                                                                                                                         cAAACTGGACAAATCATAAATAAGTCTCGTATTTTTTCGTAAGCATAAAAATATTAAAAAAGCACAAAAAGCATCTAACGGTTTTATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGt  <  2:392783/149‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                    GCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGT  >  NC_000913_3_bme_pgi/4542554‑4542838

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.