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

ACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAG  >  NC_000913/498681‑499118                                                                                                                                                                                                                                                        |                                                                                                                                                                                               cgccgctgatggtttacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCg                                                                                                                                                                                             <  1:265701‑M2/4‑1 (MQ=255) cgccgctgatggtgtacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattcccctgttctacggtcggtgcggtatgggatgaactggcacgcagtctggcgcgcaaacgtagcattccggggatGGc                                                                                                                                                                                              <  2:237152‑M2/3‑1 (MQ=255)         atggtttacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGaaaaaa                                                                                                                                                                                     <  1:417016‑M2/12‑1 (MQ=255)                                         ggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGcc                                                                                                                                                     <  1:257175‑M2/44‑1 (MQ=255)                                           cacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCgg                                                                                                                                                   <  1:414331‑M2/46‑1 (MQ=255)                                                              gcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCTAAATGATgg                                                                                                                                 <  2:274308‑M2/64‑1 (MQ=255)                                                               cccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCt                                                                                                                               >  2:82304‑M2/186‑251 (MQ=255)                                                                  gtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCgg                                                                                                                                                                     <  2:275314‑M2/28‑1 (MQ=255)                                                                  gtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCgg                                                                                                                                                                     >  1:275314‑M2/183‑210 (MQ=255)                                                                       gctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCttgt                                                                                                                       >  2:129044‑M2/178‑251 (MQ=255)                                                                          gggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCttgtt                                                                                                                      <  2:79062‑M2/75‑1 (MQ=255)                                                                                   ctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCg                                                                                                           >  1:245629‑M2/166‑251 (MQ=255)                                                                                           cgccatcactggaaagcgccgtagatgaactcctggcagagcatgcagataatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggactgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCt                                                                                                   >  1:72849‑M2/158‑251 (MQ=255)                                                                                            gccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTg                                                                                                  <  1:267839‑M2/95‑1 (MQ=255)                                                                                                        aagctccgtagatgaactcctggcagagcatgtagagcatattgtggtgctgccgctgtatccgcaattcgcctgttctacggtcggtgcggtatgggatgaactggcacgcattcgggcgcgcaaacgtagcagtccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGAGGGCTAAGCTTGTTGCCGCGTATCGCTAATGCGGAGCGGTAAAGaa                                                                                      <  2:434024‑M2/107‑1 (MQ=255)                                                                                                            gccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTg                                                                                   >  1:21133‑M2/141‑250 (MQ=255)                                                                                                                             ggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAaga                                                                   <  1:325833‑M2/126‑1 (MQ=255)                                                                                                                              gcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGcc                                                                >  2:206‑M2/123‑251 (MQ=255)                                                                                                                                             catattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtc                                                 <  1:63148‑M2/144‑1 (MQ=255)                                                                                                                                              atattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtca                                                <  1:257468‑M2/145‑1 (MQ=255)                                                                                                                                                     ggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATcc                                          <  2:121602‑M2/151‑1 (MQ=255)                                                                                                                                                        gctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCCTCCGtgg                                      >  2:393685‑M2/97‑250 (MQ=255)                                                                                                                                                             cgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGGTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGaataat                                 <  2:45978‑M2/160‑1 (MQ=255)                                                                                                                                                                   atcctcaattctccttttctacggtctgtgctgtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCTGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCTTCGCGAAGATCCTCGTCGGCGGTTTTCATCATCCGGGAATAATGCaaa                            <  1:371795‑M2/165‑1 (MQ=255)                                                                                                                                                                         aattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGt                     <  1:82304‑M2/172‑1 (MQ=255)                                                                                                                                                                            tctccagttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGc                   <  2:444174‑M2/174‑1 (MQ=255)                                                                                                                                                                                  gttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTg             <  1:393685‑M2/180‑1 (MQ=255)                                                                                                                                                                                    tctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCt          >  1:269617‑M2/69‑251 (MQ=255)                                                                                                                                                                                       acggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCt          <  1:270394‑M2/183‑1 (MQ=255)                                                                                                                                                                                       acggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTg         <  2:90919‑M2/184‑1 (MQ=255)                                                                                                                                                                                          gtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTag    <  2:54531‑M2/189‑1 (MQ=255)                                                                                                                                                                                            cggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTagag  <  2:129013‑M2/191‑1 (MQ=255)                                                                                                                                                                                            cggtgcgggatgggatgaactggcacgcattctggcgcgcaaacgtagaattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAAACTATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCCTCATCCGTGAATAATGCAAAAAGGCGTGCCGCGGGCCTGGGTag    >  1:221263‑M2/61‑249 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                               ACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAG  >  NC_000913/498681‑499118

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: ‑