BRESEQ :: Evidence
New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_012967 2772857 =NA (NA)13 (1.540) 7/44 NT 93.7% intergenic (+378/‑54) iap/ECB_RS13760 alkaline phosphatase isozyme conversion aminopeptidase/IS4 family transposase
?NC_012967 = 4524530 1 (0.100)coding (502/549 nt) fimA type 1 fimbrial major subunit FimA

GGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGT  >  NC_012967/2772827‑2772882
                              |                         
ggAACACCCTGAACGGTGGGGAAAAACCCCATAAgc                      >  1:469934/1‑36 (MQ=255)
        cagAACAGTGAGGAAAAACCCCATAAGCGCTAACtt              <  1:722775/34‑1 (MQ=17)
        cTGAACGGTGGGTAAAAACCCCATAAGCGCTAACtt              <  2:789720/36‑1 (MQ=21)
             cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAggg         >  2:868949/1‑36 (MQ=255)
               gTGGAGAAAAACCCCATAAGCGCTAACTTAAGGGtt       <  1:637595/36‑1 (MQ=14)
                    gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGt  <  2:254619/36‑1 (MQ=40)
                              |                         
GGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGT  >  NC_012967/2772827‑2772882

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