New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	895785 =	17 (0.610)	9 (0.370)	6/262	1.8	28.1%	coding (806/822 nt)	metQ_2	Methionine‑binding lipoprotein MetQ
	?	USA300TCH1516_ALE	= 909740	31 (1.260)					intergenic (‑663/‑371)	USA300HOU_RS04500/csbD_1	hypothetical protein/Stress response protein CsbD

CTGTAGAAAGTGCGAAAGATAATCCTTATGCAAACTTAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTT  >  USA300TCH1516_ALE/895638‑895912                                                                                                                                                    |                                                                                                                                ggagtgtggttgttttcaagggcaaaaaaagggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattggccttgtacttccttgataatctggctaatt                                                                                                                                <  1:247613‑M1/2‑1 (MQ=255) ggagtgtggttgttttcaagggcaaaaaaagggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaatt                                                                                                                                <  1:215521‑M1/2‑1 (MQ=255) cTGTAGAAAGTGCGAAAGATAATCCTTATGCAAACTTAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGtt                                                                                                                                >  1:196909/1‑149 (MQ=255)  gagtgtggttgttttcaagggcaaaaaaagggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaatta                                                                                                                               <  1:112342‑M1/3‑1 (MQ=255)  gagtgtggttgttttcaagggcaaaaaaagggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaatta                                                                                                                               <  2:250084‑M1/3‑1 (MQ=255)   agtgtggttgttttcaagggcaaaaaaagggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaattat                                                                                                                              <  2:151301‑M1/4‑1 (MQ=255)   agtgtggttgttttcaagggcaaaaaaagggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaattat                                                                                                                              <  1:131909‑M1/4‑1 (MQ=255)   agtgtggttgttttcaagggcaaaaaaagggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaattat                                                                                                                              <  1:15766‑M1/4‑1 (MQ=255)   agtgtggttgttttcaagggcaaaaaaagggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaattat                                                                                                                              <  2:84347‑M1/4‑1 (MQ=255)   agtgtggttgttttcaagggcaaaaaaagggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaattat                                                                                                                              <  1:10015‑M1/4‑1 (MQ=255)       aaaGTGCGAAAGATAATCCTTATGCAAACTTAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCt                                                                                                                          <  2:234799/149‑1 (MQ=255)         agggcaaaaGATAAACCATATGCAAACTTAGTTGCTGTTCAAGAAGGACATGAAGATAATAAGAATATCAAAGCATTAACTGAAGTATTACAGTCTAAAGATATTCAAGACAAAAAAAAAAAAAAATACAATGGTGCAGTTATTCCTGc                                                                                                                        <  2:266040/143‑1 (MQ=255)            gCGAAAGATAATCCTTATGCAAACTTAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTaa                                                                                                                     <  1:196826/149‑1 (MQ=255)                 aGATAATCCTTATGCAAACTTAATTGCTGTTAAAGAAGGACCTCACGATGATAAGCAAATCAAAGCCTTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAaataat                                                                                                                <  1:55461/149‑1 (MQ=255)                   aTAATCCTTATGCAAACTTAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATAc                                                                                                              <  2:92365/149‑1 (MQ=255)                               cAAACTTAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATtaaataaa                                                                                                  >  1:55223/1‑149 (MQ=255)                                     tAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATAcc                                                                                            <  2:86844/149‑1 (MQ=255)                                        ttGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCtt                                                                                         <  1:195425/149‑1 (MQ=255)                                          tgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTta                                                                                       >  1:18632‑M1/107‑149 (MQ=255)                                           ctactacagaagGACATAAAGATGATACGAAAATCAAAGCATTACTTGAACTATTAAAATCTAAAGATATTAACGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTtaa                                                                                      <  1:139979/142‑1 (MQ=255)                                                       aagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTa                                                                          >  2:115197‑M1/94‑149 (MQ=255)                                                       aagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTa                                                                          >  2:115218‑M1/94‑149 (MQ=255)                                                               gatgatAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTa                                                                  >  2:232154/1‑149 (MQ=255)                                                                       gAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAGGGGGTTATTTTtgtg                                                          <  1:221169/149‑1 (MQ=255)                                                                            tCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTAtt                                                     <  1:23266/149‑1 (MQ=255)                                                                             cAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTg                                                    >  2:176536/1‑149 (MQ=255)                                                                                                        tAAAGATATTCAAGACTTCATTAATGAAAAATATAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCa                         >  2:62787/1‑149 (MQ=255)                                                                                                        tAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCa                         >  1:176123/1‑149 (MQ=255)                                                                                                        acaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCa                         >  2:162411‑M1/45‑149 (MQ=255)                                                                                                        acaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCa                         <  1:170035‑M1/105‑1 (MQ=255)                                                                                                           atgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGta                      >  1:222223‑M1/42‑149 (MQ=255)                                                                                                           atgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGta                      >  1:189004‑M1/42‑149 (MQ=255)                                                                                                           atgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGta                      >  1:101802‑M1/42‑149 (MQ=255)                                                                                                                        ccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAAAAACGATTAAGGGGTTATTTTTGCGGTATTGAAATTCTAAACGCGCACAACACGGGCCGTATATGCAAATTatt         <  1:196202‑M1/121‑2 (MQ=255)                                                                                                                               aTGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAAtttt  >  2:160761/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                CTGTAGAAAGTGCGAAAGATAATCCTTATGCAAACTTAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTT  >  USA300TCH1516_ALE/895638‑895912

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