Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 389,707 Δ10,497 bp [hemB][ddlA] [hemB], yaiT, insF‑2, insE‑2, yaiT, iprA, ampH, sbmA, yaiW, yaiY, yaiZ, [ddlA]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 389707 400203 10497 46 [0] [0] 46 [hemB]–[ddlA] [hemB],yaiT,insF‑2,insE‑2,yaiT,iprA,ampH,sbmA,yaiW,yaiY,yaiZ,[ddlA]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3897060 (0.000)46 (1.280) 39/160 0.1 100% coding (22/975 nt) hemB porphobilinogen synthase
?NC_000913 400204 = 0 (0.000)coding (720/1095 nt) ddlA D‑alanine‑‑D‑alanine ligase A

TTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATC  >  NC_000913/400124‑400281
                                                                                |                                                                             
ggcaacaccaggtcgttaaggctaagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggt                                                                               <  1:76117‑M2/1‑1 (MQ=255)
ggcaacaccaggtcgttaaggctaagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggt                                                                               <  2:76118‑M2/1‑1 (MQ=255)
ggcaacaccaggtcgttaaggctaagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggt                                                                               >  2:76437‑M2/81‑81 (MQ=255)
 gcaacaccaggtcgttaaggctaagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggtt                                                                              <  1:76119‑M2/2‑1 (MQ=255)
  caacaccaggtcgttaaggctaagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTg                                                                             >  1:76445‑M2/79‑81 (MQ=255)
  caacaccaggtcgttaaggctaagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTg                                                                             >  1:76457‑M2/79‑81 (MQ=255)
   aacaccaggtcgttaaggctaagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGc                                                                            >  1:76436‑M2/78‑81 (MQ=255)
    acaccaggtcgttaaggctaagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCg                                                                           <  2:76120‑M2/5‑1 (MQ=255)
    acaccaggtcgttaaggctaagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCg                                                                           >  2:76454‑M2/77‑81 (MQ=255)
     caccaggtcgttaaggctaagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCg                                                                           >  1:76446‑M2/76‑80 (MQ=255)
              gttaaggctaagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGtt                                                                 >  1:76464‑M2/67‑81 (MQ=255)
               ttaaggctaagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTg                                                                <  1:76121‑M2/16‑1 (MQ=255)
               ttaaggctaagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTg                                                                >  1:76463‑M2/66‑81 (MQ=255)
                taaggctaagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGc                                                               >  2:76452‑M2/65‑81 (MQ=255)
                 aaggctaagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGcc                                                              >  1:76438‑M2/64‑81 (MQ=255)
                  aggctaagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGccc                                                             <  1:76123‑M2/19‑1 (MQ=255)
                  aggctaagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGcc                                                              <  2:76122‑M2/18‑1 (MQ=255)
                      taagtgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGaa                                                         >  1:76450‑M2/59‑81 (MQ=255)
                          tgttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTgc                                                     <  1:76124‑M2/27‑1 (MQ=255)
                           gttgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTgcg                                                    >  2:76443‑M2/54‑81 (MQ=255)
                             tgtctcttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCa                                                  >  2:76453‑M2/52‑81 (MQ=255)
                                  cttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGa                                             >  2:76455‑M2/47‑81 (MQ=255)
                                   ttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTc                                               <  1:76126‑M2/33‑1 (MQ=255)
                                   ttcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTc                                               <  2:76125‑M2/33‑1 (MQ=255)
                                    tcaaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGAtc                                           >  1:76462‑M2/45‑81 (MQ=255)
                                     caaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGAt                                            <  2:76127‑M2/36‑1 (MQ=255)
                                      aaacatagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGAtctc                                         <  2:76128‑M2/39‑1 (MQ=255)
                                         catagcgcgcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACg                                      <  1:76129‑M2/42‑1 (MQ=255)
                                                 gcagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGa                              >  1:76468‑M2/32‑81 (MQ=255)
                                                  cagcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGAt                             <  1:76130‑M2/51‑1 (MQ=255)
                                                   agcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATc                            <  1:76131‑M2/52‑1 (MQ=255)
                                                   agcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATc                            <  1:76132‑M2/52‑1 (MQ=255)
                                                   agcgcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATc                            >  2:76451‑M2/30‑81 (MQ=255)
                                                      gcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCt                         >  1:76476‑M2/27‑81 (MQ=255)
                                                      gcaggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCt                         >  1:76461‑M2/27‑81 (MQ=255)
                                                        aggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTg                       >  1:76481‑M2/25‑81 (MQ=255)
                                                        aggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTg                       >  1:76467‑M2/25‑81 (MQ=255)
                                                        aggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTg                       >  1:76473‑M2/25‑81 (MQ=255)
                                                         ggagatttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGc                      >  2:76449‑M2/24‑81 (MQ=255)
                                                              tttgcgcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAAc                 <  2:76061‑M2/63‑1 (MQ=255)
                                                                   gcaggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCa            >  2:76460‑M2/14‑81 (MQ=255)
                                                                    caggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCAc           >  2:76456‑M2/13‑81 (MQ=255)
                                                                     aggcgacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACt          >  1:76472‑M2/12‑81 (MQ=255)
                                                                         gacgaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTAt      <  1:76091‑M2/74‑1 (MQ=255)
                                                                            gaggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGAt   >  1:76459‑M2/5‑81 (MQ=255)
                                                                             aggTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATc  >  1:76474‑M2/4‑81 (MQ=255)
                                                                                |                                                                             
TTCGCGCCATCTTCGTCAATGTACTTGGTGTCGTAGGCATAGAAATCGCTGGTGAGTACGATCTCTCCACAGGTGCTGGCTTGCGGATTGTCGTTGCCCAGAACTGCGCATTCGATCTCACGACCTTTGATCCCTTGCTCAACGATCACTTTATGATC  >  NC_000913/400124‑400281

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