New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	= 4542606	26 (0.880)	4 (0.140)	4/272	3.5	13.7%	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	26 (0.940)					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.

GTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACT  >  NC_000913_3_pae_tpiA/4542758‑4543045                                                                                                                                                 |                                                                                                                                                gtgGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGc                                                                                                                                             >  1:223617/1‑149 (MQ=255)      tCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCGATGTAATTATCTATAAATAAGGATAACAATAGAATATTAAGACAACCAATAAACTGAAAAAGTTTGTCCGCGATGGTTTCCTCTATGAGTCCATATGGcccaac                                                                                                                                        >  1:176626/1‑146 (MQ=255)          cGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCCAATGt                                                                                                                                    >  1:323711/1‑149 (MQ=255)                     cATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTg                                                                                                                         <  1:214391/149‑1 (MQ=255)                           tctctcCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaa                                                                                                                   >  1:25218/1‑149 (MQ=255)                                   aaaaCCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGc                                                                                                           <  2:409639/149‑1 (MQ=255)                                           cTCATGAAATATAAACATCAAAAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTgg                                                                                                   <  2:25218/149‑1 (MQ=255)                                              aTGCAATATAAACATCTATAAGTAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAg                                                                                                >  1:394701/1‑149 (MQ=255)                                               tGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                                                                                               <  2:254780/149‑1 (MQ=255)                                                      taAACATCTATAAATAAAGATAACAATAGCATCACATGGCAGCTGTGAGACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACtt                                                                                        <  2:500720/149‑2 (MQ=255)                                                           aTCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCAGCTTTTGGAGGACAACTGTGCAGTGTTGGCAGGCCAACTCGTGTa                                                                                   >  1:21934/1‑149 (MQ=255)                                                                taaataaaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaa                                                                              >  1:107248/1‑149 (MQ=255)                                                                    taaaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaaa                                                                          <  2:223617/149‑1 (MQ=255)                                                                         aTAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTa                                                                     >  1:163182/1‑149 (MQ=255)                                                                              aaTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTCCAGaa                                                                >  2:359921/1‑149 (MQ=255)                                                                              aaTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa                                                                >  1:219126/1‑149 (MQ=255)                                                                              aaTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTAAAAAACAAAGTGTGCAGaa                                                                >  1:270024/1‑149 (MQ=255)                                                                                tAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACg                                                              >  1:447797/1‑149 (MQ=255)                                                                                         aaGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGCTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCAt                                                     <  2:2597/149‑1 (MQ=255)                                                                                           tctagcaaATAAACTGAAAAAGTTTGTCCTCGATGCTTTCCTCTATGAGTCGAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGt                                                   <  2:494914/144‑1 (MQ=255)                                                                                                 aaaTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTGACAATGGACTCAAATCTTTCCCCTTTTGGGGCTCAAGTGTGCAGTGTTCGCGCTCTAACGCGCTTACAAAACAAAGGGTACAGAACGACTGCCCATGTCGAtgt                                             >  2:386001/1‑147 (MQ=255)                                                                                                 aaaTAAACTGAAAAAGTTTGGCTGGAATGCTTTCCTCTATGAGTCAAAATTGCCCCAAATGTTTCATCTTTTGGGGGAAAATTGTGCAGTGTTGGCAGTCCAACTCGTTTACAAAACAAAGTGTAAAGAACGACTGCCCATGTCGAttt                                             <  2:194250/149‑1 (MQ=255)                                                                                                  aaTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTa                                            <  1:65565/149‑1 (MQ=255)                                                                                                               aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGtttttt                               <  1:221422/149‑1 (MQ=255)                                                                                                               aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACCCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGtttttt                               <  1:162116/149‑1 (MQ=255)                                                                                                                   tGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaa                           <  1:300706/149‑1 (MQ=255)                                                                                                                                    tctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGaaaa          >  2:341183/1‑149 (MQ=255)                                                                                                                                            tCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAACCAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGACAGCCGCACGACAGTTAAAAct  >  1:180587/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                GTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACT  >  NC_000913_3_pae_tpiA/4542758‑4543045

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.