New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 401524713 (0.500)3 (0.130) 3/118 2.1 18.4% coding (2195/2262 nt) metE 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
?NC_000913 = 4015259 15 (0.650)coding (2207/2262 nt) metE 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Rejected: Frequency below/above cutoff threshold.

GCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAA  >  NC_000913/4015186‑4015314
                                                             |                                                                   
gCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGgtttct                                                               <  1:890044‑M1/68‑7 (MQ=255)
       gCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACc                                                        >  2:554227‑M2/1‑67 (MQ=255)
       gCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACc                                                        <  1:8767‑M2/68‑2 (MQ=255)
               gCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGc                                                <  1:483026/68‑1 (MQ=255)
               gCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGc                                                <  2:515590/68‑1 (MQ=255)
                 agagCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCAc                                              >  2:910740/1‑68 (MQ=255)
                    gcgcCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTgg                                           <  2:922455/68‑1 (MQ=255)
                    gcgcCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTgg                                           <  1:642971/68‑1 (MQ=255)
                     cgcCTGTGGGTCAACCCGGACTGTGGCCTGAACACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGc                                          >  1:246202/1‑68 (MQ=255)
                            gggTCACCCCGTCCTGTGGCCTGAAACCGCGCGGCTGGCCATAAACccgcgcgttttcaggccacagt                                   >  1:500801‑M2/1‑46 (MQ=255)
                              gTCAACCCGGACTGTGGCCTGAAAACGCGCGGgtttatggcccgccgcgcgttttaaggccacagtcc                                 >  1:839946‑M1/1‑32 (MQ=255)
                              gTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGt                                 >  1:59811/1‑68 (MQ=255)
                                   cccGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAgg                            >  2:61268/1‑68 (MQ=255)
                                      ggACTGTGGCCTTAAAACGCGCGGCGGGCCATAAACccgcgcgttttcaggccacagtccgggttgac                         >  1:839946‑M2/1‑36 (MQ=255)
                                       gACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGTCGAACATGGTGCAggcggc                        <  1:249271/68‑1 (MQ=255)
                                       gACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAggcggc                        <  1:826714/68‑1 (MQ=255)
                                        aCTGTGGCCTGAAAACGCGCGGgtttatggccagccgcgcggtttcaggccacaggacggggtgaccc                       >  1:500801‑M1/1‑22 (MQ=255)
                                         cTGTGGCCTGAAAACGCGCGGCTGGCCAGAACCCCGCGCGGCACTGGCGAACATGGTGCAGGCGgcgc                      >  2:619988/1‑68 (MQ=255)
                                             ggCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGaa                  >  2:723614/1‑68 (MQ=255)
                                                         cgcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCgggg      >  2:242926/1‑68 (MQ=255)
                                                             ggtttctggccagccgcgcgttttcaggccacagtccgggttgacccacaggcgctctgccggaatgc  <  1:8767‑M1/68‑68 (MQ=255)
                                                             ggtttctggccagccgcgcgttttcaggccacagtccgggttgacccacaggcgctctgccggaatgc  >  2:554227‑M1/1‑1 (MQ=255)
                                                             |                                                                   
GCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAA  >  NC_000913/4015186‑4015314

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