Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	USA300TCH1516_ALE	2,372,484	+ATAC	21.8%	coding (1292/1335 nt)	pbuG ←	Guanine/hypoxanthine permease PbuG

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	= 2372484	28 (0.870)	8 (0.240) +ATAC	4/232	NT	21.8%	coding (1292/1335 nt)	pbuG	Guanine/hypoxanthine permease PbuG
	?	USA300TCH1516_ALE	2372485 =	29 (0.840)					coding (1291/1335 nt)	pbuG	Guanine/hypoxanthine permease PbuG

TAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  USA300TCH1516_ALE/2372346‑2372484 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTG  >  USA300TCH1516_ALE/2372485‑2372619                                                                                                                                            ||||                                                                                                                                        TAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAAT                                                                                                                                           <  2:130254/141‑1   ATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATAC                                                                                                                                         >  2:180256/1‑141     ATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGACCAAATATGAAGTACAGAATGAATAGTACCATTAATACAT                                                                                                                                       >  2:211303/1‑141                   TTTCGATTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCAGTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATACATTATCGG                                                                                                                         <  2:239939/141‑1                      ctCTTCTTAGTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATACATTATCGGATG                                                                                                                      <  2:144754/139‑1                                                                           AGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATC                                                                 <  1:379548/141‑1                                                                                      GTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATACATTATCGGATGTACTTCTTTATGTATTTTAGAGATAAGCATTGTAATTTGATAGAAAATAAATACACACGCAATA                                                      >  1:178605/1‑141                                                                                                                    acTACAGAATGAATAGTACCATTAATACATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCATTTGTAATTGAATATGATAAAGGCATC                        >  2:422524/3‑141                                                                                                                                          TAATACATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTG  >  1:354045/1‑141                                                                                                                                            ||||                                                                                                                                        TAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  USA300TCH1516_ALE/2372346‑2372484 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTG  >  USA300TCH1516_ALE/2372485‑2372619

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.