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

GGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGA  >  NC_000913/389627‑389784
                                                                               |                                                                              
ggCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGt                                                                               <  1:76117‑M1/81‑2 (MQ=255)
ggCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGt                                                                               <  2:76118‑M1/81‑2 (MQ=255)
ggCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGt                                                                               >  2:76437‑M1/1‑80 (MQ=255)
 gCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGtt                                                                              <  1:76119‑M1/81‑3 (MQ=255)
  cAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttg                                                                             >  1:76445‑M1/1‑78 (MQ=255)
  cAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttg                                                                             >  1:76457‑M1/1‑78 (MQ=255)
   aaCACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgc                                                                            >  1:76436‑M1/1‑77 (MQ=255)
    acacCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcg                                                                           <  2:76120‑M1/81‑6 (MQ=255)
    acacCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcg                                                                           >  2:76454‑M1/1‑76 (MQ=255)
     cacCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcg                                                                           >  1:76446‑M1/1‑75 (MQ=255)
              gTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgtt                                                                 >  1:76464‑M1/1‑66 (MQ=255)
               ttAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttg                                                                <  1:76121‑M1/81‑17 (MQ=255)
               ttAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttg                                                                >  1:76463‑M1/1‑65 (MQ=255)
                tAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgc                                                               >  2:76452‑M1/1‑64 (MQ=255)
                 aaGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcc                                                              >  1:76438‑M1/1‑63 (MQ=255)
                  aGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgccc                                                             <  1:76123‑M1/81‑20 (MQ=255)
                  aGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcc                                                              <  2:76122‑M1/80‑19 (MQ=255)
                      tAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaa                                                         >  1:76450‑M1/1‑58 (MQ=255)
                          tgttgtCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgc                                                     <  1:76124‑M1/81‑28 (MQ=255)
                           gttgtCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcg                                                    >  2:76443‑M1/1‑53 (MQ=255)
                             tgtCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgca                                                  >  2:76453‑M1/1‑51 (MQ=255)
                                  ctTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcga                                             >  2:76455‑M1/1‑46 (MQ=255)
                                   ttCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattc                                               <  1:76126‑M1/78‑34 (MQ=255)
                                   ttCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattc                                               <  2:76125‑M1/78‑34 (MQ=255)
                                    tCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatc                                           >  1:76462‑M1/1‑44 (MQ=255)
                                     cAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgat                                            <  2:76127‑M1/79‑37 (MQ=255)
                                      aaaCATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatctc                                         <  2:76128‑M1/81‑40 (MQ=255)
                                         cATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatctcacg                                      <  1:76129‑M1/81‑43 (MQ=255)
                                                 gcagcgcagGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttga                              >  1:76468‑M1/1‑31 (MQ=255)
                                                  cagcgcagGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttgat                             <  1:76130‑M1/81‑52 (MQ=255)
                                                   agcgcagGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttgatc                            <  1:76131‑M1/81‑53 (MQ=255)
                                                   agcgcagGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttgatc                            <  1:76132‑M1/81‑53 (MQ=255)
                                                   agcgcagGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttgatc                            >  2:76451‑M1/1‑29 (MQ=255)
                                                      gcagGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttgatccct                         >  1:76476‑M1/1‑26 (MQ=255)
                                                      gcagGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttgatccct                         >  1:76461‑M1/1‑26 (MQ=255)
                                                        agGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttgatcccttg                       >  1:76481‑M1/1‑24 (MQ=255)
                                                        agGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttgatcccttg                       >  1:76467‑M1/1‑24 (MQ=255)
                                                        agGAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttgatcccttg                       >  1:76473‑M1/1‑24 (MQ=255)
                                                         ggAGATTTGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttgatcccttgc                      >  2:76449‑M1/1‑23 (MQ=255)
                                                              tttGCGCAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttgatcccttgctcaac                 <  2:76061‑M1/81‑64 (MQ=255)
                                                                   gcAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttgatcccttgctcaacgatca            >  2:76460‑M1/1‑13 (MQ=255)
                                                                    cAGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttgatcccttgctcaacgatcac           >  2:76456‑M1/1‑12 (MQ=255)
                                                                     aGGCGACGAGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttgatcccttgctcaacgatcact          >  1:76472‑M1/1‑11 (MQ=255)
                                                                         gacgaGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttgatcccttgctcaacgatcactttat      <  1:76091‑M1/81‑75 (MQ=255)
                                                                            gaGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttgatcccttgctcaacgatcactttatgat   >  1:76459‑M1/1‑4 (MQ=255)
                                                                             aGGttgcggattgtcgttgcccagaactgcgcattcgatctcacgacctttgatcccttgctcaacgatcactttatgatc  >  1:76474‑M1/1‑3 (MQ=255)
                                                                               |                                                                              
GGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGA  >  NC_000913/389627‑389784

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: