New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? CP000730 2181134 =18 (0.460)4 (0.110) 4/262 5.3 19.1% noncoding (274/1540 nt) rrs4 16S ribosomal RNA
?CP000730 2297467 = NA (NA)noncoding (251/1540 nt) rrs5 16S ribosomal RNA
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ATCACCCTCTCAGGTCGGCTATGCATCGTTGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTAGCTCCGGTTTCCCGAAGTTATCCCAGTC  >  CP000730/2181101‑2181274
                                 |                                                                                                                                            
aTCACCCTCTCAGGTCGGCTATGCATCGTTGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCgg                                   <  1:419470/141‑1 (MQ=32)
 tCACCCTCTCAGGTCGGCTATGCATCGTTGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGt                                  >  2:401494/1‑141 (MQ=32)
   aCCCTCTAAGGTCGGCTATGCATCGTTGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTAt                                <  1:337912/141‑1 (MQ=17)
     ccTCTCAGGTCGGCTATGCATCGTTGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTa                              >  2:127142/1‑141 (MQ=32)
     ccTCTCAGGTCGGCTATGCATCGTTGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTa                              >  1:333476/1‑141 (MQ=32)
       tctcAGGTCGGCTATGCATCGTTGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTAGc                            <  1:302821/141‑1 (MQ=32)
       tctcAGGTCGGCTATGCATCGTTGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGCGTCAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTAGc                            <  1:157016/141‑1 (MQ=17)
                ggCTATGCATCGTTGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTAGCTCCGGTTTc                   >  2:126078/1‑141 (MQ=32)
                ggCTATGCATCGTTGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTAGCTCCGGTTTc                   >  2:54829/1‑141 (MQ=32)
                 gCTATGCATCGTTGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTAGCTCCGGTTTcc                  >  2:227775/1‑141 (MQ=32)
                    aTGCATCGTTGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTAGCTCCGGTTTCCCGa               >  2:171411/1‑141 (MQ=32)
                         tCGTTGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTAGCTCCGGTTTCCCGAAGTTa          >  1:407560/1‑141 (MQ=32)
                          cGTTGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTAGCTCCGGTTTCCCGAAGTTAt         >  1:142774/1‑141 (MQ=32)
                           gTTGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTAGCTCCGGTTTCCCGAAGTTATc        >  1:71709/1‑141 (MQ=32)
                            ttGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTAGCTCCGGTTTCCCGAAGTTATcc       >  1:178438/1‑141 (MQ=32)
                            ttGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTAGCTCCGGTTTCCCGAAGTTATcc       <  1:384924/141‑1 (MQ=32)
                                cTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTAGCTCCGGTTTCCCGAAGTTATCCCAGt   <  1:4280/141‑1 (MQ=32)
                                cTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTAGCTCCGGTTTCCCGAAGTTATCCCAGt   <  2:207343/141‑1 (MQ=32)
                                 ttGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTAGCTCCGGTTTCCCGAAGTTATCCCAGTc  >  2:435017/1‑141 (MQ=32)
                                 ttGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTAGCTCCGGTTTCCCGAAGTTATCCCAGTc  >  2:271806/1‑141 (MQ=32)
                                 |                                                                                                                                            
ATCACCCTCTCAGGTCGGCTATGCATCGTTGCCTTGGTAAGCCGTTACCTTACCAACTAGCTAATGCAGCGCGGATCCATCTATAAGTGACAGCAAGACCGTCTTTCACTTTTGAACCATGCGGTTCAAAATATCATCCGGTATTAGCTCCGGTTTCCCGAAGTTATCCCAGTC  >  CP000730/2181101‑2181274

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