BRESEQ :: Evidence
New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_012967 3651091 =NA (NA)8 (0.830) 7/50 NT NA intergenic (‑33/‑47) hokA/ECB_RS25190 type I toxin‑antitoxin system toxin HokA/IS3 family transposase
?NC_012967 = 3741531 NA (NA)pseudogene (380/698 nt) ECB_RS18445 IS1‑like element IS1A family transposase

CCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAA  >  NC_012967/3651061‑3651121
                              |                              
ccAGCCTCCTGTGTTGATTTCCAACGAGTATGTACt                           >  1:511609/1‑36 (MQ=255)
ccAGCCTCCTGTGTTGATTTCCAACGAGTATGTACt                           >  1:728188/1‑36 (MQ=255)
ccAGCCTCCTGTGTTGATTTCCAACGAGTATGTACt                           <  1:421149/36‑1 (MQ=255)
        ctGTGTTGATTTCCAACGAGTATGTACTGCACCCAt                   >  1:14735/1‑36 (MQ=255)
          gtgtTGATTTCCAACGAGTATGTACTGCACCCAttt                 >  1:534021/1‑36 (MQ=255)
           tgtTGATTTCCAGCGAGTATGTACTGCACCCAtttt                <  1:988650/36‑1 (MQ=37)
           tgtTGATTTCCAACGAGTATGTACTGCACCCAtttt                <  2:539508/36‑1 (MQ=255)
            gtTGATTTCCAACGAGTATGTACTGCACCCATTttg               >  2:8909/1‑36 (MQ=255)
            gtTGATTTCCAACGAGTATGTACTGCACCCATTttg               <  2:863181/36‑1 (MQ=255)
                aTTTCCAACGAGTATGTACTGCACCCATTTTGTTgg           <  1:830107/36‑1 (MQ=255)
                      aaCGAGTATGTACTGCACCCCTTTTGTTGGACGATg     >  1:386779/1‑36 (MQ=18)
                       aCGAGTATGTACTGCACCCATTTTGTTGGACGATGa    >  1:92210/1‑36 (MQ=255)
                         gAGTATGTACTGCACCCATTTTGTTGGACGATGaaa  >  1:777365/1‑36 (MQ=255)
                         gAGTATGTACTGCACCCATTTTGTTGGACGATGaaa  >  1:786603/1‑36 (MQ=255)
                              |                              
CCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAA  >  NC_012967/3651061‑3651121

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