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 28 [0] [0] 28 hemH ferrochelatase

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

ACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGG  >  NC_000913/498681‑499158
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                      
cgccgctgatgggttacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgaggtatgggatgaaatggcacgcattctggcgcgcaaacgtagcattacggggctGGCg                                                                                                                                                                                                                                     >  2:7102‑M2/248‑251 (MQ=255)
          ggtttacagccgccagcaacagcaggcgctggcccaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatacgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAat                                                                                                                                                                                                                            >  2:178711‑M2/238‑250 (MQ=255)
                          caacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacgggcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGgcg                                                                                                                                                                                                           <  1:152728‑M2/30‑1 (MQ=255)
                            acagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCt                                                                                                                                                                                                         <  1:122346‑M2/32‑1 (MQ=255)
                                  ggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGc                                                                                                                                                                                                   <  1:11939‑M2/38‑1 (MQ=255)
                                               cgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattcacctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACat                                                                                                                                                                                      <  2:174704‑M2/51‑1 (MQ=255)
                                                   taccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctccggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacttagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATc                                                                                                                                                                                   <  2:49178‑M2/54‑1 (MQ=255)
                                                    accggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATAGGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCg                                                                                                                                                                                  >  1:135148‑M2/196‑250 (MQ=255)
                                                                 gtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAAt                                                                                                                                                                    >  2:179004‑M2/183‑251 (MQ=255)
                                                                         gggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGc                                                                                                                                                            >  1:130336‑M2/175‑251 (MQ=255)
                                                                                    acggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGc                                                                                                                                                  <  2:46047‑M2/87‑1 (MQ=255)
                                                                                                                agatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGcc                                                                                                                     <  2:203683‑M2/116‑1 (MQ=255)
                                                                                                                            ggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAaga                                                                                                           <  1:31254‑M2/126‑1 (MQ=255)
                                                                                                                             gcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGCGcc                                                                                                        >  1:158178‑M2/123‑251 (MQ=255)
                                                                                                                               agagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtc                                                                                                      >  1:249443‑M2/121‑251 (MQ=255)
                                                                                                                                                      tgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCCTccggg                                                                               >  2:89210‑M2/98‑249 (MQ=255)
                                                                                                                                                                 tatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCaa                                                                     >  1:77471‑M2/87‑250 (MQ=255)
                                                                                                                                                                  atccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCaaaa                                                                   <  1:199951‑M2/166‑1 (MQ=255)
                                                                                                                                                                     cgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAgg                                                                <  1:120396‑M2/169‑1 (MQ=255)
                                                                                                                                                                                            ggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGTATAATGCAAAAAGGCGTGCCGCGTGCCTGGGTAGAGt                                         >  2:39393‑M2/60‑251 (MQ=255)
                                                                                                                                                                                            ggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGGGCCTGGGTAGAGt                                         >  1:215320‑M2/60‑251 (MQ=255)
                                                                                                                                                                                                     ggggatgaactggcacggattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAAGGAGTGCTAATCTTGGTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCTGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTGGAACTca                                <  2:231021‑M2/201‑1 (MQ=255)
                                                                                                                                                                                                         attaactggcacgcattctggcgcgcaaacgtaggcttccggggatGGGGGGAAAAAATAAGAATATATTCCGGCGCTTAATGCCACGCCGGAACATAACGGAATGATGGCTAAACTTGTTGCCGCGTTACGCTAAAGCTGGGCCGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCCGGGATAATGCCAAAAGGCGTGCAGCGTTCCTGGGTAGAGTTTGAACTCACAgg                            <  1:186969‑M2/205‑1 (MQ=255)
                                                                                                                                                                                                            aactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGgct                         <  2:77471‑M2/208‑1 (MQ=255)
                                                                                                                                                                                                                    cgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATgc                 <  2:134137‑M2/216‑1 (MQ=255)
                                                                                                                                                                                                                                cgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTTAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGGTCAGCGGGCCTGGGGAGAGTTTTAACTCACAGGGCTGCTGATGCGCCGCTAACGTc     >  2:176974‑M2/24‑251 (MQ=255)
                                                                                                                                                                                                                                 gcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCt    >  1:179255‑M2/23‑251 (MQ=255)
                                                                                                                                                                                                                                   aaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGGAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTgg  >  1:189665‑M2/21‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                      
ACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGG  >  NC_000913/498681‑499158

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 20 ≤ 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.