New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 445664740 (0.900)4 (0.090) 4/272 5.9 9.4% coding (862/1374 nt) mpl UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
?NC_000913_3_bme_pgi = 4456657 39 (0.920)coding (872/1374 nt) mpl UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_bme_pgi/4456493‑4456647
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atggcgaGCCGCTGCAATCGCCATCAGGCCATTGTGCATATTATGTTCGCCTACCAGCGACCATTTCACTTCGCCCACTTTTTCGCCATCCAGCAAAACTTCCCATTCGGAAGCATCGGTGGTCAGCTTTTTCGCCTGCCAGTGACCCTGCT  <  NC_000913_3_bme_pgi/4456657‑4456513
                                                                                                                                                                                                                                                                                                            
CGAGCAGGAGCTGGTGGGTGGGGAGGGTTAATGGCTGCGCAATGATCGGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGA                                                                                                                                                         <  1:721227/149‑1
   GCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGG                                                                                                                                                      <  1:376795/149‑1
                                  CAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGAGCCTCTGCAATCGCCATCAGGCCATTGT                                                                                                                       <  1:111865/149‑1
                                                                                             AGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGAGCCGCTGCAATCGCCAT                                                                                                                                  >  1:554094/1‑79
                                                                                             AGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGAGCCGCTGCAATCGCCAT                                                                                                                                  <  2:554094/79‑1
                                                                                                                                                  aGATGGCGAGCCGCTGCAATCGCCATCAGGCCATTGTGCATATTATGTTCGCCTACCAGCGACCATTTCACTTCGCCCACTTTTTCGCCATCCAGCAAAACTTCCCATTCGGAAGCATCGGTGGTCAGCTTTTTCGCCTGCCAGTGACC       <  2:174852/148‑1
                                                                                                                                                       cCGAGCCGCCGCAATCGCCACCAGGCCATTGTGCATATTATGGTCGCCTACCAGCGACCATTCCACTTCGCCCACTTTTTCGCCATCCAGCAAAACTTCCCATTCGGAAGCATCGGTGGTCAGCTTTTTCGCCTGCCAGTGACCCTGCT  <  1:234513/148‑1
                                                                                                                                                                                                                                                                                                            
CGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCGAAAAAGTGGGCGAAGTGAAATGGTCGCTGGTAGGCGAACATAATATGCACAATGGCCTGATGGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_bme_pgi/4456493‑4456647
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atggcgaGCCGCTGCAATCGCCATCAGGCCATTGTGCATATTATGTTCGCCTACCAGCGACCATTTCACTTCGCCCACTTTTTCGCCATCCAGCAAAACTTCCCATTCGGAAGCATCGGTGGTCAGCTTTTTCGCCTGCCAGTGACCCTGCT  <  NC_000913_3_bme_pgi/4456657‑4456513

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