New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	895785 =	32 (1.170)	8 (0.330)	4/246	2.3	25.8%	coding (806/822 nt)	metQ_2	Methionine‑binding lipoprotein MetQ
	?	USA300TCH1516_ALE	= 909740	18 (0.750)					intergenic (‑663/‑371)	USA300HOU_RS04500/csbD_1	hypothetical protein/Stress response protein CsbD

GTTGTTTTCAAGGGCAAAAAAAGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAAT  >  USA300TCH1516_ALE/909602‑909860                                                                                                                                           |                                                                                                                         gttgttTTCAAGGGCAAAAAAAGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAtt                                                                                                                        <  2:138112‑M2/141‑3 (MQ=255)  ttgttTTCAAGGGCAAAAAAAGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATATAGTTAAtta                                                                                                                       <  2:291241‑M2/141‑4 (MQ=255)                aaaaaaaGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgctaaataa                                                                                                         >  2:199286‑M2/1‑124 (MQ=255)                aaaaaaaGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgctaaataa                                                                                                         >  2:199284‑M2/1‑124 (MQ=255)                aaaaaaaGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgctaaataa                                                                                                         >  2:62711‑M2/1‑124 (MQ=255)                   aaaaGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGAGGGGGGGGGGGGGGGGTTGATAATCTGGCTAATTATTCCTGCTAAATAAATc                                                                                                      <  2:87326/141‑1 (MQ=255)                                 ttGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCa                                                                                        >  2:180403/1‑141 (MQ=255)                                          aGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGa                                                                               >  1:101077/1‑141 (MQ=255)                                               aaGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTg                                                                          <  1:64674/141‑1 (MQ=255)                                                                      ttGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAgg                                                   <  1:262586/141‑1 (MQ=255)                                                                             ttttATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATAATTACTGCTCCCTCAACCGTATAAGAGGCTTTGTCAGGTGCATTTGTTGGGCTTTTTAATTTAAAATGGTATTtat                                            >  1:197067/1‑141 (MQ=255)                                                                                tataAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATTTg                                         >  1:137538/1‑141 (MQ=255)                                                                                tataAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATg                                         >  2:49555/1‑141 (MQ=255)                                                                                 ataAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggtattgaaaa                                        <  2:211209‑M2/141‑83 (MQ=255)                                                                                           atataACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGaa                              <  2:24957/141‑1 (MQ=255)                                                                                               aaCAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggtattgaaaataaaaaagggcaca                          >  1:229940‑M2/1‑45 (MQ=255)                                                                                               aaCAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggtattgaaaataaaaaagggcaca                          >  1:243024‑M2/1‑45 (MQ=255)                                                                                               aaCAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggtattgaaaataaaaaagggcaca                          >  2:130561‑M2/1‑45 (MQ=255)                                                                                               aaCAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgattaagcggttatttttgtggtattgaaaataaaaaagggcaca                          >  1:219921‑M2/1‑45 (MQ=255)                                                                                                  aaTGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgataaataataccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggtattgaaaataaaaaagggcacaaaa                       >  2:258805‑M2/1‑42 (MQ=255)                                                                                                   aTGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTACACAGGTATTTATATGATATTTAGGAATGAGATGc                      >  1:282300/1‑141 (MQ=255)                                                                                                      ttttGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTAt                   <  2:280173/141‑1 (MQ=255)                                                                                                      ttttGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTAt                   <  2:90040/141‑1 (MQ=255)                                                                                                        ttGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTa                 >  1:14271/1‑141 (MQ=255)                                                                                                            ttGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGa             <  1:120933/141‑1 (MQ=255)                                                                                                            ttGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGa             <  1:148695/141‑1 (MQ=255)                                                                                                            ttGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGa             <  1:128941/141‑1 (MQ=255)                                                                                                               cccTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATaa          >  1:179192/1‑141 (MQ=255)                                                                                                                       cTTCCTTGATAATCTGGCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggtattgaaaataaaaaagggcacaaaaggggcagtatatgcaaattat  <  1:293756‑M2/141‑121 (MQ=255)                                                                                                                       cTTCCTTGATAATCTGGCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggtattgaaaataaaaaagggcacaaaaggggcagtatatgcaaattat  <  1:293745‑M2/141‑121 (MQ=255)                                                                                                                       cTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAAt  >  1:143094/1‑141 (MQ=255)                                                                                                                                           |                                                                                                                         GTTGTTTTCAAGGGCAAAAAAAGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAAT  >  USA300TCH1516_ALE/909602‑909860

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