New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2729630 =NA (NA)43 (0.260) 14/186 17.1 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.

CCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGAC  >  NC_000913/2729554‑2729664
                                                                                                  |            
ccTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTa              <  2:456663/99‑1 (MQ=255)
ccTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTa              <  2:2110657/99‑1 (MQ=255)
 cTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTAc             >  1:870938/1‑99 (MQ=255)
 cTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTAc             <  2:973339/99‑1 (MQ=255)
 cTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTAc             >  2:2968633/1‑99 (MQ=255)
 cTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTAc             <  2:1977837/99‑1 (MQ=255)
 cTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTAc             <  1:2471576/99‑1 (MQ=255)
 cTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTAc             <  1:1905289/99‑1 (MQ=255)
  ttGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTAcc            >  2:462450/1‑99 (MQ=255)
    gCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCtt          <  2:3798745/99‑1 (MQ=40)
    gCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCtt          >  2:1213462/1‑99 (MQ=40)
     cTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTg         >  2:599562/1‑99 (MQ=39)
      tATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGt        <  2:3741767/99‑1 (MQ=38)
       aTTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGtt       >  1:2060674/1‑99 (MQ=37)
       aTTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGtt       >  1:3317024/1‑99 (MQ=37)
       aTTCACTTTTCATCAGACAATCTGTGAGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGtt       >  1:3404132/1‑99 (MQ=37)
        ttCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTa      >  1:3774158/1‑99 (MQ=35)
        ttCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTa      <  1:16992/99‑1 (MQ=35)
         tcgcTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTCCCTTGTTAc     >  2:3918154/4‑99 (MQ=9)
         tCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCACCCGCAGGTTCCCCTCCGGTTACCTTGTTAc     >  1:1756068/1‑99 (MQ=255)
         tCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTAc     >  1:1050663/1‑99 (MQ=34)
         tCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTAc     >  1:2657447/1‑99 (MQ=34)
         tCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTAc     <  2:3244022/99‑1 (MQ=34)
           ccTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGa   <  1:2818932/98‑1 (MQ=14)
           aCTTTTCATCAGCCAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGa   <  1:2818871/99‑1 (MQ=17)
           aCTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGa   <  1:2348799/99‑1 (MQ=32)
           aCTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGa   >  1:1517243/1‑99 (MQ=32)
           aCTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGa   >  2:1406250/1‑99 (MQ=32)
           aCTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGa   >  2:1299375/1‑99 (MQ=32)
           aCTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGa   >  2:1266139/1‑99 (MQ=32)
           aCTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGa   <  1:2951388/99‑1 (MQ=32)
           aCTTTTCATCAGACAATCTGTGTGACCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCTACCGCAGGTTCCCCTACGGTTACCTTGTTACGa   <  1:4038135/99‑1 (MQ=12)
            cTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGCGGTCATCCAACCGCAGCTTCCCCTACGGTTACCTTGTTACGAc  >  1:1880474/1‑99 (MQ=255)
            cTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGAc  >  2:3801179/1‑99 (MQ=31)
                                                                                                  |            
CCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGAC  >  NC_000913/2729554‑2729664

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.