New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2729452NA (NA)17 (0.180) 7/158 NT NA intergenic (‑8/+164) gltW/rrsG tRNA‑Glu/16S ribosomal RNA of rrnG operon
?NC_000913 = 2729474 NA (NA)intergenic (‑30/+142) gltW/rrsG tRNA‑Glu/16S ribosomal RNA of rrnG operon

CGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCA  >  NC_000913/2729462‑2729572
            |                                                                                                  
cgcgTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTa             >  2:419421/1‑100 (MQ=11)
 gcgTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGTCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTAt            <  1:2199343/100‑1 (MQ=14)
   gTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTc          <  1:1293296/100‑1 (MQ=32)
    ttGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCa         <  1:1128834/100‑1 (MQ=32)
     ttccttatatTTTCGTGTAGGGTGCGCGTTCATTAATAGAAAGCGACCGGCCTTATTCTCTGAAGCGTCGCGCCCAGCGCGTAACCACCTTGCTATTCAc        <  1:1643833/94‑1 (MQ=2)
     tGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGATAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCa         <  2:779252/99‑1 (MQ=17)
      gCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACt       >  1:2032554/1‑100 (MQ=32)
       cGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtt      >  2:391113/1‑100 (MQ=32)
        gTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACttt     <  1:2289310/100‑1 (MQ=33)
          gATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc   <  1:1060931/100‑1 (MQ=35)
          gATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc   >  1:284769/1‑100 (MQ=35)
           atatTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtca  <  1:1539386/100‑1 (MQ=38)
           atatTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtca  >  1:2878954/1‑100 (MQ=38)
           atatTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtca  <  2:1555680/100‑1 (MQ=38)
            |                                                                                                  
CGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCA  >  NC_000913/2729462‑2729572

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