New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4532946 =98 (1.070)5 (0.060) 5/166 NT 5.3% coding (238/747 nt) yjhP putative methyltransferase
?NC_000913 4532970 = 86 (1.010)coding (214/747 nt) yjhP putative methyltransferase
Rejected: Frequency below cutoff threshold.

GGCTCGGGCGAGATGCTCTGCACCTGGGCCAGGGATCATGGCATTACGGGGACTGGCATCGACATGAGTTCGCTCTTCACCGCGCAGGCAAAGCGTCGCGCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4533048‑4532946
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgcggTGAAGAGCGAACTCATGTCGATGCCAGTCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGC  >  NC_000913/4532970‑4533068
                                                                                                                                                                                                          
GGCTCGGGCGAGATGCTCTGCACCTGGGCCAGGGATCATGGCATTACGGGGACTGGCATCGACATGAGTTCGCTCTTCACCGCGCAGGCAAAGCGTCGCG                                                                                                        >  1:1050943/1‑100
 GCTCGGGCGAGATGCTCTGCACCTGGGCCAGGGATCATGGCATTACGGGGACTGGCATCGACATGAGTTCGCTCTTCACCGCGCAGGCAAAGCGTCGCGC                                                                                                       >  2:2501752/1‑100
  CTCGGGCGAGATGCTCTGCACCTGGGCCAGGGATCATGGCATTACGGGGACTGGCATCGACATGAGTTCGCTCTTCACCGCGCAGGCAAAGCGTCGCGCG                                                                                                      <  1:1306158/100‑1
   TCGGGCGAGATGCTCTGCACCTGGGCCAGGGATCATGGCATTACGGGGACTGGCATCGACATGAGTTCGCTCTTCACCGCGCAGGCAAAGCGTCGCGCGG                                                                                                     <  1:646522/100‑1
   TCGGGCGAGATGCTCTGCACCTGGGCCAGGGATCATGGCATTACGGGGACTGGCATCGACATGAGTTCGCTCTTCACCGCGCAGGCAAAGCGTCGCGCGG                                                                                                     <  2:2372786/100‑1
                                                                       GCTCTTCATCGCGCAGGCAAAGCGTCGCGCGGTGAAGAGCGAACTCATGTCGATGCCAGTCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATC                                 >  1:616122/1‑100
                                                                        CTCTTCACCGCGCAGGCAAAGCGTCGCGCGGTGAAGAGCGAACTCATGTCGATGCCAGTCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCT                                >  1:301289/1‑100
                                                                         TCTTCACCGCGCAGGCAAAGCGTCGCGCGGTGAAGAGCGAACTCATGTCGATGCAAGTCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTC                               >  1:677888/1‑100
                                                                          CTTCACCGCGCAGGCAAAGCGTCGCGCGGTGAAGAGCGAACTCATGTCGATGCCAGTCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCG                              >  1:37662/1‑100
                                                                               CCGCGCAGGCAAAGCGTCGCGCGGTGAAGAGCGAACTCATGTCGATGCCAGTCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGA                         >  2:2008517/1‑100
                                                                                                 GCGCGGTGAAGAGCGAACTCATGTCGATGCCAGTCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAG       <  1:1448811/100‑1
                                                                                                 GCGCGGTGAAGAGCGAACTCATGTCGATGCCAGTCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAG       <  1:188259/100‑1
                                                                                                 GCGCGGTGAAGAGCGAACTCATGTCGATGCCAGTCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAG       <  2:1836565/100‑1
                                                                                                 GCGCGGGGAAGAGCGAACTCATGTCGATGCCAGTCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAG       <  1:2214438/100‑1
                                                                                                   GCGGTGAAGAGCGAACTCATGTCGATGCCAGTCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAA     >  1:1342321/1‑100
                                                                                                      GTGAAGAGCGAACTCATGTCGATGCCAGTCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGC  >  1:1840979/1‑100
                                                                                                                                                                                                          
GGCTCGGGCGAGATGCTCTGCACCTGGGCCAGGGATCATGGCATTACGGGGACTGGCATCGACATGAGTTCGCTCTTCACCGCGCAGGCAAAGCGTCGCGCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4533048‑4532946
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgcggTGAAGAGCGAACTCATGTCGATGCCAGTCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGC  >  NC_000913/4532970‑4533068

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.