New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1428794	NA (NA)	6 (0.130)	6/486	3.2	40.1%	noncoding (1/1196 nt)	IS5	repeat region
	?	NC_000913	= 2470044	9 (0.190)					coding (914/1332 nt)	yfdI	serotype‑specific glucosyl transferase YfdI
Rejected: Coverage evenness skew score above cutoff.

CAAAAGGTGTTGTTGAATTAAATAAGTACCATGCAACATACTTCGGTAGTTATCTTTATATGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTAATATATGGAGAATCAAACATACTAACGACTATTACTAACATAAAAGACAATATATTTAAAAACATTAGATTTATATCATTGTTATTATTTTTTATTGCTTCTATTTTTATTAGAAATAATAAAATAAAGGCATCTTTATTTGTAGTATCTC  >  NC_000913/2469810‑2470289                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      cAAAAGGTGTTGTTGAATTAAATAAGTACCATGCAACATACTTCGGTAGTTATCTTTATATGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCAt                                                                                                                                                                                                                                        <  2:49388/250‑1 (MQ=255)                              caggCAACATTCTTCGGTAGTTATCTTTATAAGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAGCCCCAACAGAATTTGGAACGGAATGTTTCGAATCTCATAAAGATGAATCGTTTTCGAATGCACTCTTGTTATGGGTTAGCAAACCAAGCACCATCTTCAAACTGCCATTTGATTATGGTGTGAt                                                                                                                                                                                                          <  2:453762/248‑1 (MQ=255)                                           tCGGTAGTTATCTTTATATGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGtata                                                                                                                                                                                               <  1:99309/249‑1 (MQ=255)                                                aGTTATCTTTATATGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTTGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGaaa                                                                                                                                                                                         >  2:94456/1‑250 (MQ=255)                                                          atatGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTTGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCa                                                                                                                                                                               <  1:94456/250‑1 (MQ=255)                                                                            ataAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGAGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTacac                                                                                                                                                            >  2:299715/1‑251 (MQ=255)                                                                                                gatgatAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTAATATATGGAGAATCaaa                                                                                                                                        <  1:111024/251‑1 (MQ=255)                                                                                                              tGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTAATATATGGAGAATCAAACATACTAACGACt                                                                                                                           >  1:261010/1‑250 (MQ=255)                                                                                                                                                                                                                                      ccttAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTAATATATGGATAATCAAACATTCTAACGACTATTACTAACATAAAAGACAATATATTTAAAAACATTAGATTTATATCAGTGTTATTATTTTTTATTGCTTCTATTTTTATTAGAAATAATAAAATAAAGGCATCTTTATTTGTAGTAgctc  <  2:17695/249‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      CAAAAGGTGTTGTTGAATTAAATAAGTACCATGCAACATACTTCGGTAGTTATCTTTATATGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTAATATATGGAGAATCAAACATACTAACGACTATTACTAACATAAAAGACAATATATTTAAAAACATTAGATTTATATCATTGTTATTATTTTTTATTGCTTCTATTTTTATTAGAAATAATAAAATAAAGGCATCTTTATTTGTAGTATCTC  >  NC_000913/2469810‑2470289

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