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

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4985210 (0.000)50 (0.590) 40/598 0.8 100% coding (467/963 nt) hemH ferrochelatase
?NC_000913 498928 = 0 (0.000)coding (874/963 nt) hemH ferrochelatase

CTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATATCGTTTATACGTGATTACGCCGATAACCACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAA  >  NC_000913/498244‑498809
                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                
cTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatata                                                                                                                                                                                                                                                                            >  1:83513‑M1/1‑278 (MQ=255)
  gTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatatt                                                                                                                                                                                                                                                                          >  1:314402‑M1/1‑276 (MQ=255)
                ggatggaAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATAAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgc                                                                                                                                                                                                                                                            <  2:206101‑M1/300‑39 (MQ=255)
                           ggCTCGCCTCTGATGGTTTACAGCCACCAGCAACAGCAGGAGCGGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGTCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaac                                                                                                                                                                                                                                                 <  2:256705‑M1/300‑50 (MQ=255)
                                       aTGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatga                                                                                                                                                                                                                                     >  2:289633‑M1/1‑239 (MQ=255)
                                          gTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatgg                                                                                                                                                                                                                                  >  1:299628‑M1/1‑236 (MQ=255)
                                                ctgcgagatgtgtataagagacaggcacaacgtttaccggagatgcccGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGTGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatc                                                                                                                                                                                                                            >  1:494319‑M1/24‑230 (MQ=255)
                                                       aGCAACAGCAGGCGCTGGCACAACGTTTACCGTAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATTGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgc                                                                                                                                                                                                                     <  2:229539‑M1/300‑78 (MQ=255)
                                                                       ggcgcgacgtttaCCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTTTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagc                                                                                                                                                                                                     <  2:8704‑M1/294‑94 (MQ=255)
                                                                       ggCACAACGTTTACCGTAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTTGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagc                                                                                                                                                                                                     <  1:165776‑M1/300‑94 (MQ=255)
                                                                       ggCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagc                                                                                                                                                                                                     <  2:297166‑M1/300‑94 (MQ=255)
                                                                             aCGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcg                                                                                                                                                                                               >  1:93293‑M1/1‑201 (MQ=255)
                                                                                      ggAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTTGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaac                                                                                                                                                                                      <  2:112988‑M1/300‑109 (MQ=255)
                                                                                      ggAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaac                                                                                                                                                                                      <  1:114976‑M1/300‑109 (MQ=255)
                                                                                        agaTGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactg                                                                                                                                                                                    <  1:639123‑M1/300‑111 (MQ=255)
                                                                                            gCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagtggtaaagaactgcgcg                                                                                                                                                                                >  1:280354‑M1/1‑186 (MQ=255)
                                                                                              ccGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCAGGGGATggcgggaaaaaatatgaatatattccggcgcttgatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgcc                                                                                                                                                                              >  2:566903‑M1/1‑184 (MQ=255)
                                                                                                               aGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgt                                                                                                                                                             >  1:252325‑M1/1‑167 (MQ=255)
                                                                                                                              tCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgctggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggtttttctca                                                                                                                                              >  2:91887‑M1/1‑152 (MQ=255)
                                                                                                                                 cTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCAGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatcc                                                                                                                                           >  2:412293‑M1/1‑149 (MQ=255)
                                                                                                                                       aGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaat                                                                                                                                     <  1:371663‑M1/300‑158 (MQ=255)
                                                                                                                                        gcgcCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtagcgaagagcctcgtcggcggttttcatcatccgtgaata                                                                                                                                    >  2:356521‑M1/1‑142 (MQ=255)
                                                                                                                                                          tGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcag                                                                                                                  >  1:101010‑M1/1‑124 (MQ=255)
                                                                                                                                                               gagCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCCTTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgc                                                                                                             >  2:509524‑M1/1‑119 (MQ=255)
                                                                                                                                                               gagCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgc                                                                                                             >  1:258060‑M1/1‑119 (MQ=255)
                                                                                                                                                               aaacatgTAGATCATATTGTGGGGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATTGGAGGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatattgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgc                                                                                                             <  2:77864‑M1/297‑182 (MQ=255)
                                                                                                                                                                 gCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcct                                                                                                           >  1:584810‑M1/1‑117 (MQ=255)
                                                                                                                                                                 gCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgntgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcct                                                                                                           >  1:485637‑M1/1‑117 (MQ=255)
                                                                                                                                                                  cATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcatg                                                                                                          >  2:254941‑M1/1‑116 (MQ=255)
                                                                                                                                                                     gTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCNGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctngtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggt                                                                                                       <  2:615460‑M1/300‑188 (MQ=255)
                                                                                                                                                                     gTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggt                                                                                                       <  1:112677‑M1/300‑188 (MQ=255)
                                                                                                                                                                       aGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtag                                                                                                     <  1:132386‑M1/300‑190 (MQ=255)
                                                                                                                                                                         aTCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtctcgctaaagctgagcggtaaagaacggagcgccgtcgcgaagagcctcgtcggnggttttcatcatccgtgaataatgcaaaaaggagtgcagcgtgcatgggtagag                                                                                                   >  2:497740‑M1/1‑109 (MQ=255)
                                                                                                                                                                               ttGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggnaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaa                                                                                             <  1:457815‑M1/300‑198 (MQ=255)
                                                                                                                                                                               ttGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaa                                                                                             >  1:338352‑M1/1‑103 (MQ=255)
                                                                                                                                                                                         gcCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggc                                                                                   <  2:607068‑M1/300‑208 (MQ=255)
                                                                                                                                                                                                   ccGCAATTCTCCTTTTCTACGGTCGGTGCGGTATGGGATGAACTTACACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcg                                                                         <  2:135143‑M1/300‑218 (MQ=255)
                                                                                                                                                                                                   ccGCAATTCTCCTGTTCTACTGTCGGTGCGGTATTGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCTGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacngatcgaaataacggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcg                                                                         <  2:161341‑M1/300‑218 (MQ=255)
                                                                                                                                                                                                            tcCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcgccgctaacg                                                                <  1:272554‑M1/300‑227 (MQ=255)
                                                                                                                                                                                                                gTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcgccgctaacgtctg                                                            >  1:585461‑M1/1‑70 (MQ=255)
                                                                                                                                                                                                                       ggtcggtGCAGTATGGGATGAACTGGCACGCATTCTGGCGCGTAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcgccgctaacgtctggtaaagc                                                     <  2:79261‑M1/300‑238 (MQ=255)
                                                                                                                                                                                                                       ggtcggnGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATNCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatccctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacaggcctgctgatgcgccgctaacgtctggtaaagc                                                     <  2:553573‑M1/300‑238 (MQ=255)
                                                                                                                                                                                                                            gtGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcgccgctaacgtctggtaaagcagacg                                                >  1:393597‑M1/1‑58 (MQ=255)
                                                                                                                                                                                                                            gtGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagaggtaaagaactnagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcgccgataacgtctggtaaagcagacg                                                >  2:161269‑M1/1‑58 (MQ=255)
                                                                                                                                                                                                                                ggTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcgccgctaacgtctggtaaagcagacggctg                                            >  1:513344‑M1/1‑54 (MQ=255)
                                                                                                                                                                                                                                   aTGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcnaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcgccgctaacgtctggtaaagcagacggctgtca                                         <  2:574516‑M1/300‑250 (MQ=255)
                                                                                                                                                                                                                                     gggATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcgccgctaacgtctggtaaagcagacggctgtcatc                                       <  1:506222‑M1/300‑252 (MQ=255)
                                                                                                                                                                                                                                         tGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcgccgctaacgtctggtaaagcagacggctgtcatccagc                                   <  1:52005‑M1/300‑256 (MQ=255)
                                                                                                                                                                                                                                          gAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcatgggtagagtttgaactcanagggctgctgatgcgccgctaacgtctggtaaagcagacggctgtcatccagca                                  >  2:496413‑M1/1‑44 (MQ=255)
                                                                                                                                                                                                                                                                          cATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcgccgctaacgtctggtaaagcagacggctgtcatccagcagcggatcgaactccgccccggcaataaaacag  <  1:55850‑M1/300‑289 (MQ=255)
                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                
CTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATATCGTTTATACGTGATTACGCCGATAACCACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAA  >  NC_000913/498244‑498809

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