New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3967987	17 (0.500)	5 (0.160)	5/228	3.5	23.4%	coding (72/1104 nt)	wecA	UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
	?	NC_000913	= 3967995	17 (0.540)					coding (80/1104 nt)	wecA	UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
Rejected: Coverage evenness skew score above cutoff.

CTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACTATATTCC  >  NC_000913/3967857‑3968127                                                                                                                                   |                                                                                                                                             cTGACTAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCAGAGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCaaaaa                                                                                                                              <  8:27707/147‑1 (MQ=255)                                    taatTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACg                                                                                         >  5:116654/1‑149 (MQ=255)                                                       gcatTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACct                                                                      >  5:6057/1‑149 (MQ=255)                                                         atTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACctct                                                                    >  2:126352/1‑149 (MQ=255)                                                               aaTTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTgg                                                              >  7:86522/1‑149 (MQ=255)                                                                       gACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTGGGGGGAttt                                                      >  3:7947/1‑149 (MQ=255)                                                                                     aTCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCg                                                                                               >  2:150209/1‑94 (MQ=255)                                                                                     aTCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCg                                                                                               <  1:150209/94‑1 (MQ=255)                                                                                          atcaGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAggg                                        <  6:6057/144‑1 (MQ=255)                                                                                               tATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATAc                                                                        <  5:4929/107‑1 (MQ=255)                                                                                               tATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATAc                                                                        >  6:4929/1‑107 (MQ=255)                                                                                                      ttATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCg                        >  5:122306/1‑148 (MQ=255)                                                                                                         ttCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCg                        >  5:59088/1‑145 (MQ=255)                                                                                                         ttCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCg                        <  6:59088/145‑1 (MQ=255)                                                                                                                  ctgtttctgtttTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTa           <  6:116654/149‑1 (MQ=255)                                                                                                                                   ccGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACTATATTcc  >  2:3846/1‑141 (MQ=255)                                                                                                                                   ccGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACTATATTcc  <  1:3846/141‑1 (MQ=255)                                                                                                                                   |                                                                                                                                             CTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACTATATTCC  >  NC_000913/3967857‑3968127

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