New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 3847603 =60 (1.070)4 (0.080) 3/96 NT 6.5% intergenic (+14/‑261) htrL/kbl hypothetical protein/glycine C‑acetyltransferase
?AP009048 3847628 = 58 (1.100)intergenic (+39/‑236) htrL/kbl hypothetical protein/glycine C‑acetyltransferase

TTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/3847675‑3847603
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑actGAGATAGGAAAAATGATGCTTTGGCATCATTTTTTCTTTTTATATCAATGCAGTAAAATCAATCATCAT  >  AP009048/3847628‑3847696
                                                                                                                                              
TTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCA                                                                         >  1:2581498/1‑71
  GATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCAC                                                                        >  1:1297543/1‑70
   ATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACT                                                                       <  1:3457835/70‑1
   ATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCAC                                                                        >  1:4261777/1‑69
                                                AGTTGGTATCAAAACGCTACCCACTGAGATAGGAAAAATGATGCTTTGGCATCATTTTTTCTTTTTATAT                          >  1:4651384/1‑70
                                                 GTTGGTATCAAAACGCTACCCACTGAGATAGGAAAAAT                                                         <  1:929170/38‑1
                                                          AAAACGCTACCAACTGAGATAGGAAAAATGATGCTTTGGCATCAT                                         <  1:158063/45‑1
                                                          AAAACGCTACCAACTGAGATAGGAAAAATGATGCTTTGGCATCAT                                         <  1:413811/45‑1
                                                                       CTGAGATAGGAAAAATGATGCTTTGGCATCATTTTTTCTTTTTATATCAATGCAGTAAAATCAATCATCAT  >  1:4210050/1‑71
                                                                       CTGAGATAGGAAAAATGATGCTTTGGCATCATTTTTTCTTTTTATATCAATGCAGTAAAATCAATCATCAT  >  1:4268904/1‑71
                                                                                                                                              
TTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/3847675‑3847603
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑actGAGATAGGAAAAATGATGCTTTGGCATCATTTTTTCTTTTTATATCAATGCAGTAAAATCAATCATCAT  >  AP009048/3847628‑3847696

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