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

GTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAG  >  NC_000913/498650‑499215
                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                               
ctgtatgcctctgtctggatggaaggtggctcgccgctgatggtttacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAatata                                                                                                                                                                                                                                                                            >  1:83513‑M2/279‑300 (MQ=255)
  gtatgcctctgtctggatggaaggtggctcgccgctgatggtttacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATAtt                                                                                                                                                                                                                                                                          >  1:314402‑M2/277‑300 (MQ=255)
                ggatggaaggtggctcgccgctgatggtttacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagataaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGc                                                                                                                                                                                                                                                            <  2:206101‑M2/38‑1 (MQ=255)
                           ggctcgcctctgatggtttacagccaccagcaacagcaggagcgggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctgtcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAAc                                                                                                                                                                                                                                                 <  2:256705‑M2/49‑1 (MQ=255)
                                       atggtttacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAatga                                                                                                                                                                                                                                     >  2:289633‑M2/240‑300 (MQ=255)
                                          gtttacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATgg                                                                                                                                                                                                                                  >  1:299628‑M2/237‑300 (MQ=255)
                                                ctgcgagatgtgtataagagacaggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgtggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATc                                                                                                                                                                                                                            >  1:494319‑M2/231‑300 (MQ=255)
                                                       agcaacagcaggcgctggcacaacgtttaccgtagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtattggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGc                                                                                                                                                                                                                     <  2:229539‑M2/77‑1 (MQ=255)
                                                                       ggcgcgacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattttggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGc                                                                                                                                                                                                     <  2:8704‑M2/93‑1 (MQ=255)
                                                                       ggcacaacgtttaccgtagatgcccgtagcgctgggaatgagctacggctcgccatcacttgaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGc                                                                                                                                                                                                     <  1:165776‑M2/93‑1 (MQ=255)
                                                                       ggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGc                                                                                                                                                                                                     <  2:297166‑M2/93‑1 (MQ=255)
                                                                             acgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCg                                                                                                                                                                                               >  1:93293‑M2/202‑300 (MQ=255)
                                                                                      ggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactccttgcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAAc                                                                                                                                                                                      <  2:112988‑M2/108‑1 (MQ=255)
                                                                                      ggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAAc                                                                                                                                                                                      <  1:114976‑M2/108‑1 (MQ=255)
                                                                                        agatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTg                                                                                                                                                                                    <  1:639123‑M2/110‑1 (MQ=255)
                                                                                            gcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGTGGTAAAGAACTGcgcg                                                                                                                                                                                >  1:280354‑M2/187‑300 (MQ=255)
                                                                                              ccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattcaggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTGATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGcc                                                                                                                                                                              >  2:566903‑M2/185‑300 (MQ=255)
                                                                                                               agctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtcgt                                                                                                                                                             >  1:252325‑M2/168‑300 (MQ=255)
                                                                                                                              tcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGCTGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCggcggtttttctca                                                                                                                                              >  2:91887‑M2/153‑295 (MQ=255)
                                                                                                                                 ctggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattcaggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATcc                                                                                                                                           >  2:412293‑M2/150‑300 (MQ=255)
                                                                                                                                       agcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGaat                                                                                                                                     <  1:371663‑M2/157‑1 (MQ=255)
                                                                                                                                        gcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTAGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGaata                                                                                                                                    >  2:356521‑M2/143‑300 (MQ=255)
                                                                                                                                                          tggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAg                                                                                                                  >  1:101010‑M2/125‑300 (MQ=255)
                                                                                                                                                               gagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgccttctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGc                                                                                                             >  2:509524‑M2/120‑300 (MQ=255)
                                                                                                                                                               gagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGc                                                                                                             >  1:258060‑M2/120‑300 (MQ=255)
                                                                                                                                                               aaacatgtagatcatattgtggggctgccgctttatccgcaattctcctgttctacggtcggtgcggtattggaggaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATTGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGc                                                                                                             <  2:77864‑M2/181‑1 (MQ=255)
                                                                                                                                                                 gcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCt                                                                                                           >  1:584810‑M2/118‑300 (MQ=255)
                                                                                                                                                                 gcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGNTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCt                                                                                                           >  1:485637‑M2/118‑300 (MQ=255)
                                                                                                                                                                  catgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCgtgcatg                                                                                                          >  2:254941‑M2/117‑297 (MQ=255)
                                                                                                                                                                     gtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcngtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTNGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGt                                                                                                       <  2:615460‑M2/187‑1 (MQ=255)
                                                                                                                                                                     gtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGt                                                                                                       <  1:112677‑M2/187‑1 (MQ=255)
                                                                                                                                                                       agatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTag                                                                                                     <  1:132386‑M2/189‑1 (MQ=255)
                                                                                                                                                                         atcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTCTCGCTAAAGCTGAGCGGTAAAGAACGGAGCGCCGTCGCGAAGAGCCTCGTCGGNGGTTTTCATCATCCGTGAATAATGCAAAAAGGAGTGCAGCGTGCATGGGTagag                                                                                                   >  2:497740‑M2/110‑300 (MQ=255)
                                                                                                                                                                               ttgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGNAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGaa                                                                                             <  1:457815‑M2/197‑1 (MQ=255)
                                                                                                                                                                               ttgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGaa                                                                                             >  1:338352‑M2/104‑300 (MQ=255)
                                                                                                                                                                                         gccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGgc                                                                                   <  2:607068‑M2/207‑1 (MQ=255)
                                                                                                                                                                                                   ccgcaattctccttttctacggtcggtgcggtatgggatgaacttacacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATgcg                                                                         <  2:135143‑M2/217‑1 (MQ=255)
                                                                                                                                                                                                   ccgcaattctcctgttctactgtcggtgcggtattggatgaactggcacgcattctggcgcgcaaacgtagcattcctgggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACNGATCGAAATAACGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATgcg                                                                         <  2:161341‑M2/217‑1 (MQ=255)
                                                                                                                                                                                                            tcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACg                                                                <  1:272554‑M2/226‑1 (MQ=255)
                                                                                                                                                                                                                gttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTg                                                            >  1:585461‑M2/71‑300 (MQ=255)
                                                                                                                                                                                                                       ggtcggtgcagtatgggatgaactggcacgcattctggcgcgtaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGc                                                     <  2:79261‑M2/237‑1 (MQ=255)
                                                                                                                                                                                                                       ggtcggngcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcatnccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCCCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGCCTGCTGATGCGCCGCTAACGTCTGGTAAAGc                                                     <  2:553573‑M2/237‑1 (MQ=255)
                                                                                                                                                                                                                            gtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACg                                                >  1:393597‑M2/59‑300 (MQ=255)
                                                                                                                                                                                                                            gtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGAGGTAAAGAACTNAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGATAACGTCTGGTAAAGCAGACg                                                >  2:161269‑M2/59‑300 (MQ=255)
                                                                                                                                                                                                                                ggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTg                                            >  1:513344‑M2/55‑300 (MQ=255)
                                                                                                                                                                                                                                   atgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCNAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCa                                         <  2:574516‑M2/249‑1 (MQ=255)
                                                                                                                                                                                                                                     gggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATc                                       <  1:506222‑M2/251‑1 (MQ=255)
                                                                                                                                                                                                                                         tgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCcagc                                   <  1:52005‑M2/255‑1 (MQ=255)
                                                                                                                                                                                                                                          gaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCATGGGTAGAGTTTGAACTCANAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCcagca                                  >  2:496413‑M2/45‑300 (MQ=255)
                                                                                                                                                                                                                                                                          cattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAg  <  1:55850‑M2/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                               
GTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAG  >  NC_000913/498650‑499215

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.