New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4389174 =123 (0.890)4 (0.030) 4/212 NT 3.6% coding (197/3324 nt) mscM mechanosensitive channel protein, miniconductance
?NC_000913 4389211 = 98 (0.760)coding (160/3324 nt) mscM mechanosensitive channel protein, miniconductance
Rejected: Frequency below/above cutoff threshold.

GGGGGCGTACGCCGCGACGGCCCCCGATAGCAAACAAATCACTCAGGAACTGGAGCAGGCAAAAGCGGCGAAACCCGCACAGCCGGAAGTCGTAGAGGCGCTCCAGTCTGCCTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4389326‑4389174
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agcgcatTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGA  >  NC_000913/4389211‑4389357
                                                                                                                                                                                                                                                                                                            
GGGGGCGTACGCCGCGACGGCCCCCGATAGCAAACAAATCACTCAGGAACTGGAGCAGGCAAAAGCGGCGAAACCCGCACAGCCGGAAGTCGTAGAGGCGCTCCAGTCTGCCTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGC                                                                                                                                                         >  2:1294194/1‑149
   GGCGTACGCCGCGACGGCCCCCGATAGCAAACAAATCACTCAGGAACTGGAGCAGGCAAAAGCGGCGAAACCCGCACAGCCGGAAGTCGTAGAGGCGCTCCAGTCTGCCTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCA                                                                                                                                                      <  4:676740/149‑1
                                            AGGAACTGGAGCAGGCAAAAGCGGCGAAACCCGCACAGCCGGAAGTCGTAGAGGCGCTCCAGTCTGCCTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGA                                                                                                                                                           <  1:1363692/103‑1
                                            AGGAACTGGAGCAGGCAAAAGCGGCGAAACCCGCACAGCCGGAAGTCGTAGAGGCGCTCCAGTCTGCCTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGA                                                                                                                                                           >  4:1190000/1‑103
                                            AGGAACTGGAGCAGGCAAAAGCGGCGAAACCCGCACAGCCGGAAGTCGTAGAGGCGCTCCAGTCTGCCTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGA                                                                                                                                                           <  3:1190000/103‑1
                                            AGGAACTGGAGCAGGCAAAAGCGGCGAAACCCGCACAGCCGGAAGTCGTAGAGGCGCTCCAGTCTGCCTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGA                                                                                                                                                           >  2:1363692/1‑103
                                                          GCAAAAGCGGCGAAACCCGCACAGCCGGAAGTCGTAGAGGCGCTCCAGTCTGCCTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGA                                                                                                                                                           >  3:494508/1‑89
                                                          GCAAAAGCGGCGAAACCCGCACAGCCGGAAGTCGTAGAGGCGCTCCAGTCTGCCTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGA                                                                                                                                                           <  4:494508/89‑1
                                                                                                               CTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCC                                                   <  2:1023182/140‑1
                                                                                                               CTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCC                                                   >  1:1023182/1‑140
                                                                                                                TTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCT                                                                             <  3:103833/113‑1
                                                                                                                TTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCT                                                                             >  4:103833/1‑113
                                                                                                                                                   GCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGA                                                                           >  3:245065/1‑80
                                                                                                                                                   GCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGA                                                                           <  4:245065/80‑1
                                                                                                                                                    CGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAG     <  3:676740/149‑1
                                                                                                                                                     GCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGTCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGT    <  2:281774/149‑1
                                                                                                                                                     GCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGT    <  2:171881/149‑1
                                                                                                                                                     GCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCAT             <  2:762518/140‑1
                                                                                                                                                     GCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCAT             >  1:762518/1‑140
                                                                                                                                                       ATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGA  <  2:294238/149‑1
                                                                                                                                                       ATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGA  <  3:440854/149‑1
                                                                                                                                                       ATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATC            >  2:1383509/1‑139
                                                                                                                                                       ATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATC            <  1:1383509/139‑1
                                                                                                                                                       ATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCC                                                                              <  4:374498/73‑1
                                                                                                                                                       ATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCC                                                                              >  3:374498/1‑73
                                                                                                                                                        TTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAA      <  3:935122/144‑1
                                                                                                                                                        TTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAA      >  4:935122/1‑144
                                                                                                                                                                                                                                                                                                            
GGGGGCGTACGCCGCGACGGCCCCCGATAGCAAACAAATCACTCAGGAACTGGAGCAGGCAAAAGCGGCGAAACCCGCACAGCCGGAAGTCGTAGAGGCGCTCCAGTCTGCCTTAAATGCGCTTGAGGAACGAAAAGGTTCCCTTGAGCGCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4389326‑4389174
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agcgcatTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGA  >  NC_000913/4389211‑4389357

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