New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3643086 =51 (0.740)3 (0.050) 3/268 10.2 5.7% coding (46/753 nt) rsmJ 16S rRNA m(2)G1516 methyltransferase, SAM‑dependent
?NC_000913 3643106 = 51 (0.790)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.

GCGGTTCAGAAGAGCCGATGGATCTGTTCAAACGCTTCCGTGGTCGTGAACCGCAGCTGGATGCGATGCTGGAGCATTACGGCATTAAGGGCTGATCATTCAGTGAAAATCTGCTTAATTGATGAAACAGGCACCGGAGACGGTGCCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3643233‑3643086
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggtgcctGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAAT  >  NC_000913/3643106‑3643246
                                                                                                                                                                                                                                                                                                 
GCGGTTCAGAAGAGCCGATGGATCTGTTCAAACGCTTCCGTGGTCGTGAACCGCAGCTGGATGCGATGCTGGAGCATTACGGCATTAAGGGCTGATCATTCAGTGAAAATCTGCTTAATTGATGAAGCAGGCACCGGAGCCGG                                                                                                                                                    >  2:960626/1‑143
       AGAAGAGCCGATGGATCTGTTAAAACGCTTAGGTGGGCGTGAACCGCAGCTGGATGCGATGCTGGAGCATTACGGCATTAAGGGCTGATCATTCAGTGAAAATCTGCTTAATGGAGGCAACAGGCATCAGAGAAGGTGCatca                                                                                                                                             >  1:797250/1‑139
                                                                                                                         ATGAAACAGGCACCGGAGACGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTA                                                                                            <  1:224497/78‑1
                                                                                                                         ATGAAACAGGCACCGGAGACGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTA                                                                                            >  2:224497/1‑78
                                                                                                                             AACAGGCACCGGAGACGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTC                       >  2:1029855/1‑143
                                                                                                                                            CGGTGCCTGTTTCATCAAGTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCT        <  2:249661/143‑1
                                                                                                                                              gagcaTGTTTCACCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGTGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCA      <  2:1063365/138‑1
                                                                                                                                                  CTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGCAT  >  2:415749/1‑143
                                                                                                                                                                                                                                                                                                 
GCGGTTCAGAAGAGCCGATGGATCTGTTCAAACGCTTCCGTGGTCGTGAACCGCAGCTGGATGCGATGCTGGAGCATTACGGCATTAAGGGCTGATCATTCAGTGAAAATCTGCTTAATTGATGAAACAGGCACCGGAGACGGTGCCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3643233‑3643086
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggtgcctGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAAT  >  NC_000913/3643106‑3643246

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