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

GCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCA  >  NC_000913/498678‑499176
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         
gctcgccgctgatggtttacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatg                                                                                                                                                                                                                                                           <  1:12873‑M2/1‑1 (MQ=255)
    gccgcgtatggtttccagccgccagcaacagcaggcgcgggcaccacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccctcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggaggaactggcacgcatgctggcgcgcaatcgtagcatgccggggatGGCgg                                                                                                                                                                                                                                                       <  2:46703‑M2/5‑1 (MQ=255)
     ccgctgatggtttacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatatagtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCggg                                                                                                                                                                                                                                                      <  1:102810‑M2/6‑1 (MQ=255)
           atcgtttacagccgccagcaacagcaggcgctggcacaacgttgaccggagatgcccgtagcgctgggaatgcgctacggctcgccattcctggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgccgccgctttatccgcaagtctcctgttctacggtcggtacggtacgggatgaactggcacgcattctcgtgcgcaaacgtagcattacggggatGGCGGGaaaaa                                                                                                                                                                                                                                                 >  2:751‑M2/240‑250 (MQ=255)
                 tacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGa                                                                                                                                                                                                                                           >  1:118224‑M2/234‑250 (MQ=255)
                 tacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGa                                                                                                                                                                                                                                           >  1:192270‑M2/234‑250 (MQ=255)
                                                  cgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACat                                                                                                                                                                                                         <  1:27252‑M2/51‑1 (MQ=255)
                                                  cgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATGCCGGCGCTTAATGCCACGCCGGAACat                                                                                                                                                                                                         <  2:252546‑M2/51‑1 (MQ=255)
                                                       accggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCg                                                                                                                                                                                                     <  1:131591‑M2/55‑1 (MQ=255)
                                                        ccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacgggcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGa                                                                                                                                                                                                    >  1:11384‑M2/195‑250 (MQ=255)
                                                                     tagcgctgggaatgagctacggctcgccatcactggaacgcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTACGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATc                                                                                                                                                                                      >  2:268346‑M2/182‑251 (MQ=255)
                                                                      agcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCt                                                                                                                                                                                     <  1:135588‑M2/71‑1 (MQ=255)
                                                                                    gctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATc                                                                                                                                                                       <  1:199724‑M2/85‑1 (MQ=255)
                                                                                       acggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCt                                                                                                                                                                    <  1:273282‑M2/88‑1 (MQ=255)
                                                                                          gctcgacatcactggaaagcgccgtagatgaactcctggccgatcatggagattatattgtggggttgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACAGATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTaaa                                                                                                                                                                 <  2:238733‑M2/91‑1 (MQ=255)
                                                                                              gccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTg                                                                                                                                                             <  1:159008‑M2/95‑1 (MQ=255)
                                                                                                catcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAg                                                                                                                                                           <  1:51452‑M2/97‑1 (MQ=255)
                                                                                                   cactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCgg                                                                                                                                                        >  1:9802‑M2/152‑251 (MQ=255)
                                                                                                                cgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAgc                                                                                                                                           <  2:194005‑M2/113‑1 (MQ=255)
                                                                                                                  tagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAgcgc                                                                                                                                         >  1:139110‑M2/137‑251 (MQ=255)
                                                                                                                       gaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCTGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTcg                                                                                                                                    <  2:142874‑M2/120‑1 (MQ=255)
                                                                                                                                                atatggtggtgctgccgctttatccgcaattctcctgttctacggtcggggcggtatgggatgagctggcaagcattccggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCTCTTAATGCCACGGCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAGGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtca                                                                                                           <  1:47093‑M2/145‑1 (MQ=255)
                                                                                                                                                                  tttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCCTCATCCGTGAATAATGc                                                                                         >  2:65488‑M2/89‑250 (MQ=255)
                                                                                                                                                                      tccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGAGAATAATGCaaaa                                                                                     >  1:76028‑M2/85‑250 (MQ=255)
                                                                                                                                                                                       ctacgatcggcgcggtatgggatgaactggcacgcattctggcgcgcacacgtagcattccggggatGGAGGGAAAAAATATGAATATATTCCGGCGCTTAATGCAACGCCGGAACATATCGAAATGCTGGCTAATCTTGTTGCCGCGTATCGCTAAAGATGAGAGGTAAAGAACTGAGAGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCCTCCGGgaagtaagcaaaaaggaggtcagagtgaaa                                                                    >  2:29234‑M2/68‑236 (MQ=255)
                                                                                                                                                                                         acggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGc                                                                                                                                                                                                                      >  2:204301‑M2/66‑103 (MQ=255)
                                                                                                                                                                                         acggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGc                                                                                                                                                                                                                      <  1:204301‑M2/38‑1 (MQ=255)
                                                                                                                                                                                         acggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTg                                                                   <  2:55307‑M2/184‑1 (MQ=255)
                                                                                                                                                                                                    ggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGGGCCTGGGGAGAGTTTGa                                                       >  2:83158‑M2/55‑250 (MQ=255)
                                                                                                                                                                                                    ggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATACGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGa                                                       >  1:16337‑M2/55‑250 (MQ=255)
                                                                                                                                                                                                         gggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTcac                                                 <  2:25076‑M2/202‑1 (MQ=255)
                                                                                                                                                                                                               aactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCg                                                                                                    >  1:191711‑M2/44‑195 (MQ=255)
                                                                                                                                                                                                               aactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGACGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCg                                                                                                    <  2:191711‑M2/152‑1 (MQ=255)
                                                                                                                                                                                                                     cacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGAt                                     <  1:11223‑M2/214‑1 (MQ=255)
                                                                                                                                                                                                                                                        atGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAgcagacggttggca  <  2:159825‑M2/249‑7 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         
GCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCA  >  NC_000913/498678‑499176

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.