New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3968570 =112 (1.320)5 (0.060) 5/522 NT 4.3% coding (655/1104 nt) wecA UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
?NC_000913 = 4499511 NA (NA)noncoding (1331/1331 nt) IS2 repeat region
Rejected: Frequency below/above cutoff threshold.

GCCAGCACGCCAATGGAAGCGAGCAGTGCTGCGGCAAGGGTAATCAGCACAAACGCCTGACGGGAAGTAAACCCGGCACGCATGATCAAATGGTGAATATGCTGACGGTCAGGAGAGAATGGGCTCATGCCTTTACGCAGGCGACGGTACATAATCGCCACCATATCCATTAGCGGAATGGCGATTATCCACAAAGCGGTAACCGGGCTGATGGGATGGGTTTTGCCCTGGGTCGTTTCGAGCAGGATCCAGATAACGGTAAAACCAATCAGCGTACTGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3968850‑3968570
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCACCGCCGTGTTCTTCGGTTCAAGTCCCAACATCCGGGCGAA  <  NC_000913/4499511‑4499225
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
GCCAGCACGCCAATGGAAGCGAGCAGTGCTGCGGCAAGGGTAATCAGCACAAACGCCTGACGGGAAGTAAACCCGGCACGCATGATCAAATGGTGAATATGCTGACGGTCAGGAGAGAATGGGCTCATGCCTTTACGCAGGCGACGGTACATAATCGCCACCATATCCATTAGCGGAATGGCGATTATCCACAAAGCGGTAACCGGGCTGATGGGATGGGTTTTGCCCTGGGTCGTTTCGAGCAGGATCCAGATAACGGTAAAACCAATCAGCGTACTGCCTGGATTTG                                                                                                                                                                                                                                                                                         >  2:350971/1‑289
                                    AGGGTAATCAGCACAAACGCCTGACGGGAAGTAAACCCGGCACGCATGATCAAATGGTGAATATGCTGACGGTCAGGAGAGAATGGGCTCATGCCTTTACGCAGGCGACGGTACATAATCGCCACCATATCCATTAGCGGAATGGCGATTATCCACAAAGCGGTAACCGGGCTGATGGGATGGGTTTTGCCCTGGGTCGTTTCGAGCAGGATCCAGATAACGGTAAAACCAATCAGCGTACTGCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACCTAACCC                                                                                                                                                                                                                                                     >  1:594328/1‑289
                                                                                        AATGGTGAATATGCTGACGGTCAGGAGAGAATGGGCTCATGCCTTTACGCAGGCGACGGTACATAATCGCCACCATATCCATTAGCGGAATGGCGATTATCCACAAAGCGGTAACCGGGCTGATGGGATGGGTTTTGCCCTGGGTCGTTTCGAGCAGGATCCAGATAACGGTAAAACCAATCAGCGTACTGCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCA                                                                                                                                                                                                 <  1:62545/289‑1
                                                                                                                                                                       CATTAGCGGAATGGCGATTATCCACAAAGCGGTAACCGGGCTGATGGGATGGGTTTTGCCCTGGGTCGTTTCGAGCAGGATCCAGATAACGGTAAAACCAATCAGCGTACTGCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCA                                                                                                                  <  1:350971/289‑1
                                                                                                                                                                                                                        TGGGTTTTGCCCTGGGTCGTTTCTAGCAGGATCCAGATAACGGTAAAACCAATCAGCGTACTGCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTCCAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTA                                                                  <  2:594328/288‑1
                                                                                                                                                                                                                                                                                       CCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCACCGCCGTGTTCTTCGGTTCAAGTCCCAACATCCGGGCGAA  <  2:273880/289‑1
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
GCCAGCACGCCAATGGAAGCGAGCAGTGCTGCGGCAAGGGTAATCAGCACAAACGCCTGACGGGAAGTAAACCCGGCACGCATGATCAAATGGTGAATATGCTGACGGTCAGGAGAGAATGGGCTCATGCCTTTACGCAGGCGACGGTACATAATCGCCACCATATCCATTAGCGGAATGGCGATTATCCACAAAGCGGTAACCGGGCTGATGGGATGGGTTTTGCCCTGGGTCGTTTCGAGCAGGATCCAGATAACGGTAAAACCAATCAGCGTACTGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3968850‑3968570
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCACCGCCGTGTTCTTCGGTTCAAGTCCCAACATCCGGGCGAA  <  NC_000913/4499511‑4499225

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ 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.