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

AGCAGTCCGGTGGCCGGCGCATCAAGCGCAGCCTGTTCATCGACGCCGCCGGCGTGCGCTTCCTCACCCGCGAAGAAGAACAGCGCCTGAGCCAGGTGCAGCTACTGAGTGACTACCTGGCGGGCAAGCGCCAGGAACTGCAAAACTGGAACGAGGCGCTGGGGCCGGTGGCGGACCTGTCGGCCAACCGCCGCAAGCTGACCAACATCGGCACCTTCCGCGCGTTTGCCCTGGCCTATCTGAAGAACCATCCCAACGTGCACCCCAACATGACC  >  NC_002947/5147719‑5147993
                                                                                                                                      |                                                                                                                                            
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)
                                                                                                                    cTGGCGGGCAAGCGCCAGGtagtcactcagtagctgcacctggctcaggcgctgttcttcttcgcgggtgaggaagcgcacgccggcggcgtcgatgaacaggctgcgcttgatgcgccggccaccggactgctgcatgccgcggtaatt           <  1:372042‑M2/150‑132 (MQ=255)
                                                                                                                             aaGCGCCAGGtagtcactcagtagctgcacctggctcaggcgctgttcttcttcgcgggtgaggaagcgcacgccggcggcgtcgatgaacaggctgcgcttgatgcgccggccaccggactgctgcatgccgcggtaattgcggaatga  >  2:372042‑M2/1‑10 (MQ=255)
                                                                                                                                      |                                                                                                                                            
AGCAGTCCGGTGGCCGGCGCATCAAGCGCAGCCTGTTCATCGACGCCGCCGGCGTGCGCTTCCTCACCCGCGAAGAAGAACAGCGCCTGAGCCAGGTGCAGCTACTGAGTGACTACCTGGCGGGCAAGCGCCAGGAACTGCAAAACTGGAACGAGGCGCTGGGGCCGGTGGCGGACCTGTCGGCCAACCGCCGCAAGCTGACCAACATCGGCACCTTCCGCGCGTTTGCCCTGGCCTATCTGAAGAACCATCCCAACGTGCACCCCAACATGACC  >  NC_002947/5147719‑5147993

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.