New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	895785 =	17 (0.890)	4 (0.240)	4/246	1.8	21.1%	coding (806/822 nt)	metQ_2	Methionine‑binding lipoprotein MetQ
	?	USA300TCH1516_ALE	= 909740	15 (0.890)					intergenic (‑663/‑371)	USA300HOU_RS04500/csbD_1	hypothetical protein/Stress response protein CsbD

AAGTGCGAAAGATAATCCTTATGCAAACTTAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTT  >  USA300TCH1516_ALE/895645‑895911                                                                                                                                             |                                                                                                                               aggTGCGAAAGATAAGCCTTATGAAAACTTAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGt                                                                                                                                <  2:194988/139‑1 (MQ=255)  aGTGCGAAAGATAATCCTTATGCAAACTTAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATGACAATCTAAAGGTGTTGGAGACACCATTAATGAAAAATACAATGGTGCAGtt                                                                                                                               >  1:135197/1‑141 (MQ=255)         aaGATAATCCTTATGCAAACTTAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTg                                                                                                                        >  1:169005/1‑141 (MQ=255)               aTCTTTATGCAAAATGAATTGCTGTTAAAGAAGTACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAaat                                                                                                                  <  2:87455/141‑1 (MQ=255)                        gggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAAt                                                                                                         <  1:133038‑M1/24‑1 (MQ=255)                              tAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATtaaata                                                                                                   <  2:21784/141‑1 (MQ=255)                                         aaaaaaGAAAATCAAGATAATAAGAAAATGAAAGAATTAATTGAAGTATTAAAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCtt                                                                                        <  2:8752/141‑1 (MQ=255)                                         aaaGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCtt                                                                                        <  1:117713/141‑1 (MQ=255)                                         aaaGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATAGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCtt                                                                                        >  2:137167/1‑141 (MQ=255)                                             aagGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTtaat                                                                                    <  1:61143/141‑1 (MQ=255)                                             aagGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTtaat                                                                                    >  1:85483/1‑141 (MQ=255)                                               ggACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTtaataa                                                                                  >  2:183729/1‑141 (MQ=255)                                                gACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATaac                                                                                 >  2:158024/1‑141 (MQ=255)                                                            atAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAggg                                                                     <  1:80468/141‑1 (MQ=255)                                                                aattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAACTAACAATACCCTTTAATAACAACGATTAAGGGGTTa                                                                 >  1:60725‑M1/78‑141 (MQ=255)                                                                         aGCATTAATTAAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTgg                                                        <  2:92950/141‑1 (MQ=255)                                                                                                        tttgtttgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGc                         >  1:64519‑M1/38‑141 (MQ=255)                                                                                                             aGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACAAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGtct                    >  2:72684/1‑139 (MQ=255)                                                                                                                             aaaTACAATGGTGCAGTTCTACCTGAGAACTACAACCAATTAAATAAAAATACCCTTTAATCACCTAGATTCAGAGGTGAGTAGTGTGGTATTGACAATAAAAAAGGGCACAAAAGGGGCAGTATAGGCAAAATATCGAAt    >  2:55588/1‑141 (MQ=255)                                                                                                                              tgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAAtt   <  2:122390‑M1/126‑1 (MQ=255)                                                                                                                               aTACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAGGGGGCAGGATATGCAAATTATCGAAttt  >  1:175419/1‑141 (MQ=255)                                                                                                                                             |                                                                                                                               AAGTGCGAAAGATAATCCTTATGCAAACTTAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTT  >  USA300TCH1516_ALE/895645‑895911

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.