New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 224957 =NA (NA)7 (0.080) 9/166 NT NA noncoding (1187/1542 nt) rrsH 16S ribosomal RNA of rrnH operon
?NC_000913 224974 = NA (NA)noncoding (1204/1542 nt) rrsH 16S ribosomal RNA of rrnH operon

AGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCA  >  NC_000913/224875‑224989
                                                                                                   |               
aGCGCAACCCTTATCTTTTGTTGTCAGCGGTCTGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtcatca                 <  1:461916/100‑1 (MQ=16)
aGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGAAAAAAACTGAAAGAAGGAGGGGATGACGTCAAGtcatca                 >  1:523851/1‑100 (MQ=11)
 gcgcAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTcatcat                >  1:312325/1‑100 (MQ=33)
   gcAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATgg              >  1:31728/1‑100 (MQ=33)
   gcAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATgg              >  2:477571/1‑100 (MQ=32)
     aaCCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGcc            <  2:1008763/100‑1 (MQ=33)
        ccTTATCTTTTGTTGCCAGCGGTCCGGTCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAATGTGGGGATGACGTCAAGTCATTATGGCCCtt         >  2:1143138/1‑100 (MQ=255)
         cTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGTCAGTGATATACTGGAGGTAGGAGGGGATGACGTCAAGTCATCATGGCCCTTa        >  2:628582/1‑100 (MQ=11)
           tATCTTTTGTTGCCAGCGGTCCGGCCGGGAATTCAAAGGTGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACg      <  1:2012596/100‑1 (MQ=12)
           tATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACg      <  1:90751/100‑1 (MQ=33)
             tCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGAc    <  1:1823190/100‑1 (MQ=32)
              cTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGAcc   >  1:2310164/1‑100 (MQ=31)
               ttttGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCa  >  2:830657/1‑100 (MQ=11)
                                                                                                   |               
AGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCA  >  NC_000913/224875‑224989

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.