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.

CTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACTATATTCCGCATGCATCTCTCTA  >  NC_000913/3967857‑3968142                                                                                                                                           |                                                                                                                                                    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)                                                                                                                                          gTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACTATATTCCGCATGCATCTCTCTa  <  6:122306/149‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                    CTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACTATATTCCGCATGCATCTCTCTA  >  NC_000913/3967857‑3968142

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.