New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 16585427 (0.740)5 (0.150) 5/226 3.2 15.5% coding (1125/2535 nt) mrcB fused glycosyl transferase and transpeptidase
?NC_000913 = 165878 29 (0.850)coding (1149/2535 nt) mrcB fused glycosyl transferase and transpeptidase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GAAGAGCTAAGCCTCGACCAGCAGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/165747‑165854
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑taatctgTTGCTGTTGCAGCAGACGCAGCACCAGATTACGTCGCTCCAGCGCCAGTTTTGGGTTACGCCACGGGTTGTAGATGGACGCCCCTTTCACCATACCGACTAACAGCGCCTGCTGGTCGAGGCTTAGCTCTTCTACCGGGCGA  <  NC_000913/165878‑165737
                                                                                                                                                                                                                                                          
GAAGAGCTAAGCCTCGACCAGCAGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGTTGCTGTTGCAGCAGACGCATCACCAG                                                                                                                     <  5:36542/135‑1
GAAGAGCTAAGCCTCGACCAGCAGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGTTGCTGTTGCAGCAGACGCATCACCAG                                                                                                                     >  6:36542/1‑135
                                     GTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGTTGCTGTTGCAGCAGACGCAGCACCAGATT                                                                                                                  <  2:28888/101‑1
                                         GGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGTTGCTGTTGCAGCAGACGCAGCACCAGATTACGTC                                                                                                             >  1:18860/1‑102
                                         GGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGTTGCTGTTGCAGCAGACGCAGCACCAGATTACGTC                                                                                                             <  2:18860/102‑1
                                                                                                     TAATCTGTTGCTGTTGCAGCAGACGCAGCACCAGATTACGTCGCTCCAGCGCCAGTTTTGGGTTACGCCACGGGTTGTAGATGGACGCCCCTTTCACCATACCGACTAACAGCGCCTGCTGGTCGAGGCTTAGCTCTTCTACCGGGCGA  >  3:38212/1‑149
                                                                                                                                                                                                                                                          
GAAGAGCTAAGCCTCGACCAGCAGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/165747‑165854
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑taatctgTTGCTGTTGCAGCAGACGCAGCACCAGATTACGTCGCTCCAGCGCCAGTTTTGGGTTACGCCACGGGTTGTAGATGGACGCCCCTTTCACCATACCGACTAACAGCGCCTGCTGGTCGAGGCTTAGCTCTTCTACCGGGCGA  <  NC_000913/165878‑165737

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