New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 4542691	35 (0.790)	8 (0.190)	8/268	4.1	19.1%	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	35 (0.840)					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.

GTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAAC  >  NC_000913_3_bme_pgi/4542843‑4543129                                                                                                                                                 |                                                                                                                                               gtgGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGc                                                                                                                                            <  1:402938/149‑1 (MQ=255)      tCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCCTGCAATATAAACATCTATAAATAAAGATAACACTAGAATATTAAGCAAACAAATAAACTAAAAAAGTTTGTCAGCGCTGCGTTCCTCTATGAGTCAAAATGGACCCaa                                                                                                                                       >  2:281406/1‑149 (MQ=255)              tCTTCCTCATTCTTCTCTCCAACAACCAACAATTACAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAATTGtcttg                                                                                                                               >  2:746416/1‑145 (MQ=255)               gttCCTCATTGTTCTCTCCAAAAACCACCTCATGCCATATGAACATCGATAAATAAAGATAACAATAGAATCTTAAGCCATCAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCAt                                                                                                                              <  2:320962/148‑1 (MQ=255)                   cTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                                                                                      <  2:77848/137‑1 (MQ=255)                   cTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                                                                                                      >  1:77848/1‑137 (MQ=255)                   cTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                                                                                          >  1:356445/1‑149 (MQ=255)                   cTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                                                                                          >  1:481667/1‑149 (MQ=255)                             tctcCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTTCGAGGGTGGTGCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaaa                                                                                                                >  1:411839/1‑149 (MQ=255)                               tcCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCCAAATGGCCACAATTGTTTCATCTTTTGGGGGAAAAca                                                                                                              >  1:143170/1‑148 (MQ=255)                                 cAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGCGGGAAAACcgt                                                                                                            >  2:696127/1‑149 (MQ=255)                                 cAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAAGGTTTCATCTTTTGGGGGAAAACtgt                                                                                                            >  2:620412/1‑149 (MQ=255)                                    aaaCCACCTCATGCCATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCa                                                                                                         <  1:370407/149‑1 (MQ=255)                                       ccaccTCATGCAATATAAACATCTATATATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATAGTTTGGGGGAAAACTGTGCAgtg                                                                                                      >  2:476946/1‑149 (MQ=255)                                                    tataAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAAc                                                                                         <  1:449414/149‑1 (MQ=255)                                                         aCATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGtt                                                                                    <  2:87350/149‑1 (MQ=255)                                                             ctctaaacaaaGTTAACAATAGAATATTAAGCCAAATAATAATCTCAAAAGGTTTGCCCGCGATGCTTTTCTCTATGCGTCCTTGTGGCCCCAAATGTTTTCTCTTCTCGGTGAAACCCGTGCTGCGTTGTCAGTCAAACTCGTTTata                                                                                <  2:128072/146‑3 (MQ=255)                                                               ataaataaAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaa                                                                              >  2:448573/1‑149 (MQ=255)                                                                         aaaaCAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTa                                                                    <  2:355179/147‑1 (MQ=255)                                                                              aaTAGAATATTAAGCCAACAAATAAGCTGAAAAAGTTTGTCCGAGATGCTTTCCTCTATGAGTCAAAATAGCCCCAAATGTTTCATCTTTTGGGGGAAAACGGTGCAGTATTGGCAGTCAAACTCGTTTAGAAAACAAAGTGTACAGaa                                                               >  1:588363/1‑149 (MQ=255)                                                                                     tatCAAGCCAACAAATAAACTGAACAAGTTTGTCCGCGATGCTTCCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGc                                                        <  1:281406/149‑1 (MQ=255)                                                                                      atTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAGGTGTACAGAACGACTGcc                                                       >  2:560571/1‑149 (MQ=255)                                                                                       ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGccc                                                      <  1:468016/149‑1 (MQ=255)                                                                                           gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGGGTACAGAACGACTGCCCATGt                                                  >  1:2420/1‑149 (MQ=255)                                                                                             caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCCTGTCg                                                >  1:586294/1‑149 (MQ=255)                                                                                             caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAAGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCg                                                >  1:507065/1‑149 (MQ=255)                                                                                               acaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAt                                              <  2:202411/149‑1 (MQ=255)                                                                                                         tGAAAAAGTTTGTCCGCGATGCTTTCCTTTATGAGTCAAAATGGACCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAg                                    <  2:613361/149‑1 (MQ=255)                                                                                                            aaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGTTTTAGAAATAGttt                                 >  2:317950/1‑149 (MQ=255)                                                                                                              aaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGttttt                               <  1:412594/149‑1 (MQ=255)                                                                                                                   tGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaa                          <  2:481667/149‑1 (MQ=255)                                                                                                                   tGTCCGCGATGCTGTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaa                          <  1:345707/149‑1 (MQ=255)                                                                                                                     tCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaagg                        >  2:76013/1‑149 (MQ=255)                                                                                                                               tttCCTCTATGAGTCAAAATGGCCCCACACCCATGTCGCGGTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCAt              >  1:574151/1‑149 (MQ=255)                                                                                                                               tttCCTCTATGAGGCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCAt              >  1:269821/1‑149 (MQ=255)                                                                                                                                          aGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAgtgt                                                                     <  1:500668/83‑1 (MQ=255)                                                                                                                                          aGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAgtgt                                                                     >  2:500668/1‑83 (MQ=255)                                                                                                                                           gTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTGGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAAc  >  1:38915/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                               GTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAAC  >  NC_000913_3_bme_pgi/4542843‑4543129

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.