Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,846,301 Δ20,968 bp [hycE][rpoS] 21 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2846301 2867268 20968 20 [0] [0] 19 [hycE]–[rpoS] [hycE],hycD,hycC,hycB,hycA,hypA,hypB,hypC,hypD,hypE,fhlA,ygbA,mutS,pphB,ygbI,ygbJ,ygbK,ygbL,ygbM,ygbN,[rpoS]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 28463000 (0.000)19 (0.680) 16/92 0.4 100% coding (172/1710 nt) hycE formate hydrogenlyase subunit HycE
?NC_000913 2867269 = 0 (0.000)coding (283/993 nt) rpoS RNA polymerase, sigma S (sigma 38) factor

CATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAC  >  NC_000913/2867224‑2867310
                                             |                                         
gtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatcatc                                       >  1:1278037‑M2/46‑50 (MQ=255)
gtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatcatc                                       >  1:1317757‑M2/46‑50 (MQ=255)
gtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatcatc                                       >  1:1005545‑M2/46‑50 (MQ=255)
    tgaccattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCcggc                                   >  1:632310‑M2/42‑50 (MQ=255)
         attcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCggcggc                                >  1:815691‑M2/37‑48 (MQ=255)
            cagtttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGagag                           >  1:1810234‑M2/34‑50 (MQ=255)
             agtttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAgg                          >  1:1585923‑M2/33‑50 (MQ=255)
               tttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCg                        <  1:1923353‑M2/20‑1 (MQ=255)
                  gcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGAc                      <  1:2058389‑M2/22‑1 (MQ=255)
                  gcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACa                     >  1:22890‑M2/28‑50 (MQ=255)
                    gttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACAtc                   >  1:2098378‑M2/26‑50 (MQ=255)
                          cgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg             >  1:1354803‑M2/20‑50 (MQ=255)
                           gttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGc            <  1:1450462‑M2/32‑1 (MQ=255)
                            ttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCa           <  1:2077518‑M2/33‑1 (MQ=255)
                                caaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgc       <  1:1806429‑M2/37‑1 (MQ=255)
                                caaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgc       >  1:442951‑M2/14‑50 (MQ=255)
                                    cagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGa   >  1:141032‑M2/10‑50 (MQ=255)
                                     agcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAc  <  1:608935‑M2/42‑1 (MQ=255)
                                     agcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAc  <  1:867091‑M2/42‑1 (MQ=255)
                                             |                                         
CATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAC  >  NC_000913/2867224‑2867310

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