New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3968663 =36 (0.960)6 (0.170) 6/220 3.3 16.0% coding (748/1104 nt) wecA UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
?NC_000913 3968701 = 30 (0.870)coding (786/1104 nt) wecA UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTACATAATCGCCACCATATCCATTAGCGGAATGGCGATTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3968704‑3968663
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggcgattatGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGC  >  NC_000913/3968701‑3968846
                                                                                                                                                                                            
GTACATAATCGCCACCATATCCATTAGCGGAATGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGA                                         >  6:83672/1‑149
 TACATAATCGCCACCATATCCATTAGCGGAATGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGAT                                        >  1:214472/1‑149
     TAATCGCCACCATATCCATTAGCGGAATGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGA                                                                                   <  5:47920/102‑1
     TAATCGCCACCATATCCATTAGCGGAATGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGA                                                                                   >  6:47920/1‑102
                       TTAGCGGAATGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTT                                 >  7:56200/1‑134
                       TTAGCGGAATGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTT                                 <  8:56200/134‑1
                                        ATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGGCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGC  <  5:83672/148‑1
                                                                                                                                                                                            
GTACATAATCGCCACCATATCCATTAGCGGAATGGCGATTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3968704‑3968663
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggcgattatGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGC  >  NC_000913/3968701‑3968846

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