New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3643087 =22 (0.760)4 (0.150) 4/236 3.2 14.9% coding (46/753 nt) rsmJ 16S rRNA m(2)G1516 methyltransferase, SAM‑dependent
?NC_000913 3643107 = 25 (0.930)coding (26/753 nt) rsmJ 16S rRNA m(2)G1516 methyltransferase, SAM‑dependent
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGGATCTGTTCAAACGCTTCCGTGGTCGTGAACCGCAGCTGGATGCGATGCTGGAGCATTACGGCATTAAGGGCTGATCATTCAGTGAAAATCTGCTTAATTGATGAAACAGGCACCGGAGACGGTGCCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3643216‑3643087
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggtgcctGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATC  >  NC_000913/3643107‑3643170
                                                                                                                                                                                                  
TGGATCTGTTCAAACGCTTCCGTGGTCGTGAACCGCAGCTGGATGCGATGCTGGAGCATTACGGCATTAAGGGCTGATCATTCAGTGAAAATCTGCTTAATTGATGAAACAGGCACCGGAGACGG                                                                       >  4:26926/1‑125
                                              GATGCTGGAGCATTACGGCATTAAGGGCTGATCATTCAGTGAAAATCTGCTTAATTGATGAAACAGGCACCGGAGACGGTG                                                                     <  5:114266/81‑1
                                              GATGCTGGAGCATTACGGCATTAAGGGCTGATCATTCAGTGAAAATCTGCTTAATTGATGAAACAGGCACCGGAGACGGTG                                                                     >  6:114266/1‑81
                                                                               ATTCAGTGAAAATCTGCTTAATTGATGAAACAGGCACCGGAGACGGT                                                                      <  5:124239/47‑1
                                                                               ATTCAGTGAAAATCTGCTTAATTGATGAAACAGGCACCGGAGACGGT                                                                      >  6:124239/1‑47
                                                                                                        TGAAACAGGCACCGGAGACGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATC  <  1:74806/90‑1
                                                                                                        TGAAACAGGCACCGGAGACGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATC  >  2:74806/1‑90
                                                                                                         GAAACAGGCACCGGAGACGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCG                 <  1:80073/74‑1
                                                                                                         GAAACAGGCACCGGAGACGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCG                 >  2:80073/1‑74
                                                                                                                                                                                                  
TGGATCTGTTCAAACGCTTCCGTGGTCGTGAACCGCAGCTGGATGCGATGCTGGAGCATTACGGCATTAAGGGCTGATCATTCAGTGAAAATCTGCTTAATTGATGAAACAGGCACCGGAGACGGTGCCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3643216‑3643087
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggtgcctGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATC  >  NC_000913/3643107‑3643170

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