BRESEQ :: Evidence
New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_012967 1733645 =0 (0.000)14 (1.580) 13/46 NT 100% coding (681/1413 nt) pykF pyruvate kinase PykF
?NC_012967 = 2775875 NA (NA)intergenic (‑45/+204) ECB_RS25125/cysH IS3 family transposase/phosphoadenosine phosphosulfate reductase

CACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGA  >  NC_012967/1733610‑1733677
                                   |                                
attccatttcatcgtccaacaaaatgggtgcagtac                                  >  2:258912‑M1/36‑36 (MQ=255)
attccatttcatcgtccaacaaaatgggtgcagtac                                  >  2:527704‑M1/36‑36 (MQ=255)
attccatttcatcgtccaacaaaatgggtgcagtac                                  <  1:730833‑M1/1‑1 (MQ=255)
actccatttcatcgtccaacaaaatgggtgcagtac                                  <  1:168861‑M1/1‑1 (MQ=255)
 ttccatttcattgtccaacaaaatgggtgcagtaca                                 >  1:961481‑M1/35‑36 (MQ=255)
 ttccatttcatcgtccaacaaaatgggtgctgtaca                                 <  2:408006‑M1/2‑1 (MQ=255)
 ttccatttcatcgtccaacaaaatgggtgcagtaca                                 >  1:798407‑M1/35‑36 (MQ=255)
 ttccatttcatcgtccaacaaaatgggtgcagtaca                                 >  2:562385‑M1/35‑36 (MQ=255)
 ttccatttcatcgtccaacaaaatgggtgcagtaca                                 <  1:725966‑M1/2‑1 (MQ=255)
 acaTCATCTCCAAAATCGAAAACCAGGAAGGCCTca                                 >  2:655813/1‑36 (MQ=255)
   ccatttcatcgtccaacaaaatgggtgcagtacaac                               >  1:518720‑M1/33‑36 (MQ=255)
         tcCAAAATCGAAAACCAGGAAGGCCTCAACAtgtac                         >  1:802586‑M1/1‑31 (MQ=255)
         catcgtccaacaaaatgggtgcagtaCAACAACTTc                         >  2:181002‑M1/27‑36 (MQ=255)
          ccAAAATCGAAAACCAGGAAGGCCTCAACAtgtact                        >  2:988009‑M1/1‑30 (MQ=255)
          atcgtccaacaaaatgggtgcagtaCAACAACTTcg                        >  1:669229‑M1/26‑36 (MQ=255)
           tcgtccaacaaaatgggtgcagtaCAACAACTTcga                       >  1:703351‑M1/25‑36 (MQ=255)
            cgtccaacaaaatgggtgtagtaCAACAACTTcgac                      >  1:428739‑M1/24‑36 (MQ=255)
              tccaacaaaatgggtgcagtaCAACAACTTcgacga                    >  1:237377‑M1/22‑36 (MQ=255)
              aaTCGAAAACCAGGAAGGCCTCAACAtgtactgacc                    >  2:429035‑M1/1‑26 (MQ=255)
               ccaacaaaatgggtgcagtaCAACAACTTCGACGaa                   <  2:14550‑M1/16‑1 (MQ=255)
               ccaacaaaatgggtgcagtaCAACAACTTCGACGaa                   <  2:649258‑M1/16‑1 (MQ=255)
               aTCGAAAACCAGGAAGGCCTCAACAtgtactgaccc                   <  1:691462‑M1/36‑12 (MQ=255)
                   aaaaCCAGGAAGGCCTCAACAtgtactgaccccaaa               <  2:323199‑M1/36‑16 (MQ=255)
                    aaagtgggtgcagtaCAACAACTTCGACGAAATCCt              >  2:762145‑M1/16‑36 (MQ=255)
                    aaaatgggtgcagtaCAACAACTTCGACGAAATCCt              <  1:817429‑M1/21‑1 (MQ=255)
                     aaatgggtgcagtaCAACAACTTCGACGAAATCCTc             >  2:589228‑M1/15‑36 (MQ=255)
                      aatgggtgcagtaCAACAACTTCGACGAAATCCTCg            <  2:968700‑M1/23‑1 (MQ=255)
                        tgggtgcagtaCAACAACTTCGACGAAATCCTCGaa          >  2:544944‑M1/12‑36 (MQ=255)
                         aggaggGCCTCAACAtgtactgaccccaaaaagttg         <  1:589780‑M1/36‑22 (MQ=255)
                         aGGAAGGTCTCAACAtgtactgacaccaaaaagttg         >  2:279697‑M1/1‑15 (MQ=255)
                          ggtgcagtaCAACAACTTCGACGAAATCCTCGAAGc        <  2:781527‑M1/27‑1 (MQ=255)
                          ggAAGGCCTCAACAtgtactgaccccaaaaagttgg        <  1:651316‑M1/36‑23 (MQ=255)
                           gtgcagtaCAACAACTTCGACGAAATCCTCGAAGcc       <  2:314777‑M1/28‑1 (MQ=255)
                              cagtaCAACAACTTCGACGAAATCCTCGAAGCctct    <  2:284243‑M1/31‑1 (MQ=255)
                               agtaCAACAACTTCGACGAAATCCTCGCAGCCTCTg   <  1:381863‑M1/32‑1 (MQ=255)
                                gaaCAACAACTTCGACGAAATCCTCGAAGCCTCTGa  >  2:556875‑M1/4‑36 (MQ=255)
                                   |                                
CACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGA  >  NC_012967/1733610‑1733677

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