Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	498,522	Δ406 bp	coding (468‑873/963 nt)	hemH →	ferrochelatase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	498522	498927	406	28 [0]	[0] 28	hemH	ferrochelatase

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

CGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATATCGTTTATACGTGATTACGCCGATAACCACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTT  >  NC_000913/498275‑498752                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        cgccgcTGATGGGTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGAGGTATGGGATGAAATGGCACGCATTCTGGCGCGCAAACGTAGCATTACGGGGCTggcg                                                                                                                                                                                                                                     >  2:7102‑M1/1‑247 (MQ=255)           ggTTTACAGCCGCCAGCAACAGCAGGCGCTGGCCCAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATACGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaat                                                                                                                                                                                                                            >  2:178711‑M1/1‑237 (MQ=255)                           cAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGGCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcg                                                                                                                                                                                                           <  1:152728‑M1/251‑31 (MQ=255)                             aCAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgct                                                                                                                                                                                                         <  1:122346‑M1/251‑33 (MQ=255)                                   ggCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgc                                                                                                                                                                                                   <  1:11939‑M1/251‑39 (MQ=255)                                                cGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCACCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacat                                                                                                                                                                                      <  2:174704‑M1/251‑52 (MQ=255)                                                    tACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTCCGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACTTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatc                                                                                                                                                                                   <  2:49178‑M1/250‑55 (MQ=255)                                                     aCCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaataggaatatattccggcgcttaatgccacgccggaacatatcg                                                                                                                                                                                  >  1:135148‑M1/1‑195 (MQ=255)                                                                  gTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaat                                                                                                                                                                    >  2:179004‑M1/1‑182 (MQ=255)                                                                          gggAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgc                                                                                                                                                            >  1:130336‑M1/1‑174 (MQ=255)                                                                                     aCGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgc                                                                                                                                                  <  2:46047‑M1/250‑88 (MQ=255)                                                                                                                 aGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgcc                                                                                                                     <  2:203683‑M1/251‑117 (MQ=255)                                                                                                                             ggCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaaga                                                                                                           <  1:31254‑M1/249‑127 (MQ=255)                                                                                                                              gCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagcgcc                                                                                                        >  1:158178‑M1/1‑122 (MQ=255)                                                                                                                                agagCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctc                                                                                                      >  1:249443‑M1/1‑120 (MQ=255)                                                                                                                                                       tgctgcCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcctccggg                                                                               >  2:89210‑M1/1‑97 (MQ=255)                                                                                                                                                                  tATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaa                                                                     >  1:77471‑M1/1‑86 (MQ=255)                                                                                                                                                                   aTCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaa                                                                   <  1:199951‑M1/251‑167 (MQ=255)                                                                                                                                                                      cGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaagg                                                                <  1:120396‑M1/251‑170 (MQ=255)                                                                                                                                                                                             ggtGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgtataatgcaaaaaggcgtgccgcgtgcctgggtagagt                                         >  2:39393‑M1/1‑59 (MQ=255)                                                                                                                                                                                             ggtGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgggcctgggtagagt                                         >  1:215320‑M1/1‑59 (MQ=255)                                                                                                                                                                                                      ggggaTGAACTGGCACGGATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaaggagtgctaatcttggtgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtctgcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagttggaactca                                <  2:231021‑M1/250‑202 (MQ=255)                                                                                                                                                                                                          attaactgGCACGCATTCTGGCGCGCAAACGTAGGCTTCCGGGGATggggggaaaaaataagaatatattccggcgcttaatgccacgccggaacataacggaatgatggctaaacttgttgccgcgttacgctaaagctgggccgtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatcccgggataatgccaaaaggcgtgcagcgttcctgggtagagtttgaactcacagg                            <  1:186969‑M1/248‑206 (MQ=255)                                                                                                                                                                                                             aaCTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggct                         <  2:77471‑M1/251‑209 (MQ=255)                                                                                                                                                                                                                     cGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgc                 <  2:134137‑M1/251‑217 (MQ=255)                                                                                                                                                                                                                                 cgcAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggttaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcggtcagcgggcctggggagagttttaactcacagggctgctgatgcgccgctaacgtc     >  2:176974‑M1/1‑23 (MQ=255)                                                                                                                                                                                                                                  gcAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcgccgctaacgtct    >  1:179255‑M1/1‑22 (MQ=255)                                                                                                                                                                                                                                    aaaCGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctggggagagtttgaactcacagggctgctgatgcgccgctaacgtctgg  >  1:189665‑M1/1‑20 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        CGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATATCGTTTATACGTGATTACGCCGATAACCACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTT  >  NC_000913/498275‑498752

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