New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3398963 =84 (1.920)4 (0.090) 3/266 NT 4.7% coding (1014/1104 nt) mreC cell wall structural complex MreBCD transmembrane component MreC
?NC_000913 3399001 = 82 (1.940)coding (976/1104 nt) mreC cell wall structural complex MreBCD transmembrane component MreC
Rejected: Frequency below/above cutoff threshold.

TGCTAATGAACGTCTGATGCAGATGATGCCGCAGGTATTGCCTTCGCCAGACGCGATGGGGCCAAAGTTACCTGAACCGGCAACGGGGATCGCTCAGCCGACTCCGCAGCAACCGGCGACAGGAAATGCAGCTACTGCGCCTGCTGCGCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3399113‑3398963
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccgGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGACGTTCATTA  >  NC_000913/3399001‑3399110
                                                                                                                                                                                                                                                                     
TGCTAATGAACGTCTGATGCAGATGCTGCCGCAGGTATTGCCTTCGCCAGACGCGATGGGGCCAAAGTTACCTGAACCGGCAACGGGGATCGCTCTGCCGACTCCGCAGCAACCGGCGAGAGGAAATGCAGCTACTGCGCCTGCTGCGC                                                                                                                  <  1:477256/149‑1
TGCTAATGAACGTCTGATGCAGATGATGCCGCAGGTATTGCCTTCGCCAGACGCGATGGGGCCAAAGTTACCTGAACCGGCAACGGGGATCGCTCAGCCGACTCCGCAGCAACCGGCGACAGGAAATGCAGCTACTGCGCCTGCTGCGC                                                                                                                  <  2:500083/149‑1
TGCTAATGAACGTCTGATGCAGATGATGCCGCAGGTATTGCCTTCGCCAGACGCGATGGGGCCAAAGTTACCTGAACCGGCAACGGGGATCGCTCAGCCGACTCCGCAGCAACCGGCGACAGGAAATGCAGCTACTGCGCCTGCTGCGC                                                                                                                  <  2:593786/149‑1
TGCTAATGAACGTCTGATGCAGATGATGCCGCAGGTATTGCCTTCGCCAGACGCGATGGGGCCAAAGTTACCTGAACCGGCAACGGGGATCGCTCAGCCGACTCCGCAGCAACCGGCGACAGGAAATGCAGCTACCGCGCCTGCTGCGC                                                                                                                  >  2:54799/1‑149
TGCTAATGAACGGCTGATGTAGATGATGCCGCAGGTATTGCTTTCGCCATACGCGATGGGGCCAAAGTTACCTGAACCGGCAACGGGGATCGCTCAGCCGACTCCGCAGCACCCGGCGACAGGAAATGCAGCTACTGCGCCTGCTGCGC                                                                                                                  <  2:334101/149‑1
 gcgAATGAACGTCTGATTCTGATGTTGCCGCTTGTATTGCCTTCGCCAGACGCCATGGGGCCACGGTGACCTGAACCGGCAACGGGGTTCCCTCAGCCGACTCCGCATCACCCGGCGACAGGAAATGCAGCTACTGCGCCTGCTGCGCC                                                                                                                 <  2:117058/146‑1
                                                                                         TCGCTCAGCCGACTCCGCAGCAACCGGCGACAGGAAATGCAGCTACTGCGCCTGCTGCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCA                         >  2:763070/1‑149
                                                                                                              AACCGGCGACAGGAAATGCAGCTACTGCGCCTGCTGCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGACGCTCAT    >  1:361524/1‑149
                                                                                                              AACCGGCGACAGGAAATGCAGCTACTGCGCCTGCTGCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGGCTGGCGAAGGCCATACCTCCGGCATCATCTGCATCAGACGTTCAT    >  1:724381/1‑149
                                                                                                                CCGGCGACAGGAAATGCAGCTACTGCGCCTGCTGCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGACGTTCATTA  >  1:664482/1‑149
                                                                                                                                                                                                                                                                     
TGCTAATGAACGTCTGATGCAGATGATGCCGCAGGTATTGCCTTCGCCAGACGCGATGGGGCCAAAGTTACCTGAACCGGCAACGGGGATCGCTCAGCCGACTCCGCAGCAACCGGCGACAGGAAATGCAGCTACTGCGCCTGCTGCGCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3399113‑3398963
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccgGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGACGTTCATTA  >  NC_000913/3399001‑3399110

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.