Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 498,522 Δ406 bp coding (468‑873/963 nt) hemH → ferrochelatase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 498522 498927 406 35 [0] [0] 35 hemH ferrochelatase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4985210 (0.000)35 (1.260) 30/488 0.1 100% coding (467/963 nt) hemH ferrochelatase
?NC_000913 498928 = 0 (0.000)coding (874/963 nt) hemH ferrochelatase

GCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATATCGTTTATACGTGATTACGCCGATAACCACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCG  >  NC_000913/498272‑498770
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
gCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATg                                                                                                                                                                                                                                                          <  1:12873‑M1/251‑2 (MQ=255)
    gccgcgtatggtttcCAGCCGCCAGCAACAGCAGGCGCGGGCACCACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCCTCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGAGGAACTGGCACGCATGCTGGCGCGCAATCGTAGCATGCCGGGGATggcgg                                                                                                                                                                                                                                                      <  2:46703‑M1/244‑6 (MQ=255)
     ccgcTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATAGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcggg                                                                                                                                                                                                                                                     <  1:102810‑M1/251‑7 (MQ=255)
           atcgtttACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTGACCGGAGATGCCCGTAGCGCTGGGAATGCGCTACGGCTCGCCATTCCTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCCGCCGCTTTATCCGCAAGTCTCCTGTTCTACGGTCGGTACGGTACGGGATGAACTGGCACGCATTCTCGTGCGCAAACGTAGCATTACGGGGATggcgggaaaaa                                                                                                                                                                                                                                                >  2:751‑M1/4‑239 (MQ=255)
                 tACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatga                                                                                                                                                                                                                                          >  1:118224‑M1/1‑233 (MQ=255)
                 tACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatga                                                                                                                                                                                                                                          >  1:192270‑M1/1‑233 (MQ=255)
                                                  cGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacat                                                                                                                                                                                                        <  1:27252‑M1/251‑52 (MQ=255)
                                                  cGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatatgccggcgcttaatgccacgccggaacat                                                                                                                                                                                                        <  2:252546‑M1/251‑52 (MQ=255)
                                                       aCCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcg                                                                                                                                                                                                    <  1:131591‑M1/250‑56 (MQ=255)
                                                        ccGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGGCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcga                                                                                                                                                                                                   >  1:11384‑M1/1‑194 (MQ=255)
                                                                     tAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAACGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattacggcgcttaatgccacgccggaacatatcgaaatgatggctaatc                                                                                                                                                                                     >  2:268346‑M1/1‑181 (MQ=255)
                                                                      aGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatct                                                                                                                                                                                    <  1:135588‑M1/251‑72 (MQ=255)
                                                                                    gCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatc                                                                                                                                                                      <  1:199724‑M1/251‑86 (MQ=255)
                                                                                       aCGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgct                                                                                                                                                                   <  1:273282‑M1/251‑89 (MQ=255)
                                                                                          gCTCGACATCACTGGAAAGCGCCGTAGATGAACTCCTGGCCGATCATGGAGATTATATTGTGGGGTTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacagatcgaaatgatggctaatcttgttgccgcgtatcgctaaa                                                                                                                                                                <  2:238733‑M1/251‑92 (MQ=255)
                                                                                              gCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctg                                                                                                                                                            <  1:159008‑M1/251‑96 (MQ=255)
                                                                                                cATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgag                                                                                                                                                          <  1:51452‑M1/251‑98 (MQ=255)
                                                                                                   cACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcgg                                                                                                                                                       >  1:9802‑M1/1‑151 (MQ=255)
                                                                                                                cGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagc                                                                                                                                          <  2:194005‑M1/251‑114 (MQ=255)
                                                                                                                  tAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgc                                                                                                                                        >  1:139110‑M1/1‑136 (MQ=255)
                                                                                                                       gAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgctggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcg                                                                                                                                   <  2:142874‑M1/251‑121 (MQ=255)
                                                                                                                                                atatGGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGGGCGGTATGGGATGAGCTGGCAAGCATTCCGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggctcttaatgccacggcggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaggctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttca                                                                                                          <  1:47093‑M1/251‑146 (MQ=255)
                                                                                                                                                                  tttATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcctcatccgtgaataatgc                                                                                         >  2:65488‑M1/1‑88 (MQ=255)
                                                                                                                                                                      tCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgagaataatgcaaaa                                                                                     >  1:76028‑M1/1‑84 (MQ=255)
                                                                                                                                                                                       cTACGATCGGCGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCACACGTAGCATTCCGGGGATggagggaaaaaatatgaatatattccggcgcttaatgcaacgccggaacatatcgaaatgctggctaatcttgttgccgcgtatcgctaaagatgagaggtaaagaactgagagccgtcgcgaagagcctcgtcggcggttttcatcctccgggaagtaagcaaaaaggaggtcagagtgaaa                                                                    >  2:29234‑M1/1‑67 (MQ=255)
                                                                                                                                                                                         aCGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctg                                                                   <  2:55307‑M1/249‑185 (MQ=255)
                                                                                                                                                                                         aCGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgc                                                                                                                                                                                                                     >  2:204301‑M1/1‑65 (MQ=255)
                                                                                                                                                                                         aCGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgc                                                                                                                                                                                                                     <  1:204301‑M1/103‑39 (MQ=255)
                                                                                                                                                                                                    ggTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgggcctggggagagtttga                                                       >  2:83158‑M1/1‑54 (MQ=255)
                                                                                                                                                                                                    ggTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatacgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttga                                                       >  1:16337‑M1/1‑54 (MQ=255)
                                                                                                                                                                                                         gggATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcac                                                 <  2:25076‑M1/251‑203 (MQ=255)
                                                                                                                                                                                                               aaCTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccg                                                                                                   >  1:191711‑M1/1‑43 (MQ=255)
                                                                                                                                                                                                               aaCTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccgacgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccg                                                                                                   <  2:191711‑M1/195‑153 (MQ=255)
                                                                                                                                                                                                                     cACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgat                                     <  1:11223‑M1/251‑215 (MQ=255)
                                                                                                                                                                                                                                                        aTggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcgccgctaacgtctggtaaagcagacggttggca  <  2:159825‑M1/251‑250 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
GCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATATCGTTTATACGTGATTACGCCGATAACCACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCG  >  NC_000913/498272‑498770

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.