New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 5147833NA (NA)7 (0.230) 7/252 2.7 100% coding (890/1302 nt) mscMB miniconductance mechanosensitive channel protein
?NC_002947 = 5147853 0 (0.000)coding (910/1302 nt) mscMB miniconductance mechanosensitive channel protein

CGCAATTACCGCGGCATGCAGCAGTCCGGTGGCCGGCGCATCAAGCGCAGCCTGTTCATCGACGCCGCCGGCGTGCGCTTCCTCACCCGCGAAGAAGAACAGCGCCTGAGCCAGGTGCAGCTACTGAGTGACTACCTGGCGGGCAAGCGCCAGGAACTGCAAAACTGGAACGAGGCGCTGGGGCCGGTGGCGGACCTGTCGGCCAACCGCCGCAAGCTGACCAACATCGGCACCTTCCGCGCGTTTGCCCTGGCCTATCTGAAGAACCATCCCAACGT  >  NC_002947/5147700‑5147977
                                                                                                                                     |                                                                                                                                                
cGCAATTACCGCGGCATGCAGCAGTCCGGTGGCCGGCGCATCAAGCGCAGCCTGTTCATCGACGCCGCCGGCGTGCGCTTCCTCACCCGCGAAGAAGAACAGCGCCTGAGCCAGGTGCAGCTACTGAGTGACTACCTGGCGGGCAAGCGC                                                                                                                                  <  1:255996‑M2/150‑1 (MQ=255)
                   agcagTCCGGTGGCCGGCGCATCAAGCGCAGCCTGTTCATCGACGCCGCCGGCGTGCGCTTCCTCACCCGCGAAGAAGAACAGCGCCTGAGCCAGGTGCAGCTACTGAGTGACTACCTGGCGCTTGCCCGCCAGGtagtcactcagtagc                                                                                                               <  2:806457‑M2/150‑16 (MQ=255)
                      agTCCGGTGGCCGGCGCATCAAGCGCAGCCTGTTCATCGACGCCGCCGGCGTGCGCTTCCTCACCCGCGAAGAAGAACAGCGCCTGAGCCAGGTGCAGCTACTGAGTGACTACCTGGCGCTTGCCCGCCAGGtagtcactcagtagctgc                                                                                                            >  1:806457‑M2/1‑132 (MQ=255)
                                                                                                                                tGACTACCTGGCGGGCAAGCGCCAGGtagtcactcagtagctgcacctggctcaggcgctgttcttcttcgcgggtgaggaagcgcacgccggcggcgtcgatgaacaggctgcgcttgatgcgccggccaccggactgctgcatgccgc  <  1:92518‑M2/150‑125 (MQ=255)
                                                                                                                                 gACTACCTGGCGGGCAAGCGCCAGGtagtcactcagtagctgcacctggctcaggcgctgttcttcttcgcgggtgaggaagc                                                                    >  1:463476‑M2/1‑25 (MQ=255)
                                                                                                                                 gACTACCTGGCGGGCAAGCGCCAGGtagtcactcagtagctgcacctggctcaggcgctgttcttcttcgcgggtgaggaagc                                                                    <  2:463476‑M2/83‑59 (MQ=255)
                                                                                                                                     |                                                                                                                                                
CGCAATTACCGCGGCATGCAGCAGTCCGGTGGCCGGCGCATCAAGCGCAGCCTGTTCATCGACGCCGCCGGCGTGCGCTTCCTCACCCGCGAAGAAGAACAGCGCCTGAGCCAGGTGCAGCTACTGAGTGACTACCTGGCGGGCAAGCGCCAGGAACTGCAAAACTGGAACGAGGCGCTGGGGCCGGTGGCGGACCTGTCGGCCAACCGCCGCAAGCTGACCAACATCGGCACCTTCCGCGCGTTTGCCCTGGCCTATCTGAAGAACCATCCCAACGT  >  NC_002947/5147700‑5147977

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