New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4583004 =88 (1.150)3 (0.040) 3/168 NT 3.8% coding (45/1590 nt) hsdM DNA methyltransferase M
?NC_000913 4583031 = 69 (0.950)coding (18/1590 nt) hsdM DNA methyltransferase M
Rejected: Frequency below cutoff threshold.

AATTTCTTCTTCTGCGGCTCGATGAATTTGGGCCGCTGCTTAATTTACGGAACTCACAATGAACAATAACGATCTGGTCGCGAAGCTGTGGAAGCTGTGCGAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4583106‑4583004
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgacCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAAT  >  NC_000913/4583031‑4583129
                                                                                                                                                                                                          
AATTTCTTCTTCTGCGGCTCGATGAATTTGGGCCGCTGCTTAATTTACGGAACTCACAATGAACAATAACGATCTGGTCGCGAAGCTGTGGAAGCTGTGC                                                                                                        <  2:1506776/100‑1
AATTTCTTCTTCTGCGGCTCGATGAATTTGGGCCGCTGCTTAATTTACGGAACTCACAATGAACAATAACGATCTGGTCGCGAAGCTGTGGAAGCTGTGC                                                                                                        >  2:465286/1‑100
 ATTTCTTCTTCTGCGGCTCGATGAATTTGGGCCGCTGCTTAATTTACGGAACTCACAATGAACAATAACGATCTGGTCGCGAAGCTGTGGAAGCTGTGCG                                                                                                       <  1:1523093/100‑1
 ATTTCTTCTTCTGCGGCTCGATGAATTTGGGCCGCTGCTTAATTTACGGAACTCACAATGAACAATAACGATCTGGTCGCGAAGCTGTGGAAGCTGTGCG                                                                                                       <  2:1106582/100‑1
   TTCTTCTTCTGCGGCTCGATGAATTTGGGCCGCTGCTTAATTTACGGAACTCACAATGAACAATAACGATCTGGTCGCGAAGCTGTGGAAGCTGTGCGAC                                                                                                     <  1:375229/100‑1
   TTCGTCTTCTGCGGCTCGATGAAGTTGGGCCGCTGCTTAATTTACGGAACTCACAATGAACAATAACGATCTGGTCGCGAAGCTGTGGAAGCTGTGCGAC                                                                                                     <  2:1588020/100‑1
                                                                    ACGATCTGGTCGCGAAGCTGTGGAAGCTGTGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAG                                    >  1:630580/1‑100
                                                                      GATCTGGTCGCGAAGCTGTGGAAGCTGTGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAA                                  >  1:1439168/1‑100
                                                                          TGGTCGCGAAGCTGTGGAAGCTGTGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAG                              >  2:338902/1‑100
                                                                                                  GCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACA      >  1:525934/1‑100
                                                                                                  GCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACA      <  2:2026819/100‑1
                                                                                                      CCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAAT  <  1:906041/100‑1
                                                                                                      CCAGATCGTTATTGTTCATTGTGAGTTCCGGAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAA   <  1:191989/99‑1
                                                                                                                                                                                                          
AATTTCTTCTTCTGCGGCTCGATGAATTTGGGCCGCTGCTTAATTTACGGAACTCACAATGAACAATAACGATCTGGTCGCGAAGCTGTGGAAGCTGTGCGAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4583106‑4583004
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgacCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAAT  >  NC_000913/4583031‑4583129

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