New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2729475 =14 (0.430)4 (0.130) 4/224 3.6 26.8% intergenic (‑31/+141) gltW/rrsG tRNA‑Glu/16S ribosomal RNA of rrnG operon
?NC_000913 = 4168363 9 (0.300)intergenic (+163/‑9) rrsB/gltT 16S ribosomal RNA of rrnB operon/tRNA‑Glu
Rejected: Coverage evenness skew score above cutoff.

GGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGA  >  NC_000913/2729413‑2729621
                                                              |                                                                                                                                                  
ggCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCTTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTa                                                              >  5:28775/1‑149 (MQ=11)
          tGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                                                     >  4:53146/1‑148 (MQ=32)
            ggCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc                                                    <  4:92578/147‑1 (MQ=33)
              ccTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatca                                                <  2:124821/149‑1 (MQ=37)
                             gACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAg                                 >  6:130442/1‑149 (MQ=37)
                               cacGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCa                               <  6:28775/149‑1 (MQ=37)
                                      aTTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGa                                              >  2:21828/1‑127 (MQ=38)
                                         aaatATCACGCAACGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAAt                                          <  1:66522/125‑1 (MQ=255)
                                         aaatATCACGCAACGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAAt                                          >  2:66522/4‑128 (MQ=255)
                                           ttATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCaaa                        >  7:19913/1‑144 (MQ=39)
                                           ttATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCaaa                        <  8:19913/144‑1 (MQ=39)
                                                aCGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttct              >  8:85659/1‑149 (MQ=39)
                                                          tGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaaggtaag    >  4:99236/1‑149 (MQ=39)
                                                            atatTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAagga  <  7:85659/149‑1 (MQ=39)
                                                              |                                                                                                                                                  
GGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGA  >  NC_000913/2729413‑2729621

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