New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP000730	= 800928	28 (0.800)	4 (0.120)	4/266	3.4	14.9%	intergenic (+62/‑289)	hisC1/USA300HOU_0748	histidinol‑phosphate transaminase/hypothetical protein
	?	CP000730	801948 =	19 (0.570)					intergenic (+189/+435)	USA300HOU_0748/dgkA1	hypothetical protein/possible diacylglycerol kinase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGTTCT  >  CP000730/801851‑802080                                                                                                  |                                                                                                                                     acaacaCTATATTGGCAGTAGTTGACTGCGGGGCCCCAATATGTAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCAt                                                                                           >  1:283534/1‑141 (MQ=255) acaacaCTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCAt                                                                                           >  2:257146/1‑141 (MQ=255) acaacaCTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCAt                                                                                           >  1:41548/1‑141 (MQ=255) acaacaCTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCAt                                                                                           >  2:287891/1‑141 (MQ=255)          tatTGGTAGTATTTTACTTCTTGGCCCCAACGTAGAGAAATTGGAACTCCAATTTCTACAGACAATGCAAGTTTGGGTGGGCCTCAACATATAGGAATGCTTTTTCTTTAGAATTTAGTATTTCTTATGCATGAGTGTAAc                                                                                  <  2:251470/141‑1 (MQ=14)          tatTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAAc                                                                                  <  1:168148/141‑1 (MQ=255)                       tGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCata                                                                     >  2:179099/1‑141 (MQ=255)                        gACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCatat                                                                    <  1:87272/141‑1 (MQ=255)                        gACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCatat                                                                    <  1:77617/141‑1 (MQ=255)                                 gCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTa                                                           >  2:364435/1‑141 (MQ=255)                                 gCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTa                                                           >  2:364458/1‑141 (MQ=255)                                 gCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATGAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTa                                                           >  2:7417/1‑141 (MQ=255)                                            gagaAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCtg                                                <  2:237904/141‑1 (MQ=255)                                                       aCGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTTTTTAGATATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGAt                                     >  2:217832/1‑141 (MQ=255)                                                           cAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAg                                 <  2:207417/141‑1 (MQ=255)                                                                    cAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTAc                        >  1:251817/1‑141 (MQ=255)                                                                        cAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATaa                    >  1:79110/1‑141 (MQ=255)                                                                               aGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAatcaatca             >  2:281339/1‑141 (MQ=255)                                                                                          gCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGTTCt  >  1:45264/1‑141 (MQ=255)                                                                                          accccAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGTTCt  <  2:189632/140‑1 (MQ=255)                                                                                                  |                                                                                                                                     ACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGTTCT  >  CP000730/801851‑802080

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