New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1691425 =NA (NA)4 (0.100)
+GATTTTGCGTTTGTCATCAGACTT		 4/234 6.3 NA noncoding (1/97 nt) RIP126 (repetitive extragenic palindromic) element; contains 2 REP sequences and 1 IHF site RIP126 (repetitive extragenic palindromic) element; contains 2 REP sequences and 1 IHF site
?NC_000913 = 4470516 NA (NA)intergenic (+195/+11) mgtA/ridA magnesium transporter/enamine/imine deaminase, reaction intermediate detoxification
Rejected: Coverage evenness skew score above cutoff.

GTAGGCCTGATAAGCGTAGCGCATCAGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1691453‑1691425
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AATCCGGGCTATCATGCCCGGATTAAGTCTGATGACAAAC  <  NC_000913/4470516‑4470477
                             ||||||||||||||||||||||||                                        
GTAGGCCTGATAAGCGTAGCGCATCAGGCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCTATCATGCCCGGATTAAGTCTGATGACAAAC  >  1:634275/1‑93
GTAGGCCTGATAAGCGTAGCGCATCAGGCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCTATCATGCCCGGATTAAGTCTGATGACAAAC  <  2:634275/93‑1
               GTAGCGCATCAGGCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCTATCATGCCCGGA                    >  1:182942/1‑60
               GTAGCGCATCAGGCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCTATCATGCCCGGA                    <  2:182942/60‑1
                             ||||||||||||||||||||||||                                        
GTAGGCCTGATAAGCGTAGCGCATCAGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1691453‑1691425
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AATCCGGGCTATCATGCCCGGATTAAGTCTGATGACAAAC  <  NC_000913/4470516‑4470477

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