New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	4542683 =	24 (0.540)	8 (0.190)	7/268	4.2	27.4%	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 =	20 (0.480)					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.

ATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAG  >  NC_000913_3_bme_pgi/4542537‑4542831                                                                                                                                                   |                                                                                                                                                     aTACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTg                                                                                                                                                    >  1:589052/1‑149 (MQ=255)                  gCCAGTACTGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGTCCCCAAACTGTCGata                                                                                                                                   >  2:352299/1‑149 (MQ=255)                         aTGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaa                                                                                                                            >  1:636070/1‑149 (MQ=255)                            ctggtCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAAGACAAGACACTTGGGGCCAAACTGTCCATATCataaataa                                                                                                                         <  1:414899/149‑1 (MQ=255)                              gctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGt                                                                                                                       >  1:269381/1‑149 (MQ=255)                               ctagTTATGTCGGATGATAGGAAAGAAATAATCACATAAACGAAAAATTAAAAAGATAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGtt                                                                                                                      <  2:382536/146‑1 (MQ=255)                                                            aatCTCATGAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATAtt                                                                                         <  1:681767/149‑1 (MQ=255)                                                              tctcATAAACGAACAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaa                                                                                       <  1:485841/149‑1 (MQ=255)                                                                 cATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGGTAAAAAAACAAATAAAAACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTGTCAATCATAAAAATATTaaaaa                                                                                    >  1:265457/1‑149 (MQ=255)                                                                  agaaaCGAAATATGAAAAAGAGAAGACGTTTGATTTAATGTAGTGATAATAAAGTTAAAAAAACAAAAAAATGCGAGACGATTGGGGCTCAACGTTACATATCGTAAATAAGTTGCGTATTTTTTCTCAAGCATAAAAATATTaaaaaa                                                                                   <  2:321965/147‑1 (MQ=255)                                                                      aCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAAGAc                                                                               >  1:482499/1‑149 (MQ=255)                                                                             aTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGc                                                                        <  1:451291/149‑1 (MQ=255)                                                                                   aagagaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATAATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTaa                                                                  >  1:388661/1‑149 (MQ=255)                                                                                            ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGat                                                         <  2:589052/149‑1 (MQ=255)                                                                                            gatttaatttaaCTAATTAATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGAGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGat                                                         <  1:300613/147‑1 (MQ=255)                                                                                                   tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTaaa                                                  <  2:636070/149‑1 (MQ=255)                                                                                                         ttATTGATCATACAGTTAAAAAAACAAATAAATACAAGACCATTGGTGCCAAACTGTCCATATCATAAATTAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTAttt                                            <  1:567810/149‑1 (MQ=255)                                                                                                         ttATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAACCATCTAACTGTTTGATATGTAAATTAttt                                            >  2:267926/1‑149 (MQ=255)                                                                                                                   tgaaGTTAAAAAAACAAAAAAAGACAAGACAATGGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTaaa                                  <  2:36741/147‑1 (MQ=255)                                                                                                                   taaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTaaa                                  >  1:269055/1‑149 (MQ=255)                                                                                                                            aaaatCAAAGTAATAAAAGACAATAGGGTCGAAACTGACCAAATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAAAAAAAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTca                         <  2:465819/149‑1 (MQ=255)                                                                                                                                     taaataCAAGACAATTGGGGCCAAACTGTCAATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTc                <  2:104467/149‑1 (MQ=255)                                                                                                                                            aaGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTCTTTCTATTGCAAATTAATTCCACATCACCCCCGCcata         >  2:108337/1‑149 (MQ=255)                                                                                                                                                 aaTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTaa    >  1:132274/1‑149 (MQ=255)                                                                                                                                                   ttGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAg  >  1:9893/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     ATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAG  >  NC_000913_3_bme_pgi/4542537‑4542831

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.