New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AE015451 5680623 =40 (1.000)6 (0.160) 4/324 NT 13.7% coding (679/720 nt) rsmE 16S rRNA m3U1498 methyltransferase
?AE015451 5680692 = 37 (0.960)coding (591/618 nt) PP_4986 Channel protein, hemolysin III family

AGTTAAAAGTCGCCCCACAACTGCTGTGCCACCGACAGCGCCAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/5680666‑5680623
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccacGAAGTGCATCAGGCTGCCAGCGATGACGAACAGGTGCCAGATGCCATGCCAGTGGCGGAAACGGCTGTCGAAGGCGAAGAAAATGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGC  >  AE015451/5680692‑5680921
                                                                                                                                                                                                                                                                                  
AGTTAAAAGTCGCCCCACAACTGCTGTGCCACCGACAGCGCCACGAAGTGCATCAGGCTGCCAGCGATGACGAACAGGTGCCAGATGCCATGCCAGTGGCGGAAACGGCTGTCGAAGGCGAAGAAAATGATGC                                                                                                                                               <  1:644514/133‑1
AGTTAAAAGTCGCCCCACAACTGCTGTGCCACCGACAGCGCCACGAAGTGCATCAGGCTGCCAGCGATGACGAACAGGTGCCAGATGCCATGCCAGTGGCGGAAACGGCTGTCGAAGGCGAAGAAAATGATGC                                                                                                                                               >  2:644513/1‑133
 GTTAAAAGTCGCCCCACAACTGCTGTGCCACCGACAGCGCCACGAAGTGCATCAGGCTGCCAGCGATGACGAACAGGTGCCAGATGCCATGCCAGTGGCGGAAACGGCTGTCGAAGGCGAAGAAAATGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTAC                                        >  1:391424/1‑235
              CCACAACTGCTGTGCCACCGACAGCGCCACGAAGTGCATCAGGCTGCCAGCGATGACGAACAGGTGCCAGATGCCATGCCAGTGGCGGAAACGGCTGTCGAAGGCGAAGAAAATGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGC                                                                                       >  1:88422/1‑175
              CCACAACTGCTGTGCCACCGACAGCGCCACGAAGTGCATCAGGCTGCCAGCGATGACGAACAGGTGCCAGATGCCATGCCAGTGGCGGAAACGGCTGTCGAAGGCGAAGAAAATGATGC                                                                                                                                               <  1:293689/119‑1
              CCACAACTGCTGTGCCACCGACAGCGCCACGAAGTGCATCAGGCTGCCAGCGATGACGAACAGGTGCCAGATGCCATGCCAGTGGCGGAAACGGCTGTCGAAGGCGAAGAAAATGATGC                                                                                                                                               >  2:293688/1‑119
                                      CGCCACGAAGTGCATCAGGCTGCCAGCGATGACGAACAGGTGCCAGATGCCATGCCAGTGGCGGAAACGGCTGTCGAAGGCGAAGAAAATGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCCGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCG                                                                          >  1:221434/1‑164
                                       GCCACGAAGTGCATCAGGCTGCCAGCGATGACGAACAGGTGCCAGATGCCATGCCAGTGGCGGAAACGGCTGTCGAAGGCGAAGAAAATGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGCAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTT                                                                      >  2:33840/1‑167
                                       GCCACGAAGTGCATCAGGCTGCCAGCGATGACGAACAGGTGCCAGATGCCATGCCAGTGGCGGAAACGGCTGTCGAAGGCGAAGAAAATGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGC  >  1:658296/1‑235
                                         CACGAAGGGCATCAGGCTGCCAGCGATGACGAACAGGTGCCAGATGCCATGCCAGTGGCGGAAACGGCTGTCGAAGGCGAAGAAAATGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACC                                       <  2:506790/196‑1
                                                                                                                                                                                                                                                                                  
AGTTAAAAGTCGCCCCACAACTGCTGTGCCACCGACAGCGCCAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/5680666‑5680623
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccacGAAGTGCATCAGGCTGCCAGCGATGACGAACAGGTGCCAGATGCCATGCCAGTGGCGGAAACGGCTGTCGAAGGCGAAGAAAATGATGCCGACGGTGTAGAACACGCCGCCGGCAGCCAGCCAGGCGAAGCCGGCCGTGCCCAGCGAGTTCAGCAGCGGCTTGACTGCCACCAGCACGATCCAGCCCATTACCGCATAGATGATGATCGACAGGATCCGCGCCTCGGAGC  >  AE015451/5680692‑5680921

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