New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	895785 =	12 (0.410)	8 (0.320)	8/226	NT	27.4%	coding (806/822 nt)	metQ_2	Methionine‑binding lipoprotein MetQ
	?	USA300TCH1516_ALE	= 909740	32 (1.270)					intergenic (‑663/‑371)	USA300HOU_RS04500/csbD_1	hypothetical protein/Stress response protein CsbD

AAGTGCGAAAGATAATCCTTATGCAAACTTAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATA  >  USA300TCH1516_ALE/895645‑895925                                                                                                                                             |                                                                                                                                             ggttgttttcaagggcaaaaaaagggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaat                                                                                                                                              >  1:116564‑M1/141‑141 (MQ=255) ggttgttttcaagggcaaaaaaagggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaat                                                                                                                                              >  1:162472‑M1/141‑141 (MQ=255)   ttgttttcaagggcaaaaaaagggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaatta                                                                                                                                            <  1:60099‑M1/3‑1 (MQ=255)    tgttttcaagggcaaaaaaagggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaattat                                                                                                                                           >  1:20006‑M1/138‑141 (MQ=255)        aaaGATAATCCTTATGCAAACTTAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTcca                                                                                                                                       >  1:371465/1‑140 (MQ=255)          caagggcaaaaaaagggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGc                                                                                                                                     >  1:567648‑M1/132‑141 (MQ=255)           aagggcaaaaaaagggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataagtgcatatgtaacaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCt                                                                                                                                    <  1:251209‑M1/11‑1 (MQ=255)           aagggcaaaaaaagggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCt                                                                                                                                    <  1:598037‑M1/11‑1 (MQ=255)            agggcaaaaaaagggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTa                                                                                                                                   >  1:122199‑M1/130‑141 (MQ=255)            aTAATCCTTATGCAAACTTAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTa                                                                                                                                   >  1:171614/1‑141 (MQ=255)              ggcaaaaaaagggcgtaattttgaaacaaagggcaagggtatgactaataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTaaa                                                                                                                                 <  1:367905‑M1/14‑1 (MQ=255)                       agggcataattttgaaacaaagggcaagggtatgaataataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCaa                                                                                                                        <  1:237336‑M1/23‑1 (MQ=255)                                aTTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATtaaataaa                                                                                                               >  1:132322/1‑141 (MQ=255)                                                gACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATaac                                                                                               >  1:620625/1‑141 (MQ=255)                                                             ataaattatcttgaaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAgggg                                                                                  <  1:275969‑M1/61‑1 (MQ=255)                                                                       aaaGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTtgg                                                                        >  1:283461/1‑140 (MQ=255)                                                                          aaaattttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGt                                                                     >  1:586371‑M1/68‑141 (MQ=255)                                                                               tttataaatgcatatataacaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTGATTTTTGTGGTATTGa                                                                <  1:639523‑M1/79‑1 (MQ=255)                                                                                  tGAAGGATTACAATCTAACGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGaaaa                                                             <  1:231478/141‑1 (MQ=255)                                                                                                aacaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGcaca                                               >  1:104848‑M1/46‑141 (MQ=255)                                                                                                  aaatgttttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACaaa                                             >  1:250043‑M1/44‑141 (MQ=255)                                                                                                   aaGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCAGaaaa                                            <  1:459305/141‑1 (MQ=255)                                                                                                         ttgattgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCa                                      <  1:387677‑M1/105‑1 (MQ=255)                                                                                                           cAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATATTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGt                                    >  1:625854/1‑141 (MQ=255)                                                                                                            attgcccttgaacttccttgataatctggctaaTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAAAGCGCAAAAGGGgtagta                                   <  1:8647‑M1/108‑1 (MQ=255)                                                                                                             aGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTTTTTTTGTGGTATTGAAAATAAAAAAGGGCACAACAGGGGCCGtaa                                  >  1:518583/1‑140 (MQ=255)                                                                                                                  tCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTGTATGCa                             >  1:197762/1‑141 (MQ=255)                                                                                                                  tCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCa                             >  1:539269/1‑141 (MQ=255)                                                                                                                          gAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCg                     >  1:519912/1‑141 (MQ=255)                                                                                                                                             TTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATa  >  1:619476‑M1/1‑141 (MQ=255)                                                                                                                                             |                                                                                                                                             AAGTGCGAAAGATAATCCTTATGCAAACTTAATTGCTGTTAAAGAAGGACATCAAGATGATAAGAAAATCAAAGCATTAATTGAAGTATTACAATCTAAAGATATTCAAGACTTCATTAATGAAAAATACAATGGTGCAGTTATTCCTGCTAAATAATACCAATTAAATAAAAATACCCTTTAATAACAACGATTAAGGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATA  >  USA300TCH1516_ALE/895645‑895925

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.