BRESEQ :: Evidence
New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_012967 1821523 =0 (0.000)10 (1.130) 9/46 NT 100% coding (577/1041 nt) ynjI DUF1266 domain‑containing protein
?NC_012967 = 2775875 NA (NA)intergenic (‑45/+204) ECB_RS25125/cysH IS3 family transposase/phosphoadenosine phosphosulfate reductase

GCTTCCTCCATGTTTGCTTTAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCAC  >  NC_012967/1821489‑1821558
                                  |                                   
gCTTCCTCCATGTTTGCTTTAAAGGTATTGGCTCCa                                    >  1:163117/1‑36 (MQ=255)
 tccatttcatcgtccaacaaaatgggtgcagtaCAt                                   <  1:190257‑M1/3‑1 (MQ=255)
 tccatttcatcgtccaacaaaatgggtgcagtaCAt                                   >  2:88995‑M1/34‑36 (MQ=255)
  ccatttcatcgtccaacaaaatgggtgcagtaCATg                                  <  1:122046‑M1/4‑1 (MQ=255)
     tttcatcgtccaacaaactgggtgcagtaCATGGTc                               <  1:699940‑M1/7‑1 (MQ=255)
     ctccATGTTTGCTTTGAAGGTATTGGCTCCATGGTg                               <  1:513012‑M1/36‑2 (MQ=255)
      ttcatcgtccaacaaaatgggtgcagtaCATGGTCg                              <  2:156332‑M1/8‑1 (MQ=255)
        catcgtccaacaaaatgggtgcagtaCATGGTCGcc                            >  1:231739‑M1/27‑36 (MQ=255)
           cgtccaacaaaatgggtgcagtaCATGGTCGCCaga                         >  2:747815‑M1/24‑36 (MQ=255)
            tttGCTTTAAAGGTATTGGCTCCATGGTgtactgtc                        >  1:532498‑M1/1‑28 (MQ=255)
             ttGCTTTAAAGGTATTGGCTCCATGGTgtactgacc                       <  1:224650‑M1/36‑10 (MQ=255)
              ccaacaaaatgggtgcagtaCATGGTCGCCagaaag                      <  2:672536‑M1/16‑1 (MQ=255)
                aacaaaatgggtgcagtaCATGGTCGCCagaaagaa                    >  2:454845‑M1/19‑36 (MQ=255)
                  caaaatgggtgcagtaCATGGTCGCCAGAAAGaaaa                  <  1:943853‑M1/20‑1 (MQ=255)
                  caaaatgggtgcagtaCATGGTCGCCAGAAAGaaaa                  <  2:639754‑M1/20‑1 (MQ=255)
                   aaaatgggtgcagtaCATGGTCGCCAGAAAGAAAAt                 >  2:862896‑M1/16‑36 (MQ=255)
                    aaaGGTATTGGCTCCATGGTgtactgaccccataaa                >  2:34333‑M1/1‑20 (MQ=255)
                       tgggtgcagtaCATGGTCGCCAGAAAGAAAATGCTc             <  1:509459‑M1/25‑1 (MQ=255)
                        gTATTGGCTCCATGGTgtactgaccccaaaaagttg            <  1:537083‑M1/36‑21 (MQ=255)
                            gcagtaCATGGTCGCCAGAAAGAAAATGCTCCATTa        >  1:621987‑M1/7‑36 (MQ=255)
                             ggCTCCATGGTgtactgaccccaaaaagttggacag       >  2:934353‑M1/1‑11 (MQ=255)
                              gCTCCATGGTgtactgaccccaaaaagttggacagt      >  1:680187‑M1/1‑10 (MQ=255)
                               gtaCATGGTCGCCAGAAAGAAAATGCTCCATTAAgg     <  1:719203‑M1/33‑1 (MQ=255)
                                 ccATGGTgtactgaccccaaaaagttggacagttaa   >  2:157284‑M1/1‑7 (MQ=255)
                                 aCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGca   <  2:384198‑M1/35‑1 (MQ=255)
                                 aCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGca   <  1:477‑M1/35‑1 (MQ=255)
                                  CATGGTCGCCAGAAAGAAAATGCTCCATTTAGGcac  <  2:696462‑M1/36‑1 (MQ=255)
                                  |                                   
GCTTCCTCCATGTTTGCTTTAAAGGTATTGGCTCCATGGTCGCCAGAAAGAAAATGCTCCATTAAGGCAC  >  NC_012967/1821489‑1821558

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