New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2238627 (0.120)7 (0.120) 5/280 8.0 50.4% noncoding (92/1542 nt) rrsH 16S ribosomal RNA of rrnH operon
?NC_000913 3428687 = NA (NA)noncoding (76/1542 nt) rrsD 16S ribosomal RNA of rrnD operon
Rejected: Coverage evenness skew score above cutoff.

GAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCA  >  NC_000913/3428546‑3428811
                                                                                                                                             |                                                                                                                            
gagGCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCTTGCTGTTTcg                                                                                                                             <  1:720438/143‑2 (MQ=12)
 ggccccGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCt                                                                                                                            <  2:1085615/140‑1 (MQ=17)
   gCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAG‑GAAACAGCTTGCTGTTTcgctga                                                                                                                         <  1:958774/143‑6 (MQ=11)
   ccccGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAG‑GAAACAGCTTGCTGTTTcgctga                                                                                                                         <  1:239031/142‑6 (MQ=11)
                                                                                                                           cgagataGCAAGCTGTTTCCTGTTCCCGTTCGACTGGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTTTTCAATTTAAAAGTTTGCTGGTCTGGTGAGGAAACTTCGTCTTGTATTACGTGTTCa  <  2:1074582/137‑1 (MQ=1)
                                                                                                                                             |                                                                                                                            
GAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGACGTTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGTTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAACAGCAAGCTGTTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCA  >  NC_000913/3428546‑3428811

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