New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 331850355 (0.990)4 (0.070) 4/248 NT 6.9% coding (534/2067 nt) treSA trehalose synthase A
?NC_002947 = 3318537 54 (1.010)coding (500/2067 nt) treSA trehalose synthase A

GCAGCAGATTGACCGCATCGCGCCCGGCGGGCACCTCGGGCAGCAGCGCCCAGTCTTCTTCGCGGATCTCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/3318433‑3318503
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctcgCCCATGGCCCCTACCCGGGGCTGTACCACATGGTCGAGATCCGCGAAGAAGACTGGGCGCTGCTGCCCGAGGTGCCCGCCGGGCGCGATGCGGTCAATCTGCTGCCAGCTCAGTGTGACGAGCTGAAGGCGCGCCATTACATCGTTG  <  NC_002947/3318537‑3318391
                                                                                                                                                                                                                          
GCAGCAGATTGACCGCATCGCGCCCGGCGGGCACCTCGGGCAGCAGCGCCCAGTCTTCTTCGCGGATCTCGCCCATGGCCCCTACCCGGGGCTGTACCACATGGTCG                                                                                                                 >  8:232866/1‑107
GCAGCAGATTGACCGCATCGCGCACGGCTGGCACCTCGGGCAGCAGCGCCCAGTCTTCTTCGCGGATCTCGCCCATGGCCCCTACCCGGGGCTGTACCACATGGTCG                                                                                                                 <  7:232866/107‑1
         TGACCGCATCGCGCCCGGCGGGCACCTCGGGCAGCAGCGCCCAGTCTTCTTCGCGGATCTCGCCCATGGCCCCTACCCGGGGCTGTACCACATGGTCGAGA                                                                                                              >  3:254163/1‑101
         TGACCGCATCGCGCCCGGCGGGCACCTCGGGCAGCAGCGCCCAGTCTTCTTCGCGGATCTCGCCCATGGCCCCTACCCGGGGCTGTACCACATGGTCGAGA                                                                                                              <  4:254163/101‑1
                                                                      GCCCATGGCCCCTACCCGGGGCTGTACCACATGGTCGAGATCCGCGAAGAAGACTGGGCGCTGCTGCCCGAGGTGCCCGCCGGGCGCGATGCGGTCAATCTGCTGCCAGCTCAGTGTGACGAGCTGAAGGCGCGCCATTACATCGTTG  <  7:33923/148‑1
                                                                      GCCCATGGCCCCTACCCGGGGCTGTACCACATGGTCGAGATCCGCGAAGAAGACTGGGCGCTGCTGCCCGAGGTGCCCGCCGGGCGCGATGCGGTCAATCTGCTGCCAGCTCAGTGTGACGAGCTGAAGGCGCGCCATTACATCGTTG  >  8:33923/1‑148
                                                                                                                                                                                                                          
GCAGCAGATTGACCGCATCGCGCCCGGCGGGCACCTCGGGCAGCAGCGCCCAGTCTTCTTCGCGGATCTCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/3318433‑3318503
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctcgCCCATGGCCCCTACCCGGGGCTGTACCACATGGTCGAGATCCGCGAAGAAGACTGGGCGCTGCTGCCCGAGGTGCCCGCCGGGCGCGATGCGGTCAATCTGCTGCCAGCTCAGTGTGACGAGCTGAAGGCGCGCCATTACATCGTTG  <  NC_002947/3318537‑3318391

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