Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 4,457,932 IS5 (–) +4 bp 5.6% intergenic (‑67/‑24) yjgA ← / → pmbA DUF615 domain‑containing protein YjgA/metalloprotease PmbA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2288919 =NA (NA)7 (0.080) 5/176 NT 7.4% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 4457935 84 (0.990)intergenic (‑70/‑24) yjgA/pmbA DUF615 domain‑containing protein YjgA/metalloprotease PmbA
* ? NC_000913 = 3653230NA (NA)3 (0.040) 3/168 NT 3.4% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 4457932 = 84 (0.990)intergenic (‑67/‑27) yjgA/pmbA DUF615 domain‑containing protein YjgA/metalloprotease PmbA
Rejected: Frequency below/above cutoff threshold.

CAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/2289012‑2288919
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTAACAGAGAACCGTTTTTCAGTGCGCATCTTAACTGGTAACATTAGCCTCTTT  <  NC_000913/4457935‑4457882
                                                                                                                                                    
CAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAC                                                   <  1:1032860/99‑1
CAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAC                                                   <  2:1773648/99‑1
 AAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACA                                                  <  2:2291887/99‑1
  AGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACAG                                                 >  2:1385042/1‑99
       TTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACAGAGAAC                                            >  2:132400/1‑99
           CGGGAAAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACAGAGAACCGTT                                        >  1:704761/1‑99
                GAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACAGAGAACCGTTTTTCA                                   <  2:357200/99‑1
                                               ATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACAGAGAACCGTTTTTCAGTGCGCATCTTAACTGGTAACATTAGCCTCT    <  1:2100270/99‑1
                                                 GACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACAGAGAACCGTTTTTCAGTGCGCATCTTAACTGGTAACATTAGCCTCTTT  <  2:1735631/99‑1
                                                                                                                                                    
CAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/2289012‑2288919
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTAACAGAGAACCGTTTTTCAGTGCGCATCTTAACTGGTAACATTAGCCTCTTT  <  NC_000913/4457935‑4457882

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