New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP000730	811555 =	42 (1.100)	4 (0.110)	3/262	6.2	10.3%	intergenic (+160/‑117)	USA300HOU_0757/USA300HOU_0758	hypothetical protein/hypothetical protein
	?	CP000730	= 1903651	30 (0.830)					intergenic (+105/+190)	arsC2/gad	arsenate reductase (glutaredoxin)/possible mannosyl‑glycoprotein endo‑beta‑N‑glucosaminidase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ATCGCAATTGAAAATTTCAAATCACGATGAGAGTGGGACAGAAATGATAAAAAAACACTAATGATTTATTATGTAGTGGTTCTTTATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCATTGGTGGGATTATCAATCAATGTGTTTTCTCTAATGCAAAAAGACGCTTTCTATCAAAACAAGATAGAAAG  >  CP000730/1903517‑1903791                                                                                                                                       |                                                                                                                                             aTCGCAATTGAAAATTTCAAATCACGATGAGAGTGGGACAGAAATGATAAAAAAACACTAATGATTTATTATGTAGTGGTTCTTTATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCAt                                                                                                                                        >  2:108559/1‑141 (MQ=255)              aTTTCAAATCACGATGAGAGTGGGACAGAAATGATAAAAAAACACTAATGATTTATTATGTAGTGGTTCTTTATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATAc                                                                                                                           >  1:99142/1‑141 (MQ=255)                     aTCACGATGAGAGTGTAAAAGAAATGATAAAAAAACACTAATGATTTATTATGTAGTGGTTCTTTATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATTCATAAGAAAttagtatttc                                                                                                                    <  2:44748/141‑10 (MQ=18)                     aTCACGATGAGAGTGGGACAGAAATGATAAAAAAACACTAATGATTTATTATGTAGTGGTTCTTTATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTc                                                                                                                    <  2:292244/141‑1 (MQ=255)                              agagTGGGAGAGAAATGATAAAAAAACACTAATGATATAATATGTAGTAGTTCTTTATCATTAGCTATAAAAAAAAAAAAAAAAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGaaaa                                                                                                           <  2:273607/141‑1 (MQ=17)                                      aCAGAAATGATAAAAAAACACTAATGATTTATTATGTAGTGGTTCTTTATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTc                                                                                                   <  1:307558/141‑1 (MQ=255)                                         gAAATGATAAAAAAACACTAATGATTTATTATGTAGTGGTTCTTTATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCttt                                                                                                <  1:340346/141‑1 (MQ=255)                                          aaaTGATAAAAAAACACTAATGATTTATTATGTAGTGGTTCTTTATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTa                                                                                               <  2:155693/141‑1 (MQ=255)                                                   aaaaaCACTAATGATTTATTATGTAGTGGTTCTTTATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTcc                                                                                      >  2:337810/1‑141 (MQ=255)                                                   aaaaaCACTAATGATTTATTATATAGTGGTTCTTTATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTcc                                                                                      >  2:404566/1‑141 (MQ=255)                                                     aaaCACTAATGATTTATTATGTAGTGGTTCTTTATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTg                                                                                    >  1:173986/1‑141 (MQ=255)                                                            aaTGATTTATTATGTAGTGGTTCTTTATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACt                                                                             <  1:235982/141‑1 (MQ=255)                                                                aTTTATTATGTAGTGGTTCTTGATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCa                                                                         <  2:203426/141‑1 (MQ=255)                                                                 tttATTATGTAGTGGTTCTTTATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCAt                                                                        <  1:125679/141‑1 (MQ=255)                                                                    attatGTAGTGGTTCTTTATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATATCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCATtgg                                                                     <  2:269224/141‑1 (MQ=255)                                                                           aGTGGTTCTTTATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCATTGGTGGGAtt                                                              >  1:247425/1‑141 (MQ=255)                                                                           aGTGGTTCTTTATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCATTGGTGGGAtt                                                              <  2:14877/141‑1 (MQ=255)                                                                                             gCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCATTGGTGGGATTATCAATCAATGTGTTTtc                                            >  1:34194/1‑141 (MQ=255)                                                                                                ataACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGACATTCACCTCTCTCAAGAACATCTTTTTCTTCATCGATGTCCTGCCCCACTTGCATTGGTGGGATTATCAATCAATGTGTTTTctct                                         <  1:351283/141‑1 (MQ=21)                                                                                                 taACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCATTGGTGGGATTATCAATCAATGTGTTTTCTCTa                                        >  2:194268/1‑141 (MQ=255)                                                                                                      aaTGTGTACTTGAAAATAGGAATACATGAGTAAAAATCAAGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTAATCGATGTCCTGCCCCACTTGCATTGGTGGGATTATCAATCAATGTGTTTTCTCTAATGCa                                   >  2:59805/1‑141 (MQ=255)                                                                                                         gtgtACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCATTGGTGGGATTATCAATCAATGTGTTTTCTCTAATGCaaaa                                >  1:69412/1‑141 (MQ=255)                                                                                                                  aaaaTAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCATTGGTGGGATTATCAATCAATGTGTTTTCTCTAATGCAAAAAGACGCttt                       >  2:166790/1‑141 (MQ=255)                                                                                                                      tAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCATTGGTGGGATTATCAATCAATGTGTTTTCTCTAATGCAAAAAGACGCTTTCTAt                   <  1:132984/141‑1 (MQ=255)                                                                                                                       aGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCATTGGTGGGATTATCAATCAATGTGTTTTCTCTAATGCAAAGAGACGCTTTCTATc                  >  2:339908/1‑141 (MQ=255)                                                                                                                        ggAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCATTGGTGGGATTATCAATCAATGTGTTTTCTCTAATGCAAAAAGACGCTTTCTATCa                 >  1:350547/1‑141 (MQ=255)                                                                                                                        ggAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCATTGGTGGGATTATCAATCAATGTGTTTTCTCTAATGCAAAAAGACGCTTTCTATCa                 <  2:348408/141‑1 (MQ=255)                                                                                                                        ggAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCATTGGTGGGATTATCAATCAATGTGTTTTCTCTAATGCAAAAAGACGCTTTCTATCa                 >  2:42397/1‑141 (MQ=255)                                                                                                                               aTGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCATTGGTGGGATTATCAATCAATGTGTTTTCTCTAATGCAAAAAGACGCTTTCTATcttgttttg          <  1:108559/141‑9 (MQ=255)                                                                                                                                   gTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCATTGGTGGGATTATCAATCAATGTGTTTTCTCTAATGCAAAAAGACGCTTTCTATcttgttttgatag      <  1:27563/141‑13 (MQ=255)                                                                                                                                       aaCGCAGGCCTAAGACATACTCGTTTCTGAAGAAAAGGTGATTCTTTATCGATGTCCTGCGCCAGTTGGATGGGTGGGATTGTCGATCAAGGTGTGTTCTGTAATGTAAAAAGACGCTGCCTCTCAAGACAGGATAGAAAg  <  2:245677/141‑1 (MQ=255)                                                                                                                                       |                                                                                                                                             ATCGCAATTGAAAATTTCAAATCACGATGAGAGTGGGACAGAAATGATAAAAAAACACTAATGATTTATTATGTAGTGGTTCTTTATCATTAGCTATAACTAATGTGTACTTAAAAATAGGAATACATGAGTAAAACTCATGCATAAGAAATACTAATTTCTAAAGAAAAAGTAATTCTTTATCGATGTCCTGCCCCACTTGCATTGGTGGGATTATCAATCAATGTGTTTTCTCTAATGCAAAAAGACGCTTTCTATCAAAACAAGATAGAAAG  >  CP000730/1903517‑1903791

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.