New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2731065 =29 (0.790)6 (0.170) 6/260 3.6 17.5% noncoding (93/1542 nt) rrsG 16S ribosomal RNA of rrnG operon
?NC_000913 3428688 = NA (NA)noncoding (75/1542 nt) rrsD 16S ribosomal RNA of rrnD operon
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAAT  >  NC_000913/3428543‑3428797
                                                                                                                                                 |                                                                                                             
cAAGAGGCCCGAGGGTCCCCCTCTTTGTTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGTCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGt                                                                                                            <  5:66126/149‑1 (MQ=17)
cAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGt                                                                                                            <  4:150479/149‑1 (MQ=32)
 aaGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGt                                                                                                            >  4:11722/1‑148 (MQ=32)
   gagGCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGTTAc                                                                                                         <  5:166581/149‑1 (MQ=32)
    agGCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGTTAcc                                                                                                        >  5:51237/1‑149 (MQ=32)
        ccGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTc                                                                                                    >  8:7367/1‑149 (MQ=32)
         cGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTc                                                                                                    >  4:138316/1‑148 (MQ=32)
                                                                                                          acgcGTCCGCCACTCGTCAGCGAAACAGCAAGGTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAAt  <  2:178202/146‑1 (MQ=1)
                                                                                                                                                 |                                                                                                             
CAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAAT  >  NC_000913/3428543‑3428797

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.