New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3802387 =	24 (0.610)	4 (0.130)	4/228	4.1	17.2%	coding (672/1020 nt)	waaO	UDP‑D‑galactose:(glucosyl)lipopolysaccharide‑ alpha‑1,3‑D‑galactosyltransferase
	?	NC_000913	3802413 =	19 (0.600)					coding (646/1020 nt)	waaO	UDP‑D‑galactose:(glucosyl)lipopolysaccharide‑ alpha‑1,3‑D‑galactosyltransferase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AAAATAGTATCGTTGGTTACTGGGTTTATAAAGCTTTCTTTGAGTTGATAATTTAAGCTGAACTGGGTGTTATATTTAATATCAGCGAAAATAAGTTTATCCGCCAGCAACATATTTAACACATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTGCAATAGCTCGTGCAGAAACCTGCTGAGCCGCCCATTGGGCAGTATTAATCAATAAAAAACCGGAGTTAAAGTAACCTT  >  NC_000913/3802271‑3802523                                                                                                                                               |                                                                                                               aaaaTAGTATCGTTGGTTACTGGGTTTATAAAGCTTTCTTTGAGTTGATAATTTAAGCTGAACTGGGTGTTATATTTAATATCAGCGAAAATAAGTTTATCCGCCAGCAACATATTTAACACATCCTGATCAGGATGTGTTAttttttt                                                                                                          <  5:40099/149‑1 (MQ=255)           cGTTGGTTACTGGGTTTATAAAGCTTTCTTTGAGTTGATAATTTAAGCTGAACTGGGTGTTATATTTAATATCAGCGAAAATAAGTTTATCCGCCAGCAACATATTTAACACATCCTAATCAGGATGTGTTATTTTTTTGATTATTTCt                                                                                                >  5:100595/1‑149 (MQ=255)                            ataAAGCTTTCTTTGAGTTGATAATTTAAGCTGAACTGGGTGTTATATTTAATATCAGCGAAAATAAGTTTATCCGCCAGCAACATATTTAACACATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTg                                                                                <  1:88004/148‑1 (MQ=255)                              aaaGCTTTCTTTGAGTTGATAATTTAAGCTGAACTGGGTGTTATATTTAATATCAGCGAAAATAAGTTTATCCGCCAGCAACATATTTAACACATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTGCaa                                                                             >  2:134665/1‑149 (MQ=255)                               aaGCTTTCTTTGAGTTGATAATTTAAGCTGAACTGGGTGTTATATTTAATATCAGCGAAAATAAGTTTATCCGCCAGCAACATATTTAACACATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTGCaa                                                                             <  5:36366/148‑1 (MQ=255)                                aGCTTTCTTTGAGTTGATAATTTAAGCTGAACTGGGTGTTATATTTAATATCAGCGAAAATAAGTTTATCCGCCAGCAACATATTTAACTCATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTGCaa                                                                             >  4:144122/1‑147 (MQ=255)                                aGCTTTCTTTGAGTTGATAATTTAAGCTGAACTGGGTGTTATATTTAATATCAGCGAAAATAAGTTTATCCGCCAGCAACATATTTAACACATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTGCAATa                                                                           <  1:116292/149‑1 (MQ=255)                                                         gCTGAACTGGGTGTTATATTTAATATCAGCGAAAATAAGTTTATCCGCCAGCAACATATTTAACACATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTGCAATAGCTCGTGCAGAAACCTGCTGAgccg                                                  >  7:69904/1‑149 (MQ=255)                                                         gCTGAACTGGGTGTTATATTTAATATCAGCGAAAATAAGTTTATCCGCCAGCAACATATTTAACACATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTGCAATAGCTCGTGCAGAAACCTGCTGAgccg                                                  >  4:144133/1‑149 (MQ=255)                                                                         tatTTAATATCAGCGAAAATAAGTTTATCCGCCAGCAACATATTTAACACATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTGCAATAGCTCGTGCAGAAACCTGCTGAGCCGCCCATTGGGCAGTAt                                   <  1:134665/148‑1 (MQ=255)                                                                          atatAATATCAGCGAAAATAAGTTTATCAGCCAGCAACATATTTAACACATCCTGATCAGGATGTGTTATATAATTGATTATTTCTGGCTCATTTAGCATTGAAATAGCTCGTGCAGAAAcctcct                                                        >  3:121687/4‑126 (MQ=255)                                                                           tttAATATCAGCGAAAATAAGTTTATCCGCCAGCAACATATTTAACACATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTGCAATAGCTCGTGCAGAAACCTGCTGAGCCGCCCATTGGGCAGTATTaa                                >  5:161678/1‑149 (MQ=255)                                                                             tAATATCAGCGAAAATAAGTTTATCCGCCAGCAACATATTTAACACATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTGCAATAGCTCGTGCAGAAACCTGCTGAGCCGCCCATTGGGCAGTATTAAt                               <  6:161678/148‑1 (MQ=255)                                                                             tAATATCAGCGAAAATAAGTTTATCCGCCAGCAACATATTTAACACATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTGCAATAGCTCGTGCAGAAACCTGCTGAGCCGCCCATTGGGCAGTATTAAt                               <  8:69904/148‑1 (MQ=255)                                                                                  tCAGCGAAAATAAGTTTATCCGCCAGCAACATATTTAACACATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTGCAATAGCTCGTGCAGAAACCTGCTGAGCCGCCCATTGGGCAGTATTAATCAATa                          >  1:86055/1‑148 (MQ=255)                                                                                      cGAAAATAAGTTTATCCGCCAGCAACATATTTAACACATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTGCAATAGCTCGTGCAGAAACCTGCTGAGCCGCCCATTGGGCAGTATTAATCAATaaaaaa                     >  7:21005/1‑149 (MQ=255)                                                                                                  tATCCGCCAGCAACATATTTAACTCATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTGCAATAGCTCGTGCAGAAACCTGCTGAGCCGCCCATTGGGCAGTATTAATCAATAAAAAACCGGAGTTaaa          <  3:144122/148‑1 (MQ=255)                                                                                                         cAGCAACATATTTAACACATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTGCAATAGCTCGTGCAGAAACCTGCTGAGCCGCCCATTGGGCAGTATTAATCAATAAAAAACCGGAGTTAAAGTAACCtt  <  2:86055/149‑1 (MQ=255)                                                                                                            cAACATATTTAACACATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTGCAATAGCTCGTGCAGAAACCTGCTGAGCCGCCCATTGGGCAGTATTAATCAATAAAAAACCGGAGTTAAAGTAACCtt  >  7:88007/1‑146 (MQ=255)                                                                                                                                               |                                                                                                               AAAATAGTATCGTTGGTTACTGGGTTTATAAAGCTTTCTTTGAGTTGATAATTTAAGCTGAACTGGGTGTTATATTTAATATCAGCGAAAATAAGTTTATCCGCCAGCAACATATTTAACACATCCTGATCAGGATGTGTTATTTTTTTGATTATTTCTGGCTCATTTAGCATTGCAATAGCTCGTGCAGAAACCTGCTGAGCCGCCCATTGGGCAGTATTAATCAATAAAAAACCGGAGTTAAAGTAACCTT  >  NC_000913/3802271‑3802523

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.