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

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 498521	0 (0.000)	33 (0.700)	33/492	0.5	100%	coding (467/963 nt)	hemH	ferrochelatase
	?	NC_000913	498928 =	0 (0.000)					coding (874/963 nt)	hemH	ferrochelatase

CGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATATCGTTTATACGTGATTACGCCGATAACCACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTT  >  NC_000913/498275‑498712                                                                                                                                                                                                                                                       |                                                                                                                                                                                                cgccgcTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcg                                                                                                                                                                                             <  1:265701‑M1/251‑5 (MQ=255) cgccgcTGATGGTGTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCCCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCAGTCTGGCGCGCAAACGTAGCATTCCGGGGATggc                                                                                                                                                                                              <  2:237152‑M1/250‑4 (MQ=255)         aTGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaa                                                                                                                                                                                     <  1:417016‑M1/251‑13 (MQ=255)                                         ggCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgcc                                                                                                                                                     <  1:257175‑M1/251‑45 (MQ=255)                                           cacaACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccgg                                                                                                                                                   <  1:414331‑M1/251‑47 (MQ=255)                                                              gCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatctaaatgatgg                                                                                                                                 <  2:274308‑M1/250‑65 (MQ=255)                                                               cccGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggct                                                                                                                               >  2:82304‑M1/1‑185 (MQ=255)                                                                  gTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccgg                                                                                                                                                                     <  2:275314‑M1/210‑29 (MQ=255)                                                                  gTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccgg                                                                                                                                                                     >  1:275314‑M1/1‑182 (MQ=255)                                                                       gcTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgt                                                                                                                       >  2:129044‑M1/1‑177 (MQ=255)                                                                          gggAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgtt                                                                                                                      <  2:79062‑M1/249‑76 (MQ=255)                                                                                   cTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcg                                                                                                           >  1:245629‑M1/1‑165 (MQ=255)                                                                                           cGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGCAGATAATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGACTGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagct                                                                                                   >  1:72849‑M1/1‑157 (MQ=255)                                                                                            gCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctg                                                                                                  <  1:267839‑M1/251‑96 (MQ=255)                                                                                                        aaGCTCCGTAGATGAACTCCTGGCAGAGCATGTAGAGCATATTGTGGTGCTGCCGCTGTATCCGCAATTCGCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCGGGCGCGCAAACGTAGCAGTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgagggctaagcttgttgccgcgtatcgctaatgcggagcggtaaagaa                                                                                      <  2:434024‑M1/251‑108 (MQ=255)                                                                                                            gcCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactg                                                                                   >  1:21133‑M1/1‑140 (MQ=255)                                                                                                                             ggCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaaga                                                                   <  1:325833‑M1/249‑127 (MQ=255)                                                                                                                              gCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcc                                                                >  2:206‑M1/1‑122 (MQ=255)                                                                                                                                             cATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttc                                                 <  1:63148‑M1/251‑145 (MQ=255)                                                                                                                                              atatTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttca                                                <  1:257468‑M1/251‑146 (MQ=255)                                                                                                                                                     ggtgCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatcc                                          <  2:121602‑M1/250‑152 (MQ=255)                                                                                                                                                        gctgcCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcctccgtgg                                      >  2:393685‑M1/1‑96 (MQ=255)                                                                                                                                                             cGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttggtgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataat                                 <  2:45978‑M1/251‑161 (MQ=255)                                                                                                                                                                   aTCCTCAATTCTCCTTTTCTACGGTCTGTGCTGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggctggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccttcgcgaagatcctcgtcggcggttttcatcatccgggaataatgcaaa                            <  1:371795‑M1/250‑166 (MQ=255)                                                                                                                                                                         aaTTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgt                     <  1:82304‑M1/251‑173 (MQ=255)                                                                                                                                                                            tctcCAGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgc                   <  2:444174‑M1/250‑175 (MQ=255)                                                                                                                                                                                  gTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtg             <  1:393685‑M1/250‑181 (MQ=255)                                                                                                                                                                                    tCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcct          >  1:269617‑M1/1‑68 (MQ=255)                                                                                                                                                                                       aCGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctg         <  2:90919‑M1/249‑185 (MQ=255)                                                                                                                                                                                       aCGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcct          <  1:270394‑M1/248‑184 (MQ=255)                                                                                                                                                                                          gtcggtGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtag    <  2:54531‑M1/251‑190 (MQ=255)                                                                                                                                                                                            cggtGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagag  <  2:129013‑M1/251‑192 (MQ=255)                                                                                                                                                                                            cggtGCGGGATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGAATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaaactatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcctcatccgtgaataatgcaaaaaggcgtgccgcgggcctgggtag    >  1:221263‑M1/1‑60 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                CGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATATCGTTTATACGTGATTACGCCGATAACCACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTT  >  NC_000913/498275‑498712

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑