New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP000730	= 2850830	25 (0.860)	3 (0.110)	3/268	3.4	11.2%	intergenic (+29/+90)	pcp/USA300HOU_2693	C15 family pyroglutamyl‑peptidase I/hypothetical protein
	?	CP000730	= 2850851	25 (0.950)					intergenic (+50/+69)	pcp/USA300HOU_2693	C15 family pyroglutamyl‑peptidase I/hypothetical protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ATGCCATTAGAAAAGATAGTTGCAGGTTTGACTGCAGCCATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACATTGTTTAGTCTTCTTTTTTTAACTTCCTAAAATTTACAAATGTATAAATAATAAATGCTATATTTAGAATGATGAGTATACTT  >  CP000730/2850706‑2850998                                                                                                                                                  |                                                                                                                                                    aTGCCATTAGAAAAGATAGTTGCAGGTTTGACTGCAGCCATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTa                                                                                                                                                  >  2:139221/1‑149 (MQ=255) aTGCCATTAGAAAAGATAGTTGCAGGTTTGACTGCAGCCATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTa                                                                                                                                                  >  2:139370/1‑149 (MQ=255)              aGATAGTTGCAGGTTTGACTGCAGCCATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGga                                                                                                                                     >  1:307274/1‑149 (MQ=255)                     tGCAGGTTTGACTGCAGCCATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTgg                                                                                                                              <  2:182082/149‑1 (MQ=255)                     tGCAGGTTTGACTGCAGCCATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTgg                                                                                                                              <  1:217455/149‑1 (MQ=255)                     tGCAGGTTTGACTGCAGCCATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTgg                                                                                                                              <  2:168405/149‑1 (MQ=255)                       cAGGTTTGACTGCAGCCATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGta                                                                                                                            >  1:99174/1‑149 (MQ=255)                              gACTGCAGCCATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGtt                                                                                                                     <  2:286870/149‑1 (MQ=255)                                   cAGCCATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACt                                                                                                                <  1:187000/149‑1 (MQ=255)                                    aGCCATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACta                                                                                                               >  1:239303/1‑149 (MQ=255)                                    aGCCATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACta                                                                                                               >  2:100264/1‑149 (MQ=255)                                    aGCCATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACta                                                                                                               <  2:296310/149‑1 (MQ=255)                                                tATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTAtt                                                                                                   >  1:165369/1‑149 (MQ=255)                                                         tgatgaGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTaaa                                                                                          >  1:115714/1‑149 (MQ=255)                                                                      cacCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACATTGTTTAg                                                                             >  1:228253/1‑149 (MQ=255)                                                                             cccTGGGCACAACGGAATAGTTATAAGTGCACAAACCAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACATTGTTTAGTcttctt                                                                      >  2:48851/1‑149 (MQ=255)                                                                                        aCGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACATTGTTTAGTCTTCTTTTTTTAACTTc                                                           >  2:224330/1‑149 (MQ=255)                                                                                              tAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACATTGTTTAGTCTTCTTTTTTTAACTTCCTaaaa                                                     >  1:293648/1‑149 (MQ=255)                                                                                                        tGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACATTGTTTAGTCTTCTTTTTTTAACTTCCTAAAATTTACAAATg                                           <  1:247429/149‑1 (MQ=255)                                                                                                        tGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACATTGTTTAGTCTTCTTTTTTTAACTTCCTAAAATTTACAAATg                                           <  1:298914/149‑1 (MQ=255)                                                                                                                caaaAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACATTGTTTAGTCTTCTTTTTTTAACTTCCTAAAATTTACAAATGTataaata                                   <  1:231429/149‑1 (MQ=255)                                                                                                                 aaaaGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACATTGTTTAGTCTTCTTTTTTTAACTTCCTAAAATTTACAAATGTataaataa                                  >  2:170301/1‑149 (MQ=255)                                                                                                                     gCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTCTTGGATGAGTGGTATACAGTCTTACTATAATGATGTATTTTAGCTAAACGACATTGTTTAGTCTTCTTTTTTTAACTTCCTAAAATTTACAAATGTATAAATAATaaa                              >  2:37965/1‑149 (MQ=255)                                                                                                                                  ttAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACATTGTTTAGTCTTCTTTTTTTAACTTCCTAAAATTTACAAATGTATAAATAATAAATGCTATATTTAGa                 >  1:94621/1‑149 (MQ=255)                                                                                                                                  ttAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACATTGTTTAGTCTTCTTTTTTTAACTTCCTAAAATTTACAAATGTATAAATAATAAATGCTATATTTAGa                 >  1:157997/1‑149 (MQ=255)                                                                                                                                          aGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACATTGTTTAGTCTTCTTTTTTTAACTTCCTAAAATTTACAAATGTATAAATAATAAATGCTATATTTAGAATGATGAg         >  1:217504/1‑149 (MQ=255)                                                                                                                                                 cTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACATTGTTTAGTCTTCTTTTTTTAACTTCCTAAAATTTACAAATGTATAAATAATAAATGCTATATTTAGAATGATGAGTAtactt  <  2:228253/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    ATGCCATTAGAAAAGATAGTTGCAGGTTTGACTGCAGCCATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACATTGTTTAGTCTTCTTTTTTTAACTTCCTAAAATTTACAAATGTATAAATAATAAATGCTATATTTAGAATGATGAGTATACTT  >  CP000730/2850706‑2850998

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