New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 22497332 (0.240)29 (0.230) 14/186 12.8 49.1% noncoding (1203/1542 nt) rrsH 16S ribosomal RNA of rrnH operon
?NC_000913 2729972 = NA (NA)noncoding (1186/1542 nt) rrsG 16S ribosomal RNA of rrnG operon
Rejected: Coverage evenness skew score above cutoff.

GAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCA  >  NC_000913/224874‑224989
                                                                                                    |                
gAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtcatc                  >  1:3329993/1‑100 (MQ=33)
 aGCGCAACCCTTCTGTTTTGTTGCCAGCGGTCCGGCCCGGAATTCATC‑CGTTACTGCGAGAGATATAGTGGAGGACGGTGGGGATGCCGTTATGtcatca                 >  2:2910228/1‑100 (MQ=2)
 aGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtcatca                 <  2:2252605/100‑1 (MQ=33)
 aGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtcatca                 <  2:1017006/100‑1 (MQ=33)
 aGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtcatca                 <  2:222629/100‑1 (MQ=33)
  gcgcAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTcatcat                <  1:1144116/100‑1 (MQ=33)
  gcgcAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTcatcat                <  1:309541/100‑1 (MQ=33)
    gcAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATgg              <  1:2996960/100‑1 (MQ=33)
    gcAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATgg              >  1:2067671/1‑100 (MQ=33)
     cAACCTTTTTTTTTTGTGGCAGGCGGTCC‑GCCGGGAACACAAAGGGAGCCAG‑CAGTGATAAACTGGAGGAAGGTGGGGATGACTTGACGtcatccc                <  2:1132019/96‑3 (MQ=2)
     cAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGc             >  1:2492910/1‑100 (MQ=33)
     cAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGc             <  1:601483/100‑1 (MQ=32)
       aCCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGccc           <  2:379450/100‑1 (MQ=32)
        cccTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCt          >  1:2410343/1‑100 (MQ=33)
        cccTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCt          >  1:1282970/1‑100 (MQ=33)
        cccTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCt          <  2:274023/100‑1 (MQ=33)
          nttATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTa        >  2:1961668/2‑100 (MQ=17)
          cTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTa        >  1:2846315/1‑100 (MQ=33)
          cTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTa        <  1:2200832/100‑1 (MQ=32)
          cTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTa        <  2:64969/100‑1 (MQ=32)
           ttATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTAc       <  2:2202419/100‑1 (MQ=32)
           ttATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTAc       <  1:39172/100‑1 (MQ=33)
           ttATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTAc       <  1:2990095/100‑1 (MQ=33)
           ttATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTAc       >  1:1872142/1‑100 (MQ=33)
            tATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACg      >  2:3043032/1‑100 (MQ=32)
             aTCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGa     >  2:916417/1‑100 (MQ=32)
               cTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGAcc   >  2:1055007/1‑100 (MQ=31)
               cTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGAcc   <  1:944175/100‑1 (MQ=31)
               cTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGAcc   >  1:885898/1‑100 (MQ=31)
               cTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGAcc   <  1:3043044/100‑1 (MQ=31)
               cTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGAcc   >  1:2999863/1‑100 (MQ=31)
                ttttGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCGGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCa  >  1:1158825/1‑100 (MQ=11)
                ttttGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCa  >  1:1689339/1‑100 (MQ=11)
                ttttGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGGGGGGATGACGTCAAGTCATCATGGCCCTTACGACCa  >  1:2886410/1‑100 (MQ=11)
                                                                                                    |                
GAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAA‑GGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCA  >  NC_000913/224874‑224989

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