New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2729630 =NA (NA)85 (0.410) 27/188 15.8 NA noncoding (1528/1542 nt) rrsG 16S ribosomal RNA of rrnG operon
?NC_000913 2729652 = NA (NA)noncoding (1506/1542 nt) rrsG 16S ribosomal RNA of rrnG operon
Rejected: Coverage evenness skew score above cutoff.

GCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACT  >  NC_000913/2729553‑2729665
                                                                                                   |             
gCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTa               >  2:1761986/1‑100 (MQ=255)
gCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTa               <  1:3496564/100‑1 (MQ=255)
 ccTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTAc              >  2:4405505/1‑100 (MQ=255)
 ccTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTAc              >  1:755252/1‑100 (MQ=255)
   ttGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCt            >  1:659193/1‑100 (MQ=255)
    tGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCtt           >  1:2493362/1‑100 (MQ=255)
     gCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTg          >  2:693783/1‑100 (MQ=40)
     gCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTg          <  2:4568326/100‑1 (MQ=40)
     gCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTg          <  1:4958881/100‑1 (MQ=40)
      cTATTCACTTTTCATCAGACAATCTGTGTGAGTACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGt         <  2:5050326/100‑1 (MQ=255)
      cTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGt         <  1:1742441/100‑1 (MQ=39)
       tATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGtt        >  1:711411/1‑100 (MQ=38)
       tATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGtt        <  1:4180907/100‑1 (MQ=38)
       tATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGtt        <  1:325432/100‑1 (MQ=38)
        aTTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTGAGGCAAGGAGGTGATCCAACCGCAGGTTCCCGTTCGCCTACCTTGTTa       <  1:4884550/100‑1 (MQ=11)
        aTTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTa       <  1:732798/100‑1 (MQ=37)
        aTTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTa       >  2:260204/1‑100 (MQ=37)
        aTTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTa       >  2:5381252/1‑100 (MQ=37)
         ttCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTAc      <  2:3363428/100‑1 (MQ=35)
          tCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACg     <  1:3583134/100‑1 (MQ=34)
          tCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACg     <  1:2501280/100‑1 (MQ=34)
           cACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGa    >  1:4344121/1‑100 (MQ=33)
           cACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGa    >  2:4441362/1‑100 (MQ=33)
           cACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGa    >  1:3012452/1‑100 (MQ=33)
           cACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGa    >  1:2034172/1‑100 (MQ=33)
            aCTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGAc   <  1:73675/100‑1 (MQ=32)
             cTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACt  >  1:3280477/1‑100 (MQ=31)
             cTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACt  >  2:2403652/1‑100 (MQ=31)
                                                                                                   |             
GCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACT  >  NC_000913/2729553‑2729665

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