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

ACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAG  >  NC_000913/2846251‑2846348
                                                 |                                                
aCGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGA                                                  >  1:1204591‑M1/1‑50 (MQ=255)
 cGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAt                                                 >  1:728952‑M1/1‑49 (MQ=255)
 cGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAt                                                 <  1:1561461‑M1/50‑2 (MQ=255)
  ggCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtc                                                <  1:671122‑M1/50‑3 (MQ=255)
  ggCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtc                                                <  1:1750172‑M1/50‑3 (MQ=255)
     gTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatc                                             >  1:1972892‑M1/1‑45 (MQ=255)
     gTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatc                                             >  1:1484830‑M1/1‑45 (MQ=255)
     gTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatc                                             >  1:1544181‑M1/1‑45 (MQ=255)
     gTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatc                                             >  1:1946796‑M1/1‑45 (MQ=255)
      tAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatcc                                            <  1:811577‑M1/50‑7 (MQ=255)
       aaTGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccg                                           >  1:2059965‑M1/1‑43 (MQ=255)
          gACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcg                                        >  1:548655‑M1/1‑40 (MQ=255)
            ccATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggc                                      >  1:1408044‑M1/1‑38 (MQ=255)
            ccATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggc                                      >  1:1190277‑M1/1‑38 (MQ=255)
               ttCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgag                                   <  1:2454119‑M1/50‑16 (MQ=255)
                 cAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagag                                 >  1:1831118‑M1/1‑33 (MQ=255)
                 cAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagag                                 >  1:1728759‑M1/1‑33 (MQ=255)
                  aGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagagg                                >  1:1274987‑M1/1‑32 (MQ=255)
                     ttGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcga                             >  1:2123640‑M1/1‑29 (MQ=255)
                         gTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcgacatc                         >  1:1610071‑M1/1‑25 (MQ=255)
                          ttCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcgacatct                        >  1:2264744‑M1/1‑24 (MQ=255)
                          ttCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcgacatct                        >  1:1304400‑M1/1‑24 (MQ=255)
                            cgtcgtTACCAAACAGCACCGAtcatccggcggcgagaggcgacatctc                       >  1:2314260‑M1/1‑22 (MQ=255)
                              tcgtTACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccac                    >  1:779439‑M1/1‑20 (MQ=255)
                                  tACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccacgcag                >  1:829875‑M1/1‑16 (MQ=255)
                                         cAGCACCGAtcatccggcggcgagaggcgacatctccacgcagtgcgcga         >  1:342226‑M1/1‑9 (MQ=255)
                                             aCCGAtcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcg     >  1:2420396‑M1/1‑5 (MQ=255)
                                              ccGAtcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcg     >  1:2298198‑M1/1‑4 (MQ=255)
                                               cGAtcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgca   <  1:1210253‑M1/50‑48 (MQ=255)
                                                gAtcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaa  <  1:671735‑M1/50‑49 (MQ=255)
                                                 |                                                
ACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAG  >  NC_000913/2846251‑2846348

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: