New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	895785 =	14 (0.680)	5 (0.280)	3/246	2.3	23.1%	coding (806/822 nt)	metQ_2	Methionine‑binding lipoprotein MetQ
	?	USA300TCH1516_ALE	= 909740	21 (1.160)					intergenic (‑663/‑371)	USA300HOU_RS04500/csbD_1	hypothetical protein/Stress response protein CsbD

GGTTGTTTTCAAGGGCAAAAAAAGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAATCCTATATGCATGCGTG  >  USA300TCH1516_ALE/909601‑909876                                                                                                                                            |                                                                                                                                         ggTTGTTTTCAAGGGCAAAAAAAGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAt                                                                                                                                         >  1:210194‑M2/1‑140 (MQ=255)                gaaagaaaggGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgctaaata                                                                                                                          >  1:54006‑M2/2‑125 (MQ=255)                  aaaaaaGGGCATAATTTTGGAACAAAGGTCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCtaaataaa                                                                                                                        <  2:64797/141‑1 (MQ=255)                         ggCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATAtt                                                                                                                 <  1:82937/141‑1 (MQ=255)                          gCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTa                                                                                                                >  1:105846/1‑141 (MQ=255)                                ttttGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGGTAAATAAGTCATATTAATTGCt                                                                                                          >  1:56324/1‑141 (MQ=255)                                 tttGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAAACTGACGAATTATTCCTGCTAAATAAATCATATTAATGGCTc                                                                                                         >  1:162453/1‑141 (MQ=255)                                    gAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGt                                                                                                      >  1:201091/1‑141 (MQ=255)                                               cAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGtt                                                                                           <  1:165183/141‑1 (MQ=255)                                                               aaaTTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTtttaatt                                                                           <  1:74395/141‑1 (MQ=255)                                                                aaTTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTtttaattt                                                                          <  1:210764/141‑1 (MQ=255)                                                                     tCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACa                                                                     <  2:125002/141‑1 (MQ=255)                                                                                     aaTGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGatat                                                     <  2:201091/141‑1 (MQ=255)                                                                                          atatatAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGTTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAgg                                                <  2:27855/141‑1 (MQ=255)                                                                                            atataACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGaa                                              <  2:162453/141‑1 (MQ=255)                                                                                                     tggtttGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCt                                     <  1:222791/138‑1 (MQ=255)                                                                                                     tGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCt                                     <  1:222585/141‑1 (MQ=255)                                                                                                          tGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgatgaaggggttatttttatagtattgaaaataaaaaagggcacaaaaggggcag                                <  2:137723‑M2/141‑107 (MQ=255)                                                                                                          tGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGATTGATATGTTATTAc                                >  2:46839/1‑141 (MQ=255)                                                                                                            aTTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggtattgaaaataaaaaaggtcacaaaaggggcagta                              <  1:155493‑M2/141‑109 (MQ=255)                                                                                                            aTTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggtattgaaaataaaaaagggcacaaaacgggcagta                              <  2:153540‑M2/141‑109 (MQ=255)                                                                                                                           ccTTGATAATCTGGCTGATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAATCCt               >  2:146658/1‑141 (MQ=255)                                                                                                                            ggtGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAATACta              >  2:222519/3‑141 (MQ=255)                                                                                                                            cTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAATCCta              >  1:89216/1‑141 (MQ=255)                                                                                                                                      tGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAATCCTATATGCATGCg    >  2:175321/1‑141 (MQ=255)                                                                                                                                        gCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggtattgaaaataaaaaagggcacaaaaggggcagtatatgcaaattatcgaattttaaacaaaa  <  1:207238‑M2/141‑137 (MQ=255)                                                                                                                                        gCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggtattgaaaataaaaaagggcacaaaaggggcagtatatgcaaattatcgaattttaaacaaaa  <  1:72218‑M2/141‑137 (MQ=255)                                                                                                                                        gCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAATCCTATATGCATGCGTg  >  1:203931/1‑141 (MQ=255)                                                                                                                                            |                                                                                                                                         GGTTGTTTTCAAGGGCAAAAAAAGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAATCCTATATGCATGCGTG  >  USA300TCH1516_ALE/909601‑909876

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.