New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4389174 =124 (1.010)4 (0.030) 4/238 NT 3.3% coding (197/3324 nt) mscM mechanosensitive channel protein, miniconductance
?NC_000913 4389211 = 119 (1.030)coding (160/3324 nt) mscM mechanosensitive channel protein, miniconductance
Rejected: Frequency below/above cutoff threshold.

ATAGCAAACAAATCACTCAGGAACTGGAGCAGGCAAAAGCGGCGAAACCCGCACAGCCGGAAGTCGTAGAGGCGCTCCAGTCTGCCTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4389300‑4389174
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agcgcatTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGAT  >  NC_000913/4389211‑4389358
                                                                                                                                                                                                                                                                                   
ATAGCAAACAAATCACTCAGGAACTGGAGCAGGCAAAAGCGGCGAAACCCGCACAGCCGGAAGTCGTAGAGGCGCTCCAGTCTGCCTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCA                                                                                                                                                       <  1:716381/126‑1
ATAGCAAACAAATCACTCAGGAACTGGAGCAGGCAAAAGCGGCGAAACCCGCACAGCCGGAAGTCGTAGAGGCGCTCCAGTCTGCCTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCA                                                                                                                                                       >  2:716381/1‑126
                                                                                     CTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTT                                                                 >  4:861595/1‑127
                                                                                     CTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTT                                                                 <  3:861595/127‑1
                                                                                      TTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGG                                                       <  4:320342/136‑1
                                                                                      TTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGG                                                       >  3:320342/1‑136
                                                                                                                        AGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAA        <  4:659354/149‑1
                                                                                                                        AGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGT                                                                                                          <  2:847276/51‑1
                                                                                                                        AGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGT                                                                                                          >  1:847276/1‑51
                                                                                                                         GCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAA       <  1:852895/149‑1
                                                                                                                         GCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAA       <  4:59546/149‑1
                                                                                                                         GCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCC                                  <  3:795634/122‑1
                                                                                                                         GCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCC                                  >  4:795634/1‑122
                                                                                                                         GCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCC                                   <  3:738209/121‑1
                                                                                                                         GCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCC                                   >  4:738209/1‑121
                                                                                                                           GCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGT     <  2:461496/149‑1
                                                                                                                           GCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGT     <  1:613046/149‑1
                                                                                                                            CATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTG    <  1:378890/149‑1
                                                                                                                             ATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGA   <  1:226614/149‑1
                                                                                                                              TTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGAT  <  2:611406/149‑1
                                                                                                                              TTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGAT  >  1:774794/1‑149
                                                                                                                              TTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGAT  >  1:456471/1‑149
                                                                                                                                                                                                                                                                                   
ATAGCAAACAAATCACTCAGGAACTGGAGCAGGCAAAAGCGGCGAAACCCGCACAGCCGGAAGTCGTAGAGGCGCTCCAGTCTGCCTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4389300‑4389174
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agcgcatTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGAT  >  NC_000913/4389211‑4389358

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