BRESEQ :: Evidence
New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_012967 2774435 =NA (NA)8 (0.870) 7/48 NT 100% intergenic (‑67/+26) ECB_RS13765/ECB_RS25125 type I toxin‑antitoxin system Hok family toxin/IS3 family transposase
?NC_012967 = 3015774 0 (0.000)coding (1195/1677 nt) ECB_RS14885 polysaccharide export protein

ATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAGCCTG  >  NC_012967/3015739‑3015809
                                   |                                   
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            >  2:411315‑M2/6‑36 (MQ=255)
                         cagtaCACTTTCTGCCCAGTCATCCACCAaagaaga            <  1:91716‑M2/31‑1 (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)
                               aCTTTtgtactgaccccaaaaagttggacagttaaa      >  2:289993‑M2/1‑5 (MQ=255)
                               aCTTTtgtactgaccccaaaaagttggacagttaaa      <  1:396217‑M2/36‑32 (MQ=255)
                               ACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAg      <  1:77937‑M2/36‑1 (MQ=255)
                                cTTTCTGCCCAGTCATCCACCAAAGAAGAagccagc     >  2:922292/1‑36 (MQ=255)
                                  ttCTGCCCAGTCATCCACCAAAGAAGAAGCCAGCCt   <  1:825625/36‑1 (MQ=255)
                                   ttgtactgaccccaaaaagttggacagttaaacacg  <  2:526166‑M2/36‑36 (MQ=255)
                                   ttgtactgaccccaaaaagttggacagttaaacacg  >  1:490564‑M2/1‑1 (MQ=255)
                                   ttgtactgaccccaaaaagttggacagttaaacacg  >  1:263649‑M2/1‑1 (MQ=255)
                                   |                                   
ATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAGCCTG  >  NC_012967/3015739‑3015809

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