New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	895785 =	25 (0.740)	5 (0.170)	5/226	NT	20.5%	coding (806/822 nt)	metQ_2	Methionine‑binding lipoprotein MetQ
	?	USA300TCH1516_ALE	= 909740	17 (0.580)					intergenic (‑663/‑371)	USA300HOU_RS04500/csbD_1	hypothetical protein/Stress response protein CsbD

GGTTGTTTTCAAGGGCAAAAAAAGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAATCCTATATGCATGCGTG  >  USA300TCH1516_ALE/909601‑909876                                                                                                                                            |                                                                                                                                         ggTTGTTTTCAAGGGCAAAAAAAGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAt                                                                                                                                         >  2:28470‑M2/1‑140 (MQ=255)   ttgttTTCAAGGGCAAAAAAAGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAtta                                                                                                                                       <  1:147641‑M2/141‑4 (MQ=255)          aaagGGCAAAAAAAGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgc                                                                                                                                >  2:88151‑M2/2‑131 (MQ=255)                aaaaaaaagGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTTATAATCTGGCTAAttattcctgctaaata                                                                                                                          >  2:336267‑M2/2‑125 (MQ=255)                    aaaaGGGCATAATTTTGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATc                                                                                                                      <  1:330298/141‑1 (MQ=255)                                  ttGAAACAAAGGGCAAGGGTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTACttcttcatgctcaataataccaattaaataaaaa                                                                                                        >  1:305079‑M2/1‑107 (MQ=255)                                                   ggTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTg                                                                                       <  1:243559/141‑1 (MQ=255)                                                    gTATGAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTgg                                                                                      <  2:218141/141‑1 (MQ=255)                                                        gAATAATAAATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCtt                                                                                  <  2:332617/141‑1 (MQ=255)                                                              taaATTATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgattaaggggt                                                                            >  2:147641‑M2/1‑79 (MQ=255)                                                               aaaTTATCTTGATAATTTTATAAATGCATATATAACAATGTTTTGATTGCCTTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTtttaatt                                                                           <  2:191920/141‑1 (MQ=255)                                                                  ttATCTTGAAAATTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTtttaatttaa                                                                        <  1:118791/141‑1 (MQ=255)                                                                           aaaTTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggta                                                               <  1:252320‑M2/141‑76 (MQ=255)                                                                             aTTTTATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATACTCTGGCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggtatt                                                             <  1:44372‑M2/141‑78 (MQ=255)                                                                                ttATAAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTACTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTatat                                                          <  1:259468/141‑1 (MQ=255)                                                                                  ataAATGCATATATAACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGa                                                        >  1:39110/1‑141 (MQ=255)                                                                                              ataACAATGTTTTGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATg                                            <  2:348830/141‑1 (MQ=255)                                                                                                        tttGATTGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTAtt                                  >  2:407834/1‑141 (MQ=255)                                                                                                         ttGATTGCTGTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTa                                 <  2:17496/141‑1 (MQ=255)                                                                                                             ttGCCCTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGa                             >  1:259991/1‑141 (MQ=255)                                                                                                                 ccTTGAACTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATaaa                         <  2:139980/141‑1 (MQ=255)                                                                                                                        cTTCCTTGATAATCTGGCTAAttattcctgctaaataataccaattaaataaaaataccctttaataacaacgattaaggggttatttttgtggtgttgcaaataaaaaggggcacaaaaggggcagtatatgcaaattat                  <  1:112524‑M2/141‑121 (MQ=255)                                                                                                                        cTTCCTTGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAAt                  >  1:195177/1‑141 (MQ=255)                                                                                                                              tGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAATCCtata            >  2:9156/1‑141 (MQ=255)                                                                                                                              tGATAATCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAATCCtata            >  1:242965/1‑141 (MQ=255)                                                                                                                                  aaTCTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAATCCTATatgca        >  1:216373/1‑141 (MQ=255)                                                                                                                                     cTGGCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGCTGGGCTTTTTAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAATCCTATatgcatgc     >  1:136612/1‑141 (MQ=255)                                                                                                                                        gCTAATTATTCCTGCTAAATAAATCATATTAATGGCTCAGTAACGTGATGTTGATGGGATCGCCAATTTAAACAGGTATTTATATGATATTTAGGAATGAGATGCTATTACGGATAAAGTAAAATCCTATATGCATGCGTg  <  1:407834/141‑1 (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.