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 30 [0] [0] 30 [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)27 (0.950) 20/92 0.2 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

AACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2846250‑2846300
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAA  >  NC_000913/2867269‑2867316
                                                                                                   
AACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACC                                                    >  1:417048/1‑49
 ACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGA                                                  >  1:1204591/1‑50
  CGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAT                                                 >  1:728952/1‑50
  CGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAT                                                 <  1:1561461/50‑1
   GGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGATC                                                <  1:671122/50‑1
   GGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGATC                                                <  1:1750172/50‑1
      GTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGATCATC                                             >  1:1484830/1‑50
      GTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGATCATC                                             >  1:1544181/1‑50
      GTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGATCATC                                             >  1:1972892/1‑50
      GTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGATCATC                                             >  1:1946796/1‑50
       TAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGATCATCC                                            <  1:811577/50‑1
        AATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGATCATCCG                                           >  1:2059965/1‑50
           GACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGATCATCCGGCG                                        >  1:548655/1‑50
             CCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGATCATCCGGCGGC                                      >  1:1408044/1‑50
             CCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGATCATCCGGCGGC                                      >  1:1190277/1‑50
                TTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGATCATCCGGCGGCGAG                                   <  1:2454119/50‑1
                  CAGTTTGCGTTCGTCGTTACCAAACAGCACCGATCATCCGGCGGCGAGAG                                 >  1:1728759/1‑50
                  CAGTTTGCGTTCGTCGTTACCAAACAGCACCGATCATCCGGCGGCGAGAG                                 >  1:1831118/1‑50
                   AGTTTGCGTTCGTCGTTACCAAACAGCACCGATCATCCGGCGGCGAGAGG                                >  1:1274987/1‑50
                      TTGCGTTCGTCGTTACCAAACAGCACCGATCATCCGGCGGCGAGAGGCGA                             >  1:2123640/1‑50
                          GTTCGTCGTTACCAAACAGCACCGATCATCCGGCGGCGAGAGGCGACATC                         >  1:1610071/1‑50
                           TTCGTCGTTACCAAACAGCACCGATCATCCGGCGGCGAGAGGCGACATCT                        >  1:2264744/1‑50
                           TTCGTCGTTACCAAACAGCACCGATCATCCGGCGGCGAGAGGCGACATCT                        >  1:1304400/1‑50
                             CGTCGTTACCAAACAGCACCGATCATCCGGCGGCGAGAGGCGACATCTC                       >  1:2314260/1‑49
                               TCGTTACCAAACAGCACCGATCATCCGGCGGCGAGAGGCGACATCTCCAC                    >  1:779439/1‑50
                                   TACCAAACAGCACCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAG                >  1:829875/1‑50
                                          CAGCACCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGA         >  1:342226/1‑50
                                              ACCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCG     >  1:2420396/1‑50
                                               CCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCG     >  1:2298198/1‑49
                                                CGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCA   <  1:1210253/50‑1
                                                 GATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAA  <  1:671735/50‑1
                                                                                                   
AACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2846250‑2846300
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAA  >  NC_000913/2867269‑2867316

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