BRESEQ :: Evidence
New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_012967 1460160 =NA (NA)9 (0.970) 7/48 NT 100% intergenic (‑115/‑55) gap/ECB_RS07290 type I glyceraldehyde‑3‑phosphate dehydrogenase/IS1 family transposase
?NC_012967 = 1604671 0 (0.000)coding (1212/2379 nt) ECB_RS07985 prophage tail fiber N‑terminal domain‑containing protein

ATTCTCACAGTGAAACGTAACGTACTGAAAACGGGGTGATGCTGCCAACTTACTGATTTAGTGTATGA  >  NC_012967/1460126‑1460193
                                  |                                 
aTTCTCACAGTGAAACGTAACGTACTGAAAACgggg                                  >  1:999801/1‑36 (MQ=255)
 ttCTCACAGTGAAACGTAACGTACTGAAAACGGGGt                                 <  1:49406/36‑1 (MQ=255)
  tctcACAGTGAAACGTAACGTACTGAAAACGGGGTg                                >  1:490362/1‑36 (MQ=255)
     cacaGTGAAACGTAACGTACTGAAAACGGGGTGAtg                             <  1:912043/36‑1 (MQ=255)
     cacaGTGAAACGTAACGTACTGAAAACGGGGTGAtg                             <  2:176128/36‑1 (MQ=255)
     cacaGTGAAACGTAACGTACTGAAAACGGGGTGAtg                             <  2:543437/36‑1 (MQ=255)
       caGTGAAACGTAACGTACTGAAAACGGGGTGAtgct                           <  2:93170/36‑1 (MQ=255)
       aaGTGAAACGTAACGTACTGAAAACGGGGTGAtgct                           <  1:278665/35‑1 (MQ=255)
        aGTGAAACGTAACGTACTGAAAACGGGGTGAtgctg                          <  2:771468/36‑1 (MQ=255)
         gTGAAACGTAACGTACTGAAAACGGGGTGAtgctgc                         >  2:899320/1‑36 (MQ=255)
              acgtaacgtaCTGAAAACGGGGTGATGCTGCCAACt                    <  1:972376/36‑1 (MQ=255)
                    cgtaCTGAAAACGGGGTGATGCTGCTAACTTTCTgt              >  1:807673/1‑35 (MQ=25)
                            aaaCGGGGTGATGCTGCCAACTTACTGATTTAgtgt      >  1:533866/1‑36 (MQ=255)
                            aaaCGGGGTGATGCTGCCAACTTACTGATTAAgtgt      >  2:218819/1‑36 (MQ=17)
                            aaaCGGGGTGATGCTGCCAACTCACTGATATAgtgt      >  1:987707/1‑36 (MQ=17)
                             agcGGGGTGATGCTGCCAACTTACTGATTTAGTGTa     <  1:222288/34‑1 (MQ=255)
                              aCGGGGTGATGCTGCCAACTTATTGATTTAGTGTat    >  1:201879/1‑36 (MQ=14)
                               cGGGGTGATGCTGCCAACTTACTGATTTAGTGTatg   <  2:82924/36‑1 (MQ=37)
                                ggggTGATGCTGCCAACTTACTGATTTAGTGTatga  >  1:949861/1‑36 (MQ=33)
                                ggggTGATGCTGCCAACTTACTGATTTAGTGTatga  >  1:586776/1‑36 (MQ=33)
                                  |                                 
ATTCTCACAGTGAAACGTAACGTACTGAAAACGGGGTGATGCTGCCAACTTACTGATTTAGTGTATGA  >  NC_012967/1460126‑1460193

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