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 40 [0] [0] 42 [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)39 (1.110) 30/92 0.0 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

CGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGC  >  NC_000913/2846252‑2846346
                                                |                                              
cGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAt                                               <  1:887496‑M1/50‑2 (MQ=255)
cGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAt                                               <  1:783193‑M1/50‑2 (MQ=255)
 ggCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtc                                              <  1:2079994‑M1/50‑3 (MQ=255)
  gCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtca                                             >  1:447129‑M1/1‑47 (MQ=255)
   cGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcat                                            >  1:537418‑M1/1‑46 (MQ=255)
   cGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcat                                            <  1:2488291‑M1/50‑5 (MQ=255)
    gTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatc                                           >  1:1793896‑M1/1‑45 (MQ=255)
    gTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatc                                           >  1:282186‑M1/1‑45 (MQ=255)
     tAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatcc                                          <  1:756394‑M1/50‑7 (MQ=255)
      aaTGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccg                                         >  1:389020‑M1/1‑43 (MQ=255)
        tGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggc                                       >  1:2610505‑M1/1‑41 (MQ=255)
           ccATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggc                                    <  1:1534112‑M1/50‑13 (MQ=255)
            cATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcg                                   <  1:1840118‑M1/50‑14 (MQ=255)
             aTTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggc                                    >  1:1922747‑M1/1‑36 (MQ=255)
             aTTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggc                                    >  1:266763‑M1/1‑36 (MQ=255)
             aTTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggc                                    >  1:1944691‑M1/1‑36 (MQ=255)
                cAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagag                               >  1:1572320‑M1/1‑33 (MQ=255)
                 aGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagagg                              >  1:2474607‑M1/1‑32 (MQ=255)
                  gTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggc                             <  1:600033‑M1/50‑20 (MQ=255)
                  gTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggc                             >  1:2709990‑M1/1‑31 (MQ=255)
                      gCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcgaca                         >  1:2695915‑M1/1‑27 (MQ=255)
                      gCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcgaca                         >  1:546096‑M1/1‑27 (MQ=255)
                      gCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcgaca                         >  1:2539862‑M1/1‑27 (MQ=255)
                           cgtcgtTACCAAACAGCACCGAtcatccggcggcgtgaggcgacatctc                     >  1:376801‑M1/1‑22 (MQ=255)
                            gtcgtTACCAAACAGCACCGAtcatccggcggcgagaggcgacatctcc                    <  1:715572‑M1/49‑29 (MQ=255)
                            gtcgtTACCAAACAGCACCGAtcatccggcggcgagaggcgacatctcc                    <  1:1668146‑M1/49‑29 (MQ=255)
                               gtTACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccacgc                <  1:1934534‑M1/50‑33 (MQ=255)
                               gtTACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccacg                 >  1:2569409‑M1/1‑18 (MQ=255)
                                ttACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccacgca               <  1:150281‑M1/50‑34 (MQ=255)
                                ttACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccacg                 >  1:207325‑M1/1‑17 (MQ=255)
                                ttACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccacg                 >  1:180162‑M1/1‑17 (MQ=255)
                                 tACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccacgcag              >  1:730162‑M1/1‑16 (MQ=255)
                                 tACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccacgcag              >  1:739314‑M1/1‑16 (MQ=255)
                                   ccAAACAGCACCGAtcatccggcggcgagaggcgacatctccacgcagtg            <  1:326021‑M1/50‑37 (MQ=255)
                                   ccAAACAGCACCGAtcatccggcggcgagaggcgacatctccacgcagtg            >  1:2155631‑M1/1‑14 (MQ=255)
                                   ccAAACAGCACCGAtcatccggcggcgagaggcgacatctccacgcagt             <  1:2375619‑M1/49‑36 (MQ=255)
                                         aGCACCGAtcatccggcggcgagaggcgacatctccacgcagtgcgcgac      <  1:466677‑M1/50‑43 (MQ=255)
                                         aGCACCGAtcatccggcggcgagaggcgacatctccacgcagtgcgcgac      >  1:202095‑M1/1‑8 (MQ=255)
                                          gCACCGAtcatccggcggcgagaggcgacatctccacgcagtgcgcgacg     <  1:2457802‑M1/50‑44 (MQ=255)
                                              cGAtcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgc  >  1:2390816‑M1/1‑3 (MQ=255)
                                                |                                              
CGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGC  >  NC_000913/2846252‑2846346

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: