BRESEQ :: Evidence
New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_012967 = 2775875NA (NA)9 (1.020) 7/46 NT 100% intergenic (‑45/+204) ECB_RS25125/cysH IS3 family transposase/phosphoadenosine phosphosulfate reductase
?NC_012967 3015769 = 0 (0.000)coding (1190/1677 nt) ECB_RS14885 polysaccharide export protein

TGAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCA  >  NC_012967/3015736‑3015804
                                 |                                   
tGAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCAc                                   >  2:972160/1‑36 (MQ=255)
 gAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACt                                  <  2:52450/36‑1 (MQ=255)
 gAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACt                                  <  1:478905/36‑1 (MQ=255)
  aaTCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACtt                                 <  1:872366/36‑1 (MQ=255)
   aTCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTT                                >  1:188746‑M2/1‑36 (MQ=255)
     ctctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtg                              <  2:328517‑M2/36‑3 (MQ=255)
      tctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgt                             <  2:372091‑M2/36‑4 (MQ=255)
       catcgtccaacaaaatgggtgcagtaCACTTTCTGc                            >  1:791289‑M2/27‑36 (MQ=255)
          gCTGCAAAAGCTGGAGGAAGCCTCACTTTtgtactg                         <  2:127834‑M2/36‑8 (MQ=255)
            tccaacaaaatgggtgcagtaCACTTTCTGCCCAGt                       <  1:312984‑M2/15‑1 (MQ=255)
             gCAAAAGCTGGAGGAAGCCTCACTTTtgtactgacc                      >  1:676952‑M2/1‑26 (MQ=255)
               aagagctGGAGGAAGCCTCACTTTtgtactgacccc                    >  1:652336‑M2/4‑24 (MQ=255)
               aaaaGCTGGATGAAGCCTCACTTTtgtactgacccc                    >  1:681680‑M2/1‑24 (MQ=255)
                acaaaattggtgcagtaCACTTTCTGCCCAGTCATc                   <  2:173679‑M2/19‑1 (MQ=255)
                  aaaatgggtgcagtaCACTTTCTGCCCAGTCATcca                 >  1:26606‑M2/16‑36 (MQ=255)
                      tgggtgcagtaCACTTTCTGCCCAGTCATCCACCaa             <  2:430658‑M2/25‑1 (MQ=255)
                       gggtgcagtaCACTTTCTGCCCAGTCATCCACCaaa            <  1:892269‑M2/26‑1 (MQ=255)
                        ggtgcagtaCACTTTCTGCCCAGTCATCCAccaacg           >  1:455081‑M2/10‑34 (MQ=255)
                         gtgcagtaCACTTTCTGCCCAGTCATCCACCAaaga          >  1:764288‑M2/9‑36 (MQ=255)
                         gtgcagtaCACTTTCTGCCCAGTCATCCACCAaaga          >  1:289823‑M2/9‑36 (MQ=255)
                          gaAGCCTCACTTTtgtactgaccccaaaacgttgga         <  2:756664‑M2/36‑24 (MQ=255)
                            cagtaCACTTTCTGCCCAGTCATCCACCAaagaaga       <  1:91716‑M2/31‑1 (MQ=255)
                            cagtaCACTTTCTGCCCAGTCATCCACCAaagaaga       >  2:411315‑M2/6‑36 (MQ=255)
                            aGCCTCACTTTtgtactgaccccaacaagttggaca       <  1:848799‑M2/36‑26 (MQ=255)
                             gCCTCACTTTtgtactgaccccaaaaagttggacag      <  1:544242‑M2/36‑27 (MQ=255)
                               cTCACTTTtgtactgaccccaaaaagttggacagtt    >  2:231592‑M2/1‑8 (MQ=255)
                                 cACTTTtgtactgaccccaaaaagttggacagttaa  <  1:577402‑M2/36‑31 (MQ=255)
                                 CACTTTCTGCCCAGTCATCCACCAAAGACGAAGCCa  >  2:796419‑M2/1‑36 (MQ=255)
                                 |                                   
TGAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCA  >  NC_012967/3015736‑3015804

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