New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3598396 =NA (NA)4 (0.090)
+GATTTTGCGTTTGTCATCAGACTT		 4/234 6.5 NA noncoding (1/99 nt) RIP259 (repetitive extragenic palindromic) element, contains 2 REP sequences and 1 IHF site RIP259 (repetitive extragenic palindromic) element, contains 2 REP sequences and 1 IHF site
?NC_000913 = 4470517 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/3598424‑3598396
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AATCCGGGCTATCATGCCCGGATTAAGTCTGATGACAAAC  <  NC_000913/4470517‑4470478
                             ||||||||||||||||||||||||                                        
GTAGGCCTGATAAGCGTAGCGCATCAGGCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCTATCATGCCCGGATTAAGTCTGATGACAAAC  >  1:634275/1‑93
GTAGGCCTGATAAGCGTAGCGCATCAGGCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCTATCATGCCCGGATTAAGTCTGATGACAAAC  <  2:634275/93‑1
               GTAGCGCATCAGGCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCTATCATGCCCGGA                    >  1:182942/1‑60
               GTAGCGCATCAGGCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCTATCATGCCCGGA                    <  2:182942/60‑1
                             ||||||||||||||||||||||||                                        
GTAGGCCTGATAAGCGTAGCGCATCAGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3598424‑3598396
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AATCCGGGCTATCATGCCCGGATTAAGTCTGATGACAAAC  <  NC_000913/4470517‑4470478

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.