New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	= 4542606	27 (0.690)	8 (0.220)	8/268	3.5	21.8%	intergenic (+57/‑425)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	NC_000913_3_pae_tpiA	= 4542902	32 (0.880)					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.

GTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACG  >  NC_000913_3_pae_tpiA/4542461‑4542752                                                                                                                                                  |                                                                                                                                                   gTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                                                                                                                 >  1:234019/1‑149 (MQ=255) gTTATACCGCCAGTAATGCTGCTAGTTACGCGATAAGATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTAATGAGAATAAAGTTAAAAAAAAAAATAAATACAAGTCAGTTGGGGCCAAAc                                                                                                                                                 >  1:676163/1‑149 (MQ=255)    ataCCTCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCattttga                                                                                                                                              >  2:353331/1‑143 (MQ=255)    ataCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGGAAGAAATAATCTCATAATCCAAAAAGTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                                                                                                              <  1:501775/149‑1 (MQ=255)                    tgctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCat                                                                                                                                          <  1:194095/137‑1 (MQ=255)                    tgctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCat                                                                                                                                          >  2:194095/1‑137 (MQ=255)                        cGTTCTGCCGTAGTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGATGTTTGATTTAACTCATTGATAATAAAGTTAAAAAAACAAATAAGTACAAGACTATTGGGTCCTAACTGTTCATATCATAAATAAGTTAc                                                                                                                          <  2:362642/149‑1 (MQ=255)                             tGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTAtt                                                                                                                     <  1:51098/149‑1 (MQ=255)                                    ttATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTAGGGCCAAAATGTCCATATCATAAATATGTTACGTATTTTTtctc                                                                                                              >  1:298414/1‑149 (MQ=255)                                       tGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGCATTATCCGGCATAAAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAg                                                                                                           <  1:501046/149‑1 (MQ=255)                                        gggAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTAGTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGc                                                                                                          >  1:474383/1‑149 (MQ=255)                                         ggcaagataTATTCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTAGGGCCAAACTTTCCCTATCATAAATAAGTTCCGTCTTTTTTCTCAAGcc                                                                                                         >  2:254925/4‑148 (MQ=255)                                          gaaagaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATCAATAAGTGACGTAGTTTTTCTCAAGCAt                                                                                                        >  1:601316/1‑149 (MQ=255)                                           aaagaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCaa                                                                                                            <  1:531463/144‑1 (MQ=255)                                           aaagaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCaa                                                                                                            >  2:531463/1‑144 (MQ=255)                                           aaagaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATa                                                                                                       >  2:52672/1‑149 (MQ=255)                                               aaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaaaa                                                                                                   <  2:278866/149‑1 (MQ=255)                                               aaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGGTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaaaa                                                                                                   <  1:515014/149‑1 (MQ=255)                                                                     ttAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAACCCACCAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCa                                                                             <  1:641115/149‑1 (MQ=255)                                                                      tAAAAAGAGAATAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCAt                                                                            <  2:190865/149‑1 (MQ=255)                                                                                                ttATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTAATAAGTACATTAttt                                                  >  2:27572/1‑149 (MQ=255)                                                                                                         ataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtctc                                                                                                              >  2:319880/1‑80 (MQ=255)                                                                                                         ataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtctc                                                                                                              <  1:319880/80‑1 (MQ=255)                                                                                                               ttAAAAAAACAAATAAATACAAGACACTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAAGTGTTTGATATGTAAATTATTTTTATTGTAaattgtt                                   >  1:52376/1‑146 (MQ=255)                                                                                                                     aaaCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTcaca                             >  2:141767/1‑149 (MQ=255)                                                                                                                                    aGACAATTGGGGCCGAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaaa                                                                                <  2:78192/83‑1 (MQ=255)                                                                                                                                    aGACAATTGGGGCCGAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaaa                                                                                >  1:78192/1‑83 (MQ=255)                                                                                                                                                ccAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAGGCTAGCg  >  1:205408/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACG  >  NC_000913_3_pae_tpiA/4542461‑4542752

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.