New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? minE = 154591946 (0.890)6 (0.130) 5/92 NT 12.2% coding (268/285 nt) yfcZ conserved hypothetical protein
?minE = 1545933 45 (0.960)coding (254/285 nt) yfcZ conserved hypothetical protein

GAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  minE/1545848‑1545919
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agctgAAATGCTGATTTTCCAGCTGGGTCTGCGTTAATTCCTCTCGGGTACGTACGCCCCTGTCAGTTGCTG  <  minE/1545933‑1545867
                                                                                                                                           
GAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCT                                                                      >  1:2625084/1‑71
         ACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAG                                                                        >  1:2051539/1‑60
                            ACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGAAATGCTGATTTTCCAGCTGGG                                               >  1:1863889/1‑66
                                       CGTACCCGAGAGGAATTAACGCAGACCCAGCTGAAATGCTGATTTTCCAGCTG                                                 <  1:3595344/53‑1
                                             CGAGAGGAATTAACGCAGACCCAGCTGAAATGCTGATTTT                                                        <  1:912659/40‑1
                                               AGAGGAATTAACGCAGACCCAGCTGAAATGCTGATTTTCC                                                      >  1:3773941/1‑40
                                                 AGGAATTAACGCAGACCCAGCTGAAATGCTGATTTTCC                                                      >  1:2654441/1‑38
                                                 AGGAATTAACGCAGACCCAGCTGAAATGCTGATTTTCC                                                      >  1:312424/1‑38
                                                                    GCTGAAATGCTGATTTTCCAGCTGGGTCTGCGTTAATTCCTCTCGGGTACGTACGCCCCTG            <  1:1065955/61‑1
                                                                       GAAATGCTGATTTTCCAGCTGGGTCTGCGTTAATTCCTCTCGGGTACGTACGCCCCTGTCAGTTGCTG  >  1:3723710/1‑68
                                                                       GAAATGCTGATTTTCCAGCTGGGTCTGCGTTAATTCCTCTCGGG                          <  1:3264375/44‑1
                                                                                                                                           
GAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  minE/1545848‑1545919
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agctgAAATGCTGATTTTCCAGCTGGGTCTGCGTTAATTCCTCTCGGGTACGTACGCCCCTGTCAGTTGCTG  <  minE/1545933‑1545867

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.