New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4377011 =56 (1.220)4 (0.090) 4/228 2.5 7.5% coding (2177/3324 nt) mscM miniconductance mechanosensitive channel MscM
?NZ_CP009273 4377031 = 46 (1.070)coding (2157/3324 nt) mscM miniconductance mechanosensitive channel MscM
Rejected: Frequency below/above cutoff threshold.

CGTTGGTTGCCATTCTGGCGTCGGCGGTGGGTTATCTGGCAACGGCACAGGCGCTGTTAGCGAGGCTTGAAACCTCGGTTGCCATCTGGTTCCTGCTACTGGTGGTTTATCACGTTATCCGCCGCTGGATGCTGATCCAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/4377152‑4377011
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atccagcgGCGGATAACGTGATAAACCACCAGTAGCAGGAACCAGATGGCAACCGAGGTTTCAAGCCTCGCTAACAGCGCCTGTGCCGTTGCCAGATAACCCACCGCCGACGCCAGAATGGCAACCAACGGCGCGCCAATCATC  >  NZ_CP009273/4377031‑4377166
                                                                                                                                                                                                                                                                                      
CGTTGGTTGCCATTCTGGCGTCGGCGGTGGGTTATCTGGCAACGGCACAGGCGCTGTTAGCGAGGCTTGAAACCTCGGTTGCCATCTGGTTCCTGCTACTGGTGGTTTATCACGTTATCCGCCGCTGGATGCTGATCCA                                                                                                                                             <  2:31338/139‑1
                                       CAACGGCACAGGCGCTGTTAGCGAGGCTTGAAACCTCGGTTGCCATCTGGTTCCTGCTACTGGTGGTTTATCACGTTATCCGCCGCTGGATGCTGATCCAGCGGCGGATA                                                                                                                                   >  1:748938/1‑110
                                       CAACGGCACAGGCGCTGTTAGCGAGGCTTGAAACCTCGGTTGCCATCTGGTTCCTGCTACTGGTGGTTTATCACGTTATCCGCCGCTGGATGCTGATCCAGCGGCGGATA                                                                                                                                   <  2:748938/110‑1
                                                                  TTGAAACCTCGGTTGCCATCTGGTTCCTGCTACTGGTGGTTTATCACGTTATCCGCCGCTGGATGCTGATCCAGCGGCGGATAACGTGATAAACCAC                                                                                                                     <  1:826324/97‑1
                                                                  TTGAAACCTCGGTTGCCATCTGGTTCCTGCTACTGGTGGTTTATCACGTTATCCGCCGCTGGATGCTGATCCAGCGGCGGATAACGTGATAAACCAC                                                                                                                     >  2:826324/1‑97
                                                                                                                                         CAGCGGCGGATAACGTGATAAACCACCAGTAGCAGGAACCAGATGGCAACCGAGGTTTCAAGCCTCGCTAACAGCGCCTGTGCCGTTGCCAGATAACCCACCGCCGACGCCAGAATGGCAACCAACGGCGCGCCAATCA    >  2:632621/1‑139
                                                                                                                                         CAGCGGCGGATAACGTGATAAACCACCAGTAGCAGGAACCAGATGGCAACCGAGGTTTCAAGCCTCGCTAACAGCGCCTGTGCCGTTGCCAGATAACCCACCGCCGACGCCAGAATGGCAACCAACGGCGCGCCAATCA    >  2:83905/1‑139
                                                                                                                                          AGCGGCGGATAACGTGATAAACCACCAGTAGCAGGAACCAGATGGCAACCGAGGTTTCAAGCCTCGCTAACAGCGCCTGTGCCGTTGCCAGATAACCCACCGCCGACGCCAGAATGGCAACCAACGGCGCGCCAATCAT   >  2:873431/1‑139
                                                                                                                                           GCGGCGGATAACGTGATAAACCACCAGTAGCAGGAACCAGATGGCAACCGAGGTTTCAAGCCTCGCTAACAGCGCCTGTGCCGTTGCCAGATAACCCACCGCCGACGCCAGAATGGCAACCAACGGCGCGCCAATCATC  <  1:873431/139‑1
                                                                                                                                                                                                                                                                                      
CGTTGGTTGCCATTCTGGCGTCGGCGGTGGGTTATCTGGCAACGGCACAGGCGCTGTTAGCGAGGCTTGAAACCTCGGTTGCCATCTGGTTCCTGCTACTGGTGGTTTATCACGTTATCCGCCGCTGGATGCTGATCCAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/4377152‑4377011
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atccagcgGCGGATAACGTGATAAACCACCAGTAGCAGGAACCAGATGGCAACCGAGGTTTCAAGCCTCGCTAACAGCGCCTGTGCCGTTGCCAGATAACCCACCGCCGACGCCAGAATGGCAACCAACGGCGCGCCAATCATC  >  NZ_CP009273/4377031‑4377166

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