New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4388385 =90 (1.180)3 (0.040) 3/262 NT 3.6% coding (986/3324 nt) mscM mechanosensitive channel protein, miniconductance
?NC_000913 4388485 = 78 (1.090)coding (886/3324 nt) mscM mechanosensitive channel protein, miniconductance
Rejected: Frequency below/above cutoff threshold.

TGGATCTCGTTGCCTCGCAACAGCGTCAGGCTGCCAGCCAGACGTTACAGGTCCGGCAGGCGTTGAATACGCTGCGTGAACAGTCGCAATGGCTGGGATCGTCCAATCTGCTCGGCGAAGCGCTGCGGGCGCAGGTGGCACGGCTGCCGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4388535‑4388385
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgccggaCCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAG  >  NC_000913/4388485‑4388632
                                                                                                                                                                                                                                                                                                           
TGGATCTCGTTGCCTCGCAACAGCGTCAGGCTGCCAGCCAGACGTTACAGGTCCGGCAGGCGTTGAATACGCTGCGTGAACAGTCGCAATGGCTGGGATCGTCCAATCTGCTCGGCGAAGCGCTGCGGGCGCAGGTGGCACGGCTGCCGG                                                                                                                                                       <  2:467957/150‑1
                                                   TCCGGCAGGCGTTGAATACGCTGCGTGAACAGTCGCAATGGCTGGGATCGTCCAATCTGCTCGGCGAAGCGCTGCGGGCGCAGGTGGCACGGCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCC                                                                                                    >  1:1063827/1‑150
                                                      GGCAGGCGTTGAATACGCTGCGTGAACAGTCGCAATGGCTGGGATCGTCGAATCTGCTCGGCGATGCGCTGCGGGCGCAGGTGGCACGGCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCTACGAGCTCCATC                                                                                                 >  1:60619/1‑150
                                                                     CGCTGCGTGAACAGTCGCAATGGCTGGGATAGTCCAATCTGCTCGGCGAAGCGCTGCGGGCGCAGGTGGCACGGCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGCTCCATCCGCTGCGCCTGTTGA                                                                                  >  1:37707/1‑150
                                                                                                                                                 GCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTG      <  1:1097103/150‑1
                                                                                                                                                   CGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCC    <  1:930144/150‑1
                                                                                                                                                     GACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCAGCGCTGTTTTCTGCCAG  <  2:261499/150‑1
                                                                                                                                                                                                                                                                                                           
TGGATCTCGTTGCCTCGCAACAGCGTCAGGCTGCCAGCCAGACGTTACAGGTCCGGCAGGCGTTGAATACGCTGCGTGAACAGTCGCAATGGCTGGGATCGTCCAATCTGCTCGGCGAAGCGCTGCGGGCGCAGGTGGCACGGCTGCCGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4388535‑4388385
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgccggaCCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAG  >  NC_000913/4388485‑4388632

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