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 92 [7] [9] 95 [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 = 28463007 (0.050)76 (0.530) 40/92 0.5 90.7% coding (172/1710 nt) hycE formate hydrogenlyase subunit HycE
?NC_000913 2867269 = 9 (0.060)coding (283/993 nt) rpoS RNA polymerase, sigma S (sigma 38) factor

TTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAA  >  NC_000913/2867220‑2867318
                                                 |                                                 
ttGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGat                                                   <  1:4360813/50‑1 (MQ=255)
cggcgtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgat                                                   >  2:2926668‑M2/50‑50 (MQ=255)
cggcgtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgat                                                   >  1:2512201‑M2/50‑50 (MQ=255)
  gcgtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatca                                                 <  1:1251599‑M2/3‑1 (MQ=255)
    gtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatcatc                                               >  2:43013‑M2/46‑50 (MQ=255)
    gtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatcatc                                               >  1:1255152‑M2/46‑50 (MQ=255)
    gtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatcatc                                               <  2:3852949‑M2/5‑1 (MQ=255)
    gtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatcatc                                               >  1:4878193‑M2/46‑50 (MQ=255)
    gtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatcatc                                               >  1:4471853‑M2/46‑50 (MQ=255)
    gtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatcatc                                               >  2:4891500‑M2/46‑50 (MQ=255)
     taatgaccattcagtttgcgttcgtcgttaccaaacagcaccgaTCATcc                                              <  1:1643005‑M2/6‑1 (MQ=255)
      aatgaccattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCcg                                             >  2:425980‑M2/44‑50 (MQ=255)
      aatgaccattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCcg                                             >  2:4199295‑M2/44‑50 (MQ=255)
      aatgaccattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCcg                                             >  1:2300496‑M2/44‑50 (MQ=255)
      aatgaccattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCcg                                             >  1:5429715‑M2/44‑50 (MQ=255)
      aatgaccattcagtttgcgttcgtcgttaccaaacagcaccgaTCATAcg                                             >  2:15297‑M2/44‑50 (MQ=255)
       atgaccattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCcgg                                            >  2:3541046‑M2/43‑50 (MQ=255)
       atgaccattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCcgg                                            >  1:1013565‑M2/43‑50 (MQ=255)
       atgaccattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCcgg                                            >  1:1799827‑M2/43‑50 (MQ=255)
         gaccattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCcggcg                                          >  2:2372625‑M2/41‑50 (MQ=255)
         cggctgGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCcggcg                                          <  2:785750/50‑1 (MQ=255)
          ggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCcggcgg                                         <  2:3013697/50‑1 (MQ=255)
          ggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCcggcgg                                         <  2:3439255/50‑1 (MQ=255)
           ccattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCggcggc                                        >  1:3253553‑M2/39‑50 (MQ=255)
            cattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCggcggc                                        >  2:4588118‑M2/38‑49 (MQ=255)
            cattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCggcggc                                        >  1:682022‑M2/38‑49 (MQ=255)
            cattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCggcggc                                        >  2:3992332‑M2/38‑49 (MQ=255)
             attcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCggcggc                                        >  2:3264853‑M2/37‑48 (MQ=255)
             attcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCggcggc                                        >  2:2127515‑M2/37‑48 (MQ=255)
              ttcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCggcggc                                        >  1:5280728‑M2/36‑47 (MQ=255)
              ttcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCggcggc                                        >  2:3076072‑M2/36‑47 (MQ=255)
              ttcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCggcggc                                        >  1:190028‑M2/36‑47 (MQ=255)
                  gtttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGc                                 >  1:3594942‑M2/32‑50 (MQ=255)
                  gtttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGc                                 <  2:2676101‑M2/19‑1 (MQ=255)
                  gtttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGc                                 >  2:4055161‑M2/32‑50 (MQ=255)
                   tttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGc                                 >  2:2589981/1‑49 (MQ=255)
                   tttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGc                                 >  1:786297‑M2/31‑49 (MQ=255)
                   tttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGc                                 >  2:1745555‑M2/31‑49 (MQ=255)
                   tttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGc                                 >  1:1212949‑M2/31‑49 (MQ=255)
                    ttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGa                               >  2:1069384/1‑50 (MQ=255)
                     tgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGa                               >  2:138264‑M2/29‑49 (MQ=255)
                      gcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACa                             >  2:3397785‑M2/28‑50 (MQ=255)
                       cgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGGGTCGACAt                            >  2:206810‑M2/27‑50 (MQ=255)
                       cgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACAt                            >  1:3857658‑M2/27‑50 (MQ=255)
                        gttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACAtc                           >  2:5644158‑M2/26‑50 (MQ=255)
                         ttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACAtct                          >  1:836991‑M2/25‑50 (MQ=255)
                          tcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACAtctc                         >  2:755258‑M2/24‑50 (MQ=255)
                           cgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACAtctc                         >  1:612893‑M2/23‑49 (MQ=255)
                           accaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACAtctc                         >  2:471073/1‑49 (MQ=255)
                            gtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCa                       >  2:3879016‑M2/22‑50 (MQ=255)
                            gtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCa                       >  1:1119385‑M2/22‑50 (MQ=255)
                            gtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCa                       >  1:2934206‑M2/22‑50 (MQ=255)
                            gtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCa                       >  1:3643015‑M2/22‑50 (MQ=255)
                            gtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCa                       >  1:4965898‑M2/22‑50 (MQ=255)
                             tcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCAc                      <  2:5213710‑M2/30‑1 (MQ=255)
                             tcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCAc                      >  1:4085180‑M2/21‑50 (MQ=255)
                             caGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCAc                      <  2:961123/50‑1 (MQ=255)
                              cgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg                     >  2:4090792‑M2/20‑50 (MQ=255)
                              cgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg                     >  1:4947749‑M2/20‑50 (MQ=255)
                              cgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg                     >  2:303865‑M2/20‑50 (MQ=255)
                               gttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg                     >  1:5740458‑M2/19‑49 (MQ=255)
                               gttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg                     >  2:5715850‑M2/19‑49 (MQ=255)
                               gttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg                     >  1:2397609‑M2/19‑49 (MQ=255)
                               gttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg                     <  1:2768286‑M2/31‑1 (MQ=255)
                               gttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg                     >  2:2792249‑M2/19‑49 (MQ=255)
                               gttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg                     <  2:3951168‑M2/31‑1 (MQ=255)
                                aCGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCa                   <  1:5518074/50‑1 (MQ=255)
                                  accaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAg                  <  1:1634790‑M2/34‑1 (MQ=255)
                                  accaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAg                  <  2:606027‑M2/34‑1 (MQ=255)
                                  accaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGt                 >  1:3908424‑M2/16‑50 (MQ=255)
                                  accaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGt                 >  2:3558654‑M2/16‑50 (MQ=255)
                                  accaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGt                 >  2:4921200‑M2/16‑50 (MQ=255)
                                   ccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTg                >  1:2125543‑M2/15‑50 (MQ=255)
                                    caaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgc               >  2:3612599‑M2/14‑50 (MQ=255)
                                    caaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgc               >  2:5269519‑M2/14‑50 (MQ=255)
                                     aaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgcg              <  1:4811838‑M2/38‑1 (MQ=255)
                                     aaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgcg              >  2:572980‑M2/13‑50 (MQ=255)
                                      aacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgcgc             >  2:2205518‑M2/12‑50 (MQ=255)
                                      aacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgcgc             >  2:1459022‑M2/12‑50 (MQ=255)
                                       acagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGcgcg            <  1:1305396‑M2/40‑1 (MQ=255)
                                        cagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGa           >  2:670793‑M2/10‑50 (MQ=255)
                                         agcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAc          >  2:153387‑M2/9‑50 (MQ=255)
                                          gcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAcg         >  2:4877491‑M2/8‑50 (MQ=255)
                                           caccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAcgc        >  1:5777679‑M2/7‑50 (MQ=255)
                                           caccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAcgc        >  1:5622572‑M2/7‑50 (MQ=255)
                                              cgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACgcgc      >  2:705628‑M2/4‑49 (MQ=255)
                                              cgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACgcgc      >  2:317554‑M2/4‑49 (MQ=255)
                                              cgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACgcgc      >  1:609009‑M2/4‑49 (MQ=255)
                                               gaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCa     >  2:5164656‑M2/3‑49 (MQ=255)
                                               gaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCa     >  1:3545289‑M2/3‑49 (MQ=255)
                                               gaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCa     >  2:3486936‑M2/3‑49 (MQ=255)
                                               gaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCa     >  1:5044710‑M2/3‑49 (MQ=255)
                                               gaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCa     >  1:115711‑M2/3‑49 (MQ=255)
                                                 tcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCaaaa  >  2:4892478/1‑50 (MQ=255)
                                                 tcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCaaaa  >  2:2110589/1‑50 (MQ=255)
                                                 |                                                 
TTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAA  >  NC_000913/2867220‑2867318

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: