New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP000730	= 940821	22 (0.720)	3 (0.100)	3/258	4.5	13.0%	intergenic (+2/+54)	comA/USA300HOU_0905	hypothetical protein/hypothetical protein
	?	CP000730	= 940872	19 (0.650)					intergenic (+53/+3)	comA/USA300HOU_0905	hypothetical protein/hypothetical protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTATGGGATATAAAAATTAAGAATGACAAAGAACAATTAATTACAGTTATGCGTGGTACAGTTGCTATTAAACCTTTAAAATAAAAGAACTGCTAGCTGAAATGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGTTTTAGAAATGCTTGTTCAATGCGTTCGGCAGCTTTACGGCCACCCATAACATTTCTACCAAATGGTCCTAATTCTAAGTCTGCAAAGCATCCTGCGAC  >  CP000730/940736‑940994                                                                                                                                         |                                                                                                                           gTATGGGATATAAAAATTAAGAATGACAAAGAACAATTAATTACAGTTATGCGTGGTACAGTTGCTATTAAACCTTTAAAATAAAAGAACTGCTAGCTGAAATGTTATGAGATATTCATAACTACGGCTAGCAGttttttt                                                                                                                        >  1:226718/1‑141 (MQ=255)               aaTTAAGAATGACAAAGAACAATTAATTACAGTTATGCGTGGTACAGTTGCTATTAAACCTTTAAAATAAAAGAACTGCTAGCTGAAATGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATAttgt                                                                                                          >  1:4292/1‑141 (MQ=255)                       aTGACAAAGAACAATTAATTACAGTTATGCGTGGTACAGTTGCTATTAAACCTTTAAAATAAAAGAACTGCTAGCTGAAATGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGttt                                                                                                  <  2:297246/141‑1 (MQ=255)                       aTGACAAAGAACAATTAATTACAGTTATGCGTGGTACAGTTGCTATTAAACCTTTAAAATAAAAGAACTGCTAGCTGAAATGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGttt                                                                                                  >  2:245635/1‑141 (MQ=255)                                     ttaattaCAGTTATGCGTGGTACAGTTGCTATTAAACCTTTAAAATAAAAGAACTGCTAGCTGAAATGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGTTTTAGAAATGCTTGtt                                                                                    <  1:45654/141‑1 (MQ=255)                                                   gCGTGGTACAGTTGCTATTAAACCTTTAAAATAAAAGAACTGCTAGCTGAAATGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGTTTTAGAAATGCTTGTTCAATGCGTTCGGCa                                                                      >  2:230267/1‑141 (MQ=255)                                                           cAGTTGCTATTAAACCTTTAAAATAAAAGAACTGCTAGCTGAAATGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGTTTTAGAAATGCTTGTTCAATGCGTTCGGCAGCTTTACg                                                              >  2:206601/1‑141 (MQ=255)                                                                gCTATTAAACCTTTAAAATAAAAGAACTGCTAGCTGAAATGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGTTTTAGAAATGCTTGTTCAATGCGTTCGGCAGCTTTACGGCCAc                                                         <  1:154207/141‑1 (MQ=255)                                                                gCTATTAAACCTTTAAAATAAAAGAACTGCTAGCTGAAATGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGTTTTAGAAATGCTTGTTCAATGCGTTCGGCAGCTTTACGGCCAc                                                         <  2:25660/141‑1 (MQ=255)                                                                                    aaaGAACTGCTAGCTGAAATGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGTTTTAGAAATGCTTGTTCAATGCGTTCGGCAGCTTTACGGCAACCCATAACATTTCTACCAAAt                                     >  1:306497/1‑141 (MQ=255)                                                                                     aaGAACTGCTAGCTGAAATGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGTTTTAGAAATGCTTGTTCAATGCGTTCGGCAGCTTTACGGCCACCCATAACATTTCTACCAAATg                                    <  1:155546/141‑1 (MQ=255)                                                                                       gAACTGCTAGCTGGAATGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATTCGTTTTATTGTTGTAGTTTTAGAAATGCTTGTTCAATGCGTTCGGCAGCTTTACGGCCACCCATAACATTTCTACCAAATGGt                                  >  2:223924/1‑141 (MQ=255)                                                                                        aaCTGCTAGCTGAAATGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGTTTTAGAAATGCTTGTTCAATGCGTTCGGCAGCTTTACGGCCACCCATAACATTTCTACCAAATGGTc                                 >  1:184666/1‑141 (MQ=255)                                                                                             cTAGCTGAAATGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGTTTTAGAAATGCTTGTTCAATGCGTTCGGCAGCTTTACGGCCACCCATAACATTTCTACCAAATGGTCCTAAt                            <  1:223693/141‑1 (MQ=255)                                                                                              tAGCTGAAATGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGTTTTAGAAATGCTTGTTCAATGCGTTCGGCAGCTTTACGGCCACCCATAACATTTCTACCAAATGGTCCTAAtt                           >  1:181530/1‑141 (MQ=255)                                                                                              tAGCTGAAATGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGTTTTAGAAATGCTTGTTCAATGCGTTCGGCAGCTTTACGGCCACCCATAACATTTCTACCAAATGGTCCTAAtt                           <  2:4292/141‑1 (MQ=255)                                                                                                    aaaTGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGTTTTAGAAATGCTTGTTCAATGCGTTCGGCAGCTTTACGGCCACCCATAACATTTCTACCAAATGGTCCTAATTCTAAGt                     >  1:154888/1‑141 (MQ=255)                                                                                                    aaaTGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGTTTTAGAAATGCTTGTTCAATGCGTTCGGCAGCTTTACGGCCACCCATAACATTTCTACCAAATGGTCCTAATTCTAAGt                     <  2:38440/141‑1 (MQ=255)                                                                                                          tATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGTTTTAGAAATGCTTGTTCAATGCGTTCGGCAGCTTTACGGCCACCCATAACATTTCTACCAAATGGTCCTAATTCTAAGTCTGCaa               <  1:278593/141‑1 (MQ=255)                                                                                                                       tAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGTTTTAGAAATGCTTGTTCAATGCGTTCGGCAGCTTTACGGCCACCCATAACATTTCTACCAAATGGTCCTAATTCTAAGTCTGCAAAGCATCCTGCGAc  >  1:72313/1‑141 (MQ=255)                                                                                                                                         |                                                                                                                           GTATGGGATATAAAAATTAAGAATGACAAAGAACAATTAATTACAGTTATGCGTGGTACAGTTGCTATTAAACCTTTAAAATAAAAGAACTGCTAGCTGAAATGTTATGAGATATTCATAACTACGGCTAGCAGTTTTTTTATGCGCTATATTGTTGTAGTTTTAGAAATGCTTGTTCAATGCGTTCGGCAGCTTTACGGCCACCCATAACATTTCTACCAAATGGTCCTAATTCTAAGTCTGCAAAGCATCCTGCGAC  >  CP000730/940736‑940994

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ 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.