New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	4532499 =	22 (0.520)	8 (0.200)	8/278	3.7	25.5%	intergenic (+49/‑433)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	CP009273	4532813 =	26 (0.650)					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.

CCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAAT  >  CP009273/4532690‑4532954                                                                                                                            |                                                                                                                                              ccAAAAACCACCTCATGCAATATAAACATCTATAAATCAAGATAACAACAGAATATTAATCCAACAAATAAACTGAAAAAGTTTGTACGAGATGGTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTCTCGGGGAAAAcgg                                                                                                                      >  2:259194/1‑147 (MQ=255) ccAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTCTGAGTCAAAATGGCCCCAAATGTTTCATCTTTCGCGGGAAAACtg                                                                                                                      >  1:223587/1‑149 (MQ=255)  cAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAgtgt                                                                                                                     >  2:60083/1‑149 (MQ=255)               aTGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCACAATGCTTCACCTTTTGGGGGAAAACTGTGCAGTGTTGGcgg                                                                                                        >  2:555081/1‑147 (MQ=255)                          aCATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCGGTCAAACTCGtt                                                                                             >  2:278366/1‑149 (MQ=255)                                 taaataaaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaa                                                                                      <  1:612445/149‑1 (MQ=255)                                   aataaaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCAGAGATGCTTTCTAATAAGATGCAAAATGGCCCCAAATGTTTCAACTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacc                                                                                    >  2:251511/1‑148 (MQ=255)                                      aaaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                                                                                 >  2:338522/1‑149 (MQ=255)                                               aaTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa                                                                        <  1:567180/149‑1 (MQ=255)                                                  aGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGa                                                                     >  2:177894/1‑149 (MQ=255)                                                  aGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGa                                                                     >  1:582914/1‑149 (MQ=255)                                                          aaGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCAt                                                             >  1:426099/1‑149 (MQ=255)                                                             ccAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTc                                                          <  1:672683/149‑1 (MQ=255)                                                             ccAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTc                                                          >  1:333995/1‑149 (MQ=255)                                                                           gAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGt                                            >  1:550200/1‑149 (MQ=255)                                                                               aaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGttttt                                        <  2:592875/149‑1 (MQ=255)                                                                                   ttGTCCTCGATGCTTTCCTCTATGAGTCAACATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaa                                    <  2:580193/149‑1 (MQ=255)                                                                                   tcgTCCTCGATGCTTTCCTCTATGAGTCAAAATGTCCCCAAATGTTTCATCTTTTGGTGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaa                                    <  2:85267/147‑1 (MQ=255)                                                                                    tgcccGCGATTCTTTCCTCTAAGAGTCAAAATGGCCCCAAATGTTTCATCATTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAACCTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaa                                   <  2:215101/146‑1 (MQ=255)                                                                                        cgcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaaggaa                               >  1:564926/1‑149 (MQ=255)                                                                                        cgcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGCACAAAACGACTGCCCATGACGATTTAGAAATAGTTTTTTgaaaggaa                               >  1:53653/1‑149 (MQ=255)                                                                                                      ctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGCAGGCAGCATGAAAAt                 >  2:109036/1‑149 (MQ=255)                                                                                                         cgAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTaa              <  1:259194/148‑1 (MQ=255)                                                                                                          gAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTaaa             <  1:177894/149‑1 (MQ=255)                                                                                                          gAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGATTTTTGAAAGGAAAGCAGCATGAAAATTaaa             >  1:554621/1‑149 (MQ=255)                                                                                                                 tttGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAAGAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTgg      >  1:656564/3‑149 (MQ=255)                                                                                                                 aaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTgg      <  2:550200/149‑1 (MQ=255)                                                                                                                     gCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAAt  <  1:538353/149‑1 (MQ=255)                                                                                                                            |                                                                                                                                              CCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAAT  >  CP009273/4532690‑4532954

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.